Starting phenix.real_space_refine on Wed Sep 25 16:08:56 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rpw_24624/09_2024/7rpw_24624.cif Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rpw_24624/09_2024/7rpw_24624.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=4.38 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rpw_24624/09_2024/7rpw_24624.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rpw_24624/09_2024/7rpw_24624.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rpw_24624/09_2024/7rpw_24624.cif" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rpw_24624/09_2024/7rpw_24624.cif" } resolution = 4.38 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= 0.001 sd= 0.101 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 6 Type Number sf(0) Gaussians Mn 4 7.51 5 P 63 5.49 5 S 25 5.16 5 C 6567 2.51 5 N 1750 2.21 5 O 2185 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 2 residue(s): 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/chem_data/mon_lib" Total number of atoms: 10594 Number of models: 1 Model: "" Number of chains: 10 Chain: "A" Number of atoms: 1953 Number of conformers: 1 Conformer: "" Number of residues, atoms: 252, 1953 Classifications: {'peptide': 252} Link IDs: {'PTRANS': 8, 'TRANS': 243} Chain: "B" Number of atoms: 1922 Number of conformers: 1 Conformer: "" Number of residues, atoms: 243, 1922 Classifications: {'peptide': 243} Link IDs: {'PTRANS': 9, 'TRANS': 233} Chain: "C" Number of atoms: 1941 Number of conformers: 1 Conformer: "" Number of residues, atoms: 245, 1941 Classifications: {'peptide': 245} Link IDs: {'PTRANS': 5, 'TRANS': 239} Chain: "X" Number of atoms: 409 Number of conformers: 1 Conformer: "" Number of residues, atoms: 20, 409 Classifications: {'DNA': 20} Link IDs: {'rna3p': 19} Chain: "Y" Number of atoms: 229 Number of conformers: 1 Conformer: "" Number of residues, atoms: 11, 229 Classifications: {'DNA': 11} Modifications used: {'p5*END': 1} Link IDs: {'rna3p': 10} Chain: "Z" Number of atoms: 634 Number of conformers: 1 Conformer: "" Number of residues, atoms: 31, 634 Classifications: {'DNA': 31} Link IDs: {'rna3p': 30} Chain: "E" Number of atoms: 3480 Number of conformers: 1 Conformer: "" Number of residues, atoms: 439, 3480 Classifications: {'peptide': 439} Incomplete info: {'truncation_to_alanine': 2} Link IDs: {'PTRANS': 13, 'TRANS': 425} Unresolved non-hydrogen bonds: 9 Unresolved non-hydrogen angles: 11 Unresolved non-hydrogen dihedrals: 7 Planarities with less than four sites: {'ARG:plan': 1, 'ASP:plan': 1} Unresolved non-hydrogen planarities: 8 Chain: "X" Number of atoms: 2 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 2 Unusual residues: {' MN': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Chain: "Y" Number of atoms: 22 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 22 Classifications: {'RNA': 1} Modifications used: {'rna2p': 1} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 3 Unresolved non-hydrogen dihedrals: 1 Chain: "E" Number of atoms: 2 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 2 Unusual residues: {' MN': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Time building chain proxies: 6.63, per 1000 atoms: 0.63 Number of scatterers: 10594 At special positions: 0 Unit cell: (117, 106.5, 123, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 6 Type Number sf(0) Mn 4 24.99 S 25 16.00 P 63 15.00 O 2185 8.00 N 1750 7.00 C 6567 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.00 Amino acid : True - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=1, symmetry=0 Number of additional bonds: simple=1, symmetry=0 Coordination: Other bonds: Time building additional restraints: 2.86 Conformation dependent library (CDL) restraints added in 1.3 seconds 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2218 Finding SS restraints... Secondary structure from input PDB file: 31 helices and 15 sheets defined 28.1% alpha, 16.5% beta 30 base pairs and 46 stacking pairs defined. Time for finding SS restraints: 3.78 Creating SS restraints... Processing helix chain 'A' and resid 8 through 21 Processing helix chain 'A' and resid 74 through 80 removed outlier: 3.712A pdb=" N ILE A 78 " --> pdb=" O SER A 74 " (cutoff:3.500A) Processing helix chain 'A' and resid 140 through 151 removed outlier: 3.864A pdb=" N ASN A 144 " --> pdb=" O GLU A 140 " (cutoff:3.500A) Processing helix chain 'A' and resid 161 through 165 Processing helix chain 'A' and resid 202 through 210 Processing helix chain 'B' and resid 10 through 20 Processing helix chain 'B' and resid 72 through 81 removed outlier: 3.501A pdb=" N VAL B 76 " --> pdb=" O LYS B 72 " (cutoff:3.500A) Processing helix chain 'B' and resid 136 through 147 removed outlier: 4.131A pdb=" N LEU B 147 " --> pdb=" O ILE B 143 " (cutoff:3.500A) Processing helix chain 'B' and resid 198 through 204 removed outlier: 3.584A pdb=" N VAL B 202 " --> pdb=" O GLY B 198 " (cutoff:3.500A) Processing helix chain 'B' and resid 205 through 207 No H-bonds generated for 'chain 'B' and resid 205 through 207' Processing helix chain 'C' and resid 9 through 20 removed outlier: 3.705A pdb=" N ASP C 13 " --> pdb=" O ARG C 9 " (cutoff:3.500A) removed outlier: 3.565A pdb=" N ILE C 14 " --> pdb=" O VAL C 10 " (cutoff:3.500A) removed outlier: 5.692A pdb=" N ILE C 15 " --> pdb=" O LEU C 11 " (cutoff:3.500A) removed outlier: 4.322A pdb=" N LEU C 18 " --> pdb=" O ILE C 14 " (cutoff:3.500A) removed outlier: 4.105A pdb=" N ARG C 20 " --> pdb=" O GLN C 16 " (cutoff:3.500A) Processing helix chain 'C' and resid 70 through 79 removed outlier: 3.662A pdb=" N ILE C 77 " --> pdb=" O TYR C 73 " (cutoff:3.500A) removed outlier: 3.751A pdb=" N LEU C 78 " --> pdb=" O LEU C 74 " (cutoff:3.500A) removed outlier: 3.918A pdb=" N LYS C 79 " --> pdb=" O MET C 75 " (cutoff:3.500A) Processing helix chain 'C' and resid 135 through 147 removed outlier: 4.070A pdb=" N LYS C 140 " --> pdb=" O SER C 136 " (cutoff:3.500A) removed outlier: 3.767A pdb=" N SER C 141 " --> pdb=" O ASP C 137 " (cutoff:3.500A) Processing helix chain 'E' and resid 2 through 12 removed outlier: 4.706A pdb=" N ILE E 6 " --> pdb=" O GLU E 2 " (cutoff:3.500A) Processing helix chain 'E' and resid 18 through 33 removed outlier: 3.804A pdb=" N THR E 23 " --> pdb=" O ARG E 19 " (cutoff:3.500A) removed outlier: 4.805A pdb=" N ALA E 24 " --> pdb=" O LEU E 20 " (cutoff:3.500A) Processing helix chain 'E' and resid 37 through 47 Processing helix chain 'E' and resid 51 through 55 removed outlier: 3.522A pdb=" N LEU E 54 " --> pdb=" O PRO E 51 " (cutoff:3.500A) Processing helix chain 'E' and resid 63 through 75 Processing helix chain 'E' and resid 77 through 88 Processing helix chain 'E' and resid 90 through 97 removed outlier: 3.840A pdb=" N ARG E 97 " --> pdb=" O GLU E 93 " (cutoff:3.500A) Processing helix chain 'E' and resid 120 through 133 removed outlier: 4.304A pdb=" N LEU E 133 " --> pdb=" O SER E 129 " (cutoff:3.500A) Processing helix chain 'E' and resid 138 through 153 removed outlier: 4.557A pdb=" N LEU E 142 " --> pdb=" O GLY E 138 " (cutoff:3.500A) Processing helix chain 'E' and resid 160 through 168 removed outlier: 4.067A pdb=" N ARG E 164 " --> pdb=" O LYS E 160 " (cutoff:3.500A) Processing helix chain 'E' and resid 175 through 189 removed outlier: 3.514A pdb=" N VAL E 179 " --> pdb=" O GLY E 175 " (cutoff:3.500A) Processing helix chain 'E' and resid 192 through 204 removed outlier: 4.544A pdb=" N ILE E 196 " --> pdb=" O SER E 192 " (cutoff:3.500A) removed outlier: 3.734A pdb=" N ILE E 197 " --> pdb=" O ALA E 193 " (cutoff:3.500A) removed outlier: 3.598A pdb=" N ARG E 204 " --> pdb=" O ALA E 200 " (cutoff:3.500A) Processing helix chain 'E' and resid 207 through 217 Processing helix chain 'E' and resid 246 through 252 Processing helix chain 'E' and resid 291 through 300 removed outlier: 4.118A pdb=" N TYR E 295 " --> pdb=" O ASP E 291 " (cutoff:3.500A) removed outlier: 4.078A pdb=" N VAL E 296 " --> pdb=" O VAL E 292 " (cutoff:3.500A) removed outlier: 4.398A pdb=" N SER E 297 " --> pdb=" O VAL E 293 " (cutoff:3.500A) removed outlier: 3.593A pdb=" N LYS E 298 " --> pdb=" O ASP E 294 " (cutoff:3.500A) Processing helix chain 'E' and resid 323 through 331 removed outlier: 4.102A pdb=" N HIS E 329 " --> pdb=" O GLN E 325 " (cutoff:3.500A) removed outlier: 4.205A pdb=" N ARG E 330 " --> pdb=" O GLU E 326 " (cutoff:3.500A) removed outlier: 3.638A pdb=" N LYS E 331 " --> pdb=" O LEU E 327 " (cutoff:3.500A) Processing helix chain 'E' and resid 366 through 375 removed outlier: 3.636A pdb=" N GLU E 373 " --> pdb=" O ARG E 369 " (cutoff:3.500A) removed outlier: 3.741A pdb=" N SER E 374 " --> pdb=" O LYS E 370 " (cutoff:3.500A) removed outlier: 3.515A pdb=" N ILE E 375 " --> pdb=" O LEU E 371 " (cutoff:3.500A) Processing helix chain 'E' and resid 392 through 407 Processing sheet with id=AA1, first strand: chain 'A' and resid 3 through 6 removed outlier: 3.646A pdb=" N ILE A 89 " --> pdb=" O TYR A 6 " (cutoff:3.500A) removed outlier: 3.854A pdb=" N THR A 92 " --> pdb=" O LYS A 99 " (cutoff:3.500A) removed outlier: 4.330A pdb=" N LYS A 99 " --> pdb=" O THR A 92 " (cutoff:3.500A) removed outlier: 3.517A pdb=" N GLU C 175 " --> pdb=" O LYS A 110 " (cutoff:3.500A) removed outlier: 3.700A pdb=" N THR A 112 " --> pdb=" O GLU C 173 " (cutoff:3.500A) removed outlier: 3.558A pdb=" N GLU C 173 " --> pdb=" O THR A 112 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 47 through 52 removed outlier: 4.116A pdb=" N ILE A 35 " --> pdb=" O ILE A 52 " (cutoff:3.500A) removed outlier: 3.785A pdb=" N ARG A 38 " --> pdb=" O ILE A 27 " (cutoff:3.500A) removed outlier: 4.554A pdb=" N ILE A 27 " --> pdb=" O ARG A 38 " (cutoff:3.500A) removed outlier: 7.819A pdb=" N LEU A 40 " --> pdb=" O SER A 25 " (cutoff:3.500A) removed outlier: 9.028A pdb=" N SER A 25 " --> pdb=" O LEU A 40 " (cutoff:3.500A) removed outlier: 3.590A pdb=" N LEU A 70 " --> pdb=" O ILE A 26 " (cutoff:3.500A) removed outlier: 4.405A pdb=" N LEU A 28 " --> pdb=" O VAL A 68 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 134 through 138 removed outlier: 3.622A pdb=" N VAL A 134 " --> pdb=" O PHE A 223 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'A' and resid 179 through 181 removed outlier: 3.603A pdb=" N ALA A 179 " --> pdb=" O ILE A 168 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'B' and resid 3 through 6 removed outlier: 3.670A pdb=" N THR B 99 " --> pdb=" O SER B 92 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'B' and resid 27 through 30 Processing sheet with id=AA7, first strand: chain 'B' and resid 37 through 40 removed outlier: 3.630A pdb=" N ASP B 49 " --> pdb=" O ASP B 238 " (cutoff:3.500A) removed outlier: 4.060A pdb=" N PHE B 51 " --> pdb=" O TYR B 236 " (cutoff:3.500A) removed outlier: 4.242A pdb=" N TYR B 236 " --> pdb=" O PHE B 51 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'B' and resid 189 through 190 removed outlier: 3.694A pdb=" N SER B 189 " --> pdb=" O LYS B 132 " (cutoff:3.500A) removed outlier: 3.515A pdb=" N LYS B 132 " --> pdb=" O SER B 189 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'B' and resid 158 through 159 Processing sheet with id=AB1, first strand: chain 'B' and resid 165 through 168 removed outlier: 3.620A pdb=" N PHE B 165 " --> pdb=" O VAL B 176 " (cutoff:3.500A) removed outlier: 3.516A pdb=" N VAL C 97 " --> pdb=" O VAL C 111 " (cutoff:3.500A) removed outlier: 3.736A pdb=" N ILE C 87 " --> pdb=" O TYR C 5 " (cutoff:3.500A) removed outlier: 3.763A pdb=" N TYR C 5 " --> pdb=" O ILE C 87 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'C' and resid 26 through 29 removed outlier: 4.749A pdb=" N VAL C 26 " --> pdb=" O VAL C 37 " (cutoff:3.500A) removed outlier: 4.387A pdb=" N VAL C 37 " --> pdb=" O VAL C 26 " (cutoff:3.500A) removed outlier: 3.670A pdb=" N LEU C 51 " --> pdb=" O VAL C 34 " (cutoff:3.500A) removed outlier: 3.530A pdb=" N ILE C 47 " --> pdb=" O ALA C 38 " (cutoff:3.500A) removed outlier: 3.535A pdb=" N LEU C 46 " --> pdb=" O LEU C 239 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'C' and resid 130 through 131 removed outlier: 3.516A pdb=" N PHE C 218 " --> pdb=" O ILE C 130 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'C' and resid 162 through 165 removed outlier: 3.890A pdb=" N SER C 192 " --> pdb=" O GLY C 157 " (cutoff:3.500A) Processing sheet with id=AB5, first strand: chain 'E' and resid 239 through 241 removed outlier: 6.364A pdb=" N GLU E 239 " --> pdb=" O LYS E 433 " (cutoff:3.500A) Processing sheet with id=AB6, first strand: chain 'E' and resid 275 through 278 removed outlier: 4.203A pdb=" N ILE E 307 " --> pdb=" O ILE E 268 " (cutoff:3.500A) removed outlier: 6.420A pdb=" N ILE E 306 " --> pdb=" O MET E 353 " (cutoff:3.500A) 311 hydrogen bonds defined for protein. 852 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 77 hydrogen bonds 154 hydrogen bond angles 0 basepair planarities 30 basepair parallelities 46 stacking parallelities Total time for adding SS restraints: 3.71 Time building geometry restraints manager: 3.31 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.21 - 1.33: 2810 1.33 - 1.45: 1960 1.45 - 1.57: 5949 1.57 - 1.69: 122 1.69 - 1.81: 50 Bond restraints: 10891 Sorted by residual: bond pdb=" C3' DA Z 22 " pdb=" O3' DA Z 22 " ideal model delta sigma weight residual 1.422 1.474 -0.052 3.00e-02 1.11e+03 3.02e+00 bond pdb=" C4' AMP Y 101 " pdb=" O4' AMP Y 101 " ideal model delta sigma weight residual 1.454 1.431 0.023 1.50e-02 4.44e+03 2.37e+00 bond pdb=" C3' DA Z 13 " pdb=" O3' DA Z 13 " ideal model delta sigma weight residual 1.422 1.464 -0.042 3.00e-02 1.11e+03 1.92e+00 bond pdb=" C1' AMP Y 101 " pdb=" N9 AMP Y 101 " ideal model delta sigma weight residual 1.465 1.483 -0.018 1.50e-02 4.44e+03 1.48e+00 bond pdb=" C2' AMP Y 101 " pdb=" O2' AMP Y 101 " ideal model delta sigma weight residual 1.412 1.429 -0.017 1.50e-02 4.44e+03 1.32e+00 ... (remaining 10886 not shown) Histogram of bond angle deviations from ideal: 0.00 - 1.25: 14372 1.25 - 2.50: 463 2.50 - 3.75: 92 3.75 - 5.00: 21 5.00 - 6.25: 6 Bond angle restraints: 14954 Sorted by residual: angle pdb=" C1' AMP Y 101 " pdb=" C2' AMP Y 101 " pdb=" C3' AMP Y 101 " ideal model delta sigma weight residual 101.50 105.38 -3.88 1.00e+00 1.00e+00 1.50e+01 angle pdb=" N GLU E 366 " pdb=" CA GLU E 366 " pdb=" C GLU E 366 " ideal model delta sigma weight residual 113.18 109.01 4.17 1.33e+00 5.65e-01 9.81e+00 angle pdb=" C3' AMP Y 101 " pdb=" C2' AMP Y 101 " pdb=" O2' AMP Y 101 " ideal model delta sigma weight residual 114.60 110.03 4.57 1.50e+00 4.44e-01 9.29e+00 angle pdb=" C2' AMP Y 101 " pdb=" C3' AMP Y 101 " pdb=" C4' AMP Y 101 " ideal model delta sigma weight residual 102.60 105.50 -2.90 1.00e+00 1.00e+00 8.43e+00 angle pdb=" C2' AMP Y 101 " pdb=" C1' AMP Y 101 " pdb=" N9 AMP Y 101 " ideal model delta sigma weight residual 114.00 109.91 4.09 1.50e+00 4.44e-01 7.44e+00 ... (remaining 14949 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 34.68: 6066 34.68 - 69.36: 385 69.36 - 104.04: 5 104.04 - 138.72: 0 138.72 - 173.39: 5 Dihedral angle restraints: 6461 sinusoidal: 3029 harmonic: 3432 Sorted by residual: dihedral pdb=" CA SER E 139 " pdb=" C SER E 139 " pdb=" N ARG E 140 " pdb=" CA ARG E 140 " ideal model delta harmonic sigma weight residual 180.00 -155.45 -24.55 0 5.00e+00 4.00e-02 2.41e+01 dihedral pdb=" CA LEU C 205 " pdb=" C LEU C 205 " pdb=" N SER C 206 " pdb=" CA SER C 206 " ideal model delta harmonic sigma weight residual 180.00 157.66 22.34 0 5.00e+00 4.00e-02 2.00e+01 dihedral pdb=" CA LEU C 11 " pdb=" C LEU C 11 " pdb=" N LYS C 12 " pdb=" CA LYS C 12 " ideal model delta harmonic sigma weight residual 180.00 159.44 20.56 0 5.00e+00 4.00e-02 1.69e+01 ... (remaining 6458 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.035: 1173 0.035 - 0.069: 372 0.069 - 0.104: 120 0.104 - 0.139: 57 0.139 - 0.173: 2 Chirality restraints: 1724 Sorted by residual: chirality pdb=" CB VAL C 244 " pdb=" CA VAL C 244 " pdb=" CG1 VAL C 244 " pdb=" CG2 VAL C 244 " both_signs ideal model delta sigma weight residual False -2.63 -2.46 -0.17 2.00e-01 2.50e+01 7.50e-01 chirality pdb=" P DG Z 27 " pdb=" OP1 DG Z 27 " pdb=" OP2 DG Z 27 " pdb=" O5' DG Z 27 " both_signs ideal model delta sigma weight residual True 2.35 -2.49 -0.14 2.00e-01 2.50e+01 5.13e-01 chirality pdb=" CA ILE E 231 " pdb=" N ILE E 231 " pdb=" C ILE E 231 " pdb=" CB ILE E 231 " both_signs ideal model delta sigma weight residual False 2.43 2.57 -0.14 2.00e-01 2.50e+01 4.77e-01 ... (remaining 1721 not shown) Planarity restraints: 1681 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C SER C 93 " 0.019 5.00e-02 4.00e+02 2.94e-02 1.39e+00 pdb=" N PRO C 94 " -0.051 5.00e-02 4.00e+02 pdb=" CA PRO C 94 " 0.015 5.00e-02 4.00e+02 pdb=" CD PRO C 94 " 0.017 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C LYS E 226 " -0.019 5.00e-02 4.00e+02 2.92e-02 1.37e+00 pdb=" N PRO E 227 " 0.051 5.00e-02 4.00e+02 pdb=" CA PRO E 227 " -0.015 5.00e-02 4.00e+02 pdb=" CD PRO E 227 " -0.016 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C GLN B 121 " -0.019 5.00e-02 4.00e+02 2.89e-02 1.33e+00 pdb=" N PRO B 122 " 0.050 5.00e-02 4.00e+02 pdb=" CA PRO B 122 " -0.015 5.00e-02 4.00e+02 pdb=" CD PRO B 122 " -0.016 5.00e-02 4.00e+02 ... (remaining 1678 not shown) Histogram of nonbonded interaction distances: 1.76 - 2.38: 62 2.38 - 3.01: 6110 3.01 - 3.64: 17020 3.64 - 4.27: 23587 4.27 - 4.90: 36833 Nonbonded interactions: 83612 Sorted by model distance: nonbonded pdb=" O GLY E 60 " pdb="MN MN X 101 " model vdw 1.756 2.320 nonbonded pdb=" OP1 DG X 22 " pdb="MN MN E 701 " model vdw 1.835 2.320 nonbonded pdb=" O3' DG X 20 " pdb="MN MN X 101 " model vdw 1.913 2.320 nonbonded pdb=" O LYS E 130 " pdb=" OG1 THR E 134 " model vdw 2.077 3.040 nonbonded pdb=" OP1 DC X 21 " pdb="MN MN X 102 " model vdw 2.146 2.320 ... (remaining 83607 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.00 ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=0.50 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 4.970 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.000 Extract box with map and model: 0.340 Check model and map are aligned: 0.070 Set scattering table: 0.110 Process input model: 29.630 Find NCS groups from input model: 0.140 Set up NCS constraints: 0.040 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.010 Load rotamer database and sin/cos tables:3.250 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 38.560 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6815 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.052 10891 Z= 0.205 Angle : 0.570 6.255 14954 Z= 0.315 Chirality : 0.042 0.173 1724 Planarity : 0.003 0.029 1681 Dihedral : 19.295 173.395 4243 Min Nonbonded Distance : 1.756 Molprobity Statistics. All-atom Clashscore : 18.05 Ramachandran Plot: Outliers : 0.09 % Allowed : 11.27 % Favored : 88.64 % Rotamer: Outliers : 9.97 % Allowed : 13.46 % Favored : 76.57 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.82 (0.24), residues: 1171 helix: -0.42 (0.30), residues: 300 sheet: -1.68 (0.31), residues: 304 loop : -2.82 (0.24), residues: 567 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.003 0.001 TRP E 50 HIS 0.005 0.001 HIS B 157 PHE 0.018 0.001 PHE E 53 TYR 0.010 0.001 TYR C 196 ARG 0.004 0.000 ARG E 265 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 256 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 103 poor density : 153 time to evaluate : 1.180 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 37 SER cc_start: 0.8094 (t) cc_final: 0.7567 (m) REVERT: A 80 SER cc_start: 0.9022 (t) cc_final: 0.8779 (p) REVERT: A 137 THR cc_start: 0.8838 (OUTLIER) cc_final: 0.8102 (t) REVERT: B 21 ASP cc_start: 0.7999 (OUTLIER) cc_final: 0.7684 (t0) REVERT: B 89 LEU cc_start: 0.5247 (OUTLIER) cc_final: 0.4884 (tp) REVERT: C 44 ILE cc_start: 0.9077 (mt) cc_final: 0.8816 (mm) REVERT: C 75 MET cc_start: 0.7674 (OUTLIER) cc_final: 0.7341 (tpp) REVERT: C 146 VAL cc_start: 0.7875 (t) cc_final: 0.6797 (t) REVERT: E 123 GLU cc_start: 0.7720 (pp20) cc_final: 0.7429 (mp0) REVERT: E 142 LEU cc_start: 0.7441 (OUTLIER) cc_final: 0.7176 (mp) REVERT: E 164 ARG cc_start: 0.7939 (mtt180) cc_final: 0.7698 (tpt170) REVERT: E 328 MET cc_start: 0.8601 (mtm) cc_final: 0.7994 (tmm) REVERT: E 360 TYR cc_start: 0.4794 (OUTLIER) cc_final: 0.4463 (t80) REVERT: E 372 LEU cc_start: 0.6464 (tp) cc_final: 0.6179 (tp) outliers start: 103 outliers final: 40 residues processed: 239 average time/residue: 0.2254 time to fit residues: 75.5004 Evaluate side-chains 149 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 46 poor density : 103 time to evaluate : 1.178 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 122 random chunks: chunk 103 optimal weight: 4.9990 chunk 92 optimal weight: 6.9990 chunk 51 optimal weight: 6.9990 chunk 31 optimal weight: 9.9990 chunk 62 optimal weight: 4.9990 chunk 49 optimal weight: 0.9990 chunk 95 optimal weight: 1.9990 chunk 37 optimal weight: 0.5980 chunk 58 optimal weight: 7.9990 chunk 71 optimal weight: 0.0170 chunk 110 optimal weight: 5.9990 overall best weight: 1.7224 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 21 ASN E 391 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6710 moved from start: 0.2512 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.050 10891 Z= 0.211 Angle : 0.642 6.264 14954 Z= 0.348 Chirality : 0.046 0.271 1724 Planarity : 0.004 0.035 1681 Dihedral : 22.083 175.539 1931 Min Nonbonded Distance : 2.008 Molprobity Statistics. All-atom Clashscore : 15.26 Ramachandran Plot: Outliers : 0.09 % Allowed : 10.25 % Favored : 89.67 % Rotamer: Outliers : 0.19 % Allowed : 3.97 % Favored : 95.84 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.71 (0.24), residues: 1171 helix: -0.50 (0.29), residues: 308 sheet: -1.48 (0.31), residues: 309 loop : -2.77 (0.24), residues: 554 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.020 0.002 TRP E 50 HIS 0.004 0.001 HIS B 157 PHE 0.023 0.002 PHE E 53 TYR 0.015 0.002 TYR B 56 ARG 0.005 0.001 ARG E 234 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 159 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 157 time to evaluate : 1.266 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 40 LEU cc_start: 0.7312 (mp) cc_final: 0.6610 (mp) REVERT: B 136 LEU cc_start: 0.7190 (tp) cc_final: 0.6789 (tp) REVERT: B 146 GLU cc_start: 0.6387 (mt-10) cc_final: 0.6087 (tt0) REVERT: B 216 GLU cc_start: 0.3679 (mm-30) cc_final: 0.3413 (mm-30) REVERT: C 1 MET cc_start: 0.7514 (tpt) cc_final: 0.7055 (tpp) REVERT: C 44 ILE cc_start: 0.8849 (mt) cc_final: 0.8590 (mm) REVERT: C 80 VAL cc_start: 0.7491 (m) cc_final: 0.5928 (m) REVERT: C 220 ASN cc_start: 0.7578 (m-40) cc_final: 0.7369 (m-40) REVERT: E 8 GLU cc_start: 0.7705 (tm-30) cc_final: 0.6784 (pt0) REVERT: E 63 GLU cc_start: 0.7878 (mt-10) cc_final: 0.7548 (mt-10) REVERT: E 133 LEU cc_start: 0.7495 (mp) cc_final: 0.7168 (mt) REVERT: E 164 ARG cc_start: 0.7935 (mtt180) cc_final: 0.7709 (mmt-90) REVERT: E 328 MET cc_start: 0.8353 (mtm) cc_final: 0.8140 (tmm) REVERT: E 360 TYR cc_start: 0.6116 (m-10) cc_final: 0.5839 (t80) REVERT: E 372 LEU cc_start: 0.6422 (tp) cc_final: 0.6157 (tp) outliers start: 2 outliers final: 0 residues processed: 158 average time/residue: 0.2325 time to fit residues: 51.6181 Evaluate side-chains 97 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 97 time to evaluate : 1.097 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 122 random chunks: chunk 61 optimal weight: 9.9990 chunk 34 optimal weight: 10.0000 chunk 92 optimal weight: 0.0470 chunk 75 optimal weight: 8.9990 chunk 30 optimal weight: 9.9990 chunk 110 optimal weight: 0.7980 chunk 119 optimal weight: 5.9990 chunk 98 optimal weight: 7.9990 chunk 37 optimal weight: 0.9990 chunk 89 optimal weight: 0.5980 chunk 109 optimal weight: 10.0000 overall best weight: 1.6882 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 240 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6728 moved from start: 0.3099 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.050 10891 Z= 0.203 Angle : 0.617 6.078 14954 Z= 0.332 Chirality : 0.045 0.178 1724 Planarity : 0.004 0.034 1681 Dihedral : 22.142 174.121 1931 Min Nonbonded Distance : 2.042 Molprobity Statistics. All-atom Clashscore : 15.16 Ramachandran Plot: Outliers : 0.09 % Allowed : 10.59 % Favored : 89.33 % Rotamer: Outliers : 0.10 % Allowed : 3.78 % Favored : 96.13 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.68 (0.24), residues: 1171 helix: -0.32 (0.30), residues: 297 sheet: -1.59 (0.31), residues: 309 loop : -2.73 (0.24), residues: 565 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.019 0.002 TRP E 50 HIS 0.003 0.001 HIS B 157 PHE 0.017 0.002 PHE E 53 TYR 0.027 0.001 TYR C 59 ARG 0.005 0.001 ARG C 20 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 126 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 125 time to evaluate : 1.143 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: B 129 PHE cc_start: 0.6455 (m-10) cc_final: 0.6178 (m-80) REVERT: B 146 GLU cc_start: 0.6320 (mt-10) cc_final: 0.5611 (tt0) REVERT: C 44 ILE cc_start: 0.8831 (mt) cc_final: 0.8525 (mm) REVERT: C 80 VAL cc_start: 0.7551 (m) cc_final: 0.6167 (m) REVERT: C 140 LYS cc_start: 0.8431 (ttmt) cc_final: 0.8099 (tmtt) REVERT: C 150 THR cc_start: 0.8565 (p) cc_final: 0.8313 (t) REVERT: E 8 GLU cc_start: 0.7675 (tm-30) cc_final: 0.6827 (pt0) REVERT: E 63 GLU cc_start: 0.7924 (mt-10) cc_final: 0.7615 (mt-10) REVERT: E 133 LEU cc_start: 0.7593 (mp) cc_final: 0.7340 (mt) REVERT: E 328 MET cc_start: 0.8547 (mtm) cc_final: 0.8087 (tmm) REVERT: E 360 TYR cc_start: 0.5573 (m-10) cc_final: 0.5361 (t80) REVERT: E 372 LEU cc_start: 0.6438 (tp) cc_final: 0.6183 (tp) outliers start: 1 outliers final: 0 residues processed: 125 average time/residue: 0.2395 time to fit residues: 42.1113 Evaluate side-chains 92 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 92 time to evaluate : 1.079 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 122 random chunks: chunk 83 optimal weight: 10.0000 chunk 57 optimal weight: 4.9990 chunk 12 optimal weight: 3.9990 chunk 52 optimal weight: 5.9990 chunk 74 optimal weight: 10.0000 chunk 111 optimal weight: 8.9990 chunk 117 optimal weight: 20.0000 chunk 58 optimal weight: 8.9990 chunk 105 optimal weight: 20.0000 chunk 31 optimal weight: 7.9990 chunk 98 optimal weight: 0.7980 overall best weight: 4.7588 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 204 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6944 moved from start: 0.3387 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.049 10891 Z= 0.382 Angle : 0.770 9.983 14954 Z= 0.409 Chirality : 0.049 0.162 1724 Planarity : 0.005 0.037 1681 Dihedral : 22.306 171.526 1931 Min Nonbonded Distance : 1.976 Molprobity Statistics. All-atom Clashscore : 22.38 Ramachandran Plot: Outliers : 0.09 % Allowed : 13.75 % Favored : 86.17 % Rotamer: Outliers : 0.19 % Allowed : 5.81 % Favored : 94.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.09 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.23 (0.23), residues: 1171 helix: -0.83 (0.29), residues: 292 sheet: -1.86 (0.29), residues: 328 loop : -3.08 (0.24), residues: 551 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.027 0.002 TRP E 50 HIS 0.006 0.002 HIS A 240 PHE 0.032 0.003 PHE E 165 TYR 0.037 0.003 TYR C 59 ARG 0.007 0.001 ARG E 164 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 107 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 105 time to evaluate : 1.151 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 6 TYR cc_start: 0.6275 (t80) cc_final: 0.5661 (t80) REVERT: A 229 MET cc_start: 0.7728 (mmm) cc_final: 0.7103 (mmm) REVERT: A 241 MET cc_start: 0.5147 (mtp) cc_final: 0.4936 (mtp) REVERT: B 216 GLU cc_start: 0.4211 (mm-30) cc_final: 0.3801 (mm-30) REVERT: C 20 ARG cc_start: 0.8105 (mtt-85) cc_final: 0.7895 (mtt-85) REVERT: C 44 ILE cc_start: 0.8766 (mt) cc_final: 0.8524 (mm) REVERT: C 75 MET cc_start: 0.8399 (mmt) cc_final: 0.8131 (tpp) REVERT: C 140 LYS cc_start: 0.8539 (ttmt) cc_final: 0.8290 (tmtt) REVERT: C 186 LEU cc_start: 0.8465 (tp) cc_final: 0.8107 (tp) REVERT: C 213 TYR cc_start: 0.8214 (m-80) cc_final: 0.7628 (m-80) REVERT: C 230 MET cc_start: 0.4683 (mmt) cc_final: 0.3817 (mmt) REVERT: E 63 GLU cc_start: 0.7774 (mt-10) cc_final: 0.7553 (mt-10) REVERT: E 133 LEU cc_start: 0.7591 (mp) cc_final: 0.7279 (mt) REVERT: E 328 MET cc_start: 0.8574 (mtm) cc_final: 0.8156 (tmm) REVERT: E 360 TYR cc_start: 0.6316 (m-10) cc_final: 0.5138 (t80) REVERT: E 372 LEU cc_start: 0.6221 (tp) cc_final: 0.5844 (tp) outliers start: 2 outliers final: 1 residues processed: 107 average time/residue: 0.2311 time to fit residues: 35.1191 Evaluate side-chains 89 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 88 time to evaluate : 1.135 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 122 random chunks: chunk 66 optimal weight: 4.9990 chunk 1 optimal weight: 0.9990 chunk 87 optimal weight: 6.9990 chunk 48 optimal weight: 0.9990 chunk 100 optimal weight: 7.9990 chunk 81 optimal weight: 8.9990 chunk 0 optimal weight: 7.9990 chunk 60 optimal weight: 7.9990 chunk 105 optimal weight: 7.9990 chunk 29 optimal weight: 0.9980 chunk 39 optimal weight: 8.9990 overall best weight: 2.9988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: E 387 HIS E 392 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6863 moved from start: 0.3640 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.047 10891 Z= 0.271 Angle : 0.681 7.525 14954 Z= 0.365 Chirality : 0.047 0.174 1724 Planarity : 0.004 0.037 1681 Dihedral : 22.319 172.504 1931 Min Nonbonded Distance : 2.019 Molprobity Statistics. All-atom Clashscore : 18.87 Ramachandran Plot: Outliers : 0.09 % Allowed : 11.96 % Favored : 87.96 % Rotamer: Outliers : 0.00 % Allowed : 2.81 % Favored : 97.19 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.18 (0.23), residues: 1171 helix: -0.88 (0.28), residues: 299 sheet: -1.76 (0.31), residues: 301 loop : -3.04 (0.23), residues: 571 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.039 0.003 TRP E 50 HIS 0.007 0.002 HIS B 157 PHE 0.019 0.002 PHE C 228 TYR 0.030 0.002 TYR C 59 ARG 0.007 0.001 ARG E 204 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 109 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 109 time to evaluate : 1.044 Fit side-chains revert: symmetry clash revert: symmetry clash REVERT: A 6 TYR cc_start: 0.5813 (t80) cc_final: 0.4999 (t80) REVERT: A 229 MET cc_start: 0.7803 (mmm) cc_final: 0.7167 (mmm) REVERT: B 136 LEU cc_start: 0.7090 (tp) cc_final: 0.6804 (tp) REVERT: B 216 GLU cc_start: 0.4147 (mm-30) cc_final: 0.3912 (mm-30) REVERT: C 44 ILE cc_start: 0.8667 (mt) cc_final: 0.8388 (mm) REVERT: C 140 LYS cc_start: 0.8410 (ttmt) cc_final: 0.8190 (tmtt) REVERT: C 186 LEU cc_start: 0.8398 (tp) cc_final: 0.8113 (tp) REVERT: E 63 GLU cc_start: 0.7712 (mt-10) cc_final: 0.7467 (mt-10) REVERT: E 328 MET cc_start: 0.8551 (mtm) cc_final: 0.7980 (tmm) REVERT: E 360 TYR cc_start: 0.5808 (m-10) cc_final: 0.5021 (t80) REVERT: E 372 LEU cc_start: 0.6189 (tp) cc_final: 0.5827 (tp) outliers start: 0 outliers final: 0 residues processed: 109 average time/residue: 0.2235 time to fit residues: 35.0172 Evaluate side-chains 85 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 85 time to evaluate : 1.093 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 122 random chunks: chunk 106 optimal weight: 10.0000 chunk 23 optimal weight: 10.0000 chunk 69 optimal weight: 2.9990 chunk 29 optimal weight: 1.9990 chunk 118 optimal weight: 10.0000 chunk 97 optimal weight: 5.9990 chunk 54 optimal weight: 0.0980 chunk 9 optimal weight: 9.9990 chunk 39 optimal weight: 5.9990 chunk 61 optimal weight: 7.9990 chunk 113 optimal weight: 7.9990 overall best weight: 3.4188 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6898 moved from start: 0.3887 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.046 10891 Z= 0.296 Angle : 0.702 8.588 14954 Z= 0.374 Chirality : 0.047 0.182 1724 Planarity : 0.004 0.038 1681 Dihedral : 22.325 171.864 1931 Min Nonbonded Distance : 2.014 Molprobity Statistics. All-atom Clashscore : 19.88 Ramachandran Plot: Outliers : 0.09 % Allowed : 13.75 % Favored : 86.17 % Rotamer: Outliers : 0.10 % Allowed : 2.81 % Favored : 97.10 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.09 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.35 (0.23), residues: 1171 helix: -1.08 (0.28), residues: 301 sheet: -1.88 (0.30), residues: 310 loop : -3.10 (0.23), residues: 560 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.035 0.004 TRP E 50 HIS 0.010 0.001 HIS B 157 PHE 0.021 0.002 PHE C 228 TYR 0.035 0.002 TYR C 59 ARG 0.005 0.000 ARG C 20 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 100 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 99 time to evaluate : 1.229 Fit side-chains revert: symmetry clash revert: symmetry clash REVERT: A 6 TYR cc_start: 0.6317 (t80) cc_final: 0.5381 (t80) REVERT: A 229 MET cc_start: 0.8216 (mmm) cc_final: 0.7885 (mmm) REVERT: B 36 ARG cc_start: 0.6959 (mmt180) cc_final: 0.6599 (mmt90) REVERT: B 136 LEU cc_start: 0.7154 (tp) cc_final: 0.6862 (tp) REVERT: B 208 MET cc_start: 0.8947 (mmm) cc_final: 0.8739 (tpt) REVERT: B 216 GLU cc_start: 0.4217 (mm-30) cc_final: 0.3955 (mm-30) REVERT: C 1 MET cc_start: 0.7407 (tpp) cc_final: 0.7182 (tpp) REVERT: C 44 ILE cc_start: 0.8666 (mt) cc_final: 0.8402 (mm) REVERT: C 75 MET cc_start: 0.8128 (tpp) cc_final: 0.7908 (tpp) REVERT: C 140 LYS cc_start: 0.8460 (ttmt) cc_final: 0.8132 (tmtt) REVERT: C 186 LEU cc_start: 0.8394 (tp) cc_final: 0.8145 (tp) REVERT: C 213 TYR cc_start: 0.7772 (m-80) cc_final: 0.7448 (m-80) REVERT: E 63 GLU cc_start: 0.7746 (mt-10) cc_final: 0.7496 (mt-10) REVERT: E 328 MET cc_start: 0.8472 (mtm) cc_final: 0.7966 (tmm) REVERT: E 360 TYR cc_start: 0.5410 (m-10) cc_final: 0.4486 (t80) outliers start: 1 outliers final: 1 residues processed: 100 average time/residue: 0.2220 time to fit residues: 32.2494 Evaluate side-chains 81 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 80 time to evaluate : 1.068 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 122 random chunks: chunk 13 optimal weight: 0.6980 chunk 67 optimal weight: 9.9990 chunk 86 optimal weight: 10.0000 chunk 66 optimal weight: 7.9990 chunk 99 optimal weight: 2.9990 chunk 65 optimal weight: 0.9980 chunk 117 optimal weight: 0.9990 chunk 73 optimal weight: 7.9990 chunk 71 optimal weight: 7.9990 chunk 54 optimal weight: 8.9990 chunk 72 optimal weight: 0.8980 overall best weight: 1.3184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 157 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6760 moved from start: 0.4053 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.049 10891 Z= 0.186 Angle : 0.613 6.055 14954 Z= 0.329 Chirality : 0.044 0.162 1724 Planarity : 0.004 0.038 1681 Dihedral : 22.250 173.319 1931 Min Nonbonded Distance : 2.118 Molprobity Statistics. All-atom Clashscore : 13.81 Ramachandran Plot: Outliers : 0.09 % Allowed : 9.65 % Favored : 90.26 % Rotamer: Outliers : 0.00 % Allowed : 2.03 % Favored : 97.97 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.96 (0.23), residues: 1171 helix: -0.71 (0.29), residues: 296 sheet: -1.54 (0.31), residues: 300 loop : -2.94 (0.23), residues: 575 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.026 0.002 TRP E 50 HIS 0.004 0.001 HIS B 157 PHE 0.023 0.001 PHE C 228 TYR 0.023 0.001 TYR C 59 ARG 0.005 0.000 ARG A 38 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 107 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 107 time to evaluate : 1.195 Fit side-chains revert: symmetry clash REVERT: A 6 TYR cc_start: 0.5852 (t80) cc_final: 0.5123 (t80) REVERT: A 229 MET cc_start: 0.7747 (mmm) cc_final: 0.7288 (mmm) REVERT: B 36 ARG cc_start: 0.7134 (mmt180) cc_final: 0.6738 (mmt90) REVERT: B 126 ASN cc_start: 0.7459 (t0) cc_final: 0.7204 (t0) REVERT: C 44 ILE cc_start: 0.8774 (mt) cc_final: 0.8517 (mm) REVERT: C 75 MET cc_start: 0.7655 (tpp) cc_final: 0.7427 (tpp) REVERT: C 186 LEU cc_start: 0.8378 (tp) cc_final: 0.8119 (tp) REVERT: C 213 TYR cc_start: 0.7784 (m-80) cc_final: 0.7481 (m-80) REVERT: E 63 GLU cc_start: 0.7639 (mt-10) cc_final: 0.7427 (mt-10) REVERT: E 142 LEU cc_start: 0.7592 (mp) cc_final: 0.7379 (tt) REVERT: E 150 LEU cc_start: 0.8977 (mm) cc_final: 0.8755 (tp) REVERT: E 251 MET cc_start: 0.6450 (ttt) cc_final: 0.5982 (ttt) REVERT: E 328 MET cc_start: 0.8485 (mtm) cc_final: 0.8028 (tmm) outliers start: 0 outliers final: 0 residues processed: 107 average time/residue: 0.2101 time to fit residues: 32.8951 Evaluate side-chains 81 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 81 time to evaluate : 1.191 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 122 random chunks: chunk 46 optimal weight: 10.0000 chunk 70 optimal weight: 5.9990 chunk 35 optimal weight: 1.9990 chunk 23 optimal weight: 3.9990 chunk 22 optimal weight: 5.9990 chunk 74 optimal weight: 10.0000 chunk 80 optimal weight: 0.6980 chunk 58 optimal weight: 0.9990 chunk 10 optimal weight: 9.9990 chunk 92 optimal weight: 0.9990 chunk 106 optimal weight: 5.9990 overall best weight: 1.7388 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 157 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6778 moved from start: 0.4183 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.048 10891 Z= 0.203 Angle : 0.623 6.311 14954 Z= 0.332 Chirality : 0.044 0.163 1724 Planarity : 0.004 0.038 1681 Dihedral : 22.197 173.671 1931 Min Nonbonded Distance : 2.095 Molprobity Statistics. All-atom Clashscore : 14.68 Ramachandran Plot: Outliers : 0.09 % Allowed : 11.53 % Favored : 88.39 % Rotamer: Outliers : 0.00 % Allowed : 1.26 % Favored : 98.74 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.97 (0.23), residues: 1171 helix: -0.81 (0.29), residues: 305 sheet: -1.47 (0.31), residues: 294 loop : -2.95 (0.23), residues: 572 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.023 0.002 TRP E 50 HIS 0.006 0.001 HIS B 157 PHE 0.019 0.002 PHE C 228 TYR 0.025 0.001 TYR C 59 ARG 0.003 0.000 ARG B 210 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 106 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 106 time to evaluate : 1.128 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 6 TYR cc_start: 0.5912 (t80) cc_final: 0.5059 (t80) REVERT: A 229 MET cc_start: 0.8019 (mmm) cc_final: 0.7705 (mmm) REVERT: B 121 GLN cc_start: 0.8489 (tm-30) cc_final: 0.7810 (tm-30) REVERT: B 126 ASN cc_start: 0.7550 (t0) cc_final: 0.7325 (t0) REVERT: C 44 ILE cc_start: 0.8806 (mt) cc_final: 0.8532 (mm) REVERT: C 186 LEU cc_start: 0.8302 (tp) cc_final: 0.8088 (tp) REVERT: C 213 TYR cc_start: 0.7857 (m-80) cc_final: 0.7498 (m-80) REVERT: E 63 GLU cc_start: 0.7615 (mt-10) cc_final: 0.7400 (mt-10) REVERT: E 294 ASP cc_start: 0.8298 (m-30) cc_final: 0.8077 (p0) REVERT: E 328 MET cc_start: 0.8488 (mtm) cc_final: 0.8019 (tmm) outliers start: 0 outliers final: 0 residues processed: 106 average time/residue: 0.2071 time to fit residues: 32.4143 Evaluate side-chains 84 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 84 time to evaluate : 1.067 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 122 random chunks: chunk 112 optimal weight: 0.0870 chunk 102 optimal weight: 9.9990 chunk 109 optimal weight: 3.9990 chunk 65 optimal weight: 9.9990 chunk 47 optimal weight: 7.9990 chunk 86 optimal weight: 2.9990 chunk 33 optimal weight: 5.9990 chunk 99 optimal weight: 3.9990 chunk 103 optimal weight: 10.0000 chunk 71 optimal weight: 8.9990 chunk 115 optimal weight: 2.9990 overall best weight: 2.8166 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 157 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6847 moved from start: 0.4247 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.045 10891 Z= 0.258 Angle : 0.669 7.778 14954 Z= 0.356 Chirality : 0.046 0.156 1724 Planarity : 0.004 0.035 1681 Dihedral : 22.243 173.318 1931 Min Nonbonded Distance : 2.052 Molprobity Statistics. All-atom Clashscore : 16.89 Ramachandran Plot: Outliers : 0.09 % Allowed : 12.72 % Favored : 87.19 % Rotamer: Outliers : 0.00 % Allowed : 0.39 % Favored : 99.61 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.14 (0.23), residues: 1171 helix: -0.99 (0.29), residues: 305 sheet: -1.64 (0.31), residues: 299 loop : -2.99 (0.23), residues: 567 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.046 0.004 TRP E 50 HIS 0.007 0.001 HIS B 157 PHE 0.021 0.002 PHE C 228 TYR 0.034 0.002 TYR C 59 ARG 0.003 0.000 ARG E 164 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 101 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 101 time to evaluate : 1.099 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 6 TYR cc_start: 0.6217 (t80) cc_final: 0.5354 (t80) REVERT: B 121 GLN cc_start: 0.8297 (tm-30) cc_final: 0.7806 (tm-30) REVERT: C 1 MET cc_start: 0.7468 (tpp) cc_final: 0.7162 (tpp) REVERT: C 44 ILE cc_start: 0.8692 (mt) cc_final: 0.8450 (mm) REVERT: C 95 ASP cc_start: 0.6284 (m-30) cc_final: 0.6081 (m-30) REVERT: C 186 LEU cc_start: 0.8341 (tp) cc_final: 0.8123 (tp) REVERT: C 213 TYR cc_start: 0.7758 (m-80) cc_final: 0.7422 (m-80) REVERT: E 25 LEU cc_start: 0.9094 (mm) cc_final: 0.8839 (pp) REVERT: E 328 MET cc_start: 0.8510 (mtm) cc_final: 0.8013 (tmm) outliers start: 0 outliers final: 0 residues processed: 101 average time/residue: 0.2103 time to fit residues: 31.3572 Evaluate side-chains 83 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 83 time to evaluate : 1.239 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 122 random chunks: chunk 70 optimal weight: 2.9990 chunk 54 optimal weight: 0.9980 chunk 80 optimal weight: 4.9990 chunk 121 optimal weight: 8.9990 chunk 111 optimal weight: 7.9990 chunk 96 optimal weight: 4.9990 chunk 10 optimal weight: 5.9990 chunk 74 optimal weight: 8.9990 chunk 59 optimal weight: 0.9990 chunk 76 optimal weight: 10.0000 chunk 103 optimal weight: 8.9990 overall best weight: 2.9988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: E 21 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6863 moved from start: 0.4388 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.045 10891 Z= 0.276 Angle : 0.681 7.675 14954 Z= 0.362 Chirality : 0.046 0.185 1724 Planarity : 0.004 0.035 1681 Dihedral : 22.265 173.083 1931 Min Nonbonded Distance : 2.065 Molprobity Statistics. All-atom Clashscore : 18.39 Ramachandran Plot: Outliers : 0.09 % Allowed : 12.72 % Favored : 87.19 % Rotamer: Outliers : 0.00 % Allowed : 0.39 % Favored : 99.61 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.16 (0.23), residues: 1171 helix: -1.00 (0.28), residues: 301 sheet: -1.68 (0.30), residues: 304 loop : -3.00 (0.23), residues: 566 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.045 0.003 TRP E 50 HIS 0.005 0.001 HIS B 157 PHE 0.020 0.002 PHE C 228 TYR 0.033 0.002 TYR C 59 ARG 0.005 0.001 ARG C 20 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2342 Ramachandran restraints generated. 1171 Oldfield, 0 Emsley, 1171 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 95 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 95 time to evaluate : 1.097 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 6 TYR cc_start: 0.6040 (t80) cc_final: 0.5125 (t80) REVERT: B 2 MET cc_start: 0.6933 (tmm) cc_final: 0.6697 (tmm) REVERT: B 121 GLN cc_start: 0.8433 (tm-30) cc_final: 0.7924 (tm-30) REVERT: C 1 MET cc_start: 0.7270 (tpp) cc_final: 0.6987 (tpp) REVERT: C 44 ILE cc_start: 0.8701 (mt) cc_final: 0.8462 (mm) REVERT: C 186 LEU cc_start: 0.8373 (tp) cc_final: 0.8161 (tp) REVERT: C 213 TYR cc_start: 0.7788 (m-80) cc_final: 0.7466 (m-80) REVERT: E 25 LEU cc_start: 0.9072 (mm) cc_final: 0.8775 (pp) REVERT: E 142 LEU cc_start: 0.7674 (tt) cc_final: 0.7436 (tt) REVERT: E 328 MET cc_start: 0.8574 (mtm) cc_final: 0.8074 (tmm) outliers start: 0 outliers final: 0 residues processed: 95 average time/residue: 0.2071 time to fit residues: 29.5736 Evaluate side-chains 81 residues out of total 1035 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 81 time to evaluate : 0.958 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 122 random chunks: chunk 29 optimal weight: 0.6980 chunk 89 optimal weight: 7.9990 chunk 14 optimal weight: 5.9990 chunk 26 optimal weight: 1.9990 chunk 96 optimal weight: 8.9990 chunk 40 optimal weight: 7.9990 chunk 99 optimal weight: 10.0000 chunk 12 optimal weight: 3.9990 chunk 17 optimal weight: 0.5980 chunk 85 optimal weight: 0.8980 chunk 5 optimal weight: 9.9990 overall best weight: 1.6384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4214 r_free = 0.4214 target = 0.111551 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 34)----------------| | r_work = 0.3907 r_free = 0.3907 target = 0.093612 restraints weight = 46153.703| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 27)----------------| | r_work = 0.3901 r_free = 0.3901 target = 0.093807 restraints weight = 41236.349| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 27)----------------| | r_work = 0.3912 r_free = 0.3912 target = 0.094498 restraints weight = 37614.788| |-----------------------------------------------------------------------------| r_work (final): 0.3889 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3889 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3889 r_free = 0.3889 target_work(ls_wunit_k1) = 0.093 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3887 r_free = 0.3887 target_work(ls_wunit_k1) = 0.093 | | occupancies: max = 1.00 min = 0.00 number of occupancies < 0.1: 5 | |-----------------------------------------------------------------------------| r_final: 0.3887 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6825 moved from start: 0.4514 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.047 10891 Z= 0.201 Angle : 0.624 5.898 14954 Z= 0.332 Chirality : 0.045 0.173 1724 Planarity : 0.004 0.037 1681 Dihedral : 22.211 173.669 1931 Min Nonbonded Distance : 2.123 Molprobity Statistics. All-atom Clashscore : 15.02 Ramachandran Plot: Outliers : 0.09 % Allowed : 10.50 % Favored : 89.41 % Rotamer: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.97 (0.23), residues: 1171 helix: -0.84 (0.29), residues: 305 sheet: -1.47 (0.32), residues: 286 loop : -2.92 (0.23), residues: 580 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.036 0.003 TRP E 50 HIS 0.005 0.001 HIS B 157 PHE 0.014 0.001 PHE C 228 TYR 0.027 0.001 TYR C 59 ARG 0.006 0.000 ARG C 20 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 1880.50 seconds wall clock time: 34 minutes 44.31 seconds (2084.31 seconds total)