Starting phenix.real_space_refine on Mon Dec 11 09:13:59 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h1t_34431/12_2023/8h1t_34431.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h1t_34431/12_2023/8h1t_34431.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.0 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h1t_34431/12_2023/8h1t_34431.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h1t_34431/12_2023/8h1t_34431.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h1t_34431/12_2023/8h1t_34431.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h1t_34431/12_2023/8h1t_34431.pdb" } resolution = 3.0 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.013 sd= 1.080 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: 5 Type Number sf(0) Gaussians P 334 5.49 5 S 28 5.16 5 C 9881 2.51 5 N 3214 2.21 5 O 3908 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped None Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/chem_data/mon_lib" Total number of atoms: 17365 Number of models: 1 Model: "" Number of chains: 14 Chain: "A" Number of atoms: 790 Number of conformers: 1 Conformer: "" Number of residues, atoms: 96, 790 Classifications: {'peptide': 96} Link IDs: {'PTRANS': 3, 'TRANS': 92} Chain: "B" Number of atoms: 641 Number of conformers: 1 Conformer: "" Number of residues, atoms: 80, 641 Classifications: {'peptide': 80} Link IDs: {'PTRANS': 1, 'TRANS': 78} Chain: "C" Number of atoms: 858 Number of conformers: 1 Conformer: "" Number of residues, atoms: 111, 858 Classifications: {'peptide': 111} Link IDs: {'PTRANS': 5, 'TRANS': 105} Chain: "D" Number of atoms: 737 Number of conformers: 1 Conformer: "" Number of residues, atoms: 94, 737 Classifications: {'peptide': 94} Link IDs: {'PTRANS': 2, 'TRANS': 91} Chain: "E" Number of atoms: 807 Number of conformers: 1 Conformer: "" Number of residues, atoms: 98, 807 Classifications: {'peptide': 98} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 3, 'TRANS': 94} Chain: "F" Number of atoms: 653 Number of conformers: 1 Conformer: "" Number of residues, atoms: 82, 653 Classifications: {'peptide': 82} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 1, 'TRANS': 80} Chain: "G" Number of atoms: 838 Number of conformers: 1 Conformer: "" Number of residues, atoms: 109, 838 Classifications: {'peptide': 109} Link IDs: {'PTRANS': 5, 'TRANS': 103} Chain: "H" Number of atoms: 737 Number of conformers: 1 Conformer: "" Number of residues, atoms: 94, 737 Classifications: {'peptide': 94} Link IDs: {'PTRANS': 2, 'TRANS': 91} Chain: "I" Number of atoms: 3398 Number of conformers: 1 Conformer: "" Number of residues, atoms: 167, 3398 Classifications: {'DNA': 167} Link IDs: {'rna3p': 166} Chain: "J" Number of atoms: 3449 Number of conformers: 1 Conformer: "" Number of residues, atoms: 167, 3449 Classifications: {'DNA': 167} Link IDs: {'rna3p': 166} Chain: "K" Number of atoms: 520 Number of conformers: 1 Conformer: "" Number of residues, atoms: 73, 520 Classifications: {'peptide': 73} Link IDs: {'PTRANS': 2, 'TRANS': 70} Chain: "L" Number of atoms: 2433 Number of conformers: 1 Conformer: "" Number of residues, atoms: 301, 2433 Classifications: {'peptide': 301} Link IDs: {'PTRANS': 15, 'TRANS': 285} Chain breaks: 2 Chain: "M" Number of atoms: 601 Number of conformers: 1 Conformer: "" Number of residues, atoms: 76, 601 Classifications: {'peptide': 76} Link IDs: {'PTRANS': 3, 'TRANS': 72} Chain: "N" Number of atoms: 903 Number of conformers: 1 Conformer: "" Number of residues, atoms: 109, 903 Classifications: {'peptide': 109} Link IDs: {'PTRANS': 3, 'TRANS': 105} Time building chain proxies: 9.32, per 1000 atoms: 0.54 Number of scatterers: 17365 At special positions: 0 Unit cell: (121, 118, 159, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) S 28 16.00 P 334 15.00 O 3908 8.00 N 3214 7.00 C 9881 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.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: 9.17 Conformation dependent library (CDL) restraints added in 2.1 seconds 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2468 Finding SS restraints... Secondary structure from input PDB file: 58 helices and 14 sheets defined 61.6% alpha, 5.7% beta 166 base pairs and 273 stacking pairs defined. Time for finding SS restraints: 7.92 Creating SS restraints... Processing helix chain 'A' and resid 44 through 56 Processing helix chain 'A' and resid 63 through 79 removed outlier: 3.982A pdb=" N ASP A 77 " --> pdb=" O GLU A 73 " (cutoff:3.500A) removed outlier: 4.115A pdb=" N PHE A 78 " --> pdb=" O ILE A 74 " (cutoff:3.500A) Processing helix chain 'A' and resid 85 through 114 Processing helix chain 'A' and resid 120 through 132 Processing helix chain 'B' and resid 25 through 29 Processing helix chain 'B' and resid 30 through 42 Processing helix chain 'B' and resid 48 through 76 removed outlier: 4.078A pdb=" N GLU B 52 " --> pdb=" O GLY B 48 " (cutoff:3.500A) removed outlier: 3.808A pdb=" N ARG B 55 " --> pdb=" O TYR B 51 " (cutoff:3.500A) removed outlier: 3.911A pdb=" N GLU B 63 " --> pdb=" O LYS B 59 " (cutoff:3.500A) removed outlier: 3.893A pdb=" N ARG B 67 " --> pdb=" O GLU B 63 " (cutoff:3.500A) removed outlier: 4.667A pdb=" N ASP B 68 " --> pdb=" O ASN B 64 " (cutoff:3.500A) Processing helix chain 'B' and resid 82 through 93 removed outlier: 3.538A pdb=" N GLN B 93 " --> pdb=" O ALA B 89 " (cutoff:3.500A) Processing helix chain 'C' and resid 16 through 21 Processing helix chain 'C' and resid 26 through 38 Processing helix chain 'C' and resid 45 through 73 removed outlier: 4.471A pdb=" N VAL C 49 " --> pdb=" O ALA C 45 " (cutoff:3.500A) Processing helix chain 'C' and resid 79 through 90 Processing helix chain 'C' and resid 90 through 97 Processing helix chain 'D' and resid 37 through 49 Processing helix chain 'D' and resid 58 through 84 Processing helix chain 'D' and resid 90 through 102 Processing helix chain 'D' and resid 104 through 124 removed outlier: 3.925A pdb=" N LYS D 108 " --> pdb=" O GLY D 104 " (cutoff:3.500A) Processing helix chain 'E' and resid 45 through 56 Processing helix chain 'E' and resid 63 through 75 Processing helix chain 'E' and resid 85 through 114 Processing helix chain 'E' and resid 120 through 131 removed outlier: 3.516A pdb=" N ILE E 124 " --> pdb=" O MET E 120 " (cutoff:3.500A) Processing helix chain 'F' and resid 25 through 29 Processing helix chain 'F' and resid 30 through 41 Processing helix chain 'F' and resid 49 through 76 removed outlier: 3.651A pdb=" N GLU F 53 " --> pdb=" O LEU F 49 " (cutoff:3.500A) Processing helix chain 'F' and resid 82 through 93 removed outlier: 3.512A pdb=" N GLN F 93 " --> pdb=" O ALA F 89 " (cutoff:3.500A) Processing helix chain 'G' and resid 16 through 22 Processing helix chain 'G' and resid 27 through 37 Processing helix chain 'G' and resid 45 through 73 removed outlier: 4.347A pdb=" N VAL G 49 " --> pdb=" O ALA G 45 " (cutoff:3.500A) Processing helix chain 'G' and resid 79 through 90 removed outlier: 3.526A pdb=" N ASP G 90 " --> pdb=" O ALA G 86 " (cutoff:3.500A) Processing helix chain 'G' and resid 90 through 97 Processing helix chain 'H' and resid 37 through 49 Processing helix chain 'H' and resid 55 through 84 Processing helix chain 'H' and resid 90 through 102 removed outlier: 3.559A pdb=" N ILE H 94 " --> pdb=" O THR H 90 " (cutoff:3.500A) Processing helix chain 'H' and resid 103 through 124 removed outlier: 3.658A pdb=" N ALA H 107 " --> pdb=" O PRO H 103 " (cutoff:3.500A) Processing helix chain 'K' and resid 38 through 49 removed outlier: 4.154A pdb=" N LEU K 42 " --> pdb=" O PRO K 38 " (cutoff:3.500A) Processing helix chain 'K' and resid 58 through 69 Processing helix chain 'K' and resid 75 through 85 Processing helix chain 'L' and resid 11 through 23 Processing helix chain 'L' and resid 54 through 60 Processing helix chain 'L' and resid 74 through 80 Processing helix chain 'L' and resid 90 through 102 Processing helix chain 'L' and resid 109 through 120 Processing helix chain 'L' and resid 123 through 133 Processing helix chain 'L' and resid 134 through 144 removed outlier: 3.616A pdb=" N HIS L 144 " --> pdb=" O ALA L 140 " (cutoff:3.500A) Processing helix chain 'L' and resid 147 through 152 Processing helix chain 'L' and resid 202 through 218 Processing helix chain 'L' and resid 236 through 243 removed outlier: 3.990A pdb=" N LYS L 240 " --> pdb=" O ASP L 236 " (cutoff:3.500A) Processing helix chain 'L' and resid 648 through 668 removed outlier: 3.545A pdb=" N THR L 668 " --> pdb=" O ASP L 664 " (cutoff:3.500A) Processing helix chain 'L' and resid 671 through 686 Processing helix chain 'L' and resid 687 through 695 removed outlier: 3.625A pdb=" N LEU L 691 " --> pdb=" O MET L 687 " (cutoff:3.500A) removed outlier: 3.707A pdb=" N GLU L 693 " --> pdb=" O ALA L 689 " (cutoff:3.500A) Processing helix chain 'M' and resid 22 through 35 removed outlier: 3.548A pdb=" N VAL M 26 " --> pdb=" O THR M 22 " (cutoff:3.500A) Processing helix chain 'M' and resid 37 through 39 No H-bonds generated for 'chain 'M' and resid 37 through 39' Processing helix chain 'M' and resid 55 through 59 Processing helix chain 'N' and resid 263 through 268 Processing helix chain 'N' and resid 269 through 275 Processing helix chain 'N' and resid 277 through 286 removed outlier: 3.732A pdb=" N PHE N 286 " --> pdb=" O GLN N 282 " (cutoff:3.500A) Processing helix chain 'N' and resid 310 through 327 Processing helix chain 'N' and resid 330 through 352 removed outlier: 4.396A pdb=" N ARG N 336 " --> pdb=" O GLU N 332 " (cutoff:3.500A) removed outlier: 3.620A pdb=" N GLU N 344 " --> pdb=" O GLU N 340 " (cutoff:3.500A) removed outlier: 4.279A pdb=" N LYS N 345 " --> pdb=" O MET N 341 " (cutoff:3.500A) removed outlier: 3.839A pdb=" N LYS N 346 " --> pdb=" O GLU N 342 " (cutoff:3.500A) removed outlier: 3.626A pdb=" N VAL N 347 " --> pdb=" O LYS N 343 " (cutoff:3.500A) removed outlier: 3.582A pdb=" N GLN N 349 " --> pdb=" O LYS N 345 " (cutoff:3.500A) removed outlier: 5.350A pdb=" N TRP N 350 " --> pdb=" O LYS N 346 " (cutoff:3.500A) removed outlier: 4.378A pdb=" N LYS N 351 " --> pdb=" O VAL N 347 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 83 through 84 removed outlier: 7.038A pdb=" N ARG A 83 " --> pdb=" O VAL B 81 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA1 Processing sheet with id=AA2, first strand: chain 'B' and resid 96 through 98 removed outlier: 6.450A pdb=" N THR B 96 " --> pdb=" O THR G 101 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA2 Processing sheet with id=AA3, first strand: chain 'C' and resid 42 through 43 removed outlier: 7.464A pdb=" N ARG C 42 " --> pdb=" O ILE D 89 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA3 Processing sheet with id=AA4, first strand: chain 'C' and resid 77 through 78 Processing sheet with id=AA5, first strand: chain 'C' and resid 100 through 102 removed outlier: 6.577A pdb=" N THR C 101 " --> pdb=" O TYR F 98 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA5 Processing sheet with id=AA6, first strand: chain 'E' and resid 83 through 84 removed outlier: 7.026A pdb=" N ARG E 83 " --> pdb=" O VAL F 81 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA6 Processing sheet with id=AA7, first strand: chain 'E' and resid 118 through 119 Processing sheet with id=AA8, first strand: chain 'G' and resid 42 through 43 removed outlier: 7.223A pdb=" N ARG G 42 " --> pdb=" O ILE H 89 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA8 Processing sheet with id=AA9, first strand: chain 'G' and resid 77 through 78 Processing sheet with id=AB1, first strand: chain 'K' and resid 56 through 57 Processing sheet with id=AB2, first strand: chain 'L' and resid 8 through 9 removed outlier: 3.770A pdb=" N LEU L 8 " --> pdb=" O GLY M 75 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'L' and resid 28 through 30 removed outlier: 4.375A pdb=" N ALA L 167 " --> pdb=" O TRP L 52 " (cutoff:3.500A) removed outlier: 3.580A pdb=" N ILE L 191 " --> pdb=" O GLU L 182 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'M' and resid 12 through 15 Processing sheet with id=AB5, first strand: chain 'N' and resid 293 through 295 598 hydrogen bonds defined for protein. 1749 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 429 hydrogen bonds 858 hydrogen bond angles 0 basepair planarities 166 basepair parallelities 273 stacking parallelities Total time for adding SS restraints: 7.59 Time building geometry restraints manager: 9.64 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.33: 3600 1.33 - 1.45: 5397 1.45 - 1.57: 8649 1.57 - 1.69: 666 1.69 - 1.81: 48 Bond restraints: 18360 Sorted by residual: bond pdb=" N ARG M 74 " pdb=" CA ARG M 74 " ideal model delta sigma weight residual 1.452 1.490 -0.037 1.25e-02 6.40e+03 9.00e+00 bond pdb=" N GLY M 75 " pdb=" CA GLY M 75 " ideal model delta sigma weight residual 1.447 1.475 -0.028 1.03e-02 9.43e+03 7.23e+00 bond pdb=" N GLU L 7 " pdb=" CA GLU L 7 " ideal model delta sigma weight residual 1.456 1.487 -0.031 1.23e-02 6.61e+03 6.55e+00 bond pdb=" N LYS A 122 " pdb=" CA LYS A 122 " ideal model delta sigma weight residual 1.459 1.488 -0.029 1.21e-02 6.83e+03 5.77e+00 bond pdb=" N LEU M 73 " pdb=" CA LEU M 73 " ideal model delta sigma weight residual 1.456 1.490 -0.034 1.50e-02 4.44e+03 5.18e+00 ... (remaining 18355 not shown) Histogram of bond angle deviations from ideal: 98.36 - 105.48: 1592 105.48 - 112.60: 10177 112.60 - 119.72: 5950 119.72 - 126.84: 7556 126.84 - 133.96: 934 Bond angle restraints: 26209 Sorted by residual: angle pdb=" C GLY M 75 " pdb=" N GLY M 76 " pdb=" CA GLY M 76 " ideal model delta sigma weight residual 121.70 133.28 -11.58 1.80e+00 3.09e-01 4.14e+01 angle pdb=" O GLY M 75 " pdb=" C GLY M 75 " pdb=" N GLY M 76 " ideal model delta sigma weight residual 123.92 127.46 -3.54 9.60e-01 1.09e+00 1.36e+01 angle pdb=" N PHE L 660 " pdb=" CA PHE L 660 " pdb=" C PHE L 660 " ideal model delta sigma weight residual 111.36 107.39 3.97 1.09e+00 8.42e-01 1.33e+01 angle pdb=" C1' DG J 159 " pdb=" N9 DG J 159 " pdb=" C4 DG J 159 " ideal model delta sigma weight residual 127.00 122.04 4.96 1.50e+00 4.44e-01 1.10e+01 angle pdb=" CA ARG M 74 " pdb=" C ARG M 74 " pdb=" O ARG M 74 " ideal model delta sigma weight residual 121.82 118.26 3.56 1.21e+00 6.83e-01 8.68e+00 ... (remaining 26204 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 27.39: 8216 27.39 - 54.78: 1664 54.78 - 82.17: 308 82.17 - 109.56: 4 109.56 - 136.96: 2 Dihedral angle restraints: 10194 sinusoidal: 6377 harmonic: 3817 Sorted by residual: dihedral pdb=" C4' DC J 119 " pdb=" C3' DC J 119 " pdb=" O3' DC J 119 " pdb=" P DT J 120 " ideal model delta sinusoidal sigma weight residual 220.00 83.04 136.96 1 3.50e+01 8.16e-04 1.36e+01 dihedral pdb=" C4' DG J 27 " pdb=" C3' DG J 27 " pdb=" O3' DG J 27 " pdb=" P DT J 28 " ideal model delta sinusoidal sigma weight residual 220.00 87.80 132.20 1 3.50e+01 8.16e-04 1.31e+01 dihedral pdb=" CA PRO H 103 " pdb=" C PRO H 103 " pdb=" N GLY H 104 " pdb=" CA GLY H 104 " ideal model delta harmonic sigma weight residual 180.00 163.41 16.59 0 5.00e+00 4.00e-02 1.10e+01 ... (remaining 10191 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.034: 2194 0.034 - 0.068: 558 0.068 - 0.103: 144 0.103 - 0.137: 54 0.137 - 0.171: 10 Chirality restraints: 2960 Sorted by residual: chirality pdb=" C3' DG J 159 " pdb=" C4' DG J 159 " pdb=" O3' DG J 159 " pdb=" C2' DG J 159 " both_signs ideal model delta sigma weight residual False -2.66 -2.83 0.17 2.00e-01 2.50e+01 7.32e-01 chirality pdb=" C3' DC J 7 " pdb=" C4' DC J 7 " pdb=" O3' DC J 7 " pdb=" C2' DC J 7 " both_signs ideal model delta sigma weight residual False -2.66 -2.50 -0.16 2.00e-01 2.50e+01 6.72e-01 chirality pdb=" C3' DG J 10 " pdb=" C4' DG J 10 " pdb=" O3' DG J 10 " pdb=" C2' DG J 10 " both_signs ideal model delta sigma weight residual False -2.66 -2.82 0.16 2.00e-01 2.50e+01 6.28e-01 ... (remaining 2957 not shown) Planarity restraints: 2176 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA LYS L 659 " 0.010 2.00e-02 2.50e+03 1.95e-02 3.79e+00 pdb=" C LYS L 659 " -0.034 2.00e-02 2.50e+03 pdb=" O LYS L 659 " 0.013 2.00e-02 2.50e+03 pdb=" N PHE L 660 " 0.011 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C VAL L 234 " -0.027 5.00e-02 4.00e+02 4.09e-02 2.68e+00 pdb=" N PRO L 235 " 0.071 5.00e-02 4.00e+02 pdb=" CA PRO L 235 " -0.021 5.00e-02 4.00e+02 pdb=" CD PRO L 235 " -0.023 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C1' DC I 165 " 0.013 2.00e-02 2.50e+03 9.02e-03 1.83e+00 pdb=" N1 DC I 165 " -0.020 2.00e-02 2.50e+03 pdb=" C2 DC I 165 " 0.011 2.00e-02 2.50e+03 pdb=" O2 DC I 165 " -0.007 2.00e-02 2.50e+03 pdb=" N3 DC I 165 " -0.000 2.00e-02 2.50e+03 pdb=" C4 DC I 165 " 0.002 2.00e-02 2.50e+03 pdb=" N4 DC I 165 " 0.002 2.00e-02 2.50e+03 pdb=" C5 DC I 165 " -0.000 2.00e-02 2.50e+03 pdb=" C6 DC I 165 " -0.001 2.00e-02 2.50e+03 ... (remaining 2173 not shown) Histogram of nonbonded interaction distances: 1.34 - 2.05: 1 2.05 - 2.76: 2413 2.76 - 3.47: 24272 3.47 - 4.19: 47582 4.19 - 4.90: 72594 Nonbonded interactions: 146862 Sorted by model distance: nonbonded pdb=" NZ LYS C 119 " pdb=" C GLY M 76 " model vdw 1.336 3.350 nonbonded pdb=" CD1 LEU M 73 " pdb=" O ARG M 74 " model vdw 2.164 3.460 nonbonded pdb=" OH TYR E 99 " pdb=" OE1 GLU E 133 " model vdw 2.247 2.440 nonbonded pdb=" NZ LYS C 119 " pdb=" O GLY M 76 " model vdw 2.259 2.520 nonbonded pdb=" OG1 THR D 90 " pdb=" OE1 GLU D 93 " model vdw 2.278 2.440 ... (remaining 146857 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = chain 'A' selection = (chain 'E' and resid 38 through 133) } ncs_group { reference = chain 'B' selection = (chain 'F' and resid 22 through 101) } ncs_group { reference = (chain 'C' and resid 11 through 118) selection = (chain 'G' and resid 11 through 118) } ncs_group { reference = chain 'D' selection = chain 'H' } 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.00 max=1.00 mean=0.96 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: 2.250 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.020 Extract box with map and model: 9.390 Check model and map are aligned: 0.260 Set scattering table: 0.180 Process input model: 55.800 Find NCS groups from input model: 0.640 Set up NCS constraints: 0.080 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.480 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 80.100 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6614 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.037 18360 Z= 0.149 Angle : 0.532 11.575 26209 Z= 0.325 Chirality : 0.037 0.171 2960 Planarity : 0.003 0.041 2176 Dihedral : 24.612 136.956 7726 Min Nonbonded Distance : 1.336 Molprobity Statistics. All-atom Clashscore : 10.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.24 % Favored : 97.76 % Rotamer: Outliers : 0.00 % Allowed : 4.63 % Favored : 95.37 % 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.26 (0.24), residues: 1295 helix: 2.55 (0.19), residues: 759 sheet: -0.00 (0.61), residues: 70 loop : -0.02 (0.29), residues: 466 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP N 320 HIS 0.003 0.000 HIS L 169 PHE 0.005 0.001 PHE C 25 TYR 0.006 0.001 TYR C 50 ARG 0.002 0.000 ARG C 81 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 351 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 351 time to evaluate : 1.471 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 1 residues processed: 351 average time/residue: 1.4706 time to fit residues: 567.2125 Evaluate side-chains 221 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 220 time to evaluate : 1.535 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 1 residues processed: 1 average time/residue: 0.2977 time to fit residues: 2.3005 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 159 random chunks: chunk 134 optimal weight: 0.6980 chunk 120 optimal weight: 0.7980 chunk 66 optimal weight: 0.9990 chunk 41 optimal weight: 2.9990 chunk 81 optimal weight: 3.9990 chunk 64 optimal weight: 0.3980 chunk 124 optimal weight: 0.9990 chunk 48 optimal weight: 4.9990 chunk 75 optimal weight: 6.9990 chunk 92 optimal weight: 5.9990 chunk 144 optimal weight: 1.9990 overall best weight: 0.7784 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 93 GLN B 27 GLN B 93 GLN C 38 ASN E 39 HIS E 93 GLN ** F 25 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** L 84 HIS ** L 89 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** L 193 HIS M 62 GLN ** N 281 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6775 moved from start: 0.1910 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.042 18360 Z= 0.213 Angle : 0.601 15.353 26209 Z= 0.335 Chirality : 0.036 0.198 2960 Planarity : 0.004 0.045 2176 Dihedral : 29.001 136.399 5082 Min Nonbonded Distance : 2.161 Molprobity Statistics. All-atom Clashscore : 10.38 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.63 % Favored : 97.37 % Rotamer: Outliers : 3.74 % Allowed : 16.21 % Favored : 80.05 % 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.16 (0.24), residues: 1295 helix: 2.41 (0.19), residues: 785 sheet: -0.16 (0.61), residues: 71 loop : -0.08 (0.29), residues: 439 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.009 0.002 TRP L 5 HIS 0.005 0.001 HIS F 75 PHE 0.025 0.002 PHE L 118 TYR 0.012 0.001 TYR E 54 ARG 0.009 0.001 ARG N 302 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 271 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 42 poor density : 229 time to evaluate : 1.483 Fit side-chains revert: symmetry clash outliers start: 42 outliers final: 20 residues processed: 251 average time/residue: 1.3428 time to fit residues: 374.6563 Evaluate side-chains 221 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 20 poor density : 201 time to evaluate : 1.308 Switching outliers to nearest non-outliers outliers start: 20 outliers final: 10 residues processed: 10 average time/residue: 0.3056 time to fit residues: 6.1854 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 159 random chunks: chunk 80 optimal weight: 3.9990 chunk 44 optimal weight: 0.9990 chunk 120 optimal weight: 5.9990 chunk 98 optimal weight: 2.9990 chunk 39 optimal weight: 1.9990 chunk 144 optimal weight: 1.9990 chunk 156 optimal weight: 1.9990 chunk 128 optimal weight: 3.9990 chunk 143 optimal weight: 2.9990 chunk 49 optimal weight: 0.7980 chunk 116 optimal weight: 1.9990 overall best weight: 1.5588 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 25 ASN B 27 GLN C 24 GLN D 47 GLN D 109 HIS ** E 39 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** F 25 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** H 95 GLN L 84 HIS ** L 133 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** M 62 GLN N 281 GLN ** N 310 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6891 moved from start: 0.2809 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.049 18360 Z= 0.303 Angle : 0.632 24.653 26209 Z= 0.348 Chirality : 0.037 0.288 2960 Planarity : 0.004 0.058 2176 Dihedral : 29.300 135.419 5082 Min Nonbonded Distance : 2.176 Molprobity Statistics. All-atom Clashscore : 11.36 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.01 % Favored : 96.99 % Rotamer: Outliers : 5.43 % Allowed : 17.63 % Favored : 76.94 % 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: 1.77 (0.23), residues: 1295 helix: 2.09 (0.18), residues: 793 sheet: -0.15 (0.63), residues: 71 loop : -0.30 (0.28), residues: 431 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.002 TRP L 196 HIS 0.012 0.001 HIS N 331 PHE 0.026 0.002 PHE L 122 TYR 0.017 0.002 TYR A 54 ARG 0.007 0.001 ARG L 179 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 282 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 61 poor density : 221 time to evaluate : 1.528 Fit side-chains outliers start: 61 outliers final: 30 residues processed: 256 average time/residue: 1.2370 time to fit residues: 354.9961 Evaluate side-chains 233 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 30 poor density : 203 time to evaluate : 1.563 Switching outliers to nearest non-outliers outliers start: 30 outliers final: 16 residues processed: 16 average time/residue: 0.4941 time to fit residues: 12.4644 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 159 random chunks: chunk 142 optimal weight: 2.9990 chunk 108 optimal weight: 20.0000 chunk 75 optimal weight: 6.9990 chunk 16 optimal weight: 0.6980 chunk 69 optimal weight: 0.7980 chunk 97 optimal weight: 3.9990 chunk 145 optimal weight: 0.0980 chunk 153 optimal weight: 0.8980 chunk 137 optimal weight: 0.6980 chunk 41 optimal weight: 0.0040 chunk 127 optimal weight: 1.9990 overall best weight: 0.4592 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 108 ASN A 125 GLN B 27 GLN C 24 GLN D 47 GLN ** E 125 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** H 95 GLN L 84 HIS ** L 133 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** N 281 GLN ** N 310 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 8 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6835 moved from start: 0.2883 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.039 18360 Z= 0.177 Angle : 0.581 18.036 26209 Z= 0.322 Chirality : 0.035 0.212 2960 Planarity : 0.004 0.050 2176 Dihedral : 29.227 136.292 5082 Min Nonbonded Distance : 2.176 Molprobity Statistics. All-atom Clashscore : 10.15 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.55 % Favored : 97.45 % Rotamer: Outliers : 4.01 % Allowed : 19.86 % Favored : 76.14 % 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: 1.93 (0.23), residues: 1295 helix: 2.22 (0.18), residues: 791 sheet: -0.17 (0.63), residues: 71 loop : -0.21 (0.29), residues: 433 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.001 TRP L 196 HIS 0.004 0.001 HIS F 75 PHE 0.020 0.001 PHE L 118 TYR 0.011 0.001 TYR L 173 ARG 0.009 0.000 ARG M 72 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 263 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 45 poor density : 218 time to evaluate : 1.513 Fit side-chains revert: symmetry clash outliers start: 45 outliers final: 25 residues processed: 244 average time/residue: 1.2240 time to fit residues: 335.9864 Evaluate side-chains 222 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 25 poor density : 197 time to evaluate : 1.455 Switching outliers to nearest non-outliers outliers start: 25 outliers final: 16 residues processed: 10 average time/residue: 0.4741 time to fit residues: 8.1713 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 159 random chunks: chunk 87 optimal weight: 10.0000 chunk 2 optimal weight: 0.2980 chunk 114 optimal weight: 0.4980 chunk 63 optimal weight: 2.9990 chunk 131 optimal weight: 0.7980 chunk 106 optimal weight: 9.9990 chunk 0 optimal weight: 5.9990 chunk 78 optimal weight: 4.9990 chunk 137 optimal weight: 1.9990 chunk 38 optimal weight: 1.9990 chunk 51 optimal weight: 0.0070 overall best weight: 0.7200 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 125 GLN B 27 GLN D 95 GLN E 39 HIS E 125 GLN ** L 133 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** N 310 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6855 moved from start: 0.3050 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.041 18360 Z= 0.201 Angle : 0.573 14.902 26209 Z= 0.322 Chirality : 0.034 0.190 2960 Planarity : 0.004 0.048 2176 Dihedral : 29.190 135.668 5082 Min Nonbonded Distance : 2.169 Molprobity Statistics. All-atom Clashscore : 10.12 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.93 % Favored : 97.07 % Rotamer: Outliers : 3.92 % Allowed : 20.84 % Favored : 75.24 % 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.00 (0.24), residues: 1295 helix: 2.27 (0.19), residues: 791 sheet: -0.15 (0.64), residues: 71 loop : -0.16 (0.29), residues: 433 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.001 TRP L 196 HIS 0.004 0.001 HIS F 75 PHE 0.024 0.001 PHE L 118 TYR 0.011 0.001 TYR L 173 ARG 0.011 0.000 ARG M 72 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 256 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 44 poor density : 212 time to evaluate : 1.403 Fit side-chains outliers start: 44 outliers final: 24 residues processed: 239 average time/residue: 1.2355 time to fit residues: 330.0351 Evaluate side-chains 228 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 24 poor density : 204 time to evaluate : 1.293 Switching outliers to nearest non-outliers outliers start: 24 outliers final: 14 residues processed: 11 average time/residue: 0.3207 time to fit residues: 6.7614 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 159 random chunks: chunk 138 optimal weight: 0.7980 chunk 30 optimal weight: 0.9990 chunk 90 optimal weight: 10.0000 chunk 37 optimal weight: 3.9990 chunk 153 optimal weight: 0.9980 chunk 127 optimal weight: 0.6980 chunk 71 optimal weight: 0.9990 chunk 12 optimal weight: 0.9990 chunk 50 optimal weight: 0.8980 chunk 80 optimal weight: 3.9990 chunk 148 optimal weight: 2.9990 overall best weight: 0.8782 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 125 GLN B 27 GLN D 63 ASN E 125 GLN ** L 133 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** L 229 ASN ** N 310 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6870 moved from start: 0.3226 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.042 18360 Z= 0.215 Angle : 0.590 19.764 26209 Z= 0.327 Chirality : 0.035 0.179 2960 Planarity : 0.004 0.045 2176 Dihedral : 29.219 135.766 5082 Min Nonbonded Distance : 2.160 Molprobity Statistics. All-atom Clashscore : 10.12 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.93 % Favored : 97.07 % Rotamer: Outliers : 3.92 % Allowed : 21.02 % Favored : 75.07 % 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: 1.95 (0.24), residues: 1295 helix: 2.27 (0.18), residues: 788 sheet: -0.21 (0.60), residues: 81 loop : -0.22 (0.29), residues: 426 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.009 0.001 TRP L 196 HIS 0.004 0.001 HIS E 39 PHE 0.022 0.002 PHE L 118 TYR 0.013 0.001 TYR F 88 ARG 0.010 0.000 ARG M 72 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 264 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 44 poor density : 220 time to evaluate : 1.535 Fit side-chains outliers start: 44 outliers final: 27 residues processed: 244 average time/residue: 1.2953 time to fit residues: 356.2095 Evaluate side-chains 233 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 27 poor density : 206 time to evaluate : 1.548 Switching outliers to nearest non-outliers outliers start: 27 outliers final: 16 residues processed: 11 average time/residue: 0.5148 time to fit residues: 9.1088 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 159 random chunks: chunk 17 optimal weight: 9.9990 chunk 87 optimal weight: 8.9990 chunk 112 optimal weight: 0.8980 chunk 129 optimal weight: 2.9990 chunk 85 optimal weight: 4.9990 chunk 153 optimal weight: 0.9990 chunk 95 optimal weight: 4.9990 chunk 93 optimal weight: 3.9990 chunk 70 optimal weight: 0.9990 chunk 94 optimal weight: 4.9990 chunk 61 optimal weight: 1.9990 overall best weight: 1.5788 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 125 GLN B 27 GLN E 125 GLN ** L 133 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** L 229 ASN L 695 ASN M 62 GLN ** N 310 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 6 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6932 moved from start: 0.3681 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.049 18360 Z= 0.307 Angle : 0.662 30.637 26209 Z= 0.354 Chirality : 0.037 0.187 2960 Planarity : 0.005 0.055 2176 Dihedral : 29.381 134.691 5082 Min Nonbonded Distance : 2.119 Molprobity Statistics. All-atom Clashscore : 11.55 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.78 % Favored : 96.22 % Rotamer: Outliers : 4.45 % Allowed : 20.75 % Favored : 74.80 % 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: 1.57 (0.23), residues: 1295 helix: 1.96 (0.18), residues: 796 sheet: -0.27 (0.62), residues: 81 loop : -0.45 (0.29), residues: 418 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.002 TRP L 196 HIS 0.008 0.001 HIS L 144 PHE 0.037 0.002 PHE L 122 TYR 0.011 0.002 TYR A 54 ARG 0.011 0.001 ARG M 72 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 269 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 50 poor density : 219 time to evaluate : 1.489 Fit side-chains outliers start: 50 outliers final: 30 residues processed: 248 average time/residue: 1.2722 time to fit residues: 352.9157 Evaluate side-chains 234 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 30 poor density : 204 time to evaluate : 1.500 Switching outliers to nearest non-outliers outliers start: 30 outliers final: 19 residues processed: 11 average time/residue: 0.2035 time to fit residues: 5.7025 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 159 random chunks: chunk 91 optimal weight: 7.9990 chunk 46 optimal weight: 2.9990 chunk 30 optimal weight: 0.9990 chunk 29 optimal weight: 0.8980 chunk 97 optimal weight: 3.9990 chunk 104 optimal weight: 5.9990 chunk 75 optimal weight: 6.9990 chunk 14 optimal weight: 1.9990 chunk 120 optimal weight: 3.9990 chunk 139 optimal weight: 0.9980 chunk 146 optimal weight: 0.8980 overall best weight: 1.1584 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 125 GLN B 27 GLN B 93 GLN D 63 ASN D 95 GLN E 125 GLN ** L 133 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** L 229 ASN L 694 GLN L 710 HIS ** N 310 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6914 moved from start: 0.3787 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.049 18360 Z= 0.259 Angle : 0.649 24.522 26209 Z= 0.349 Chirality : 0.036 0.194 2960 Planarity : 0.004 0.057 2176 Dihedral : 29.379 135.220 5082 Min Nonbonded Distance : 2.091 Molprobity Statistics. All-atom Clashscore : 11.04 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.17 % Favored : 96.83 % Rotamer: Outliers : 3.83 % Allowed : 22.71 % Favored : 73.46 % 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: 1.61 (0.23), residues: 1295 helix: 1.98 (0.18), residues: 796 sheet: -0.30 (0.61), residues: 81 loop : -0.40 (0.29), residues: 418 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.002 TRP L 196 HIS 0.006 0.001 HIS L 710 PHE 0.020 0.002 PHE L 118 TYR 0.015 0.002 TYR H 83 ARG 0.016 0.001 ARG C 71 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 257 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 43 poor density : 214 time to evaluate : 1.639 Fit side-chains outliers start: 43 outliers final: 27 residues processed: 242 average time/residue: 1.3465 time to fit residues: 363.8584 Evaluate side-chains 237 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 27 poor density : 210 time to evaluate : 1.568 Switching outliers to nearest non-outliers outliers start: 27 outliers final: 20 residues processed: 7 average time/residue: 0.2106 time to fit residues: 4.4254 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 159 random chunks: chunk 133 optimal weight: 2.9990 chunk 142 optimal weight: 2.9990 chunk 146 optimal weight: 2.9990 chunk 85 optimal weight: 4.9990 chunk 62 optimal weight: 1.9990 chunk 112 optimal weight: 0.8980 chunk 43 optimal weight: 0.6980 chunk 129 optimal weight: 0.0870 chunk 135 optimal weight: 0.0370 chunk 93 optimal weight: 4.9990 chunk 151 optimal weight: 0.2980 overall best weight: 0.4036 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 125 GLN B 27 GLN D 49 HIS E 125 GLN G 73 ASN ** L 133 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** L 702 GLN ** N 281 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** N 310 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 6 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6864 moved from start: 0.3740 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.060 18360 Z= 0.191 Angle : 0.634 18.813 26209 Z= 0.342 Chirality : 0.035 0.246 2960 Planarity : 0.004 0.036 2176 Dihedral : 29.241 135.948 5082 Min Nonbonded Distance : 2.133 Molprobity Statistics. All-atom Clashscore : 11.04 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.09 % Favored : 96.91 % Rotamer: Outliers : 3.03 % Allowed : 23.86 % Favored : 73.11 % 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: 1.83 (0.24), residues: 1295 helix: 2.15 (0.18), residues: 787 sheet: -0.31 (0.63), residues: 73 loop : -0.21 (0.29), residues: 435 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.009 0.001 TRP L 196 HIS 0.006 0.001 HIS D 49 PHE 0.021 0.002 PHE L 122 TYR 0.027 0.002 TYR H 83 ARG 0.009 0.000 ARG C 71 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 260 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 34 poor density : 226 time to evaluate : 1.495 Fit side-chains outliers start: 34 outliers final: 23 residues processed: 251 average time/residue: 1.2509 time to fit residues: 352.8524 Evaluate side-chains 232 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 23 poor density : 209 time to evaluate : 1.506 Switching outliers to nearest non-outliers outliers start: 23 outliers final: 16 residues processed: 7 average time/residue: 0.4910 time to fit residues: 6.2901 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 159 random chunks: chunk 92 optimal weight: 5.9990 chunk 71 optimal weight: 2.9990 chunk 104 optimal weight: 4.9990 chunk 158 optimal weight: 5.9990 chunk 145 optimal weight: 0.3980 chunk 126 optimal weight: 0.6980 chunk 13 optimal weight: 0.1980 chunk 97 optimal weight: 3.9990 chunk 77 optimal weight: 4.9990 chunk 100 optimal weight: 3.9990 chunk 134 optimal weight: 0.6980 overall best weight: 0.9982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 125 GLN B 27 GLN ** C 73 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 49 HIS D 63 ASN E 125 GLN G 73 ASN ** L 133 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** L 710 HIS N 281 GLN ** N 294 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** N 310 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 8 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6894 moved from start: 0.3818 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.081 18360 Z= 0.246 Angle : 0.653 22.319 26209 Z= 0.353 Chirality : 0.036 0.250 2960 Planarity : 0.004 0.054 2176 Dihedral : 29.259 135.388 5082 Min Nonbonded Distance : 2.116 Molprobity Statistics. All-atom Clashscore : 11.46 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.40 % Favored : 96.60 % Rotamer: Outliers : 2.94 % Allowed : 24.13 % Favored : 72.93 % 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: 1.71 (0.23), residues: 1295 helix: 2.05 (0.18), residues: 793 sheet: -0.29 (0.63), residues: 73 loop : -0.31 (0.29), residues: 429 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.002 TRP L 202 HIS 0.004 0.001 HIS L 144 PHE 0.025 0.002 PHE L 122 TYR 0.019 0.002 TYR D 83 ARG 0.013 0.001 ARG C 71 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2590 Ramachandran restraints generated. 1295 Oldfield, 0 Emsley, 1295 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 237 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 33 poor density : 204 time to evaluate : 1.452 Fit side-chains outliers start: 33 outliers final: 21 residues processed: 229 average time/residue: 1.2825 time to fit residues: 328.8073 Evaluate side-chains 218 residues out of total 1123 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 21 poor density : 197 time to evaluate : 1.511 Switching outliers to nearest non-outliers outliers start: 21 outliers final: 18 residues processed: 3 average time/residue: 0.1590 time to fit residues: 2.5738 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 159 random chunks: chunk 38 optimal weight: 0.7980 chunk 116 optimal weight: 0.6980 chunk 18 optimal weight: 0.8980 chunk 35 optimal weight: 0.8980 chunk 126 optimal weight: 0.8980 chunk 52 optimal weight: 0.9980 chunk 129 optimal weight: 0.7980 chunk 15 optimal weight: 2.9990 chunk 23 optimal weight: 2.9990 chunk 110 optimal weight: 10.0000 chunk 7 optimal weight: 0.0370 overall best weight: 0.6458 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 125 GLN B 27 GLN D 49 HIS E 125 GLN G 73 ASN ** L 133 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** L 694 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** L 710 HIS N 281 GLN ** N 294 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** N 310 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 7 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4117 r_free = 0.4117 target = 0.159652 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 33)----------------| | r_work = 0.3845 r_free = 0.3845 target = 0.136804 restraints weight = 43406.310| |-----------------------------------------------------------------------------| r_work (start): 0.3802 rms_B_bonded: 2.47 r_work: 0.3657 rms_B_bonded: 3.65 restraints_weight: 0.5000 r_work (final): 0.3657 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3665 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3665 r_free = 0.3665 target_work(ls_wunit_k1) = 0.128 | | occupancies: max = 1.00 min = 0.00 number of occupancies < 0.1: 495 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3656 r_free = 0.3656 target_work(ls_wunit_k1) = 0.127 | | occupancies: max = 1.00 min = 0.00 number of occupancies < 0.1: 561 | |-----------------------------------------------------------------------------| r_final: 0.3656 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7268 moved from start: 0.3838 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.069 18360 Z= 0.212 Angle : 0.644 21.139 26209 Z= 0.347 Chirality : 0.036 0.251 2960 Planarity : 0.004 0.052 2176 Dihedral : 29.235 135.989 5082 Min Nonbonded Distance : 2.167 Molprobity Statistics. All-atom Clashscore : 11.95 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.40 % Favored : 96.60 % Rotamer: Outliers : 2.49 % Allowed : 24.93 % Favored : 72.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: 1.74 (0.23), residues: 1295 helix: 2.07 (0.18), residues: 791 sheet: -0.25 (0.63), residues: 73 loop : -0.27 (0.29), residues: 431 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.022 0.002 TRP L 202 HIS 0.004 0.001 HIS E 39 PHE 0.026 0.002 PHE L 118 TYR 0.018 0.002 TYR D 83 ARG 0.013 0.000 ARG C 71 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 6065.00 seconds wall clock time: 108 minutes 9.92 seconds (6489.92 seconds total)