RNADiscrepancy 0.0.2__py3-none-any.whl

RNADiscrepancy/measures.py

@@ -0,0 +1,171 @@
+import math
+
+
+import numpy as np
+from scipy.spatial.transform import Rotation
+
+from RNADiscrepancy import references
+from RNADiscrepancy.rna import Atom, Nucleotide
+
+
+def c1prime_distance_difference(pair1, pair2):
+    """Returns the difference in distance between the C1' atoms (connecting sugar to residue) in the two pairs
+    always positive
+    """
+
+
+    (nt1_left, nt1_right), (nt2_left, nt2_right) = pair1, pair2
+    dist_c1 = abs(nt1_left["C1'"].distance(nt1_right["C1'"]) - nt2_left["C1'"].distance(nt2_right["C2'"]))
+
+    return dist_c1
+
+
+def c1prime_translation(pair1, pair2):
+    """
+    Here we have to do things in regards to the _left and _right side symmetricaly
+
+    1'. Use N1/N9 of their reference as (0, 0)
+    2'. Rotate to align the first base to the reference (and apply rotation to the pair)
+    """
+
+    (nt1_left, nt1_right), (nt2_left, nt2_right) = pair1, pair2
+
+    ref_nt1_left = references.nucleotide[nt1_left.name]
+    ref_nt1_right = references.nucleotide[nt1_right.name]
+    ref_nt2_left = references.nucleotide[nt2_left.name]
+    ref_nt2_right = references.nucleotide[nt2_right.name]
+
+    #Retrieve position of nucleotide + ref for alignment
+    nt1_left_pos = []
+    nt1_left_pos_ref = []
+    for atom_name, atom in nt1_left.items():
+        try:
+            nt1_left_pos_ref.append(ref_nt1_left[atom_name].position)
+            nt1_left_pos.append(atom.position)
+        except KeyError:
+            continue
+
+    nt1_right_pos = []
+    nt1_right_pos_ref = []
+    for atom_name, atom in nt1_right.items():
+        try:
+            nt1_right_pos_ref.append(ref_nt1_right[atom_name].position)
+            nt1_right_pos.append(atom.position)
+        except KeyError:
+            continue
+
+
+    nt2_left_pos = []
+    nt2_left_pos_ref = []
+    for atom_name, atom in nt2_left.items():
+        try:
+            nt2_left_pos_ref.append(ref_nt2_left[atom_name].position)
+            nt2_left_pos.append(atom.position)
+        except KeyError:
+            continue
+
+    nt2_right_pos = []
+    nt2_right_pos_ref = []
+    for atom_name, atom in nt2_right.items():
+        try:
+            nt2_right_pos_ref.append(ref_nt2_right[atom_name].position)
+            nt2_right_pos.append(atom.position)
+        except KeyError:
+            continue
+
+ 
+
+    nt1_left_rotation = best_rotation(nt1_left_pos_ref, np.array(ref_nt1_left.ribose_fulcrum.position), 
+                                      nt1_left_pos, np.array(nt1_left.ribose_fulcrum.position))[0]
+    nt1_right_rotation = best_rotation(nt1_right_pos_ref, np.array(ref_nt1_right.ribose_fulcrum.position), 
+                                      nt1_right_pos, np.array(nt1_right.ribose_fulcrum.position))[0]
+    nt2_left_rotation = best_rotation(nt2_left_pos_ref, np.array(ref_nt2_left.ribose_fulcrum.position), 
+                                      nt2_left_pos, np.array(nt2_left.ribose_fulcrum.position))[0]
+    nt2_right_rotation = best_rotation(nt2_right_pos_ref, np.array(ref_nt2_right.ribose_fulcrum.position), 
+                                      nt2_right_pos, np.array(nt2_right.ribose_fulcrum.position))[0]
+
+
+    c1p_1right = nt1_right["C1'"].position
+    c1p_1right = (c1p_1right - np.array(nt1_left.ribose_fulcrum.position)) @ nt1_left_rotation
+    c1p_2right = nt2_right["C1'"].position
+    c1p_2right = (c1p_2right - np.array(nt2_left.ribose_fulcrum.position)) @ nt2_left_rotation
+    c1p_t1 = math.sqrt((c1p_1right[0] - c1p_2right[0])**2 + (c1p_1right[1] - c1p_2right[1])**2)
+
+
+    c1p_1left = nt1_left["C1'"].position
+    c1p_1left = (c1p_1left - np.array(nt1_right.ribose_fulcrum.position)) @ nt1_right_rotation
+    c1p_2left = nt2_left["C1'"].position
+    c1p_2left = (c1p_2left - np.array(nt2_right.ribose_fulcrum.position)) @ nt2_right_rotation
+    c1p_t2 = math.sqrt((c1p_1left[0] - c1p_2left[0])**2 + (c1p_1left[1] - c1p_2left[1])**2)
+
+
+    bond_1right = [nt1_right["C1'"].position, nt1_right.ribose_fulcrum.position]
+    bond_1right = (bond_1right - np.array(nt1_left.ribose_fulcrum.position)) @ nt1_left_rotation
+    bond_2right = [nt2_right["C1'"].position, nt2_left.ribose_fulcrum.position]
+    bond_2right = (bond_2right - np.array(nt2_left.ribose_fulcrum.position)) @ nt2_left_rotation
+    bond_1right[:, 2:] = 0
+    bond_2right[:, 2:] = 0
+    output = best_rotation(bond_1right[:2, :], 0, bond_2right[:2, :], 0)
+    theta1 = math.acos(output[0][0][0])
+
+    bond_1left = [nt1_left["C1'"].position, nt1_left.ribose_fulcrum.position]
+    bond_1left = (bond_1left - np.array(nt1_right.ribose_fulcrum.position)) @ nt1_right_rotation
+    bond_2left = [nt2_left["C1'"].position, nt2_left.ribose_fulcrum.position]
+    bond_2left = (bond_2left - np.array(nt2_right.ribose_fulcrum.position)) @ nt2_right_rotation
+    bond_1left[:, 2:] = 0
+    bond_2left[:, 2:] = 0
+    output = best_rotation(bond_1left, 0, bond_2left, 0)
+    theta2 = math.acos(output[0][0][0])
+
+
+    return c1p_t1, c1p_t2, theta1, theta2
+
+
+
+
+def isodiscrepancy(pair1, pair2):
+    """
+    Return IDI between pairs (nt1_left, nt1_right) and (nt2_left, nt2_right)
+    it assumes the orientation is to align nt1_left with nt2_left  / nt1_right with nt2_right
+
+    we follow https://doi.org/10.1093/nar/gkp011
+
+    1) delta c1' distances
+    2) align on one N1/N9 and see other C1' translation distance
+    3) rotation in plane to orient same second base in relation to first
+    """
+
+    c1prime_dist_delta = c1prime_distance_difference(pair1, pair2)
+    c1prime_t1, c1prime_t2, theta1, theta2 = c1prime_translation(pair1, pair2)
+
+    return math.sqrt(c1prime_dist_delta**2 
+                     + (c1prime_t1**2 + c1prime_t2**2)/2 
+                     + ((2 * theta1)**2 + (2 * theta2)**2)/2
+                     )
+
+def best_rotation(A, center_A, B, center_B):
+    A_centered = A - center_A
+    B_centered = B - center_B
+    rotation, _ = Rotation.align_vectors(A_centered, B_centered)
+    rotation = rotation.as_matrix()
+    B2A = B_centered @ rotation + center_A
+    difference = A_centered - B2A + center_A
+    sum_square_error = float(np.sum(difference ** 2))
+    RMSD = float(np.sqrt(sum_square_error / len(A)))
+    return rotation, B2A, center_A, center_B, sum_square_error, RMSD
+
+
+def best_mean_centered_rotation(A, B):
+    """Least-squares rigid-body fit transform of two numpy array
+    The center of mass of the nucleotide is at (0, 0, 0)
+    returns 
+        rotation :  matrix
+        B2A : B moved to over A
+        mean_A, 
+        mean_B,
+        sum_square_error,
+        RMSD
+    """
+    mean_A = A.mean(axis=0)
+    mean_B = B.mean(axis=0)
+    return best_rotation(A, mean_A, B, mean_B)

RNADiscrepancy/references.py

@@ -0,0 +1,66 @@
+from RNADiscrepancy import rna
+
+"""
+Position of idealized nucleotides heavy atoms with center at (0,0,0)
+It follows FR3D standard orientation with the Watson-Crick edge in 
+the positive x and y quadrant.
+"""
+nucleotide = {}
+
+ref_heavy_atoms_positions = {
+    'A' : {
+        'N9' : ([-1.110515, -1.823319, 0.0] ),
+        'C4' : ([0.007975, -1.020192, 0.0] ),
+        'N3' : ([1.298514, -1.383716, 0.0] ),
+        'N1' : ([1.754329, 1.001072, 0.0] ),
+        'C6' : ([0.453129, 1.320188, 0.0] ),
+        'N6' : ([0.092589, 2.623715, 0.0] ),
+        'C8' : ([-2.200426, -0.989732, 0.0] ),
+        'C5' : ([-0.504186, 0.282479, 0.0] ),
+        'C2' : ([2.092069, -0.307357, 0.0] ),
+        'N7' : ([-1.883477, 0.296861, 0.0] ),
+        "C1'" :([ -1.110515,  (-1.823319 - 1.48),   0.000000])
+    },
+    'C' : {
+        'N1' : ([-0.380579, -1.484583, 0.0] ),
+        'C2' : ([0.908756, -0.88533, 0.0] ),
+        'O2' : ([1.888871, -1.605366, 0.0] ),
+        'N3' : ([0.931558, 0.493192, 0.0] ),
+        'C4' : ([-0.197209, 1.170043, 0.0] ),
+        'N4' : ([-0.099223, 2.524619, 0.0] ),
+        'C6' : ([-1.542632, -0.786228, 0.0] ),
+        'C5' : ([-1.509542, 0.573654, 0.0] ),
+        "C1'" : ( [ -0.380579,  (-1.484583-1.48),   0.000000] )
+    },
+    'G' : {
+        'N9' : ([-1.45668, -1.711888, 0.0] ),
+        'C4' : ([-0.339093, -0.921183, 0.0] ),
+        'N3' : ([0.947138, -1.370892, 0.0] ),
+        'N1' : ([1.440727, 0.946157, 0.0] ),
+        'C6' : ([0.110442, 1.47918, 0.0] ),
+        'O6' : ([-0.059694, 2.682325, 0.0] ),
+        'C8' : ([-2.54532, -0.870005, 0.0] ),
+        'C5' : ([-0.831466, 0.385014, 0.0] ),
+        'C2' : ([1.802732, -0.378447, 0.0] ),
+        'N7' : ([-2.208072, 0.408244, 0.0] ),
+        'N2' : ([3.139287, -0.648506, 0.0] ),
+        "C1'" : ( [ -1.456680,  (-1.711888-1.48),   0.000000])
+    },
+    'U' : {
+        'N1' : ([-0.32642, -1.514422, 0.0] ),
+        'C2' : ([0.933866, -0.925171, 0.0] ),
+        'O2' : ([1.96607, -1.562816, 0.0] ),
+        'N3' : ([0.868896, 0.45947, 0.0] ),
+        'C4' : ([-0.266417, 1.293935, 0.0] ),
+        'O4' : ([-0.149264, 2.505638, 0.0] ),
+        'C6' : ([-1.504181, -0.804727, 0.0] ),
+        'C5' : ([-1.52255, 0.548093, 0.0] ),
+        "C1'" : ([ -0.326420,  (-1.514422-1.48),   0.000000])
+    }
+}
+
+for nt_name, v in ref_heavy_atoms_positions.items():
+    atoms = []
+    for atom_name, pos in v.items():
+       atoms.append(rna.Atom(atom_name, atom_name[0], *pos) )
+    nucleotide[nt_name] = rna.Nucleotide(nt_name, atoms)

RNADiscrepancy/rna.py

@@ -0,0 +1,57 @@
+import math
+
+class Atom:
+    def __init__(self, name, element, x, y, z):
+        self.name = name
+        self.element = element
+        self.x = x
+        self.y = y
+        self.z = z
+
+    @property
+    def position(self):
+        return self.x, self.y, self.z
+
+    def distance(self, other):
+        x, y, z = self.position
+        x2, y2, z2 = other.position
+        return math.sqrt((x - x2)**2 + (y - y2)**2 + (z - z2)**2)
+
+
+    def __eq__(self, other):
+        if isinstance(other, self.__class__):
+            return self.name == other.name
+        else:
+            return False
+
+    def __neq__(self, other):
+        return not self.__eq__(other)
+
+    def __str__(self):
+        return f"Atom class: {self.name}, {self.element}, {self.x}, {self.y}, {self.z}"
+
+    def __repr__(self):
+        return f"Atom class: {self.name}, {self.element}, {self.x}, {self.y}, {self.z}"
+
+
+
+class Nucleotide(dict):
+    def __init__(self, name, atoms):
+        self.name = name
+        for atom in atoms:
+            self[atom.name] = atom
+
+        if name not in ('A', 'C', 'G', 'U'):
+            raise Exception(f"WTF this should be ACGU not {name}")
+
+        #Set fulcrum
+        if name in ('C', 'U'):
+            try:
+                self.ribose_fulcrum = self['N1']
+            except KeyError:
+                self.ribose_fulcrum = None
+        elif name in ('A', 'G'):
+            try:
+                self.ribose_fulcrum = self['N9']
+            except KeyError:
+                self.ribose_fulcrum = None

RNADiscrepancy/test.py

@@ -0,0 +1,61 @@
+from collections import defaultdict
+import pickle
+
+
+import measures
+import rna
+
+
+def read_samples(path):
+    with open(path, 'rb') as f:
+        data = pickle.load(f)
+    loaded_data = defaultdict(list)
+
+    for (lw, nt1, nt2), values in data.items():
+        if not (nt1 in ('A', 'C', 'G', 'U') and nt2 in ('A', 'C', 'G', 'U')):
+            continue
+        for instance in values:
+            atoms = []
+            for name, element, x, y, z in zip(*(instance[0][k] for k in ('atoms_name', 'atoms_element', 'atoms_x', 'atoms_y', 'atoms_z'))):
+                atoms.append(rna.Atom(name, element, x, y, z))
+            nt1_nuc = rna.Nucleotide(nt1, atoms)
+
+            atoms = []
+            for name, element, x, y, z in zip(*(instance[1][k] for k in ('atoms_name', 'atoms_element', 'atoms_x', 'atoms_y', 'atoms_z'))):
+                atoms.append(rna.Atom(name, element, x, y, z))
+            nt2_nuc = rna.Nucleotide(nt2, atoms)
+
+            if nt1_nuc.ribose_fulcrum is None or nt2_nuc.ribose_fulcrum is None:
+                continue
+
+            loaded_data[lw, nt1, nt2].append((nt1_nuc, nt2_nuc))
+
+    return loaded_data
+
+
+def compute_isodiscrepancy(pair1, pair2):
+    return measures.isodiscrepancy(pair1, pair2)
+
+def main():
+
+    path_samples = '../../test/10samples.pickle'
+    data = read_samples(path_samples)
+    
+    print(data['CWW', 'U', 'A'][0][0].name, data['CWW', 'C', 'G'][0][1].name)
+    print(compute_isodiscrepancy(data['CWW', 'U', 'A'][0], data['CWW', 'C', 'G'][0]))
+    print("\n")
+    print(compute_isodiscrepancy(data['CWW', 'A', 'U'][0], data['TWW', 'G', 'G'][0]))
+
+    U, C = data['TWW', 'U', 'C'][0]
+    for atom in C:
+        print(f'"{C[atom].name}"', end=",")
+    print()
+    for atom in C:
+        print(f'"{C[atom].element}"', end=",")
+    print()
+    for atom in C:
+        print(f"{C[atom].position}", end=",")
+
+
+if __name__ == '__main__':
+    main()

RNADiscrepancy/utilities.py

@@ -0,0 +1,7 @@
+from RNADiscrepancy import rna
+
+def rawdata_format_nucleotide(nucleotide_name, list_atoms_names, list_atoms_elements, list_positions):
+    atoms = []
+    for name, element, pos in zip(list_atoms_names, list_atoms_elements, list_positions):
+        atoms.append(rna.Atom(name, element, *pos))
+    return rna.Nucleotide(nucleotide_name, atoms)

rnadiscrepancy-0.0.2.dist-info/METADATA

@@ -0,0 +1,63 @@
+Metadata-Version: 2.4
+Name: RNADiscrepancy
+Version: 0.0.2
+Summary: Compute discrepancy between 3D RNA (sub)structures
+Project-URL: Homepage, https://codeberg.org/vreinharz/RNADiscrepancy
+Project-URL: Issues, https://codeberg.org/vreinharz/RNADiscrepancy/issues
+Author-email: Vladimir Reinharz <reinharz.vladimir@uqam.ca>
+License-Expression: GPL-3.0
+License-File: LICENSE
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python :: 3
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+
+#RNADiscrepancy
+
+isostericity
+RMSD
+maybe more?
+
+
+```python
+>>> from RNADiscrepancy.utilities import rawdata_format_nucleotide
+>>> from RNADiscrepancy.measures import isodiscrepancy
+>>> cwwAU_A_list_atoms = ["P","OP1","OP2","O5'","C5'","C4'","O4'","C3'","O3'","C2'","O2'","C1'","N9","C8","N7","C5","C6","N6","N1","C2","N3","C4"]
+>>> cwwAU_A_list_atoms_elements = ["P","O","O","O","C","C","O","C","O","C","O","C","N","C","N","C","C","N","N","C","N","C"]
+>>> cwwAU_A_list_positions = [(-46.6870002746582, -80.67400360107422, -42.457000732421875),(-45.249000549316406, -80.74099731445312, -42.09299850463867),(-47.303001403808594, -79.36100006103516, -42.77299880981445),(-46.92900085449219, -81.6449966430664, -43.694000244140625),(-47.27000045776367, -83.01000213623047, -43.49399948120117),(-47.70500183105469, -83.62999725341797, -44.79600143432617),(-49.13199996948242, -83.92400360107422, -44.75299835205078),(-47.60100173950195, -82.73999786376953, -46.02299880981445),(-46.290000915527344, -82.58599853515625, -46.52899932861328),(-48.512001037597656, -83.47599792480469, -46.97999954223633),(-47.96799850463867, -84.72100067138672, -47.375),(-49.6870002746582, -83.72000122070312, -46.04600143432617),(-50.590999603271484, -82.5719985961914, -46.012001037597656),(-50.654998779296875, -81.52100372314453, -45.125),(-51.59299850463867, -80.6449966430664, -45.40800094604492),(-52.19300079345703, -81.16400146484375, -46.55400085449219),(-53.26100158691406, -80.71800231933594, -47.35499954223633),(-53.959999084472656, -79.60900115966797, -47.10100173950195),(-53.59199905395508, -81.46099853515625, -48.4379997253418),(-52.902000427246094, -82.58399963378906, -48.68299865722656),(-51.887001037597656, -83.11399841308594, -48.000999450683594),(-51.58000183105469, -82.34600067138672, -46.9370002746582)]
+>>> cwwAU_U_list_atoms = ["P","OP1","OP2","O5'","C5'","C4'","O4'","C3'","O3'","C2'","O2'","C1'","N1","C2","O2","N3","C4","O4","C5","C6"]
+>>> cwwAU_U_list_atoms_elements = ["P","O","O","O","C","C","O","C","O","C","O","C","N","C","O","N","C","O","C","C"]
+>>> cwwAU_U_list_positions = [(-62.04399871826172, -78.197998046875, -52.75299835205078),(-63.35599899291992, -78.33399963378906, -53.4370002746582),(-61.986000061035156, -77.62999725341797, -51.382999420166016),(-61.347999572753906, -79.62899780273438, -52.72700119018555),(-61.56800079345703, -80.56199645996094, -53.777000427246094),(-60.268001556396484, -81.21900177001953, -54.17399978637695),(-59.19599914550781, -80.2249984741211, -54.268001556396484),(-59.689998626708984, -82.20500183105469, -53.178001403808594),(-60.340999603271484, -83.45899963378906, -53.125),(-58.26900100708008, -82.29499816894531, -53.6879997253418),(-58.19900131225586, -82.8759994506836, -54.974998474121094),(-57.983001708984375, -80.80400085449219, -53.79199981689453),(-57.66400146484375, -80.25, -52.470001220703125),(-56.564998626708984, -80.7750015258789, -51.81100082397461),(-55.832000732421875, -81.61399841308594, -52.30500030517578),(-56.35599899291992, -80.28099822998047, -50.54899978637695),(-57.10900115966797, -79.3290023803711, -49.89699935913086),(-56.82600021362305, -79.02300262451172, -48.737998962402344),(-58.20600128173828, -78.81199645996094, -50.659000396728516),(-58.439998626708984, -79.2750015258789, -51.88999938964844)]
+>>> twwUC_U_list_atoms = ["P","OP1","OP2","O5'","C5'","C4'","O4'","C3'","O3'","C2'","O2'","C1'","N1","C2","O2","N3","C4","O4","C5","C6"]
+>>> twwUC_U_list_atoms_elements = ["P","O","O","O","C","C","O","C","O","C","O","C","N","C","O","N","C","O","C","C"]
+>>> twwUC_U_list_positions = [(52.62300109863281, 163.0240020751953, 82.34400177001953),(51.7239990234375, 162.08799743652344, 81.61499786376953),(53.42900085449219, 164.01199340820312, 81.5770034790039),(51.78300094604492, 163.79200744628906, 83.45999908447266),(50.75199890136719, 163.1179962158203, 84.20899963378906),(50.0260009765625, 164.10000610351562, 85.0989990234375),(50.957000732421875, 164.63400268554688, 86.0719985961914),(49.45600128173828, 165.31300354003906, 84.38600158691406),(48.15599822998047, 165.03599548339844, 83.86499786376953),(49.45199966430664, 166.3730010986328, 85.48100280761719),(48.358001708984375, 166.2570037841797, 86.36199951171875),(50.731998443603516, 166.02000427246094, 86.24299621582031),(51.94499969482422, 166.73500061035156, 85.80999755859375),(51.97999954223633, 168.10499572753906, 85.95700073242188),(51.03900146484375, 168.75, 86.38800048828125),(53.1619987487793, 168.69700622558594, 85.58200073242188),(54.284000396728516, 168.07400512695312, 85.08399963378906),(55.29899978637695, 168.73300170898438, 84.87999725341797),(54.16299819946289, 166.66700744628906, 84.93399810791016),(53.02899932861328, 166.0590057373047, 85.28900146484375)]
+>>> twwUC_C_list_atoms = ["P","OP1","OP2","O5'","C5'","C4'","O4'","C3'","O3'","C2'","O2'","C1'","N1","C2","O2","N3","C4","N4","C5","C6"]
+>>> twwUC_C_list_atoms_elements = ["P","O","O","O","C","C","O","C","O","C","O","C","N","C","O","N","C","N","C","C"]
+>>> twwUC_C_list_positions = [(51.45399856567383, 178.2270050048828, 85.64199829101562),(51.52799987792969, 179.64500427246094, 85.23200225830078),(50.196998596191406, 177.45799255371094, 85.41799926757812),(52.66400146484375, 177.4499969482422, 84.95800018310547),(54.012001037597656, 177.9250030517578, 85.10099792480469),(54.97700119018555, 176.7689971923828, 85.08399963378906),(54.652000427246094, 175.8520050048828, 86.15599822998047),(54.93899917602539, 175.89599609375, 83.8489990234375),(55.71099853515625, 176.46200561523438, 82.81400299072266),(55.54800033569336, 174.60000610351562, 84.3550033569336),(56.957000732421875, 174.67799377441406, 84.4469985961914),(54.93899917602539, 174.52000427246094, 85.75399780273438),(53.6879997253418, 173.73199462890625, 85.80500030517578),(53.76300048828125, 172.32899475097656, 85.83999633789062),(54.87699890136719, 171.78500366210938, 85.81199645996094),(52.61399841308594, 171.60699462890625, 85.89900207519531),(51.43299865722656, 172.22999572753906, 85.91899871826172),(50.332000732421875, 171.48800659179688, 85.98400115966797),(51.32899856567383, 173.6479949951172, 85.87699890136719),(52.46799850463867, 174.35299682617188, 85.82099914550781)]
+>>> cww_AU_A = rawdata_format_nucleotide('A', cwwAU_A_list_atoms, cwwAU_A_list_atoms_elements, cwwAU_A_list_positions)
+>>> cww_AU_U = rawdata_format_nucleotide('U', cwwAU_U_list_atoms, cwwAU_U_list_atoms_elements, cwwAU_U_list_positions)
+>>> tww_UC_U = rawdata_format_nucleotide('U', twwUC_U_list_atoms, twwUC_U_list_atoms_elements, twwUC_U_list_positions)
+>>> tww_UC_C = rawdata_format_nucleotide('C', twwUC_C_list_atoms, twwUC_C_list_atoms_elements, twwUC_C_list_positions)
+>>> >>> from RNADiscrepancy.utilities import rawdata_format_nucleotide
+>>> from RNADiscrepancy.measures import isodiscrepancy
+>>> cwwAU_A_list_atoms = ["P","OP1","OP2","O5'","C5'","C4'","O4'","C3'","O3'","C2'","O2'","C1'","N9","C8","N7","C5","C6","N6","N1","C2","N3","C4"]
+>>> cwwAU_A_list_atoms_elements = ["P","O","O","O","C","C","O","C","O","C","O","C","N","C","N","C","C","N","N","C","N","C"]
+>>> cwwAU_A_list_positions = [(-46.6870002746582, -80.67400360107422, -42.457000732421875),(-45.249000549316406, -80.74099731445312, -42.09299850463867),(-47.303001403808594, -79.36100006103516, -42.77299880981445),(-46.92900085449219, -81.6449966430664, -43.694000244140625),(-47.27000045776367, -83.01000213623047, -43.49399948120117),(-47.70500183105469, -83.62999725341797, -44.79600143432617),(-49.13199996948242, -83.92400360107422, -44.75299835205078),(-47.60100173950195, -82.73999786376953, -46.02299880981445),(-46.290000915527344, -82.58599853515625, -46.52899932861328),(-48.512001037597656, -83.47599792480469, -46.97999954223633),(-47.96799850463867, -84.72100067138672, -47.375),(-49.6870002746582, -83.72000122070312, -46.04600143432617),(-50.590999603271484, -82.5719985961914, -46.012001037597656),(-50.654998779296875, -81.52100372314453, -45.125),(-51.59299850463867, -80.6449966430664, -45.40800094604492),(-52.19300079345703, -81.16400146484375, -46.55400085449219),(-53.26100158691406, -80.71800231933594, -47.35499954223633),(-53.959999084472656, -79.60900115966797, -47.10100173950195),(-53.59199905395508, -81.46099853515625, -48.4379997253418),(-52.902000427246094, -82.58399963378906, -48.68299865722656),(-51.887001037597656, -83.11399841308594, -48.000999450683594),(-51.58000183105469, -82.34600067138672, -46.9370002746582)]
+>>> cwwAU_U_list_atoms = ["P","OP1","OP2","O5'","C5'","C4'","O4'","C3'","O3'","C2'","O2'","C1'","N1","C2","O2","N3","C4","O4","C5","C6"]
+>>> cwwAU_U_list_atoms_elements = ["P","O","O","O","C","C","O","C","O","C","O","C","N","C","O","N","C","O","C","C"]
+>>> cwwAU_U_list_positions = [(-62.04399871826172, -78.197998046875, -52.75299835205078),(-63.35599899291992, -78.33399963378906, -53.4370002746582),(-61.986000061035156, -77.62999725341797, -51.382999420166016),(-61.347999572753906, -79.62899780273438, -52.72700119018555),(-61.56800079345703, -80.56199645996094, -53.777000427246094),(-60.268001556396484, -81.21900177001953, -54.17399978637695),(-59.19599914550781, -80.2249984741211, -54.268001556396484),(-59.689998626708984, -82.20500183105469, -53.178001403808594),(-60.340999603271484, -83.45899963378906, -53.125),(-58.26900100708008, -82.29499816894531, -53.6879997253418),(-58.19900131225586, -82.8759994506836, -54.974998474121094),(-57.983001708984375, -80.80400085449219, -53.79199981689453),(-57.66400146484375, -80.25, -52.470001220703125),(-56.564998626708984, -80.7750015258789, -51.81100082397461),(-55.832000732421875, -81.61399841308594, -52.30500030517578),(-56.35599899291992, -80.28099822998047, -50.54899978637695),(-57.10900115966797, -79.3290023803711, -49.89699935913086),(-56.82600021362305, -79.02300262451172, -48.737998962402344),(-58.20600128173828, -78.81199645996094, -50.659000396728516),(-58.439998626708984, -79.2750015258789, -51.88999938964844)]
+>>> twwUC_U_list_atoms = ["P","OP1","OP2","O5'","C5'","C4'","O4'","C3'","O3'","C2'","O2'","C1'","N1","C2","O2","N3","C4","O4","C5","C6"]
+>>> twwUC_U_list_atoms_elements = ["P","O","O","O","C","C","O","C","O","C","O","C","N","C","O","N","C","O","C","C"]
+>>> twwUC_U_list_positions = [(52.62300109863281, 163.0240020751953, 82.34400177001953),(51.7239990234375, 162.08799743652344, 81.61499786376953),(53.42900085449219, 164.01199340820312, 81.5770034790039),(51.78300094604492, 163.79200744628906, 83.45999908447266),(50.75199890136719, 163.1179962158203, 84.20899963378906),(50.0260009765625, 164.10000610351562, 85.0989990234375),(50.957000732421875, 164.63400268554688, 86.0719985961914),(49.45600128173828, 165.31300354003906, 84.38600158691406),(48.15599822998047, 165.03599548339844, 83.86499786376953),(49.45199966430664, 166.3730010986328, 85.48100280761719),(48.358001708984375, 166.2570037841797, 86.36199951171875),(50.731998443603516, 166.02000427246094, 86.24299621582031),(51.94499969482422, 166.73500061035156, 85.80999755859375),(51.97999954223633, 168.10499572753906, 85.95700073242188),(51.03900146484375, 168.75, 86.38800048828125),(53.1619987487793, 168.69700622558594, 85.58200073242188),(54.284000396728516, 168.07400512695312, 85.08399963378906),(55.29899978637695, 168.73300170898438, 84.87999725341797),(54.16299819946289, 166.66700744628906, 84.93399810791016),(53.02899932861328, 166.0590057373047, 85.28900146484375)]
+>>> twwUC_C_list_atoms = ["P","OP1","OP2","O5'","C5'","C4'","O4'","C3'","O3'","C2'","O2'","C1'","N1","C2","O2","N3","C4","N4","C5","C6"]
+>>> twwUC_C_list_atoms_elements = ["P","O","O","O","C","C","O","C","O","C","O","C","N","C","O","N","C","N","C","C"]
+>>> twwUC_C_list_positions = [(51.45399856567383, 178.2270050048828, 85.64199829101562),(51.52799987792969, 179.64500427246094, 85.23200225830078),(50.196998596191406, 177.45799255371094, 85.41799926757812),(52.66400146484375, 177.4499969482422, 84.95800018310547),(54.012001037597656, 177.9250030517578, 85.10099792480469),(54.97700119018555, 176.7689971923828, 85.08399963378906),(54.652000427246094, 175.8520050048828, 86.15599822998047),(54.93899917602539, 175.89599609375, 83.8489990234375),(55.71099853515625, 176.46200561523438, 82.81400299072266),(55.54800033569336, 174.60000610351562, 84.3550033569336),(56.957000732421875, 174.67799377441406, 84.4469985961914),(54.93899917602539, 174.52000427246094, 85.75399780273438),(53.6879997253418, 173.73199462890625, 85.80500030517578),(53.76300048828125, 172.32899475097656, 85.83999633789062),(54.87699890136719, 171.78500366210938, 85.81199645996094),(52.61399841308594, 171.60699462890625, 85.89900207519531),(51.43299865722656, 172.22999572753906, 85.91899871826172),(50.332000732421875, 171.48800659179688, 85.98400115966797),(51.32899856567383, 173.6479949951172, 85.87699890136719),(52.46799850463867, 174.35299682617188, 85.82099914550781)]
+>>> cww_AU_A = rawdata_format_nucleotide('A', cwwAU_A_list_atoms, cwwAU_A_list_atoms_elements, cwwAU_A_list_positions)
+>>> cww_AU_U = rawdata_format_nucleotide('U', cwwAU_U_list_atoms, cwwAU_U_list_atoms_elements, cwwAU_U_list_positions)
+>>> tww_UC_U = rawdata_format_nucleotide('U', twwUC_U_list_atoms, twwUC_U_list_atoms_elements, twwUC_U_list_positions)
+>>> tww_UC_C = rawdata_format_nucleotide('C', twwUC_C_list_atoms, twwUC_C_list_atoms_elements, twwUC_C_list_positions)
+>>> cww_AU = (cww_AU_A, cww_AU_U)
+>>> tww_UC = (tww_UC_U, tww_UC_C)
+>>> isodiscrepancy(cww_AU, tww_UC)
+15.096495858783285
+```

rnadiscrepancy-0.0.2.dist-info/RECORD

@@ -0,0 +1,11 @@
+RNADiscrepancy/.test.py.swp,sha256=O0hTheO5Qz2CKDgdfs9ps0rNNW2MVBlJB2h1sP8_JUw,12288
+RNADiscrepancy/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+RNADiscrepancy/measures.py,sha256=7L5awwoMGEF1qaLo7aXWk9llWbQNfVza13vEVUtandQ,6569
+RNADiscrepancy/references.py,sha256=WX8psZMoWil7Hxphpu1lmO1u2k0Yh2BHc-WxeuAn5O4,2505
+RNADiscrepancy/rna.py,sha256=7KxwlxlE9DWsefZPLJEmk3fFLghBDcHKSyzU-iGpYys,1508
+RNADiscrepancy/test.py,sha256=wr0Oi21za2FUrSgtBbB6T31rv3iMp8CaJp5Wp9i-vkA,1891
+RNADiscrepancy/utilities.py,sha256=F6dzmQUGKz_725_pPD5VY8gtvv6dWh1-CX-CTHzWSYQ,342
+rnadiscrepancy-0.0.2.dist-info/METADATA,sha256=BFpbWOWhcnlbt4_wQS5JDvNGSmo4D8DECK01a7vC8DQ,13951
+rnadiscrepancy-0.0.2.dist-info/WHEEL,sha256=QccIxa26bgl1E6uMy58deGWi-0aeIkkangHcxk2kWfw,87
+rnadiscrepancy-0.0.2.dist-info/licenses/LICENSE,sha256=dP7Glyxd-54WGXrTuhOs_HelxCoX5jwPCluWFnvMHUE,749
+rnadiscrepancy-0.0.2.dist-info/RECORD,,

rnadiscrepancy-0.0.2.dist-info/WHEEL

@@ -0,0 +1,4 @@
+Wheel-Version: 1.0
+Generator: hatchling 1.29.0
+Root-Is-Purelib: true
+Tag: py3-none-any

rnadiscrepancy-0.0.2.dist-info/licenses/LICENSE

@@ -0,0 +1,17 @@
+	RNADiscrepancy - Computing RNA discrepancy values
+
+    Copyright (C) 2026  Vladimir Reinharz
+
+    This program is free software: you can redistribute it and/or modify
+    it under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    This program is distributed in the hope that it will be useful,
+    but WITHOUT ANY WARRANTY; without even the implied warranty of
+    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+    GNU General Public License for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with this program.  If not, see <https://www.gnu.org/licenses/>.
+