Цель - знакомство с возможностями гомологичного моделирования комплекса белка с лигандом
In [2]:
import sys
import modeller
import _modeller
import modeller.automodel
import nglview
import ipywidgets
from IPython.display import Image
In [3]:
env=modeller.environ()
env.io.hetatm=True
MODELLER 9.19, 2017/07/19, r11078
PROTEIN STRUCTURE MODELLING BY SATISFACTION OF SPATIAL RESTRAINTS
Copyright(c) 1989-2017 Andrej Sali
All Rights Reserved
Written by A. Sali
with help from
B. Webb, M.S. Madhusudhan, M-Y. Shen, G.Q. Dong,
M.A. Marti-Renom, N. Eswar, F. Alber, M. Topf, B. Oliva,
A. Fiser, R. Sanchez, B. Yerkovich, A. Badretdinov,
F. Melo, J.P. Overington, E. Feyfant
University of California, San Francisco, USA
Rockefeller University, New York, USA
Harvard University, Cambridge, USA
Imperial Cancer Research Fund, London, UK
Birkbeck College, University of London, London, UK
Kind, OS, HostName, Kernel, Processor: 4, Linux shadbox 3.2.0-29-generic x86_64
Date and time of compilation : 2017/07/19 14:40:43
MODELLER executable type : x86_64-intel8
Job starting time (YY/MM/DD HH:MM:SS): 2017/12/29 09:01:33
! Образец (лизоцим форели)
wget http://www.pdb.org/pdb/files/1lmp.pdb
wget http://www.uniprot.org/uniprot/P11941.fasta
! Целевой комплекс (лизоцим курицы)
In [20]:
lmp = nglview.show_structure_file('1lmp.pdb')
lmp
In [23]:
env = modeller.environ()
env.io.hetatm = True
alignm = modeller.alignment(env)
alignm.append(file='P00698.fasta', align_codes='all', alignment_format='FASTA')
mdl = modeller.model(env, file='1lmp.pdb', model_segment=('FIRST:'+'A', 'LAST:'+'A'))
alignm.append_model(mdl, atom_files='1lmp.pdb', align_codes='1lmp')
alignm[0].code = 'bad_protein'
alignm.salign()
alignm.write(file='bad_protein.ali', alignment_format='PIR')
s = alignm[0]
pdb = alignm[1]
print s.code, pdb.code
a = modeller.automodel.automodel(env, alnfile='bad_protein.ali', knowns= pdb.code , sequence = s.code )
a.name='mod'+s.code
a.starting_model = 1
a.ending_model = 1
a.make()
read_pd_459W> Residue type NAG not recognized. 'automodel' model building
will treat this residue as a rigid body.
To use real parameters, add the residue type to ${LIB}/restyp.lib,
its topology to ${LIB}/top_*.lib, and suitable forcefield
parameters to ${LIB}/par.lib.
rdpdb___459W> Residue type NDG not recognized. 'automodel' model building
will treat this residue as a rigid body.
To use real parameters, add the residue type to ${LIB}/restyp.lib,
its topology to ${LIB}/top_*.lib, and suitable forcefield
parameters to ${LIB}/par.lib.
SALIGN_____> adding the next group to the alignment; iteration 1
bad_protein 1lmp
fndatmi_285W> Only 129 residues out of 132 contain atoms of type CA
(This is usually caused by non-standard residues, such
as ligands, or by PDB files with missing atoms.)
fndatmi_285W> Only 129 residues out of 132 contain atoms of type CA
(This is usually caused by non-standard residues, such
as ligands, or by PDB files with missing atoms.)
check_ali___> Checking the sequence-structure alignment.
Implied intrachain target CA(i)-CA(i+1) distances longer than 8.0 angstroms:
ALN_POS TMPL RID1 RID2 NAM1 NAM2 DIST
----------------------------------------------
END OF TABLE
read_to_681_> topology.submodel read from topology file: 3
fndatmi_285W> Only 129 residues out of 132 contain atoms of type CA
(This is usually caused by non-standard residues, such
as ligands, or by PDB files with missing atoms.)
patch_s_522_> Number of disulfides patched in MODEL: 4
The following 1 residues contain 6-membered rings with poor geometries
after transfer from templates. Rebuilding rings from internal coordinates:
<Residue 52 (type PHE)>
mdtrsr__446W> A potential that relies on one protein is used, yet you have at
least one known structure available. MDT, not library, potential is used.
iup2crm_280W> No topology library in memory or assigning a BLK residue.
Default CHARMM atom type assigned: C1 --> CT2
This message is written only for the first such atom.
0 atoms in HETATM/BLK residues constrained
to protein atoms within 2.30 angstroms
and protein CA atoms within 10.00 angstroms
0 atoms in residues without defined topology
constrained to be rigid bodies
condens_443_> Restraints marked for deletion were removed.
Total number of restraints before, now: 12801 11802
iupac_m_397W> Atoms were not swapped because of the uncertainty of how to handle the H atom.
>> ENERGY; Differences between the model's features and restraints:
Number of all residues in MODEL : 147
Number of all, selected real atoms : 1135 1135
Number of all, selected pseudo atoms : 0 0
Number of all static, selected restraints : 11802 11802
COVALENT_CYS : F
NONBONDED_SEL_ATOMS : 1
Number of non-bonded pairs (excluding 1-2,1-3,1-4): 2392
Dynamic pairs routine : 2, NATM x NATM cell sorting
Atomic shift for contacts update (UPDATE_DYNAMIC) : 0.390
LENNARD_JONES_SWITCH : 6.500 7.500
COULOMB_JONES_SWITCH : 6.500 7.500
RESIDUE_SPAN_RANGE : 0 99999
NLOGN_USE : 15
CONTACT_SHELL : 4.000
DYNAMIC_PAIRS,_SPHERE,_COULOMB,_LENNARD,_MODELLER : T T F F F
SPHERE_STDV : 0.050
RADII_FACTOR : 0.820
Current energy : 758.1279
Summary of the restraint violations:
NUM ... number of restraints.
NUMVI ... number of restraints with RVIOL > VIOL_REPORT_CUT[i].
RVIOL ... relative difference from the best value.
NUMVP ... number of restraints with -Ln(pdf) > VIOL_REPORT_CUT2[i].
RMS_1 ... RMS(feature, minimally_violated_basis_restraint, NUMB).
RMS_2 ... RMS(feature, best_value, NUMB).
MOL.PDF ... scaled contribution to -Ln(Molecular pdf).
# RESTRAINT_GROUP NUM NUMVI NUMVP RMS_1 RMS_2 MOL.PDF S_i
------------------------------------------------------------------------------------------------------
1 Bond length potential : 1157 0 0 0.006 0.006 11.178 1.000
2 Bond angle potential : 1566 0 4 2.055 2.055 133.89 1.000
3 Stereochemical cosine torsion poten: 738 0 27 48.906 48.906 272.75 1.000
4 Stereochemical improper torsion pot: 468 0 0 1.349 1.349 19.379 1.000
5 Soft-sphere overlap restraints : 2392 0 0 0.001 0.001 0.40765 1.000
6 Lennard-Jones 6-12 potential : 0 0 0 0.000 0.000 0.0000 1.000
7 Coulomb point-point electrostatic p: 0 0 0 0.000 0.000 0.0000 1.000
8 H-bonding potential : 0 0 0 0.000 0.000 0.0000 1.000
9 Distance restraints 1 (CA-CA) : 2404 0 0 0.106 0.106 31.526 1.000
10 Distance restraints 2 (N-O) : 2562 0 0 0.146 0.146 64.685 1.000
11 Mainchain Phi dihedral restraints : 0 0 0 0.000 0.000 0.0000 1.000
12 Mainchain Psi dihedral restraints : 0 0 0 0.000 0.000 0.0000 1.000
13 Mainchain Omega dihedral restraints: 146 0 2 4.280 4.280 31.540 1.000
14 Sidechain Chi_1 dihedral restraints: 120 0 0 70.984 70.984 18.141 1.000
15 Sidechain Chi_2 dihedral restraints: 86 0 0 79.560 79.560 43.141 1.000
16 Sidechain Chi_3 dihedral restraints: 29 0 0 85.395 85.395 21.413 1.000
17 Sidechain Chi_4 dihedral restraints: 18 0 0 92.041 92.041 12.366 1.000
18 Disulfide distance restraints : 4 0 0 0.008 0.008 0.39303E-01 1.000
19 Disulfide angle restraints : 8 0 0 1.786 1.786 0.56371 1.000
20 Disulfide dihedral angle restraints: 4 0 0 28.199 28.199 3.1127 1.000
21 Lower bound distance restraints : 0 0 0 0.000 0.000 0.0000 1.000
22 Upper bound distance restraints : 0 0 0 0.000 0.000 0.0000 1.000
23 Distance restraints 3 (SDCH-MNCH) : 1609 0 0 0.342 0.342 35.233 1.000
24 Sidechain Chi_5 dihedral restraints: 0 0 0 0.000 0.000 0.0000 1.000
25 Phi/Psi pair of dihedral restraints: 145 16 11 24.118 58.325 22.832 1.000
26 Distance restraints 4 (SDCH-SDCH) : 738 0 0 0.572 0.572 35.929 1.000
27 Distance restraints 5 (X-Y) : 0 0 0 0.000 0.000 0.0000 1.000
28 NMR distance restraints 6 (X-Y) : 0 0 0 0.000 0.000 0.0000 1.000
29 NMR distance restraints 7 (X-Y) : 0 0 0 0.000 0.000 0.0000 1.000
30 Minimal distance restraints : 0 0 0 0.000 0.000 0.0000 1.000
31 Non-bonded restraints : 0 0 0 0.000 0.000 0.0000 1.000
32 Atomic accessibility restraints : 0 0 0 0.000 0.000 0.0000 1.000
33 Atomic density restraints : 0 0 0 0.000 0.000 0.0000 1.000
34 Absolute position restraints : 0 0 0 0.000 0.000 0.0000 1.000
35 Dihedral angle difference restraint: 0 0 0 0.000 0.000 0.0000 1.000
36 GBSA implicit solvent potential : 0 0 0 0.000 0.000 0.0000 1.000
37 EM density fitting potential : 0 0 0 0.000 0.000 0.0000 1.000
38 SAXS restraints : 0 0 0 0.000 0.000 0.0000 1.000
39 Symmetry restraints : 0 0 0 0.000 0.000 0.0000 1.000
# Heavy relative violation of each residue is written to: bad_protein.V99990001
# The profile is NOT normalized by the number of restraints.
# The profiles are smoothed over a window of residues: 1
# The sum of all numbers in the file: 14786.3799
List of the violated restraints:
A restraint is violated when the relative difference
from the best value (RVIOL) is larger than CUTOFF.
ICSR ... index of a restraint in the current set.
RESNO ... residue numbers of the first two atoms.
ATM ... IUPAC atom names of the first two atoms.
FEAT ... the value of the feature in the model.
restr ... the mean of the basis restraint with the smallest
difference from the model (local minimum).
viol ... difference from the local minimum.
rviol ... relative difference from the local minimum.
RESTR ... the best value (global minimum).
VIOL ... difference from the best value.
RVIOL ... relative difference from the best value.
-------------------------------------------------------------------------------------------------
Feature 25 : Phi/Psi pair of dihedral restraints
List of the RVIOL violations larger than : 6.5000
# ICSR RESNO1/2 ATM1/2 INDATM1/2 FEAT restr viol rviol RESTR VIOL RVIOL
1 3946 1M 2R C N 7 9 -84.53 -72.10 19.29 1.25 -63.00 163.67 23.59
1 2R 2R N CA 9 10 156.66 141.90 -41.10
2 3947 2R 3S C N 18 20 -80.81 -72.40 27.91 1.38 -64.10 146.93 11.72
2 3S 3S N CA 20 21 179.02 152.40 -35.00
3 3948 3S 4L C N 24 26 -84.60 -108.50 29.44 1.30 -63.50 157.94 20.81
3 4L 4L N CA 26 27 115.32 132.50 -41.20
4 3949 4L 5L C N 32 34 -83.67 -70.70 24.41 1.51 -63.50 157.82 23.12
4 5L 5L N CA 34 35 162.28 141.60 -41.20
5 3950 5L 6I C N 40 42 74.01 -120.60 165.39 11.54 -120.60 165.39 11.54
5 6I 6I N CA 42 43 128.51 130.30 130.30
6 3951 6I 7L C N 48 50 -119.98 -108.50 12.12 0.63 -63.50 178.95 22.33
6 7L 7L N CA 50 51 128.61 132.50 -41.20
7 3953 8V 9L C N 63 65 -93.07 -108.50 23.57 1.08 -63.50 158.67 20.53
7 9L 9L N CA 65 66 114.69 132.50 -41.20
8 3955 10C 11F C N 77 79 -75.26 -71.40 15.65 1.35 -63.20 170.26 22.84
8 11F 11F N CA 79 80 125.53 140.70 -44.30
9 3956 11F 12L C N 88 90 -122.58 -108.50 17.17 0.94 -63.50 174.20 21.63
9 12L 12L N CA 90 91 122.68 132.50 -41.20
10 3958 13P 14L C N 103 105 -82.71 -108.50 47.79 2.25 -63.50 134.85 17.71
10 14L 14L N CA 105 106 92.27 132.50 -41.20
11 3959 14L 15A C N 111 113 -147.79 -134.00 14.64 0.68 -62.50 -163.50 36.77
11 15A 15A N CA 113 114 142.08 147.00 -40.90
12 3960 15A 16A C N 116 118 161.70 -134.00 65.38 1.75 -62.50 -149.92 41.47
12 16A 16A N CA 118 119 158.82 147.00 -40.90
13 3961 16A 17L C N 121 123 -105.80 -108.50 36.32 1.93 -63.50 143.85 18.04
13 17L 17L N CA 123 124 96.29 132.50 -41.20
14 3962 17L 18G C N 129 131 -79.63 -80.20 18.87 0.71 82.20 -141.55 7.23
14 18G 18G N CA 131 132 155.24 174.10 8.50
15 3963 18G 19K C N 133 135 -115.88 -118.00 8.60 0.38 -62.90 179.56 21.05
15 19K 19K N CA 135 136 130.76 139.10 -40.80
16 3998 53E 54S C N 410 412 -134.73 -64.10 78.92 8.20 -64.10 78.92 8.20
16 54S 54S N CA 412 413 0.21 -35.00 -35.00
report______> Distribution of short non-bonded contacts:
DISTANCE1: 0.00 2.10 2.20 2.30 2.40 2.50 2.60 2.70 2.80 2.90 3.00 3.10 3.20 3.30 3.40
DISTANCE2: 2.10 2.20 2.30 2.40 2.50 2.60 2.70 2.80 2.90 3.00 3.10 3.20 3.30 3.40 3.50
FREQUENCY: 0 0 0 0 0 6 6 40 80 109 108 147 158 194 181
<< end of ENERGY.
>> Summary of successfully produced models:
Filename molpdf
----------------------------------------
bad_protein.B99990001.pdb 758.12793
In [24]:
bad_protein = nglview.show_structure_file('bad_protein.B99990001.pdb')
bad_protein
In [25]:
Image('bad_protein.png')
Out[25]:
In [26]:
## Получаем список остаков
l1 = alignm[1].residues[129:132]
l1
Out[26]:
[Residue 130:A (type NAG), Residue 131:A (type NAG), Residue 132:A (type NDG)]
In [28]:
# Припишем точки, как обозначение лиганда
stri = ''
for i in alignm[0].residues:
stri += i.code
stri += '...'
# Добавляем последовательность из строки
alignm.append_sequence(stri)
## Выбираем объект для моделирования
s = alignm[2]
pdb = alignm[1]
s.code = 'good'
alignm.salign()
alignm.write(file='alignment.ali', alignment_format='PIR')
# создадим другое окружение, чтобы точно все работало
env1=modeller.environ()
env1.io.atom_files_directory = ['.', '../atom_files']
env1.io.hetatm=True
## Создаем новый объект automodel
a = modeller.automodel.automodel(env1, alnfile='alignment.ali', knowns= pdb.code , sequence = s.code )
a.name='mod'+s.code
a.starting_model = 1
a.ending_model = 1
a.make()
SALIGN_____> adding the next group to the alignment; iteration 1
SALIGN_____> adding the next group to the alignment; iteration 2
SALIGN_____> adding the next group to the alignment; iteration 3
read_pd_459W> Residue type NAG not recognized. 'automodel' model building
will treat this residue as a rigid body.
To use real parameters, add the residue type to ${LIB}/restyp.lib,
its topology to ${LIB}/top_*.lib, and suitable forcefield
parameters to ${LIB}/par.lib.
rdpdb___459W> Residue type NDG not recognized. 'automodel' model building
will treat this residue as a rigid body.
To use real parameters, add the residue type to ${LIB}/restyp.lib,
its topology to ${LIB}/top_*.lib, and suitable forcefield
parameters to ${LIB}/par.lib.
fndatmi_285W> Only 129 residues out of 132 contain atoms of type CA
(This is usually caused by non-standard residues, such
as ligands, or by PDB files with missing atoms.)
fndatmi_285W> Only 129 residues out of 132 contain atoms of type CA
(This is usually caused by non-standard residues, such
as ligands, or by PDB files with missing atoms.)
check_ali___> Checking the sequence-structure alignment.
Implied intrachain target CA(i)-CA(i+1) distances longer than 8.0 angstroms:
ALN_POS TMPL RID1 RID2 NAM1 NAM2 DIST
----------------------------------------------
END OF TABLE
read_to_681_> topology.submodel read from topology file: 3
getf_______W> RTF restraint not found in the atoms list:
residue type, indices: 10 147
atom names : C +N
atom indices : 1133 0
getf_______W> RTF restraint not found in the atoms list:
residue type, indices: 10 147
atom names : C CA +N O
atom indices : 1133 1128 0 1134
fndatmi_285W> Only 129 residues out of 132 contain atoms of type CA
(This is usually caused by non-standard residues, such
as ligands, or by PDB files with missing atoms.)
patch_s_522_> Number of disulfides patched in MODEL: 4
The following 1 residues contain 6-membered rings with poor geometries
after transfer from templates. Rebuilding rings from internal coordinates:
<Residue 52 (type PHE)>
mdtrsr__446W> A potential that relies on one protein is used, yet you have at
least one known structure available. MDT, not library, potential is used.
iup2crm_280W> No topology library in memory or assigning a BLK residue.
Default CHARMM atom type assigned: C1 --> CT2
This message is written only for the first such atom.
43 atoms in HETATM/BLK residues constrained
to protein atoms within 2.30 angstroms
and protein CA atoms within 10.00 angstroms
43 atoms in residues without defined topology
constrained to be rigid bodies
condens_443_> Restraints marked for deletion were removed.
Total number of restraints before, now: 14202 13203
iupac_m_397W> Atoms were not swapped because of the uncertainty of how to handle the H atom.
>> ENERGY; Differences between the model's features and restraints:
Number of all residues in MODEL : 150
Number of all, selected real atoms : 1178 1178
Number of all, selected pseudo atoms : 0 0
Number of all static, selected restraints : 13203 13203
COVALENT_CYS : F
NONBONDED_SEL_ATOMS : 1
Number of non-bonded pairs (excluding 1-2,1-3,1-4): 2515
Dynamic pairs routine : 2, NATM x NATM cell sorting
Atomic shift for contacts update (UPDATE_DYNAMIC) : 0.390
LENNARD_JONES_SWITCH : 6.500 7.500
COULOMB_JONES_SWITCH : 6.500 7.500
RESIDUE_SPAN_RANGE : 0 99999
NLOGN_USE : 15
CONTACT_SHELL : 4.000
DYNAMIC_PAIRS,_SPHERE,_COULOMB,_LENNARD,_MODELLER : T T F F F
SPHERE_STDV : 0.050
RADII_FACTOR : 0.820
Current energy : 875.7300
Summary of the restraint violations:
NUM ... number of restraints.
NUMVI ... number of restraints with RVIOL > VIOL_REPORT_CUT[i].
RVIOL ... relative difference from the best value.
NUMVP ... number of restraints with -Ln(pdf) > VIOL_REPORT_CUT2[i].
RMS_1 ... RMS(feature, minimally_violated_basis_restraint, NUMB).
RMS_2 ... RMS(feature, best_value, NUMB).
MOL.PDF ... scaled contribution to -Ln(Molecular pdf).
# RESTRAINT_GROUP NUM NUMVI NUMVP RMS_1 RMS_2 MOL.PDF S_i
------------------------------------------------------------------------------------------------------
1 Bond length potential : 1157 0 0 0.006 0.006 12.507 1.000
2 Bond angle potential : 1566 0 5 2.164 2.164 147.37 1.000
3 Stereochemical cosine torsion poten: 738 0 27 48.745 48.745 273.75 1.000
4 Stereochemical improper torsion pot: 468 0 0 1.459 1.459 22.626 1.000
5 Soft-sphere overlap restraints : 2515 2 2 0.008 0.008 17.511 1.000
6 Lennard-Jones 6-12 potential : 0 0 0 0.000 0.000 0.0000 1.000
7 Coulomb point-point electrostatic p: 0 0 0 0.000 0.000 0.0000 1.000
8 H-bonding potential : 0 0 0 0.000 0.000 0.0000 1.000
9 Distance restraints 1 (CA-CA) : 2404 0 0 0.099 0.099 29.696 1.000
10 Distance restraints 2 (N-O) : 2562 0 1 0.151 0.151 73.454 1.000
11 Mainchain Phi dihedral restraints : 0 0 0 0.000 0.000 0.0000 1.000
12 Mainchain Psi dihedral restraints : 0 0 0 0.000 0.000 0.0000 1.000
13 Mainchain Omega dihedral restraints: 146 0 3 4.568 4.568 35.931 1.000
14 Sidechain Chi_1 dihedral restraints: 120 0 1 77.686 77.686 33.797 1.000
15 Sidechain Chi_2 dihedral restraints: 86 0 1 77.941 77.941 39.949 1.000
16 Sidechain Chi_3 dihedral restraints: 29 0 0 86.812 86.812 16.266 1.000
17 Sidechain Chi_4 dihedral restraints: 18 0 0 100.544 100.544 12.622 1.000
18 Disulfide distance restraints : 4 0 0 0.007 0.007 0.32430E-01 1.000
19 Disulfide angle restraints : 8 0 0 1.713 1.713 0.51852 1.000
20 Disulfide dihedral angle restraints: 4 0 0 28.841 28.841 3.0360 1.000
21 Lower bound distance restraints : 0 0 0 0.000 0.000 0.0000 1.000
22 Upper bound distance restraints : 0 0 0 0.000 0.000 0.0000 1.000
23 Distance restraints 3 (SDCH-MNCH) : 1609 0 0 0.393 0.393 45.453 1.000
24 Sidechain Chi_5 dihedral restraints: 0 0 0 0.000 0.000 0.0000 1.000
25 Phi/Psi pair of dihedral restraints: 145 16 10 20.756 58.719 34.771 1.000
26 Distance restraints 4 (SDCH-SDCH) : 738 0 1 0.779 0.779 64.998 1.000
27 Distance restraints 5 (X-Y) : 1401 0 0 0.029 0.029 11.440 1.000
28 NMR distance restraints 6 (X-Y) : 0 0 0 0.000 0.000 0.0000 1.000
29 NMR distance restraints 7 (X-Y) : 0 0 0 0.000 0.000 0.0000 1.000
30 Minimal distance restraints : 0 0 0 0.000 0.000 0.0000 1.000
31 Non-bonded restraints : 0 0 0 0.000 0.000 0.0000 1.000
32 Atomic accessibility restraints : 0 0 0 0.000 0.000 0.0000 1.000
33 Atomic density restraints : 0 0 0 0.000 0.000 0.0000 1.000
34 Absolute position restraints : 0 0 0 0.000 0.000 0.0000 1.000
35 Dihedral angle difference restraint: 0 0 0 0.000 0.000 0.0000 1.000
36 GBSA implicit solvent potential : 0 0 0 0.000 0.000 0.0000 1.000
37 EM density fitting potential : 0 0 0 0.000 0.000 0.0000 1.000
38 SAXS restraints : 0 0 0 0.000 0.000 0.0000 1.000
39 Symmetry restraints : 0 0 0 0.000 0.000 0.0000 1.000
# Heavy relative violation of each residue is written to: good.V99990001
# The profile is NOT normalized by the number of restraints.
# The profiles are smoothed over a window of residues: 1
# The sum of all numbers in the file: 16185.3486
List of the violated restraints:
A restraint is violated when the relative difference
from the best value (RVIOL) is larger than CUTOFF.
ICSR ... index of a restraint in the current set.
RESNO ... residue numbers of the first two atoms.
ATM ... IUPAC atom names of the first two atoms.
FEAT ... the value of the feature in the model.
restr ... the mean of the basis restraint with the smallest
difference from the model (local minimum).
viol ... difference from the local minimum.
rviol ... relative difference from the local minimum.
RESTR ... the best value (global minimum).
VIOL ... difference from the best value.
RVIOL ... relative difference from the best value.
-------------------------------------------------------------------------------------------------
Feature 25 : Phi/Psi pair of dihedral restraints
List of the RVIOL violations larger than : 6.5000
# ICSR RESNO1/2 ATM1/2 INDATM1/2 FEAT restr viol rviol RESTR VIOL RVIOL
1 3946 1M 2R C N 7 9 -73.10 -72.10 6.20 0.45 -63.00 171.18 23.81
1 2R 2R N CA 9 10 148.02 141.90 -41.10
2 3947 2R 3S C N 18 20 -142.89 -136.60 23.08 1.02 -64.10 170.85 18.03
2 3S 3S N CA 20 21 173.40 151.20 -35.00
3 3948 3S 4L C N 24 26 -65.27 -70.70 9.84 0.61 -63.50 174.60 24.14
3 4L 4L N CA 26 27 133.39 141.60 -41.20
4 3949 4L 5L C N 32 34 -113.87 -108.50 18.07 0.91 -63.50 176.39 27.60
4 5L 5L N CA 34 35 149.75 132.50 -41.20
5 3951 6I 7L C N 48 50 -63.06 -70.70 14.24 0.88 -63.50 170.79 23.74
5 7L 7L N CA 50 51 129.59 141.60 -41.20
6 3953 8V 9L C N 63 65 -122.71 -108.50 24.28 1.38 -63.50 165.00 20.42
6 9L 9L N CA 65 66 112.81 132.50 -41.20
7 3955 10C 11F C N 77 79 -156.89 -124.20 32.84 1.19 -63.20 -166.51 32.81
7 11F 11F N CA 79 80 146.41 143.30 -44.30
8 3956 11F 12L C N 88 90 -117.31 -108.50 13.71 0.78 -63.50 171.85 21.45
8 12L 12L N CA 90 91 122.00 132.50 -41.20
9 3958 13P 14L C N 103 105 -93.45 -108.50 32.34 1.56 -63.50 148.14 19.07
9 14L 14L N CA 105 106 103.88 132.50 -41.20
10 3959 14L 15A C N 111 113 -179.68 -134.00 52.05 1.20 -62.50 -171.87 36.95
10 15A 15A N CA 113 114 171.93 147.00 -40.90
11 3960 15A 16A C N 116 118 -137.74 -134.00 4.13 0.09 -62.50 -173.78 34.55
11 16A 16A N CA 118 119 148.76 147.00 -40.90
12 3961 16A 17L C N 121 123 -121.00 -108.50 13.06 0.59 -63.50 -173.44 23.32
12 17L 17L N CA 123 124 136.28 132.50 -41.20
13 3962 17L 18G C N 129 131 146.83 -167.20 46.31 0.76 82.20 -176.58 11.94
13 18G 18G N CA 131 132 -179.84 174.60 8.50
14 3963 18G 19K C N 133 135 -121.17 -118.00 18.90 0.87 -62.90 171.66 25.77
14 19K 19K N CA 135 136 157.73 139.10 -40.80
15 3998 53E 54S C N 410 412 -137.71 -64.10 83.09 8.52 -64.10 83.09 8.52
15 54S 54S N CA 412 413 3.55 -35.00 -35.00
16 4064 119D 120G C N 909 911 49.59 82.20 85.39 5.66 -62.40 115.74 19.29
16 120G 120G N CA 911 912 -70.42 8.50 -41.20
report______> Distribution of short non-bonded contacts:
DISTANCE1: 0.00 2.10 2.20 2.30 2.40 2.50 2.60 2.70 2.80 2.90 3.00 3.10 3.20 3.30 3.40
DISTANCE2: 2.10 2.20 2.30 2.40 2.50 2.60 2.70 2.80 2.90 3.00 3.10 3.20 3.30 3.40 3.50
FREQUENCY: 0 0 0 0 1 9 14 57 69 129 120 143 179 188 198
<< end of ENERGY.
>> Summary of successfully produced models:
Filename molpdf
----------------------------------------
good.B99990001.pdb 875.73004
In [29]:
good_protein = nglview.show_structure_file('good.B99990001.pdb')
good_protein
In [30]:
Image('good_protein.png')
Out[30]:
