B ‰°bn ã@s:dZddlmZddlmZddlmZGdd„dƒZdS)zHMap residues of two structures to each other based on a FASTA alignment.é)ÚSCOPData)Ú Selection)Úis_aac@s2eZdZdZd dd„Zdd„Zdd „Zd d „Zd S)ÚStructureAlignmentzGClass to align two structures based on an alignment of their sequences.récCs| ¡}t |d¡}t |d¡}d} d} i} i} g} xÊt|ƒD]¾}|dd…|f}||}||}|dkr˜x|| }| d} t|ƒrlPqlW| ||¡nd}|dkrÒx|| }| d} t|ƒr¦Pq¦W| ||¡nd}|râ|| |<|rî|| |<|  ||f¡q>W| |_| |_| |_ dS)zÚInitialize. Attributes: - fasta_align - Alignment object - m1, m2 - two models - si, sj - the sequences in the Alignment object that correspond to the structures ÚRrNú-r) Zget_alignment_lengthrZunfold_entitiesÚrangerÚ_test_equivalenceÚappendÚmap12Úmap21Úduos)ÚselfZ fasta_alignZm1Zm2ZsiZsjÚlengthZrl1Zrl2Zp1Zp2r r rÚiÚcolumnÚaa1Zaa2Úr1Zr2©rú9lib/python3.7/site-packages/Bio/PDB/StructureAlignment.pyÚ__init__sH   zStructureAlignment.__init__cCs"| ¡}tj|}||kst‚dS)z6Test if aa in sequence fits aa in structure (PRIVATE).N)Z get_resnamerZprotein_letters_3to1ÚAssertionError)rrrZresnamerrrr Ps z$StructureAlignment._test_equivalencecCs |j|jfS)z­Map residues between the structures. Return two dictionaries that map a residue in one structure to the equivealent residue in the other structure. )r r )rrrrÚget_mapsVszStructureAlignment.get_mapsccs*x$tdt|jƒƒD]}|j|VqWdS)z*Create an iterator over all residue pairs.rN)r Úlenr)rrrrrÚ get_iterator^szStructureAlignment.get_iteratorN)rr)Ú__name__Ú __module__Ú __qualname__Ú__doc__rr rrrrrrrs  <rN)rZBio.DatarZBio.PDBrZBio.PDB.PolypeptiderrrrrrÚs