From 7f1ecf44728de038ca9d4f200fa881516e0e64b3 Mon Sep 17 00:00:00 2001 From: Thomas Forest Date: Wed, 18 Sep 2019 18:00:59 +0200 Subject: [PATCH] new ladders print --- src/__pycache__/atom.cpython-37.pyc | Bin 0 -> 11068 bytes src/__pycache__/geom.cpython-37.pyc | Bin 0 -> 2211 bytes src/__pycache__/pdb.cpython-37.pyc | Bin 0 -> 3248 bytes src/__pycache__/structure.cpython-37.pyc | Bin 0 -> 1436 bytes src/atom.py | 133 +++++------------------ src/dssp.py | 15 ++- src/structure.py | 8 +- 7 files changed, 41 insertions(+), 115 deletions(-) create mode 100644 src/__pycache__/atom.cpython-37.pyc create mode 100644 src/__pycache__/geom.cpython-37.pyc create mode 100644 src/__pycache__/pdb.cpython-37.pyc create mode 100644 src/__pycache__/structure.cpython-37.pyc diff --git a/src/__pycache__/atom.cpython-37.pyc b/src/__pycache__/atom.cpython-37.pyc new file mode 100644 index 0000000000000000000000000000000000000000..cf2b4a735c9fcaa901a814d851ca0bbe6d4e9616 GIT binary patch literal 11068 zcmds7TWlQHd7d*fd*wwGMN_0D%UjuT!_<{5+DQtRk*&K_)3hQfmbJ08UF{i>OD=cG zGeenTciY4?TDL#~wGX`*hARa{fabwZc_@k^E&9~7K+vb6m?sB)(rb~2LMY&V-+yL@ zyGtrH`q*8}Kj%JY{`>!*IXyB`bnv@%_5GWdE;`PCu+sf`D7=m<{18ApS~Z+9?`}g~ zb1pdA)BXpJ_P5=#i&{o!QOly{>734g;FNt`&_&!cdPooBp4B6I6!)Cor^j&5>-~Bh z_kup4CvY$7r}QN5L;9ee!hKjz>!)!a(TDV5+()&#>dZXz2)l6>W>l=cQe(0Eb5VF5 zS9lm8aLP(MWmhZhVyb61i!EiaL^k$MMXl8uXYj{+2G*Y^%1h8Z@Ku_?l@g2T`Tg+Az`a(A17|*0HsF zaZ%U9sKQ#v*9fx!j`GzJwY85gdTgh^mS-P+ymPbi&hsd=@x~RW?EP6{(@0h^AUmZ0cMJEjvK7hbzG1_v$`6@ zc@u;cTdFNp>rHkS*%k)tmFDeuwB8H@6V+QyytWPu-^}C)xc@jSs10ZGrogsQf_(&I z1p5ib2@U|nLo&GDtVTg7h_+yeefdhI-mFKJ3Txp6fa7M>&cCSm-x?8YcE>J&(3r^W zS1NI_Qdw>3jRwoZmCE`?wPAbAFjg~@1SD?{;KK~1j8{lB<=uQ%T571IKho)ST%iVd zGG|i&A#7IovWf8J65-1y!dFOyub2qmP$GQ8iSUi+gP`A3oVyZ)b-fYPxHtZ|Jzi&U zh3^2=$q3Hod|RD$Zl3LUWSV8C?JcW0LNXHCASPrV=Rt<$hh(xWN@SHB?s?qv!gPi8 z;p@;Yc(_^sw}?oDmYSqtoB>$-dM%jo zg+IzgGG16|Kw8YRT+?U6AVEl}XwpQ0YFOp*9~HPYmSE;smk&}sLiCV#1VwjH zp!T8Wp>iZ$NE6b7v>+Ww1LA*iqT>?(RD4u=RCxLseMCR2OZr#zQT?2LUVllypnp|= z8DoiR_W>mKA*2xqn^`D(r{+%qT)6;nmDQ`4%iftASIgdo`EL=R?ZW(xvUl+k+AhtX zD|_?j&@s;;Z=SnW_AXyU+vO{7mA$Lyu9UrN7YVLlm1{Rpdi&ITtWL%1N~{)Qbv9OK zVs$Q7-;C9TSY3?OrC4FW8?kybR_9}NIaXI=buCucV|61|Z$n#WK#O2AA%&2!oS{s= zY;Pw9$?y#D_dKqU;^U|>{D-|GYSi7D>Jng(3sOQ*&Qf0eAu3uySwpEmbXIa^T)S6c zy?}DRM_rt;PT;hBq4 z7I{8`+tErEpE2KC$;CRTLRvdE#3Iq1mwi5B$msSvK{m6Bv?Q6B- znWe~wjl!sDXF0P?qKb9jnk$x)--RG5RU3^`v=o%04bu!ucT8)wR4s*TL9KqLUMm4j zRCCADDb{d#*!3g6!y+#HKKuLgGg$I}KS_<_j=5X=yO>VLr^qIs^P#f9k5EVv z?Y80mtYpk5Z8rk{Z!4{)9H8-&$bH~New2x_KpxEavKaH_{I;?!9M{{r?RLj%b+E5z z_X^wYwo4+n;c~6`=1+0=7eN}_GrU6)+B+3jZb6)b*gbRV(Rpm$yaW&@IRDPw*JlcX za&wY>bZy`}28$*DVn^-pA)K{rd3J_a3g9D ziTm)(Wr$7z9=4~8(=Z_IEGRqB2D9XKd~z1N_1vP{$**6gmLvy78EDHP%=h)Mz>#TZ z8~#H(PrrT6>FGz09hLOMxS<~^0;&VTjMaYDokzbzkK&>j?^_UK^oXiR)Tc;wsySnS zt!k=_ilA&ynBTyt*snuCV}F@2b33}X7Qh*6Ms@Qh*WeIynP7q7Er722wYuONY^L%R zwHs$oo4Q^Mpe`tKhN2zkZcAGjXKf)KgI;Jsb0bu;Es6A>tH6;ol+k$DhH6RofKeQx z?oAf0Vsv;1fMFK@MRgdkViR2zIxaX}q<6C`y7H<<5_^b-Yin9(^Cnc^>Zcl|Rd-+RDw{Bym2rT;U zVRusX+fW&oVvvTO%qti)D5T2SV$d_q6njm<9uzdxjm2KtCu-5vJ$+E~jD-AUcg4jT zK9UF?F0?L!0v~&DeBiQZA9&8#kU4Sd;!GG-kv#BVfXr3w=sfHe-4DaiD;}~uAlQj> z5UWaaV>R(Ug#86W<>5xFR&69qTOO&k5KY&TYMgCUHA3n^G#1I`w(Ctvh7D6c2+Vn*&d`?ubJr`n~}&1KFj+OM}m%F1?p(WQ!p1f zR@@GEtzPzOt#Bqc$b?i?<=mOKD%5S=!mPQ(sRIfmC(c`~^%L}^{zM*#9Hx+qpd-hy zuNS3Jt(D@FOqzNVdpdd_rH-EX59T9U1p3r;1jVI;Q^ z*@qYIn-6f$tSIw=&Vqf}rpym;&!O)JI^W4O5$AS%U1V=toV}SLYas;j<{WXu%+RAB z13Xs0R-ga;^UoivU$+(mexJF3j+xAk^B0eRk82fOUkbFTHg;&Vi|6A)ZK;kS_2`}u z0ZF*}VsnrbBP2yS@L%vW$r?K1;87Qkg4ZWqXk=Vjk8#h|Q@w;q=ceQ#pQPmZTT}*g zh-;>_+wq#G>B4pOZeyRbssn`eK|+X8 zMP~(lpJ9=d3^0B!$v{GY1Rw*;lLMwOnYl1CY7TM>R91$`8=G%q4#Bz!v?%7c2&g|z znc%kxXvqy<)GQP1h^&+(L-#Z~xUhu6FxNx)PvrOKNAf50`|@MBKV>+kMx61-+$5)q zIXUWy->)Ex;`_s!|E~xM3F>Swz z#_a3@UOM}B&D(3dZIkC)5qR=G@NE06o;~zoVxYZmJol5#mv{JIpUddLzJDJ+``*(v zOI)`XKYR4=IsZxXzVN5CJjrPLcM zNJ}?fzzUu(kpd4ZCeqRrl%?;AvP&%#9SB~cmVSl;NsnZlxA_lk4=Tr|-u{W{IV9d} zwlr%t!F+HU7LgI%JkN|KLeOsCUE9Z8zwM#*fknF8&LBij+a4+<3h3)dd#ukgFh#<{ zH$O>cQ@zN#R%8(LGaZ<=X!m0>LUP8F{K1 z@Dd{=7)I3K;ehk9$@?WyQa?!<7d7z|(<0cRv~Skfuui~;ZqL^n-}QR4@ACbIF`@f@ zh~GUEq26cSQ#pDS<9W>U0j0uu2nubh&W^EHn^Qa=yRUAJcKrKl6B$?3oSXmMVX6!5P z8L8o%&8~%UEX!I@$S*<8xcV<>3cyLfpp@B>(~NGhJDH!h_08p+oUvmMN{akeg^RMbW>`C7cxE^taTyiyj^eKtd!WoB0mnfIoKE!PPsu1YZQW$6kv4pb-q7OFk~n zQGRQ!)wPpNGAWnMI$K8=%7jZn5QWF)B;F5CqR<_B{{OZOvw_vkU4njp|FcAawA{C;!bWH_ReY~h742Gv!b6ESdXbm{0R@chDJF>3gwu=jcBWn zP=HNToA=w)v>!t6zAvRMgwdcGqcT9~q;1zs(1Gu5p?;QRMSm;iK-MAUb~KY0NeL3I zA?1)7oF8)4q>m-^Z5NcYJ^@)NLbvB(-wOYKWyK~2eiQgN6@pU)q=~5Qt$ZSwvkS-q z(%OhZjx%ef-i%!voJCwJISV`#VJP;*p{JQ2xLhptb-N8u+z@KMrtn{llJDl+29 z)??MM{R2WrtCNJF?ds3_PLsZDg^zd+bo8RXSw_$s3zW*+=1qcBo zq~d2za;lv*lbd&tD#{||#!>I!SP5nLX|8#fC1ea>q2T-F4mxwNN=y}@Z&Uh^;T-7m z(Kjc3v>xfa^*QkKBymA|vFHIK$a@OiDTRjK|08Bn;HDG;-ge*B^=3Mh9P%B-MblY=izoq4q1R;mCuix55jzVw7R0#8Bypj+` zzD+RSVU;QuqWE2w%LF72K0UCmvn1CmFUL7sc%@sgbTRK^Y@Ct4a?uW9mcLxEW#(Xi z3B8g*3B5igy-HZHpI+lQKH=v&JZ#}L<(dyMrtf1WiZWG<*bh2yJRD)9P6rh(7C#or znxkY?IqiabEvH>cnDJPZAi$E6r@4b0%9odyGk?I+DT01Ejt_{&s`S0XKk^ksuEA=I z4$L26=H6S~LeCSodW;d#6s(VIJ!g)|sn5=kXMJS zQcNGKlaJM$`9l<+h=SMAxhD$7EehmtCy~0J5(A=pfdR(Jhe0yjUGqIYVxQ#T;{wj_ zHfop`$QN>wa@Y;!!h*)h(E-Sdy2<&{8QdNbdg;Nz(lwOP~YrV6vw|GC883&Bk z7s8*xAso5$*Q_{zIOW10zy*FYYddMGO6+Lfyf+_j-h1=jZ;zYJ2En&k|MBIYE+Kzu zVev6ApP|dAAe?YICcSE;J;oWgej+`K+uXrubC=gJI^5%Rj4p5R6^u2$%9|KHzQ)%v z);WDn+8dY9NKQh1gZ^+cvCto*%SRvu<@63IaOWMeMP>{#GcqV-3$kToP{?+{3ahY( z*32m!=(kSfEQz`b zG?3jVSvX0eRLO2GvR9GnhuK7^t`uQ+7-dN(pO))qQK&M}KNVTh&oilgktVQo!?Xq! z9_B)Akp8qy4!Uu;14fYxTDi;Frgj{TY3{b>;dI6SLT9wl0h~K1^hH7|YqQc$0yQda z`K?fp=vir}nMkhO5;g2ASLUYbHEp5PZWH;g zP_U?=hzV0jvU57KFj9=D8->a_r=a#ZZINvAQ_pdz2gqlk-O?3)n! zNd)}j9sdhpl0S&EXH&>&=h>6vhL7)Vz=K>x+*g^O=AA%@VA?h-R~psm4XuMneHTTk zf8fs((i{;{E+PqxjME{O!YoYzh)|wqF@Uqfoi2T*H)jV^|8o(C{?RxLVjgwO$f&~i%?JD)!{X$YS+cw^{8wbhCTd&h?xiA<;LtC^Y zQXY=Fwjdu0yG0NESo9F|Qb5kFe@}0`_9UQ(o{Jv(dn3vAx+zL&9%tU0d32sy<18_n{N<2ewVO%gQ)oMIT2qE6^Za-%J|4 zP8h3ny(^QWQ2C9dC)Ii?8|z*1`C9ihD;;WXtT4+Vn=BB=DCr|5}V137P8cC;{^kUI3KHK?p_wg2SvwVsTYlQvc#;1GRJCoHx58Epd8Ak>$ zKcSFhCIV42WgcJ3oKC0$5vU-uC|)`)V?u^xMcU7?0cDKSH$WJ-c-biPDmMN`GxNMZ ze=`i@Wmm~$V`T+RN8z9F^^gsiKc@VUogs72%vb#6&xnh$$zJjO-e7FHrw1uZlE6#@Q|F z6_?&1cT^(IATR7?w!)4|En~tjzs4`OutS*7j`+|~GXq0wqTIEY*ttR9p`E;?|I4q? zPkO%me|k{l6YG!fAy@zFtK%;YOnlEc=WjBx)I!|kKU765pKG6ovja|6X@FR%>lF8_ zjlG>GUVeLPZ)5Mt=o93`s1B0x?U614O3YD-2<6J79KA{89+9sT zStIf;NX^ddKqVbN5?MiaH263&O#N}M(F~%P#6{wz;jtfoo>iha4W){b81N&)tkg&n zDWW(~VS4yo{rv66{m<%QC|s4efe>!gXW>;=OjRIN z`Z7|@ESD-7bu()N6Eic4RhBzwCxOcHK^G`2vRo94@I=xEN;-7pH$gmmcdRkXbtNsD zffN3f6E#}{gQRlxiYn}MeQk(m1CC>Wft-gVvxu4T=r&*!2r%>lib?tigz-Fn7JgM; z;fvfstKe7U^SaIRHw@s3S>?<4y~QnqdU*4Sv24y8j?tgHi7Q(B6|B%%MNPvB<%`0} zO3A0sPP5)$+^X*ZN*`};J^t9;{>A3*&WGUCxYSuw0HC#+=EC zJ2(>=%F;8LDQD~8y&@lIAnxzs0&pA=7#3DJapfpBm2_emJaaAFA7GUaEQUr3KwAXl6;>E%QPa{21H_dl5H~lxGl94q z=FS0Z?&5&^;yF8RKG)VC**N58>d3d#PSWn14~`pNA&%6zYWjI9M>m#B*dniRm5BXU zsY+xh*1Adhf{I|@{h|9}|9$s1du^QU_esqvO`3d zdcWp|Z4}}EQupGPzWtiis6)v%q8LlcHk1qDj{F9Z`$TF)uBuJ0Vpv;}5ApJ<2T>fQ zO}`!V;zl!y4>PmT5jc61KzMmvlfs7ba!nrviJ|u1GvU5B6RA)Ff_zlJ3Oy?gei?eD zF>d$-E0S<8GXV%$j(mnv@=<=9jM4S4GyK{>LIU|cB4l|je_ixhaiSM^{-5eBmKIdRm=68){@C*$C1Ox;K6F|j$T8d zt0JO}rlyPNZFB}~h+ESo$K3&KXk5}Q!bH!yanj%_nv-;bXwa0ni8Dzji2MP_{1RbT|-o1e)?obH7Sa S<~$tR^dm}M?U7?ZpYb2k)&)NR literal 0 HcmV?d00001 diff --git a/src/__pycache__/structure.cpython-37.pyc b/src/__pycache__/structure.cpython-37.pyc new file mode 100644 index 0000000000000000000000000000000000000000..ee79b270cecd5469277e267eddedf020446bea25 GIT binary patch literal 1436 zcmbW1L2uJA6vyqvNz-)Q#)NOz|1Z zyOHMs>0lO~vbRiyN?b7|E`wZ_zBBK*A$-} zXWa3p1>C*`?8$LZfXKIAQhdqd}G5n2Bll@qBVVuk@rLJGGNH^P;VE-Pp;(%rA5!i+O9Bi8aN< z_M09omM)R9tjk79#IQaKzS`U5bp2)mr7a9RXc;fCFQ4fI_qXkLScM9 zH_QOP(E+ij-y`8(-x`=WZB1Pu6dJ}2t6dTut98(DNoU*<3*C?a*RokIwBOW2%R+-> z%oqj~|4tWmjHp{oG9L5Q7<=F3VTLqga7G3z|6lqR&0?*`n*;{>)}_Zsrc`ITK-y>+ zH^o0A(HBqGg|9=h0f7jmSPRAX8R>|;zM=r39165>>4$F$n^8QzJ|0}8W`fWHVqXNeQlv^mcsI^@MN59&oO p&h)(M&KM`;Ng6D^k)h8OX}aWP_Yd$2*T)}e!s!XN#csTt{shF0{bc|E literal 0 HcmV?d00001 diff --git a/src/atom.py b/src/atom.py index d3c7388..c3db1b2 100644 --- a/src/atom.py +++ b/src/atom.py @@ -243,7 +243,6 @@ def get_bridges(residues): bridges = {} bridge = {} strongest_bridge = {} - selected_bridges = [] for i in range(1,len(residues)-4): E_min = 0 for j in range(i+2,len(residues)-1): @@ -293,47 +292,23 @@ def get_bridges(residues): if(bridge['res1']+bridge['res2']1): - bridges[strongest_bridge['i']] = (Bridge(strongest_bridge['btype'], - strongest_bridge['ipos'], - strongest_bridge['jpos'], - [strongest_bridge['i'], - strongest_bridge['j']], - selected_bridges[-2]['btype'], - selected_bridges[-2]['ipos'], - selected_bridges[-2]['jpos'], - [selected_bridges[-2]['i'], - selected_bridges[-2]['j']])) + + bridges[strongest_bridge['i']] = (Bridge(strongest_bridge['btype'], + strongest_bridge['ipos'], + strongest_bridge['jpos'], + [strongest_bridge['i'], + strongest_bridge['j']])) - bridges[strongest_bridge['j']] = (Bridge(strongest_bridge['btype'], - strongest_bridge['ipos'], - strongest_bridge['jpos'], - [strongest_bridge['i'], - strongest_bridge['j']], - selected_bridges[-2]['btype'], - selected_bridges[-2]['ipos'], - selected_bridges[-2]['jpos'], - [selected_bridges[-2]['i'], - selected_bridges[-2]['j']])) - else: - bridges[strongest_bridge['i']] = (Bridge(strongest_bridge['btype'], - strongest_bridge['ipos'], - strongest_bridge['jpos'], - [strongest_bridge['i'], - strongest_bridge['j']])) - - bridges[strongest_bridge['j']] = (Bridge(strongest_bridge['btype'], - strongest_bridge['ipos'], - strongest_bridge['jpos'], - [strongest_bridge['i'], - strongest_bridge['j']])) + bridges[strongest_bridge['j']] = (Bridge(strongest_bridge['btype'], + strongest_bridge['ipos'], + strongest_bridge['jpos'], + [strongest_bridge['i'], + strongest_bridge['j']])) if(len(bridges)>0): return(bridges) else: @@ -366,12 +341,21 @@ def connected_ladders(ladd_1, ladd_2): links = [] for bridge in ladd_1['bridges']: if bridge.res_partner in res_list(ladd_2): - # return([ladd_1['i'], ladd_1['j'], bridge.i, bridge.j, - # ladd_2['i'], ladd_2['j']]) return ladd_2 return False +def connected_ladders2(ladd_1, ladd_2): + links = [] + for bridge in ladd_1['bridges']: + if bridge.res_partner in res_list(ladd_2): + return([ladd_1['i'], ladd_1['j'], bridge.i, bridge.j, + ladd_2['i'], ladd_2['j']]) + #return ladd_2 + + return False + + def get_sheets(ladders): """ Bridges between ladders. @@ -413,80 +397,19 @@ def get_sheets2(ladders): sheets = {} for ladder in ladders: for ladd2 in ladders: - if connected_ladders(ladders[ladder], ladders[ladd2]): - pass - print("ladder",ladders[ladder]['i'], ladders[ladder]['j'],"bridge",bridge.i, bridge.j, - "ladder 2",ladders[ladd2]['i'], ladders[ladd2]['j']) + if connected_ladders(ladders[ladder], ladders[ladd2]): + bridge_i = connected_ladders2(ladders[ladder], ladders[ladd2])[2] + bridge_j = connected_ladders2(ladders[ladder], ladders[ladd2])[3] + print("ladder",ladders[ladder]['i'], ladders[ladder]['j'],"bridge",bridge_i, bridge_j, + "ladder 2",ladders[ladd2]['i'], ladders[ladd2]['j']) def res_list(ladder): # TODO : method in ladder class l=[] - for i in range(ladder['start'], ladder['end']): + for i in range(ladder['i'], ladder['j']): l.append(i) return(l) - -def get_bends(residues): - bends = {} - for i in range(2,len(residues)-2): - angle = math.degrees(vector_angles(vectors_substr(position_vector(residues[i].atoms["CA"].coords), - position_vector(residues[i-2].atoms["CA"].coords)), - vectors_substr(position_vector(residues[i+2].atoms["CA"].coords), - position_vector(residues[i].atoms["CA"].coords)))) - if(angle>70): - print("angle", residues[i].resid, angle) - bends[residues[i].resid] = angle - return(bends) - -def get_chirality(residues): - - for i in range(1,len(residues)-2): - chirality = {} - angle = calc_dihedral(residues[i-1].atoms["CA"].coords, - residues[i].atoms["CA"].coords, - residues[i+1].atoms["CA"].coords, - residues[i+2].atoms["CA"].coords) - - if(angle>0 and angle<=180): - sign="+" - print("chirality", residues[i].resid, "+", angle) - - if(angle<=0 and angle > -180): - sign="-" - print("chirality", residues[i].resid, "-", angle) - - chirality[residues[i].resid] = sign - - return chirality - - -def get_phi_psi(residues): - for i in range(len(residues)): - if(i==0): - phi = 360.0 - else: - phi = calc_dihedral(residues[i-1].atoms["C"].coords, - residues[i].atoms["N"].coords, - residues[i].atoms["CA"].coords, - residues[i].atoms["C"].coords) - if(i==len(residues)-1): - psi = 360.0 - else: - psi = calc_dihedral(residues[i].atoms["N"].coords, - residues[i].atoms["CA"].coords, - residues[i].atoms["C"].coords, - residues[i+1].atoms["N"].coords) - - print("ANGLES", i, phi, psi) - -def get_TCO(res1, res2): - CO_res1 = vector_from_pos(res1.atoms["C"].coords, - res1.atoms["O"].coords) - CO_res2 = vector_from_pos(res2.atoms["C"].coords, - res2.atoms["O"].coords) - angle = vector_angles(CO_res1, CO_res2) - return(math.cos(angle)) - def build_turns_patterns(residues): turns_3 = {} diff --git a/src/dssp.py b/src/dssp.py index 148f747..3f1c9a2 100644 --- a/src/dssp.py +++ b/src/dssp.py @@ -79,7 +79,6 @@ for i,res in enumerate(residues): turn_3 = print_turn_pattern(residues, res, turns[0]) turn_4 = print_turn_pattern(residues, res, turns[1]) turn_5 = print_turn_pattern(residues, res, turns[2]) - helix_3 = print_helix_pattern(residues, res, helix[0]) helix_4 = print_helix_pattern(residues, res, helix[1]) helix_5 = print_helix_pattern(residues, res, helix[2]) @@ -91,5 +90,15 @@ for i,res in enumerate(residues): round(alpha, 1), round(phi, 1), round(psi,1), round(x_ca,1), round(y_ca,1), round(z_ca,1)) -for bridge in bridges: - print(bridges[bridge].i, bridges[bridge].j) +for ladd1 in ladders: + ladd_1 = ladders[ladd1] + for ladd2 in ladders: + ladd_2 = ladders[ladd2] + for bridge in ladd_1['bridges']: + if bridge.j in res_list(ladd_2): + print("ladder", ladd_1['i'],"-",ladd_1['j'], "|", bridge.i, "...", bridge.j, "| ladder", ladd_2['i'], ladd_2['j']) + +# for ind, sheet in sheets.items(): +# print(ind, sheet) +# for ladder in sheet: +# print(ladder['start'], ladder['end']) diff --git a/src/structure.py b/src/structure.py index ab6d585..1c206ad 100644 --- a/src/structure.py +++ b/src/structure.py @@ -14,19 +14,13 @@ class Turn(Structure): class Bridge(Structure): - def __init__(self, bridge_type, res_num, res_partner, indices, - second_bridge_type=None, second_resnum=None, - second_res_partner=None, second_indices=None): + def __init__(self, bridge_type, res_num, res_partner, indices): self.bridge_type = bridge_type self.res_num = res_num self.res_partner = res_partner Structure.res = res_num self.i = indices[0] self.j = indices[1] - self.second_bridge_type = second_bridge_type - self.second_resnum = second_resnum - self.second_res_partner = second_res_partner - self.second_indices = second_indices class Helix(Structure):