Starting phenix.real_space_refine (version: dev) on Sat Feb 18 20:41:52 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6n9u_0379/02_2023/6n9u_0379_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6n9u_0379/02_2023/6n9u_0379.map Found restraint, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6n9u_0379/02_2023/6n9u_0379_ligands.cif Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.7 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { restraint_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6n9u_0379/02_2023/6n9u_0379_ligands.cif" default_restraint = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6n9u_0379/02_2023/6n9u_0379_ligands.cif" real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6n9u_0379/02_2023/6n9u_0379.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6n9u_0379/02_2023/6n9u_0379.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6n9u_0379/02_2023/6n9u_0379_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6n9u_0379/02_2023/6n9u_0379_updated.pdb" } resolution = 3.7 write_initial_geo_file = False refinement { macro_cycles = 10 } 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 ------------------------------------------------------------------------------- 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.000 sd= 0.003 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 Process input model Symmetric amino acids flipped Residue "E ARG 124": "NH1" <-> "NH2" Residue "E ARG 229": "NH1" <-> "NH2" Residue "E ARG 248": "NH1" <-> "NH2" Residue "F ARG 213": "NH1" <-> "NH2" Residue "H ARG 85": "NH1" <-> "NH2" Residue "H ARG 96": "NH1" <-> "NH2" Residue "H PHE 232": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "H ARG 255": "NH1" <-> "NH2" Residue "H ARG 429": "NH1" <-> "NH2" Residue "H ARG 547": "NH1" <-> "NH2" Residue "H ARG 566": "NH1" <-> "NH2" Residue "H ARG 591": "NH1" <-> "NH2" Residue "H ARG 599": "NH1" <-> "NH2" Residue "H TYR 648": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "H ARG 687": "NH1" <-> "NH2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/chem_data/mon_lib" Total number of atoms: 8822 Number of models: 1 Model: "" Number of chains: 7 Chain: "E" Number of atoms: 1740 Number of conformers: 1 Conformer: "" Number of residues, atoms: 222, 1740 Classifications: {'peptide': 222} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 5, 'TRANS': 216} Chain breaks: 2 Unresolved non-hydrogen bonds: 3 Unresolved non-hydrogen angles: 3 Unresolved non-hydrogen dihedrals: 3 Chain: "F" Number of atoms: 1538 Number of conformers: 1 Conformer: "" Number of residues, atoms: 196, 1538 Classifications: {'peptide': 196} Link IDs: {'PTRANS': 4, 'TRANS': 191} Chain: "H" Number of atoms: 5144 Number of conformers: 1 Conformer: "" Number of residues, atoms: 645, 5144 Classifications: {'peptide': 645} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 6} Link IDs: {'PCIS': 1, 'PTRANS': 21, 'TRANS': 622} Chain breaks: 2 Unresolved non-hydrogen bonds: 22 Unresolved non-hydrogen angles: 28 Unresolved non-hydrogen dihedrals: 16 Planarities with less than four sites: {'GLN:plan1': 1, 'GLU:plan': 3, 'ASP:plan': 2} Unresolved non-hydrogen planarities: 22 Chain: "P" Number of atoms: 122 Number of conformers: 1 Conformer: "" Number of residues, atoms: 6, 122 Unusual residues: {'DOC': 1} Classifications: {'RNA': 5, 'undetermined': 1} Modifications used: {'5*END': 1, 'rna3p_pur': 3, 'rna3p_pyr': 2} Link IDs: {'rna3p': 4, None: 1} Not linked: pdbres=" G P 5 " pdbres="DOC P 6 " Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 3 Chain: "T" Number of atoms: 246 Number of conformers: 1 Conformer: "" Number of residues, atoms: 12, 246 Classifications: {'DNA': 12} Link IDs: {'rna3p': 11} Chain: "E" Number of atoms: 1 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 1 Unusual residues: {' ZN': 1} Classifications: {'undetermined': 1} Chain: "H" Number of atoms: 31 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 31 Unusual residues: {' MG': 2, 'TTP': 1} Classifications: {'undetermined': 3} Link IDs: {None: 2} List of CYS excluded from plausible disulfide bonds: (reason: may participate in coordination) ATOM 69 SG CYS E 17 59.453 36.079 46.143 1.00 84.80 S ATOM 91 SG CYS E 20 58.126 33.933 44.451 1.00 85.19 S ATOM 204 SG CYS E 36 56.118 35.082 46.849 1.00 66.69 S ATOM 229 SG CYS E 39 58.279 32.837 47.885 1.00 69.94 S Time building chain proxies: 5.85, per 1000 atoms: 0.66 Number of scatterers: 8822 At special positions: 0 Unit cell: (86.86, 100.62, 121.26, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 7 Type Number sf(0) Zn 1 29.99 S 47 16.00 P 20 15.00 Mg 2 11.99 O 1694 8.00 N 1531 7.00 C 5527 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=1, symmetry=0 Simple disulfide: pdb=" SG CYS F 235 " - pdb=" SG CYS F 241 " 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=1, symmetry=0 Number of additional bonds: simple=1, symmetry=0 Coordination: Other bonds: Time building additional restraints: 3.94 Conformation dependent library (CDL) restraints added in 1.6 seconds Dynamic metal coordination Zn2+ tetrahedral coordination pdb=" ZN E 601 " pdb="ZN ZN E 601 " - pdb=" SG CYS E 17 " pdb="ZN ZN E 601 " - pdb=" SG CYS E 39 " pdb="ZN ZN E 601 " - pdb=" SG CYS E 36 " pdb="ZN ZN E 601 " - pdb=" SG CYS E 20 " Number of angles added : 6 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1950 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 33 helices and 11 sheets defined 30.9% alpha, 10.6% beta 5 base pairs and 10 stacking pairs defined. Time for finding SS restraints: 2.95 Creating SS restraints... Processing helix chain 'E' and resid 88 through 94 Processing helix chain 'E' and resid 101 through 103 No H-bonds generated for 'chain 'E' and resid 101 through 103' Processing helix chain 'E' and resid 160 through 169 Processing helix chain 'E' and resid 188 through 191 No H-bonds generated for 'chain 'E' and resid 188 through 191' Processing helix chain 'E' and resid 196 through 198 No H-bonds generated for 'chain 'E' and resid 196 through 198' Processing helix chain 'E' and resid 238 through 244 removed outlier: 3.744A pdb=" N ASN E 244 " --> pdb=" O GLU E 240 " (cutoff:3.500A) Processing helix chain 'E' and resid 246 through 255 removed outlier: 4.490A pdb=" N GLU E 252 " --> pdb=" O GLU E 249 " (cutoff:3.500A) removed outlier: 3.601A pdb=" N VAL E 254 " --> pdb=" O MET E 251 " (cutoff:3.500A) removed outlier: 3.928A pdb=" N TRP E 255 " --> pdb=" O GLU E 252 " (cutoff:3.500A) Processing helix chain 'F' and resid 88 through 94 Processing helix chain 'F' and resid 160 through 168 removed outlier: 3.601A pdb=" N VAL F 165 " --> pdb=" O ASP F 161 " (cutoff:3.500A) Processing helix chain 'F' and resid 185 through 192 removed outlier: 3.647A pdb=" N ALA F 191 " --> pdb=" O LYS F 187 " (cutoff:3.500A) removed outlier: 4.065A pdb=" N ASN F 192 " --> pdb=" O THR F 188 " (cutoff:3.500A) Processing helix chain 'F' and resid 210 through 219 Processing helix chain 'F' and resid 238 through 243 Processing helix chain 'F' and resid 250 through 253 No H-bonds generated for 'chain 'F' and resid 250 through 253' Processing helix chain 'H' and resid 39 through 51 Processing helix chain 'H' and resid 65 through 75 removed outlier: 3.527A pdb=" N LYS H 74 " --> pdb=" O THR H 70 " (cutoff:3.500A) Processing helix chain 'H' and resid 95 through 98 No H-bonds generated for 'chain 'H' and resid 95 through 98' Processing helix chain 'H' and resid 125 through 131 Processing helix chain 'H' and resid 139 through 150 Processing helix chain 'H' and resid 165 through 185 removed outlier: 3.856A pdb=" N LYS H 184 " --> pdb=" O ALA H 180 " (cutoff:3.500A) Processing helix chain 'H' and resid 213 through 230 removed outlier: 3.635A pdb=" N GLU H 217 " --> pdb=" O ALA H 213 " (cutoff:3.500A) removed outlier: 3.941A pdb=" N HIS H 218 " --> pdb=" O VAL H 214 " (cutoff:3.500A) Processing helix chain 'H' and resid 237 through 246 removed outlier: 3.917A pdb=" N LEU H 242 " --> pdb=" O ALA H 238 " (cutoff:3.500A) Processing helix chain 'H' and resid 255 through 258 No H-bonds generated for 'chain 'H' and resid 255 through 258' Processing helix chain 'H' and resid 341 through 347 Processing helix chain 'H' and resid 377 through 399 removed outlier: 3.926A pdb=" N ALA H 381 " --> pdb=" O PRO H 377 " (cutoff:3.500A) removed outlier: 4.509A pdb=" N ASP H 384 " --> pdb=" O GLN H 380 " (cutoff:3.500A) removed outlier: 3.729A pdb=" N LEU H 385 " --> pdb=" O ALA H 381 " (cutoff:3.500A) Processing helix chain 'H' and resid 406 through 408 No H-bonds generated for 'chain 'H' and resid 406 through 408' Processing helix chain 'H' and resid 449 through 452 No H-bonds generated for 'chain 'H' and resid 449 through 452' Processing helix chain 'H' and resid 479 through 487 removed outlier: 3.642A pdb=" N ALA H 485 " --> pdb=" O LEU H 481 " (cutoff:3.500A) Processing helix chain 'H' and resid 505 through 513 Processing helix chain 'H' and resid 518 through 529 removed outlier: 3.722A pdb=" N TYR H 526 " --> pdb=" O LYS H 522 " (cutoff:3.500A) Processing helix chain 'H' and resid 545 through 556 removed outlier: 4.054A pdb=" N LYS H 552 " --> pdb=" O GLY H 548 " (cutoff:3.500A) Processing helix chain 'H' and resid 562 through 573 Processing helix chain 'H' and resid 609 through 636 removed outlier: 4.130A pdb=" N SER H 616 " --> pdb=" O THR H 612 " (cutoff:3.500A) removed outlier: 3.876A pdb=" N ALA H 617 " --> pdb=" O LEU H 613 " (cutoff:3.500A) removed outlier: 4.142A pdb=" N LEU H 620 " --> pdb=" O SER H 616 " (cutoff:3.500A) removed outlier: 4.164A pdb=" N ILE H 621 " --> pdb=" O ALA H 617 " (cutoff:3.500A) removed outlier: 3.867A pdb=" N CYS H 622 " --> pdb=" O GLY H 618 " (cutoff:3.500A) removed outlier: 3.785A pdb=" N ILE H 627 " --> pdb=" O LYS H 623 " (cutoff:3.500A) removed outlier: 4.016A pdb=" N GLU H 630 " --> pdb=" O ILE H 626 " (cutoff:3.500A) removed outlier: 4.400A pdb=" N LEU H 633 " --> pdb=" O THR H 629 " (cutoff:3.500A) removed outlier: 3.750A pdb=" N GLU H 635 " --> pdb=" O GLU H 631 " (cutoff:3.500A) Processing helix chain 'H' and resid 663 through 683 removed outlier: 4.577A pdb=" N GLN H 674 " --> pdb=" O ILE H 670 " (cutoff:3.500A) removed outlier: 3.736A pdb=" N GLU H 675 " --> pdb=" O GLU H 671 " (cutoff:3.500A) Processing sheet with id= A, first strand: chain 'E' and resid 34 through 36 removed outlier: 3.667A pdb=" N TYR E 13 " --> pdb=" O LEU E 28 " (cutoff:3.500A) Processing sheet with id= B, first strand: chain 'E' and resid 96 through 99 removed outlier: 3.740A pdb=" N SER E 120 " --> pdb=" O TYR E 111 " (cutoff:3.500A) Processing sheet with id= C, first strand: chain 'E' and resid 175 through 177 removed outlier: 7.986A pdb=" N VAL E 176 " --> pdb=" O LYS E 152 " (cutoff:3.500A) removed outlier: 5.831A pdb=" N VAL E 154 " --> pdb=" O VAL E 176 " (cutoff:3.500A) removed outlier: 6.845A pdb=" N GLN E 201 " --> pdb=" O ILE E 153 " (cutoff:3.500A) removed outlier: 7.829A pdb=" N VAL E 155 " --> pdb=" O GLN E 201 " (cutoff:3.500A) removed outlier: 6.545A pdb=" N ILE E 203 " --> pdb=" O VAL E 155 " (cutoff:3.500A) Processing sheet with id= D, first strand: chain 'F' and resid 97 through 100 removed outlier: 4.018A pdb=" N SER F 120 " --> pdb=" O TYR F 111 " (cutoff:3.500A) Processing sheet with id= E, first strand: chain 'F' and resid 153 through 155 Processing sheet with id= F, first strand: chain 'H' and resid 32 through 34 removed outlier: 4.188A pdb=" N GLY H 21 " --> pdb=" O TYR H 34 " (cutoff:3.500A) removed outlier: 6.808A pdb=" N GLU H 7 " --> pdb=" O HIS H 19 " (cutoff:3.500A) removed outlier: 4.408A pdb=" N GLY H 21 " --> pdb=" O ASP H 5 " (cutoff:3.500A) removed outlier: 6.572A pdb=" N ASP H 5 " --> pdb=" O GLY H 21 " (cutoff:3.500A) removed outlier: 4.621A pdb=" N ILE H 23 " --> pdb=" O VAL H 3 " (cutoff:3.500A) removed outlier: 6.564A pdb=" N VAL H 3 " --> pdb=" O ILE H 23 " (cutoff:3.500A) Processing sheet with id= G, first strand: chain 'H' and resid 264 through 267 Processing sheet with id= H, first strand: chain 'H' and resid 419 through 421 Processing sheet with id= I, first strand: chain 'H' and resid 692 through 697 removed outlier: 3.685A pdb=" N ALA H 650 " --> pdb=" O GLN H 657 " (cutoff:3.500A) removed outlier: 6.032A pdb=" N GLY H 659 " --> pdb=" O TYR H 648 " (cutoff:3.500A) removed outlier: 5.401A pdb=" N TYR H 648 " --> pdb=" O GLY H 659 " (cutoff:3.500A) Processing sheet with id= J, first strand: chain 'H' and resid 592 through 594 Processing sheet with id= K, first strand: chain 'H' and resid 574 through 579 removed outlier: 4.182A pdb=" N GLU H 575 " --> pdb=" O LYS H 587 " (cutoff:3.500A) removed outlier: 7.503A pdb=" N LYS H 587 " --> pdb=" O GLU H 575 " (cutoff:3.500A) removed outlier: 4.868A pdb=" N SER H 577 " --> pdb=" O GLN H 585 " (cutoff:3.500A) removed outlier: 6.323A pdb=" N GLN H 585 " --> pdb=" O SER H 577 " (cutoff:3.500A) No H-bonds generated for sheet with id= K 225 hydrogen bonds defined for protein. 636 hydrogen bond angles defined for protein. Warning, Cannot make NA restraints for OD residue Restraints generated for nucleic acids: 13 hydrogen bonds 26 hydrogen bond angles 0 basepair planarities 5 basepair parallelities 9 stacking parallelities Total time for adding SS restraints: 2.46 Time building geometry restraints manager: 4.00 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.18 - 1.31: 1466 1.31 - 1.44: 2556 1.44 - 1.56: 4926 1.56 - 1.69: 37 1.69 - 1.82: 72 Bond restraints: 9057 Sorted by residual: bond pdb=" C4' DOC P 6 " pdb=" C3' DOC P 6 " ideal model delta sigma weight residual 1.334 1.516 -0.182 2.00e-02 2.50e+03 8.29e+01 bond pdb=" C4' DOC P 6 " pdb=" O4' DOC P 6 " ideal model delta sigma weight residual 1.609 1.442 0.167 2.00e-02 2.50e+03 6.99e+01 bond pdb=" O4' DOC P 6 " pdb=" C1' DOC P 6 " ideal model delta sigma weight residual 1.268 1.416 -0.148 2.00e-02 2.50e+03 5.50e+01 bond pdb=" C4 DOC P 6 " pdb=" N4 DOC P 6 " ideal model delta sigma weight residual 1.494 1.348 0.146 2.00e-02 2.50e+03 5.36e+01 bond pdb=" C2 DOC P 6 " pdb=" N3 DOC P 6 " ideal model delta sigma weight residual 1.492 1.374 0.118 2.00e-02 2.50e+03 3.50e+01 ... (remaining 9052 not shown) Histogram of bond angle deviations from ideal: 99.51 - 106.41: 309 106.41 - 113.31: 4743 113.31 - 120.21: 3526 120.21 - 127.12: 3586 127.12 - 134.02: 159 Bond angle restraints: 12323 Sorted by residual: angle pdb=" N ILE H 464 " pdb=" CA ILE H 464 " pdb=" C ILE H 464 " ideal model delta sigma weight residual 113.71 109.54 4.17 9.50e-01 1.11e+00 1.93e+01 angle pdb=" N GLY F 134 " pdb=" CA GLY F 134 " pdb=" C GLY F 134 " ideal model delta sigma weight residual 110.29 115.14 -4.85 1.28e+00 6.10e-01 1.44e+01 angle pdb=" C ALA H 513 " pdb=" N GLU H 514 " pdb=" CA GLU H 514 " ideal model delta sigma weight residual 122.08 127.17 -5.09 1.47e+00 4.63e-01 1.20e+01 angle pdb=" N PHE E 206 " pdb=" CA PHE E 206 " pdb=" C PHE E 206 " ideal model delta sigma weight residual 110.80 117.99 -7.19 2.13e+00 2.20e-01 1.14e+01 angle pdb=" CA GLU H 514 " pdb=" CB GLU H 514 " pdb=" CG GLU H 514 " ideal model delta sigma weight residual 114.10 120.32 -6.22 2.00e+00 2.50e-01 9.67e+00 ... (remaining 12318 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 33.27: 5167 33.27 - 66.54: 91 66.54 - 99.81: 9 99.81 - 133.07: 0 133.07 - 166.34: 1 Dihedral angle restraints: 5268 sinusoidal: 2221 harmonic: 3047 Sorted by residual: dihedral pdb=" C2 TTP H 801 " pdb=" C1' TTP H 801 " pdb=" N1 TTP H 801 " pdb=" O4' TTP H 801 " ideal model delta sinusoidal sigma weight residual 301.68 135.34 166.34 1 2.00e+01 2.50e-03 4.73e+01 dihedral pdb=" CA TYR H 64 " pdb=" C TYR H 64 " pdb=" N ASP H 65 " pdb=" CA ASP H 65 " ideal model delta harmonic sigma weight residual -180.00 -153.26 -26.74 0 5.00e+00 4.00e-02 2.86e+01 dihedral pdb=" C5' TTP H 801 " pdb=" O5' TTP H 801 " pdb=" PA TTP H 801 " pdb=" O3A TTP H 801 " ideal model delta sinusoidal sigma weight residual 179.97 84.34 95.63 1 2.00e+01 2.50e-03 2.64e+01 ... (remaining 5265 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.062: 1017 0.062 - 0.124: 250 0.124 - 0.187: 34 0.187 - 0.249: 5 0.249 - 0.311: 1 Chirality restraints: 1307 Sorted by residual: chirality pdb=" C4' TTP H 801 " pdb=" C3' TTP H 801 " pdb=" C5' TTP H 801 " pdb=" O4' TTP H 801 " both_signs ideal model delta sigma weight residual False -2.50 -2.81 0.31 2.00e-01 2.50e+01 2.42e+00 chirality pdb=" CB ILE H 505 " pdb=" CA ILE H 505 " pdb=" CG1 ILE H 505 " pdb=" CG2 ILE H 505 " both_signs ideal model delta sigma weight residual False 2.64 2.41 0.23 2.00e-01 2.50e+01 1.34e+00 chirality pdb=" C1' DOC P 6 " pdb=" O4' DOC P 6 " pdb=" C2' DOC P 6 " pdb=" N1 DOC P 6 " both_signs ideal model delta sigma weight residual False 2.33 2.55 -0.21 2.00e-01 2.50e+01 1.14e+00 ... (remaining 1304 not shown) Planarity restraints: 1526 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB TRP H 700 " 0.024 2.00e-02 2.50e+03 2.01e-02 1.01e+01 pdb=" CG TRP H 700 " -0.053 2.00e-02 2.50e+03 pdb=" CD1 TRP H 700 " 0.024 2.00e-02 2.50e+03 pdb=" CD2 TRP H 700 " -0.002 2.00e-02 2.50e+03 pdb=" NE1 TRP H 700 " -0.002 2.00e-02 2.50e+03 pdb=" CE2 TRP H 700 " -0.002 2.00e-02 2.50e+03 pdb=" CE3 TRP H 700 " -0.000 2.00e-02 2.50e+03 pdb=" CZ2 TRP H 700 " 0.005 2.00e-02 2.50e+03 pdb=" CZ3 TRP H 700 " 0.006 2.00e-02 2.50e+03 pdb=" CH2 TRP H 700 " 0.000 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TRP F 97 " 0.017 2.00e-02 2.50e+03 1.84e-02 8.46e+00 pdb=" CG TRP F 97 " -0.050 2.00e-02 2.50e+03 pdb=" CD1 TRP F 97 " 0.022 2.00e-02 2.50e+03 pdb=" CD2 TRP F 97 " -0.000 2.00e-02 2.50e+03 pdb=" NE1 TRP F 97 " 0.003 2.00e-02 2.50e+03 pdb=" CE2 TRP F 97 " 0.000 2.00e-02 2.50e+03 pdb=" CE3 TRP F 97 " 0.005 2.00e-02 2.50e+03 pdb=" CZ2 TRP F 97 " -0.005 2.00e-02 2.50e+03 pdb=" CZ3 TRP F 97 " 0.006 2.00e-02 2.50e+03 pdb=" CH2 TRP F 97 " 0.001 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TRP E 147 " 0.020 2.00e-02 2.50e+03 1.79e-02 7.97e+00 pdb=" CG TRP E 147 " -0.047 2.00e-02 2.50e+03 pdb=" CD1 TRP E 147 " 0.023 2.00e-02 2.50e+03 pdb=" CD2 TRP E 147 " -0.002 2.00e-02 2.50e+03 pdb=" NE1 TRP E 147 " -0.001 2.00e-02 2.50e+03 pdb=" CE2 TRP E 147 " -0.003 2.00e-02 2.50e+03 pdb=" CE3 TRP E 147 " -0.001 2.00e-02 2.50e+03 pdb=" CZ2 TRP E 147 " 0.000 2.00e-02 2.50e+03 pdb=" CZ3 TRP E 147 " 0.007 2.00e-02 2.50e+03 pdb=" CH2 TRP E 147 " 0.003 2.00e-02 2.50e+03 ... (remaining 1523 not shown) Histogram of nonbonded interaction distances: 2.21 - 2.74: 792 2.74 - 3.28: 7736 3.28 - 3.82: 13959 3.82 - 4.36: 16053 4.36 - 4.90: 27841 Nonbonded interactions: 66381 Sorted by model distance: nonbonded pdb=" OD2 ASP H 654 " pdb="MG MG H 802 " model vdw 2.205 2.170 nonbonded pdb=" O GLY F 158 " pdb=" OG SER F 177 " model vdw 2.254 2.440 nonbonded pdb=" O ALA H 476 " pdb="MG MG H 802 " model vdw 2.274 2.170 nonbonded pdb=" O2B TTP H 801 " pdb="MG MG H 802 " model vdw 2.283 2.170 nonbonded pdb=" OD2 ASP H 197 " pdb=" OG1 THR H 199 " model vdw 2.293 2.440 ... (remaining 66376 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 Set scattering table Set to: electron Number of scattering types: 7 Type Number sf(0) Gaussians Zn 1 6.06 5 P 20 5.49 5 Mg 2 5.21 5 S 47 5.16 5 C 5527 2.51 5 N 1531 2.21 5 O 1694 1.98 5 sf(0) = scattering factor at diffraction angle 0. 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 group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 6.550 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: 13.350 Check model and map are aligned: 0.140 Process input model: 29.680 Find NCS groups from input model: 0.260 Set up NCS constraints: 0.030 Set refine NCS operators: 0.000 Set scattering table: 0.100 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:2.650 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 52.770 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8243 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.182 9057 Z= 0.453 Angle : 0.837 9.427 12323 Z= 0.450 Chirality : 0.055 0.311 1307 Planarity : 0.006 0.059 1526 Dihedral : 12.721 166.343 3315 Min Nonbonded Distance : 2.205 Molprobity Statistics. All-atom Clashscore : 3.66 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.34 % Favored : 94.66 % Rotamer Outliers : 0.34 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % 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: -4.67 (0.20), residues: 1049 helix: -3.21 (0.18), residues: 422 sheet: -2.67 (0.40), residues: 142 loop : -3.04 (0.23), residues: 485 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 240 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 237 time to evaluate : 1.027 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 3 outliers final: 0 residues processed: 237 average time/residue: 0.2934 time to fit residues: 88.2596 Evaluate side-chains 131 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 131 time to evaluate : 1.005 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.3490 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 89 optimal weight: 0.9990 chunk 80 optimal weight: 0.8980 chunk 44 optimal weight: 1.9990 chunk 27 optimal weight: 0.9980 chunk 54 optimal weight: 1.9990 chunk 42 optimal weight: 0.5980 chunk 83 optimal weight: 0.6980 chunk 32 optimal weight: 1.9990 chunk 50 optimal weight: 0.9990 chunk 61 optimal weight: 1.9990 chunk 96 optimal weight: 0.4980 overall best weight: 0.7380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: E 70 ASN E 129 ASN E 201 GLN E 244 ASN E 253 GLN E 256 ASN F 136 HIS H 150 GLN ** H 347 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** H 439 GLN ** H 471 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** H 520 ASN H 615 GLN ** H 657 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 11 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8340 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.035 9057 Z= 0.236 Angle : 0.588 10.413 12323 Z= 0.309 Chirality : 0.044 0.188 1307 Planarity : 0.005 0.054 1526 Dihedral : 11.663 173.109 1332 Min Nonbonded Distance : 1.921 Molprobity Statistics. All-atom Clashscore : 6.80 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.67 % Favored : 95.33 % Rotamer Outliers : 2.17 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.54 (0.23), residues: 1049 helix: -1.92 (0.23), residues: 424 sheet: -2.46 (0.40), residues: 153 loop : -2.64 (0.25), residues: 472 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 168 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 19 poor density : 149 time to evaluate : 0.907 Fit side-chains outliers start: 19 outliers final: 14 residues processed: 163 average time/residue: 0.2612 time to fit residues: 55.9187 Evaluate side-chains 131 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 117 time to evaluate : 0.953 Switching outliers to nearest non-outliers outliers start: 14 outliers final: 0 residues processed: 14 average time/residue: 0.0816 time to fit residues: 3.4403 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 53 optimal weight: 2.9990 chunk 29 optimal weight: 2.9990 chunk 80 optimal weight: 1.9990 chunk 65 optimal weight: 1.9990 chunk 26 optimal weight: 0.6980 chunk 96 optimal weight: 0.6980 chunk 104 optimal weight: 0.7980 chunk 85 optimal weight: 0.9990 chunk 95 optimal weight: 1.9990 chunk 32 optimal weight: 1.9990 chunk 77 optimal weight: 0.8980 overall best weight: 0.8182 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: H 101 ASN H 171 ASN ** H 471 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 657 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8374 moved from start: 0.2771 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.048 9057 Z= 0.252 Angle : 0.567 8.912 12323 Z= 0.299 Chirality : 0.044 0.171 1307 Planarity : 0.004 0.060 1526 Dihedral : 11.581 172.121 1332 Min Nonbonded Distance : 1.920 Molprobity Statistics. All-atom Clashscore : 6.97 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.86 % Favored : 95.14 % Rotamer Outliers : 2.29 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % 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.95 (0.24), residues: 1049 helix: -1.33 (0.25), residues: 419 sheet: -2.06 (0.42), residues: 146 loop : -2.46 (0.25), residues: 484 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 153 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 20 poor density : 133 time to evaluate : 1.039 Fit side-chains outliers start: 20 outliers final: 10 residues processed: 140 average time/residue: 0.2463 time to fit residues: 46.5470 Evaluate side-chains 124 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 114 time to evaluate : 1.035 Switching outliers to nearest non-outliers outliers start: 10 outliers final: 0 residues processed: 10 average time/residue: 0.1539 time to fit residues: 3.7634 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 95 optimal weight: 0.7980 chunk 72 optimal weight: 3.9990 chunk 50 optimal weight: 1.9990 chunk 10 optimal weight: 2.9990 chunk 46 optimal weight: 1.9990 chunk 64 optimal weight: 1.9990 chunk 96 optimal weight: 0.7980 chunk 102 optimal weight: 0.7980 chunk 91 optimal weight: 3.9990 chunk 27 optimal weight: 3.9990 chunk 85 optimal weight: 0.9980 overall best weight: 1.0782 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** H 171 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 471 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 657 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8412 moved from start: 0.3157 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.047 9057 Z= 0.311 Angle : 0.605 8.654 12323 Z= 0.315 Chirality : 0.045 0.174 1307 Planarity : 0.005 0.065 1526 Dihedral : 11.618 173.320 1332 Min Nonbonded Distance : 1.895 Molprobity Statistics. All-atom Clashscore : 7.79 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.34 % Favored : 94.66 % Rotamer Outliers : 2.40 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % 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.81 (0.24), residues: 1049 helix: -1.19 (0.25), residues: 426 sheet: -2.07 (0.43), residues: 148 loop : -2.37 (0.26), residues: 475 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 142 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 21 poor density : 121 time to evaluate : 0.974 Fit side-chains outliers start: 21 outliers final: 15 residues processed: 133 average time/residue: 0.2283 time to fit residues: 41.4917 Evaluate side-chains 129 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 114 time to evaluate : 0.926 Switching outliers to nearest non-outliers outliers start: 15 outliers final: 0 residues processed: 15 average time/residue: 0.0813 time to fit residues: 3.6174 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 58 optimal weight: 0.9990 chunk 1 optimal weight: 3.9990 chunk 76 optimal weight: 1.9990 chunk 42 optimal weight: 0.5980 chunk 87 optimal weight: 0.0010 chunk 70 optimal weight: 3.9990 chunk 0 optimal weight: 0.9990 chunk 52 optimal weight: 2.9990 chunk 91 optimal weight: 0.0870 chunk 25 optimal weight: 0.6980 chunk 34 optimal weight: 0.0170 overall best weight: 0.2802 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: E 19 ASN ** F 239 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 171 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** H 398 GLN ** H 471 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 657 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8320 moved from start: 0.3314 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.033 9057 Z= 0.144 Angle : 0.520 9.955 12323 Z= 0.266 Chirality : 0.041 0.161 1307 Planarity : 0.004 0.067 1526 Dihedral : 11.176 173.078 1332 Min Nonbonded Distance : 1.964 Molprobity Statistics. All-atom Clashscore : 6.62 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.15 % Favored : 96.85 % Rotamer Outliers : 0.80 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % 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.25), residues: 1049 helix: -0.67 (0.26), residues: 421 sheet: -1.86 (0.41), residues: 164 loop : -2.07 (0.27), residues: 464 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 126 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 119 time to evaluate : 1.048 Fit side-chains outliers start: 7 outliers final: 2 residues processed: 125 average time/residue: 0.2558 time to fit residues: 43.7239 Evaluate side-chains 108 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 106 time to evaluate : 0.920 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0825 time to fit residues: 1.6258 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 92 optimal weight: 0.1980 chunk 20 optimal weight: 2.9990 chunk 60 optimal weight: 1.9990 chunk 25 optimal weight: 0.7980 chunk 102 optimal weight: 0.7980 chunk 85 optimal weight: 0.0980 chunk 47 optimal weight: 0.9990 chunk 8 optimal weight: 0.4980 chunk 33 optimal weight: 2.9990 chunk 53 optimal weight: 0.6980 chunk 98 optimal weight: 0.9980 overall best weight: 0.4580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: H 150 GLN ** H 171 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 471 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 657 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8331 moved from start: 0.3453 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.050 9057 Z= 0.172 Angle : 0.518 8.590 12323 Z= 0.265 Chirality : 0.041 0.161 1307 Planarity : 0.004 0.067 1526 Dihedral : 11.045 172.818 1332 Min Nonbonded Distance : 1.923 Molprobity Statistics. All-atom Clashscore : 6.33 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.72 % Favored : 96.28 % Rotamer Outliers : 0.69 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % 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.98 (0.26), residues: 1049 helix: -0.43 (0.26), residues: 421 sheet: -1.74 (0.43), residues: 152 loop : -1.92 (0.27), residues: 476 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 129 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 123 time to evaluate : 1.037 Fit side-chains outliers start: 6 outliers final: 1 residues processed: 127 average time/residue: 0.2382 time to fit residues: 41.7162 Evaluate side-chains 109 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 108 time to evaluate : 1.009 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0826 time to fit residues: 1.5759 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 11 optimal weight: 1.9990 chunk 58 optimal weight: 0.6980 chunk 74 optimal weight: 0.5980 chunk 86 optimal weight: 0.8980 chunk 57 optimal weight: 0.6980 chunk 102 optimal weight: 0.5980 chunk 63 optimal weight: 0.5980 chunk 62 optimal weight: 0.6980 chunk 47 optimal weight: 0.9990 chunk 40 optimal weight: 0.0000 chunk 61 optimal weight: 1.9990 overall best weight: 0.4984 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** H 171 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 471 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 657 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8330 moved from start: 0.3572 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.047 9057 Z= 0.182 Angle : 0.526 8.692 12323 Z= 0.270 Chirality : 0.042 0.157 1307 Planarity : 0.004 0.065 1526 Dihedral : 11.056 176.688 1332 Min Nonbonded Distance : 1.919 Molprobity Statistics. All-atom Clashscore : 6.22 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.34 % Favored : 96.66 % Rotamer Outliers : 0.92 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % 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.79 (0.26), residues: 1049 helix: -0.26 (0.26), residues: 420 sheet: -1.63 (0.43), residues: 152 loop : -1.84 (0.27), residues: 477 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 123 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 115 time to evaluate : 1.068 Fit side-chains outliers start: 8 outliers final: 4 residues processed: 120 average time/residue: 0.2593 time to fit residues: 41.9095 Evaluate side-chains 111 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 107 time to evaluate : 0.997 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.0846 time to fit residues: 2.0365 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 30 optimal weight: 0.9990 chunk 20 optimal weight: 0.6980 chunk 19 optimal weight: 0.8980 chunk 64 optimal weight: 0.9980 chunk 69 optimal weight: 1.9990 chunk 50 optimal weight: 1.9990 chunk 9 optimal weight: 2.9990 chunk 80 optimal weight: 0.9990 chunk 92 optimal weight: 0.7980 chunk 97 optimal weight: 0.8980 chunk 89 optimal weight: 0.8980 overall best weight: 0.8380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: E 19 ASN ** H 171 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 471 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 657 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8386 moved from start: 0.3653 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.054 9057 Z= 0.262 Angle : 0.584 9.730 12323 Z= 0.297 Chirality : 0.044 0.152 1307 Planarity : 0.004 0.064 1526 Dihedral : 11.266 179.447 1332 Min Nonbonded Distance : 1.897 Molprobity Statistics. All-atom Clashscore : 6.68 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.15 % Favored : 94.85 % Rotamer Outliers : 1.49 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % 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.89 (0.26), residues: 1049 helix: -0.36 (0.26), residues: 421 sheet: -1.56 (0.45), residues: 141 loop : -1.93 (0.27), residues: 487 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 125 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 112 time to evaluate : 0.949 Fit side-chains outliers start: 13 outliers final: 8 residues processed: 122 average time/residue: 0.2426 time to fit residues: 39.5697 Evaluate side-chains 117 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 109 time to evaluate : 0.921 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 0 residues processed: 8 average time/residue: 0.1188 time to fit residues: 2.7587 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 95 optimal weight: 0.0570 chunk 97 optimal weight: 0.8980 chunk 57 optimal weight: 1.9990 chunk 41 optimal weight: 0.0970 chunk 74 optimal weight: 0.9990 chunk 29 optimal weight: 0.7980 chunk 86 optimal weight: 0.2980 chunk 90 optimal weight: 0.8980 chunk 94 optimal weight: 0.9980 chunk 62 optimal weight: 0.5980 chunk 100 optimal weight: 0.7980 overall best weight: 0.3696 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: E 19 ASN ** H 171 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 471 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 657 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8322 moved from start: 0.3750 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.043 9057 Z= 0.159 Angle : 0.539 9.698 12323 Z= 0.271 Chirality : 0.042 0.170 1307 Planarity : 0.004 0.063 1526 Dihedral : 10.993 179.207 1332 Min Nonbonded Distance : 1.938 Molprobity Statistics. All-atom Clashscore : 6.80 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.53 % Favored : 96.47 % Rotamer Outliers : 0.46 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % 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.54 (0.26), residues: 1049 helix: -0.02 (0.27), residues: 415 sheet: -1.54 (0.43), residues: 152 loop : -1.70 (0.28), residues: 482 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 121 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 117 time to evaluate : 0.920 Fit side-chains outliers start: 4 outliers final: 3 residues processed: 120 average time/residue: 0.2504 time to fit residues: 40.4593 Evaluate side-chains 111 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 108 time to evaluate : 1.083 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0887 time to fit residues: 1.9609 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 61 optimal weight: 2.9990 chunk 47 optimal weight: 0.7980 chunk 69 optimal weight: 1.9990 chunk 105 optimal weight: 2.9990 chunk 97 optimal weight: 0.9990 chunk 84 optimal weight: 0.8980 chunk 8 optimal weight: 0.9980 chunk 64 optimal weight: 0.0570 chunk 51 optimal weight: 0.8980 chunk 66 optimal weight: 1.9990 chunk 89 optimal weight: 0.8980 overall best weight: 0.7098 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: E 169 GLN E 239 ASN ** H 171 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 471 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 657 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8361 moved from start: 0.3785 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.047 9057 Z= 0.228 Angle : 0.568 9.850 12323 Z= 0.288 Chirality : 0.043 0.176 1307 Planarity : 0.004 0.063 1526 Dihedral : 11.044 178.687 1332 Min Nonbonded Distance : 1.917 Molprobity Statistics. All-atom Clashscore : 7.15 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.15 % Favored : 94.85 % Rotamer Outliers : 0.34 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % 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.26), residues: 1049 helix: -0.09 (0.27), residues: 414 sheet: -1.52 (0.43), residues: 152 loop : -1.74 (0.28), residues: 483 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2098 Ramachandran restraints generated. 1049 Oldfield, 0 Emsley, 1049 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 112 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 109 time to evaluate : 0.991 Fit side-chains outliers start: 3 outliers final: 0 residues processed: 111 average time/residue: 0.2677 time to fit residues: 40.2187 Evaluate side-chains 108 residues out of total 881 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 108 time to evaluate : 1.148 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.6364 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 106 random chunks: chunk 25 optimal weight: 1.9990 chunk 77 optimal weight: 0.0870 chunk 12 optimal weight: 0.7980 chunk 23 optimal weight: 3.9990 chunk 84 optimal weight: 0.6980 chunk 35 optimal weight: 1.9990 chunk 86 optimal weight: 0.9990 chunk 10 optimal weight: 3.9990 chunk 15 optimal weight: 0.4980 chunk 73 optimal weight: 2.9990 chunk 4 optimal weight: 0.9990 overall best weight: 0.6160 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** H 171 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 471 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 657 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3756 r_free = 0.3756 target = 0.154125 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 70)----------------| | r_work = 0.3572 r_free = 0.3572 target = 0.139047 restraints weight = 12451.522| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 34)----------------| | r_work = 0.3590 r_free = 0.3590 target = 0.140661 restraints weight = 8745.680| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 47)----------------| | r_work = 0.3606 r_free = 0.3606 target = 0.141888 restraints weight = 6824.112| |-----------------------------------------------------------------------------| r_work (final): 0.3661 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8311 moved from start: 0.3840 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.045 9057 Z= 0.206 Angle : 0.555 10.059 12323 Z= 0.281 Chirality : 0.043 0.169 1307 Planarity : 0.004 0.062 1526 Dihedral : 10.986 178.640 1332 Min Nonbonded Distance : 1.944 Molprobity Statistics. All-atom Clashscore : 7.26 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.10 % Favored : 95.90 % Rotamer Outliers : 0.11 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.23 % 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.26), residues: 1049 helix: -0.08 (0.27), residues: 416 sheet: -1.40 (0.44), residues: 147 loop : -1.74 (0.27), residues: 486 =============================================================================== Job complete usr+sys time: 1932.89 seconds wall clock time: 36 minutes 4.15 seconds (2164.15 seconds total)