Starting phenix.real_space_refine (version: 1.21rc1) on Wed Oct 4 07:47:08 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wy0_32883/10_2023/7wy0_32883_neut.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wy0_32883/10_2023/7wy0_32883.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.83 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wy0_32883/10_2023/7wy0_32883.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wy0_32883/10_2023/7wy0_32883.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wy0_32883/10_2023/7wy0_32883_neut.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wy0_32883/10_2023/7wy0_32883_neut.pdb" } resolution = 2.83 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.020 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: 4 Type Number sf(0) Gaussians S 56 5.16 5 C 5591 2.51 5 N 1481 2.21 5 O 1595 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A TYR 210": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 32": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5097/modules/chem_data/mon_lib" Total number of atoms: 8723 Number of models: 1 Model: "" Number of chains: 5 Chain: "A" Number of atoms: 1732 Number of conformers: 1 Conformer: "" Number of residues, atoms: 209, 1732 Classifications: {'peptide': 209} Link IDs: {'PTRANS': 3, 'TRANS': 205} Chain breaks: 4 Chain: "B" Number of atoms: 2600 Number of conformers: 1 Conformer: "" Number of residues, atoms: 338, 2600 Classifications: {'peptide': 338} Link IDs: {'PTRANS': 5, 'TRANS': 332} Chain: "C" Number of atoms: 1776 Number of conformers: 1 Conformer: "" Number of residues, atoms: 231, 1776 Classifications: {'peptide': 231} Link IDs: {'PCIS': 1, 'PTRANS': 9, 'TRANS': 220} Chain breaks: 1 Chain: "D" Number of atoms: 436 Number of conformers: 1 Conformer: "" Number of residues, atoms: 57, 436 Classifications: {'peptide': 57} Link IDs: {'PTRANS': 4, 'TRANS': 52} Chain: "R" Number of atoms: 2179 Number of conformers: 1 Conformer: "" Number of residues, atoms: 277, 2179 Classifications: {'peptide': 277} Link IDs: {'PTRANS': 10, 'TRANS': 266} Chain breaks: 1 Time building chain proxies: 5.09, per 1000 atoms: 0.58 Number of scatterers: 8723 At special positions: 0 Unit cell: (91.8, 119.88, 127.44, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 4 Type Number sf(0) S 56 16.00 O 1595 8.00 N 1481 7.00 C 5591 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=2, symmetry=0 Simple disulfide: pdb=" SG CYS C 161 " - pdb=" SG CYS C 231 " distance=2.04 Simple disulfide: pdb=" SG CYS R 657 " - pdb=" SG CYS R 733 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 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=0, symmetry=0 Time building additional restraints: 3.33 Conformation dependent library (CDL) restraints added in 1.4 seconds 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2078 Finding SS restraints... Secondary structure from input PDB file: 30 helices and 13 sheets defined 36.8% alpha, 23.7% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.87 Creating SS restraints... Processing helix chain 'A' and resid 14 through 40 removed outlier: 3.860A pdb=" N THR A 35 " --> pdb=" O SER A 31 " (cutoff:3.500A) removed outlier: 4.426A pdb=" N TYR A 36 " --> pdb=" O ARG A 32 " (cutoff:3.500A) Processing helix chain 'A' and resid 50 through 66 Processing helix chain 'A' and resid 236 through 241 removed outlier: 4.060A pdb=" N SER A 240 " --> pdb=" O CYS A 237 " (cutoff:3.500A) Processing helix chain 'A' and resid 267 through 278 Processing helix chain 'A' and resid 279 through 283 removed outlier: 3.555A pdb=" N SER A 283 " --> pdb=" O ARG A 280 " (cutoff:3.500A) Processing helix chain 'A' and resid 293 through 303 Processing helix chain 'A' and resid 320 through 332 removed outlier: 3.764A pdb=" N GLN A 324 " --> pdb=" O LEU A 320 " (cutoff:3.500A) Processing helix chain 'A' and resid 355 through 375 removed outlier: 3.674A pdb=" N ARG A 361 " --> pdb=" O ARG A 357 " (cutoff:3.500A) removed outlier: 4.055A pdb=" N ASP A 362 " --> pdb=" O LEU A 358 " (cutoff:3.500A) Processing helix chain 'B' and resid 4 through 26 removed outlier: 3.734A pdb=" N ALA B 26 " --> pdb=" O ARG B 22 " (cutoff:3.500A) Processing helix chain 'B' and resid 29 through 34 removed outlier: 3.624A pdb=" N ILE B 33 " --> pdb=" O THR B 29 " (cutoff:3.500A) Processing helix chain 'C' and resid 28 through 32 Processing helix chain 'C' and resid 87 through 91 Processing helix chain 'D' and resid 7 through 24 removed outlier: 3.532A pdb=" N ASN D 24 " --> pdb=" O LYS D 20 " (cutoff:3.500A) Processing helix chain 'D' and resid 29 through 44 removed outlier: 3.595A pdb=" N ALA D 33 " --> pdb=" O LYS D 29 " (cutoff:3.500A) Processing helix chain 'R' and resid 568 through 573 Processing helix chain 'R' and resid 581 through 610 removed outlier: 4.319A pdb=" N TRP R 587 " --> pdb=" O PRO R 583 " (cutoff:3.500A) Processing helix chain 'R' and resid 611 through 614 Processing helix chain 'R' and resid 617 through 647 Processing helix chain 'R' and resid 658 through 687 Processing helix chain 'R' and resid 694 through 707 Processing helix chain 'R' and resid 707 through 722 Processing helix chain 'R' and resid 723 through 725 No H-bonds generated for 'chain 'R' and resid 723 through 725' Processing helix chain 'R' and resid 742 through 746 removed outlier: 3.501A pdb=" N LEU R 745 " --> pdb=" O LYS R 742 " (cutoff:3.500A) Processing helix chain 'R' and resid 747 through 770 Processing helix chain 'R' and resid 785 through 800 Proline residue: R 798 - end of helix Processing helix chain 'R' and resid 802 through 805 Processing helix chain 'R' and resid 806 through 813 Processing helix chain 'R' and resid 819 through 828 Processing helix chain 'R' and resid 829 through 838 Processing helix chain 'R' and resid 842 through 854 Processing sheet with id=AA1, first strand: chain 'A' and resid 210 through 213 removed outlier: 6.277A pdb=" N VAL A 41 " --> pdb=" O LYS A 220 " (cutoff:3.500A) removed outlier: 6.502A pdb=" N LYS A 42 " --> pdb=" O SER A 243 " (cutoff:3.500A) removed outlier: 7.369A pdb=" N LEU A 245 " --> pdb=" O LYS A 42 " (cutoff:3.500A) removed outlier: 6.275A pdb=" N LEU A 44 " --> pdb=" O LEU A 245 " (cutoff:3.500A) removed outlier: 7.397A pdb=" N LEU A 247 " --> pdb=" O LEU A 44 " (cutoff:3.500A) removed outlier: 6.674A pdb=" N LEU A 46 " --> pdb=" O LEU A 247 " (cutoff:3.500A) removed outlier: 8.730A pdb=" N SER A 249 " --> pdb=" O LEU A 46 " (cutoff:3.500A) removed outlier: 6.390A pdb=" N ILE A 287 " --> pdb=" O TYR A 344 " (cutoff:3.500A) removed outlier: 8.064A pdb=" N HIS A 346 " --> pdb=" O ILE A 287 " (cutoff:3.500A) removed outlier: 6.714A pdb=" N LEU A 289 " --> pdb=" O HIS A 346 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'B' and resid 47 through 52 removed outlier: 4.116A pdb=" N PHE B 335 " --> pdb=" O SER B 331 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'B' and resid 58 through 63 removed outlier: 3.717A pdb=" N LYS B 78 " --> pdb=" O SER B 74 " (cutoff:3.500A) removed outlier: 6.299A pdb=" N LEU B 79 " --> pdb=" O ALA B 92 " (cutoff:3.500A) removed outlier: 4.518A pdb=" N ALA B 92 " --> pdb=" O LEU B 79 " (cutoff:3.500A) removed outlier: 6.421A pdb=" N ILE B 81 " --> pdb=" O VAL B 90 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'B' and resid 100 through 105 removed outlier: 3.621A pdb=" N GLY B 115 " --> pdb=" O THR B 102 " (cutoff:3.500A) removed outlier: 6.008A pdb=" N ASN B 125 " --> pdb=" O VAL B 135 " (cutoff:3.500A) removed outlier: 6.025A pdb=" N VAL B 135 " --> pdb=" O ASN B 125 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'B' and resid 146 through 153 removed outlier: 3.656A pdb=" N CYS B 148 " --> pdb=" O SER B 160 " (cutoff:3.500A) removed outlier: 6.720A pdb=" N GLN B 156 " --> pdb=" O LEU B 152 " (cutoff:3.500A) removed outlier: 3.744A pdb=" N THR B 165 " --> pdb=" O SER B 161 " (cutoff:3.500A) removed outlier: 6.365A pdb=" N ASP B 170 " --> pdb=" O GLN B 176 " (cutoff:3.500A) removed outlier: 6.070A pdb=" N GLN B 176 " --> pdb=" O ASP B 170 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'B' and resid 189 through 192 removed outlier: 6.610A pdb=" N ASP B 212 " --> pdb=" O CYS B 218 " (cutoff:3.500A) removed outlier: 6.335A pdb=" N CYS B 218 " --> pdb=" O ASP B 212 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 229 through 234 removed outlier: 6.416A pdb=" N CYS B 250 " --> pdb=" O THR B 263 " (cutoff:3.500A) removed outlier: 4.351A pdb=" N THR B 263 " --> pdb=" O CYS B 250 " (cutoff:3.500A) removed outlier: 6.401A pdb=" N LEU B 252 " --> pdb=" O LEU B 261 " (cutoff:3.500A) removed outlier: 3.553A pdb=" N ASP B 254 " --> pdb=" O GLN B 259 " (cutoff:3.500A) removed outlier: 3.641A pdb=" N GLN B 259 " --> pdb=" O ASP B 254 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'B' and resid 273 through 278 removed outlier: 5.948A pdb=" N ASP B 298 " --> pdb=" O ARG B 304 " (cutoff:3.500A) removed outlier: 6.138A pdb=" N ARG B 304 " --> pdb=" O ASP B 298 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'C' and resid 3 through 7 removed outlier: 5.974A pdb=" N GLY C 16 " --> pdb=" O SER C 85 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'C' and resid 10 through 12 removed outlier: 6.322A pdb=" N GLY C 10 " --> pdb=" O THR C 118 " (cutoff:3.500A) removed outlier: 6.809A pdb=" N MET C 34 " --> pdb=" O TYR C 50 " (cutoff:3.500A) removed outlier: 4.496A pdb=" N TYR C 50 " --> pdb=" O MET C 34 " (cutoff:3.500A) removed outlier: 6.351A pdb=" N TRP C 36 " --> pdb=" O VAL C 48 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'C' and resid 10 through 12 removed outlier: 6.322A pdb=" N GLY C 10 " --> pdb=" O THR C 118 " (cutoff:3.500A) removed outlier: 3.963A pdb=" N PHE C 110 " --> pdb=" O ARG C 98 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'C' and resid 142 through 143 Processing sheet with id=AB4, first strand: chain 'C' and resid 148 through 150 removed outlier: 3.679A pdb=" N GLU C 248 " --> pdb=" O VAL C 149 " (cutoff:3.500A) removed outlier: 6.788A pdb=" N LEU C 176 " --> pdb=" O TYR C 192 " (cutoff:3.500A) removed outlier: 5.008A pdb=" N TYR C 192 " --> pdb=" O LEU C 176 " (cutoff:3.500A) removed outlier: 6.490A pdb=" N TRP C 178 " --> pdb=" O LEU C 190 " (cutoff:3.500A) 445 hydrogen bonds defined for protein. 1233 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 0 stacking parallelities Total time for adding SS restraints: 3.03 Time building geometry restraints manager: 3.58 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.34: 2779 1.34 - 1.46: 2061 1.46 - 1.58: 3983 1.58 - 1.70: 0 1.70 - 1.82: 82 Bond restraints: 8905 Sorted by residual: bond pdb=" C ARG R 624 " pdb=" N ILE R 625 " ideal model delta sigma weight residual 1.334 1.287 0.047 1.25e-02 6.40e+03 1.42e+01 bond pdb=" N ASP C 139 " pdb=" CA ASP C 139 " ideal model delta sigma weight residual 1.456 1.488 -0.032 1.23e-02 6.61e+03 6.95e+00 bond pdb=" N ASN R 816 " pdb=" CA ASN R 816 " ideal model delta sigma weight residual 1.457 1.488 -0.030 1.29e-02 6.01e+03 5.57e+00 bond pdb=" C TYR C 192 " pdb=" N ARG C 193 " ideal model delta sigma weight residual 1.333 1.299 0.034 1.45e-02 4.76e+03 5.47e+00 bond pdb=" CA SER R 814 " pdb=" CB SER R 814 " ideal model delta sigma weight residual 1.530 1.496 0.034 1.57e-02 4.06e+03 4.70e+00 ... (remaining 8900 not shown) Histogram of bond angle deviations from ideal: 98.32 - 105.48: 138 105.48 - 112.64: 4759 112.64 - 119.80: 2864 119.80 - 126.95: 4188 126.95 - 134.11: 112 Bond angle restraints: 12061 Sorted by residual: angle pdb=" N GLY C 56 " pdb=" CA GLY C 56 " pdb=" C GLY C 56 " ideal model delta sigma weight residual 110.21 116.71 -6.50 9.10e-01 1.21e+00 5.11e+01 angle pdb=" N ALA B 203 " pdb=" CA ALA B 203 " pdb=" C ALA B 203 " ideal model delta sigma weight residual 110.41 117.45 -7.04 1.18e+00 7.18e-01 3.56e+01 angle pdb=" N GLU A 236 " pdb=" CA GLU A 236 " pdb=" C GLU A 236 " ideal model delta sigma weight residual 114.75 108.66 6.09 1.26e+00 6.30e-01 2.34e+01 angle pdb=" N LYS A 232 " pdb=" CA LYS A 232 " pdb=" C LYS A 232 " ideal model delta sigma weight residual 112.72 118.75 -6.03 1.28e+00 6.10e-01 2.22e+01 angle pdb=" C ASP B 290 " pdb=" CA ASP B 290 " pdb=" CB ASP B 290 " ideal model delta sigma weight residual 110.56 102.61 7.95 1.69e+00 3.50e-01 2.22e+01 ... (remaining 12056 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.62: 4891 17.62 - 35.24: 328 35.24 - 52.86: 41 52.86 - 70.49: 3 70.49 - 88.11: 5 Dihedral angle restraints: 5268 sinusoidal: 2057 harmonic: 3211 Sorted by residual: dihedral pdb=" CA ILE A 351 " pdb=" C ILE A 351 " pdb=" N ASN A 352 " pdb=" CA ASN A 352 " ideal model delta harmonic sigma weight residual 180.00 160.02 19.98 0 5.00e+00 4.00e-02 1.60e+01 dihedral pdb=" CA GLU A 15 " pdb=" C GLU A 15 " pdb=" N ASP A 16 " pdb=" CA ASP A 16 " ideal model delta harmonic sigma weight residual 180.00 160.05 19.95 0 5.00e+00 4.00e-02 1.59e+01 dihedral pdb=" CA PHE B 234 " pdb=" C PHE B 234 " pdb=" N PHE B 235 " pdb=" CA PHE B 235 " ideal model delta harmonic sigma weight residual 180.00 160.32 19.68 0 5.00e+00 4.00e-02 1.55e+01 ... (remaining 5265 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.146: 1372 0.146 - 0.292: 17 0.292 - 0.438: 2 0.438 - 0.583: 0 0.583 - 0.729: 1 Chirality restraints: 1392 Sorted by residual: chirality pdb=" CG LEU C 176 " pdb=" CB LEU C 176 " pdb=" CD1 LEU C 176 " pdb=" CD2 LEU C 176 " both_signs ideal model delta sigma weight residual False -2.59 -1.86 -0.73 2.00e-01 2.50e+01 1.33e+01 chirality pdb=" CA CYS R 657 " pdb=" N CYS R 657 " pdb=" C CYS R 657 " pdb=" CB CYS R 657 " both_signs ideal model delta sigma weight residual False 2.51 2.18 0.33 2.00e-01 2.50e+01 2.77e+00 chirality pdb=" CA CYS B 204 " pdb=" N CYS B 204 " pdb=" C CYS B 204 " pdb=" CB CYS B 204 " both_signs ideal model delta sigma weight residual False 2.51 2.18 0.33 2.00e-01 2.50e+01 2.66e+00 ... (remaining 1389 not shown) Planarity restraints: 1506 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA ALA A 14 " -0.016 2.00e-02 2.50e+03 3.20e-02 1.02e+01 pdb=" C ALA A 14 " 0.055 2.00e-02 2.50e+03 pdb=" O ALA A 14 " -0.021 2.00e-02 2.50e+03 pdb=" N GLU A 15 " -0.019 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE C 214 " -0.029 2.00e-02 2.50e+03 1.90e-02 6.33e+00 pdb=" CG PHE C 214 " 0.036 2.00e-02 2.50e+03 pdb=" CD1 PHE C 214 " 0.008 2.00e-02 2.50e+03 pdb=" CD2 PHE C 214 " 0.006 2.00e-02 2.50e+03 pdb=" CE1 PHE C 214 " -0.003 2.00e-02 2.50e+03 pdb=" CE2 PHE C 214 " -0.001 2.00e-02 2.50e+03 pdb=" CZ PHE C 214 " -0.016 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C PHE B 235 " -0.040 5.00e-02 4.00e+02 6.10e-02 5.96e+00 pdb=" N PRO B 236 " 0.105 5.00e-02 4.00e+02 pdb=" CA PRO B 236 " -0.031 5.00e-02 4.00e+02 pdb=" CD PRO B 236 " -0.034 5.00e-02 4.00e+02 ... (remaining 1503 not shown) Histogram of nonbonded interaction distances: 1.69 - 2.33: 15 2.33 - 2.97: 4498 2.97 - 3.62: 13150 3.62 - 4.26: 20622 4.26 - 4.90: 34950 Nonbonded interactions: 73235 Sorted by model distance: nonbonded pdb=" CB LEU C 176 " pdb=" SD MET C 194 " model vdw 1.690 3.800 nonbonded pdb=" OG1 THR B 274 " pdb=" O VAL B 315 " model vdw 2.089 2.440 nonbonded pdb=" O ILE B 58 " pdb=" OG SER B 316 " model vdw 2.147 2.440 nonbonded pdb=" OE2 GLU A 236 " pdb=" NH1 ARG B 314 " model vdw 2.157 2.520 nonbonded pdb=" OG SER B 281 " pdb=" OD2 ASP D 48 " model vdw 2.201 2.440 ... (remaining 73230 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== 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=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic 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: 1.160 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.010 Extract box with map and model: 3.830 Check model and map are aligned: 0.110 Set scattering table: 0.080 Process input model: 26.460 Find NCS groups from input model: 0.260 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.000 Load rotamer database and sin/cos tables:11.200 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 43.150 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6920 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.007 0.064 8905 Z= 0.478 Angle : 0.869 10.467 12061 Z= 0.501 Chirality : 0.058 0.729 1392 Planarity : 0.005 0.061 1506 Dihedral : 11.976 88.107 3184 Min Nonbonded Distance : 1.690 Molprobity Statistics. All-atom Clashscore : 7.73 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.24 % Favored : 93.76 % Rotamer: Outliers : 0.00 % Allowed : 2.39 % Favored : 97.61 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: -1.00 (0.25), residues: 1090 helix: 0.90 (0.28), residues: 368 sheet: -1.03 (0.29), residues: 276 loop : -1.71 (0.28), residues: 446 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 184 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 184 time to evaluate : 0.892 Fit side-chains revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 184 average time/residue: 1.1477 time to fit residues: 226.5788 Evaluate side-chains 162 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 162 time to evaluate : 0.995 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.4117 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 108 random chunks: chunk 91 optimal weight: 0.3980 chunk 81 optimal weight: 0.6980 chunk 45 optimal weight: 0.0770 chunk 27 optimal weight: 0.9990 chunk 55 optimal weight: 0.5980 chunk 43 optimal weight: 0.0970 chunk 84 optimal weight: 0.5980 chunk 32 optimal weight: 0.6980 chunk 51 optimal weight: 0.6980 chunk 63 optimal weight: 0.6980 chunk 98 optimal weight: 0.8980 overall best weight: 0.3536 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 318 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** B 75 GLN ** B 220 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** C 39 GLN C 173 ASN C 181 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6872 moved from start: 0.1281 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.067 8905 Z= 0.211 Angle : 0.625 10.065 12061 Z= 0.323 Chirality : 0.043 0.154 1392 Planarity : 0.004 0.043 1506 Dihedral : 5.098 24.233 1189 Min Nonbonded Distance : 2.133 Molprobity Statistics. All-atom Clashscore : 11.62 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.50 % Favored : 95.50 % Rotamer: Outliers : 1.87 % Allowed : 11.20 % Favored : 86.93 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: -0.29 (0.26), residues: 1090 helix: 1.84 (0.28), residues: 362 sheet: -0.93 (0.28), residues: 281 loop : -1.38 (0.30), residues: 447 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 196 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 178 time to evaluate : 0.979 Fit side-chains outliers start: 18 outliers final: 9 residues processed: 186 average time/residue: 1.0451 time to fit residues: 209.3895 Evaluate side-chains 177 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 168 time to evaluate : 1.049 Switching outliers to nearest non-outliers outliers start: 9 outliers final: 3 residues processed: 6 average time/residue: 0.1024 time to fit residues: 2.4775 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 108 random chunks: chunk 54 optimal weight: 3.9990 chunk 30 optimal weight: 1.9990 chunk 81 optimal weight: 1.9990 chunk 66 optimal weight: 2.9990 chunk 27 optimal weight: 1.9990 chunk 98 optimal weight: 2.9990 chunk 106 optimal weight: 1.9990 chunk 87 optimal weight: 0.9990 chunk 97 optimal weight: 0.6980 chunk 33 optimal weight: 1.9990 chunk 78 optimal weight: 4.9990 overall best weight: 1.5388 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 155 ASN ** B 220 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 293 ASN ** B 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** C 173 ASN ** C 233 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6970 moved from start: 0.1478 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.008 0.086 8905 Z= 0.526 Angle : 0.741 14.670 12061 Z= 0.382 Chirality : 0.049 0.471 1392 Planarity : 0.005 0.039 1506 Dihedral : 5.519 30.345 1189 Min Nonbonded Distance : 2.075 Molprobity Statistics. All-atom Clashscore : 15.28 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.06 % Favored : 93.94 % Rotamer: Outliers : 4.25 % Allowed : 13.28 % Favored : 82.47 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: -0.62 (0.26), residues: 1090 helix: 1.45 (0.27), residues: 365 sheet: -1.10 (0.28), residues: 288 loop : -1.48 (0.30), residues: 437 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 220 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 41 poor density : 179 time to evaluate : 0.988 Fit side-chains outliers start: 41 outliers final: 26 residues processed: 195 average time/residue: 1.0149 time to fit residues: 213.6084 Evaluate side-chains 205 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 26 poor density : 179 time to evaluate : 0.935 Switching outliers to nearest non-outliers outliers start: 26 outliers final: 8 residues processed: 18 average time/residue: 0.3638 time to fit residues: 9.2465 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 108 random chunks: chunk 97 optimal weight: 2.9990 chunk 73 optimal weight: 0.9980 chunk 50 optimal weight: 1.9990 chunk 10 optimal weight: 0.7980 chunk 46 optimal weight: 2.9990 chunk 65 optimal weight: 0.4980 chunk 98 optimal weight: 1.9990 chunk 104 optimal weight: 0.0570 chunk 51 optimal weight: 0.9980 chunk 93 optimal weight: 0.5980 chunk 28 optimal weight: 0.5980 overall best weight: 0.5098 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 175 GLN ** B 220 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 293 ASN ** B 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 144 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 173 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6892 moved from start: 0.1834 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.058 8905 Z= 0.237 Angle : 0.606 11.529 12061 Z= 0.313 Chirality : 0.044 0.417 1392 Planarity : 0.004 0.040 1506 Dihedral : 4.983 25.306 1189 Min Nonbonded Distance : 2.104 Molprobity Statistics. All-atom Clashscore : 12.31 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.59 % Favored : 95.41 % Rotamer: Outliers : 3.73 % Allowed : 15.04 % Favored : 81.22 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: -0.22 (0.26), residues: 1090 helix: 1.88 (0.28), residues: 366 sheet: -0.86 (0.29), residues: 265 loop : -1.37 (0.29), residues: 459 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 219 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 36 poor density : 183 time to evaluate : 1.054 Fit side-chains outliers start: 36 outliers final: 18 residues processed: 200 average time/residue: 0.9963 time to fit residues: 215.5301 Evaluate side-chains 191 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 173 time to evaluate : 0.911 Switching outliers to nearest non-outliers outliers start: 18 outliers final: 10 residues processed: 9 average time/residue: 0.4934 time to fit residues: 6.4646 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 108 random chunks: chunk 86 optimal weight: 0.9980 chunk 59 optimal weight: 0.0060 chunk 1 optimal weight: 0.4980 chunk 77 optimal weight: 2.9990 chunk 43 optimal weight: 0.7980 chunk 89 optimal weight: 1.9990 chunk 72 optimal weight: 0.9980 chunk 0 optimal weight: 1.9990 chunk 53 optimal weight: 0.9980 chunk 93 optimal weight: 1.9990 chunk 26 optimal weight: 0.6980 overall best weight: 0.5996 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 155 ASN ** B 220 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 239 ASN ** B 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 144 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6895 moved from start: 0.1922 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.058 8905 Z= 0.259 Angle : 0.619 11.405 12061 Z= 0.318 Chirality : 0.046 0.430 1392 Planarity : 0.004 0.036 1506 Dihedral : 4.949 29.288 1189 Min Nonbonded Distance : 2.113 Molprobity Statistics. All-atom Clashscore : 12.76 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.95 % Favored : 95.05 % Rotamer: Outliers : 3.63 % Allowed : 16.08 % Favored : 80.29 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: -0.14 (0.26), residues: 1090 helix: 1.98 (0.28), residues: 367 sheet: -0.84 (0.29), residues: 265 loop : -1.35 (0.29), residues: 458 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 213 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 35 poor density : 178 time to evaluate : 0.856 Fit side-chains outliers start: 35 outliers final: 20 residues processed: 195 average time/residue: 0.9954 time to fit residues: 209.4408 Evaluate side-chains 195 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 20 poor density : 175 time to evaluate : 0.936 Switching outliers to nearest non-outliers outliers start: 20 outliers final: 10 residues processed: 11 average time/residue: 0.2948 time to fit residues: 5.3238 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 108 random chunks: chunk 35 optimal weight: 0.6980 chunk 94 optimal weight: 0.6980 chunk 20 optimal weight: 4.9990 chunk 61 optimal weight: 0.8980 chunk 25 optimal weight: 0.6980 chunk 104 optimal weight: 0.5980 chunk 86 optimal weight: 0.9990 chunk 48 optimal weight: 0.3980 chunk 8 optimal weight: 0.9990 chunk 34 optimal weight: 0.5980 chunk 54 optimal weight: 0.9990 overall best weight: 0.5980 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 220 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 144 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6891 moved from start: 0.2081 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.077 8905 Z= 0.262 Angle : 0.604 11.123 12061 Z= 0.312 Chirality : 0.045 0.495 1392 Planarity : 0.004 0.038 1506 Dihedral : 4.886 27.260 1189 Min Nonbonded Distance : 2.124 Molprobity Statistics. All-atom Clashscore : 12.88 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.59 % Favored : 95.41 % Rotamer: Outliers : 3.11 % Allowed : 17.01 % Favored : 79.88 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: -0.10 (0.26), residues: 1090 helix: 2.01 (0.27), residues: 369 sheet: -0.87 (0.29), residues: 268 loop : -1.29 (0.30), residues: 453 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 212 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 30 poor density : 182 time to evaluate : 0.893 Fit side-chains outliers start: 30 outliers final: 16 residues processed: 196 average time/residue: 1.0007 time to fit residues: 211.7721 Evaluate side-chains 191 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 16 poor density : 175 time to evaluate : 0.867 Switching outliers to nearest non-outliers outliers start: 16 outliers final: 9 residues processed: 7 average time/residue: 0.1297 time to fit residues: 2.7545 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 108 random chunks: chunk 100 optimal weight: 0.9990 chunk 11 optimal weight: 0.5980 chunk 59 optimal weight: 0.3980 chunk 76 optimal weight: 3.9990 chunk 87 optimal weight: 0.9990 chunk 58 optimal weight: 0.0670 chunk 104 optimal weight: 0.9990 chunk 65 optimal weight: 0.6980 chunk 63 optimal weight: 0.7980 chunk 48 optimal weight: 0.6980 chunk 64 optimal weight: 2.9990 overall best weight: 0.4918 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 220 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6877 moved from start: 0.2194 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.050 8905 Z= 0.228 Angle : 0.609 11.150 12061 Z= 0.310 Chirality : 0.043 0.328 1392 Planarity : 0.004 0.042 1506 Dihedral : 4.770 25.320 1189 Min Nonbonded Distance : 2.129 Molprobity Statistics. All-atom Clashscore : 12.65 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.14 % Favored : 94.86 % Rotamer: Outliers : 3.01 % Allowed : 18.15 % Favored : 78.84 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: -0.02 (0.26), residues: 1090 helix: 2.10 (0.27), residues: 368 sheet: -0.82 (0.29), residues: 270 loop : -1.27 (0.30), residues: 452 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 199 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 29 poor density : 170 time to evaluate : 0.954 Fit side-chains outliers start: 29 outliers final: 17 residues processed: 186 average time/residue: 0.9503 time to fit residues: 191.5878 Evaluate side-chains 183 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 17 poor density : 166 time to evaluate : 0.870 Switching outliers to nearest non-outliers outliers start: 17 outliers final: 10 residues processed: 7 average time/residue: 0.4806 time to fit residues: 5.1055 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 108 random chunks: chunk 41 optimal weight: 0.0060 chunk 62 optimal weight: 1.9990 chunk 31 optimal weight: 0.7980 chunk 20 optimal weight: 4.9990 chunk 66 optimal weight: 2.9990 chunk 70 optimal weight: 2.9990 chunk 51 optimal weight: 0.9980 chunk 9 optimal weight: 1.9990 chunk 81 optimal weight: 0.9980 chunk 94 optimal weight: 0.6980 chunk 99 optimal weight: 0.4980 overall best weight: 0.5996 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 220 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6883 moved from start: 0.2291 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.050 8905 Z= 0.257 Angle : 0.616 11.133 12061 Z= 0.314 Chirality : 0.045 0.549 1392 Planarity : 0.004 0.037 1506 Dihedral : 4.783 25.286 1189 Min Nonbonded Distance : 2.120 Molprobity Statistics. All-atom Clashscore : 12.93 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.14 % Favored : 94.86 % Rotamer: Outliers : 2.90 % Allowed : 19.19 % Favored : 77.90 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: 0.01 (0.26), residues: 1090 helix: 2.18 (0.27), residues: 365 sheet: -0.84 (0.29), residues: 268 loop : -1.25 (0.30), residues: 457 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 198 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 28 poor density : 170 time to evaluate : 1.048 Fit side-chains outliers start: 28 outliers final: 17 residues processed: 186 average time/residue: 0.9374 time to fit residues: 189.3402 Evaluate side-chains 191 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 17 poor density : 174 time to evaluate : 0.925 Switching outliers to nearest non-outliers outliers start: 17 outliers final: 10 residues processed: 7 average time/residue: 0.2666 time to fit residues: 3.7001 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 108 random chunks: chunk 90 optimal weight: 0.5980 chunk 97 optimal weight: 0.4980 chunk 99 optimal weight: 0.2980 chunk 58 optimal weight: 0.7980 chunk 42 optimal weight: 2.9990 chunk 76 optimal weight: 0.9980 chunk 29 optimal weight: 0.6980 chunk 87 optimal weight: 0.6980 chunk 91 optimal weight: 0.0980 chunk 96 optimal weight: 0.9980 chunk 63 optimal weight: 0.7980 overall best weight: 0.4380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6865 moved from start: 0.2413 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.049 8905 Z= 0.215 Angle : 0.602 10.322 12061 Z= 0.306 Chirality : 0.044 0.442 1392 Planarity : 0.003 0.035 1506 Dihedral : 4.683 25.310 1189 Min Nonbonded Distance : 2.129 Molprobity Statistics. All-atom Clashscore : 12.08 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.77 % Favored : 95.23 % Rotamer: Outliers : 2.28 % Allowed : 20.02 % Favored : 77.70 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: 0.06 (0.26), residues: 1090 helix: 2.21 (0.27), residues: 365 sheet: -0.84 (0.29), residues: 268 loop : -1.18 (0.30), residues: 457 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 198 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 22 poor density : 176 time to evaluate : 0.958 Fit side-chains outliers start: 22 outliers final: 14 residues processed: 190 average time/residue: 1.0376 time to fit residues: 212.4801 Evaluate side-chains 187 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 173 time to evaluate : 0.880 Switching outliers to nearest non-outliers outliers start: 14 outliers final: 11 residues processed: 3 average time/residue: 0.6924 time to fit residues: 3.5730 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 108 random chunks: chunk 102 optimal weight: 3.9990 chunk 62 optimal weight: 1.9990 chunk 48 optimal weight: 0.9980 chunk 71 optimal weight: 0.0980 chunk 107 optimal weight: 0.7980 chunk 99 optimal weight: 0.2980 chunk 85 optimal weight: 0.7980 chunk 8 optimal weight: 0.7980 chunk 66 optimal weight: 0.6980 chunk 52 optimal weight: 0.5980 chunk 68 optimal weight: 0.0670 overall best weight: 0.3518 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6844 moved from start: 0.2559 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.045 8905 Z= 0.197 Angle : 0.613 11.458 12061 Z= 0.310 Chirality : 0.043 0.376 1392 Planarity : 0.004 0.057 1506 Dihedral : 4.606 23.762 1189 Min Nonbonded Distance : 2.130 Molprobity Statistics. All-atom Clashscore : 12.08 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.95 % Favored : 95.05 % Rotamer: Outliers : 1.76 % Allowed : 21.16 % Favored : 77.07 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: 0.20 (0.26), residues: 1090 helix: 2.39 (0.27), residues: 366 sheet: -0.80 (0.28), residues: 270 loop : -1.15 (0.30), residues: 454 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2180 Ramachandran restraints generated. 1090 Oldfield, 0 Emsley, 1090 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 188 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 17 poor density : 171 time to evaluate : 0.936 Fit side-chains outliers start: 17 outliers final: 14 residues processed: 184 average time/residue: 1.0623 time to fit residues: 210.4593 Evaluate side-chains 181 residues out of total 964 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 167 time to evaluate : 0.856 Switching outliers to nearest non-outliers outliers start: 14 outliers final: 12 residues processed: 2 average time/residue: 1.0913 time to fit residues: 3.5642 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 108 random chunks: chunk 91 optimal weight: 0.3980 chunk 26 optimal weight: 0.9980 chunk 79 optimal weight: 0.6980 chunk 12 optimal weight: 0.9990 chunk 23 optimal weight: 1.9990 chunk 85 optimal weight: 2.9990 chunk 35 optimal weight: 1.9990 chunk 88 optimal weight: 0.5980 chunk 10 optimal weight: 0.9990 chunk 15 optimal weight: 4.9990 chunk 75 optimal weight: 0.0870 overall best weight: 0.5558 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** R 695 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4119 r_free = 0.4119 target = 0.163092 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 40)----------------| | r_work = 0.3851 r_free = 0.3851 target = 0.141351 restraints weight = 12565.662| |-----------------------------------------------------------------------------| r_work (start): 0.3852 rms_B_bonded: 2.10 r_work: 0.3740 rms_B_bonded: 2.76 restraints_weight: 0.5000 r_work: 0.3588 rms_B_bonded: 4.78 restraints_weight: 0.2500 r_work (final): 0.3588 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7237 moved from start: 0.2540 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.051 8905 Z= 0.251 Angle : 0.635 10.540 12061 Z= 0.321 Chirality : 0.044 0.363 1392 Planarity : 0.004 0.058 1506 Dihedral : 4.693 24.140 1189 Min Nonbonded Distance : 2.125 Molprobity Statistics. All-atom Clashscore : 12.48 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.32 % Favored : 94.68 % Rotamer: Outliers : 1.87 % Allowed : 21.78 % Favored : 76.35 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.12 % 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: 0.15 (0.26), residues: 1090 helix: 2.38 (0.27), residues: 366 sheet: -0.82 (0.28), residues: 270 loop : -1.23 (0.30), residues: 454 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 3836.13 seconds wall clock time: 68 minutes 59.52 seconds (4139.52 seconds total)