Starting phenix.real_space_refine (version: dev) on Mon Dec 12 16:32:49 2022 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7kve_23048/12_2022/7kve_23048.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7kve_23048/12_2022/7kve_23048.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.3 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7kve_23048/12_2022/7kve_23048.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7kve_23048/12_2022/7kve_23048.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7kve_23048/12_2022/7kve_23048.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7kve_23048/12_2022/7kve_23048.pdb" } resolution = 3.3 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.067 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 Convert atoms to be neutral Process input model Symmetric amino acids flipped Residue "B ARG 18": "NH1" <-> "NH2" Residue "B ARG 90": "NH1" <-> "NH2" Residue "B ARG 118": "NH1" <-> "NH2" Residue "B ARG 306": "NH1" <-> "NH2" Residue "B ARG 313": "NH1" <-> "NH2" Residue "B ARG 316": "NH1" <-> "NH2" Residue "B ARG 317": "NH1" <-> "NH2" Residue "B ARG 321": "NH1" <-> "NH2" Residue "B ARG 348": "NH1" <-> "NH2" Residue "B ARG 396": "NH1" <-> "NH2" Residue "B ARG 439": "NH1" <-> "NH2" Residue "B ARG 445": "NH1" <-> "NH2" Residue "B ARG 475": "NH1" <-> "NH2" Residue "B ARG 501": "NH1" <-> "NH2" Residue "B ARG 505": "NH1" <-> "NH2" Residue "B ARG 506": "NH1" <-> "NH2" Residue "B ARG 510": "NH1" <-> "NH2" Residue "B ARG 547": "NH1" <-> "NH2" Residue "B PHE 603": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 608": "NH1" <-> "NH2" Residue "B ARG 619": "NH1" <-> "NH2" Residue "B ARG 643": "NH1" <-> "NH2" Residue "B ARG 652": "NH1" <-> "NH2" Residue "B GLU 662": "OE1" <-> "OE2" Residue "B ARG 679": "NH1" <-> "NH2" Residue "B ARG 684": "NH1" <-> "NH2" Residue "B ARG 701": "NH1" <-> "NH2" Residue "B ARG 709": "NH1" <-> "NH2" Residue "B ARG 712": "NH1" <-> "NH2" Residue "B ARG 1545": "NH1" <-> "NH2" Residue "B ARG 1551": "NH1" <-> "NH2" Residue "B ARG 1571": "NH1" <-> "NH2" Residue "B GLU 1572": "OE1" <-> "OE2" Residue "B GLU 1576": "OE1" <-> "OE2" Residue "B ARG 1603": "NH1" <-> "NH2" Residue "B ARG 1606": "NH1" <-> "NH2" Residue "B ARG 1631": "NH1" <-> "NH2" Residue "B GLU 1691": "OE1" <-> "OE2" Residue "B ARG 1698": "NH1" <-> "NH2" Residue "B GLU 1750": "OE1" <-> "OE2" Residue "B ARG 1761": "NH1" <-> "NH2" Residue "B ARG 1765": "NH1" <-> "NH2" Residue "B ARG 1799": "NH1" <-> "NH2" Residue "B ARG 1865": "NH1" <-> "NH2" Residue "B ARG 1877": "NH1" <-> "NH2" Residue "B ARG 1880": "NH1" <-> "NH2" Residue "B ARG 1907": "NH1" <-> "NH2" Residue "B GLU 1923": "OE1" <-> "OE2" Residue "B ARG 1985": "NH1" <-> "NH2" Residue "B TYR 2012": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 2019": "NH1" <-> "NH2" Residue "B GLU 2034": "OE1" <-> "OE2" Residue "B ARG 2072": "NH1" <-> "NH2" Residue "B ARG 2080": "NH1" <-> "NH2" Residue "B ARG 2140": "NH1" <-> "NH2" Residue "B ARG 2171": "NH1" <-> "NH2" Time to flip residues: 0.04s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/chem_data/mon_lib" Total number of atoms: 11142 Number of models: 1 Model: "" Number of chains: 1 Chain: "B" Number of atoms: 11142 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1374, 11142 Classifications: {'peptide': 1374} Modifications used: {'COO': 1} Link IDs: {'CIS': 1, 'PCIS': 1, 'PTRANS': 66, 'TRANS': 1305} Chain breaks: 1 Time building chain proxies: 6.44, per 1000 atoms: 0.58 Number of scatterers: 11142 At special positions: 0 Unit cell: (100.1, 106.7, 141.9, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 4 Type Number sf(0) S 55 16.00 O 2108 8.00 N 1905 7.00 C 7074 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=6, symmetry=0 Simple disulfide: pdb=" SG CYS B 139 " - pdb=" SG CYS B 165 " distance=2.04 Simple disulfide: pdb=" SG CYS B 220 " - pdb=" SG CYS B 301 " distance=2.03 Simple disulfide: pdb=" SG CYS B 472 " - pdb=" SG CYS B 498 " distance=2.03 Simple disulfide: pdb=" SG CYS B1697 " - pdb=" SG CYS B1723 " distance=2.03 Simple disulfide: pdb=" SG CYS B1879 " - pdb=" SG CYS B2033 " distance=2.04 Simple disulfide: pdb=" SG CYS B2038 " - pdb=" SG CYS B2193 " distance=2.04 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: 4.11 Conformation dependent library (CDL) restraints added in 1.8 seconds 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2588 Finding SS restraints... Secondary structure from input PDB file: 17 helices and 21 sheets defined 6.8% alpha, 23.7% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.83 Creating SS restraints... Processing helix chain 'B' and resid 23 through 28 removed outlier: 4.475A pdb=" N LEU B 28 " --> pdb=" O SER B 24 " (cutoff:3.500A) Processing helix chain 'B' and resid 42 through 46 removed outlier: 3.904A pdb=" N LYS B 46 " --> pdb=" O PRO B 43 " (cutoff:3.500A) Processing helix chain 'B' and resid 191 through 195 Processing helix chain 'B' and resid 344 through 352 removed outlier: 3.948A pdb=" N ARG B 348 " --> pdb=" O ASP B 344 " (cutoff:3.500A) removed outlier: 3.809A pdb=" N SER B 349 " --> pdb=" O LYS B 345 " (cutoff:3.500A) Processing helix chain 'B' and resid 482 through 490 Processing helix chain 'B' and resid 526 through 528 No H-bonds generated for 'chain 'B' and resid 526 through 528' Processing helix chain 'B' and resid 549 through 556 removed outlier: 3.779A pdb=" N ASN B 556 " --> pdb=" O PHE B 552 " (cutoff:3.500A) Processing helix chain 'B' and resid 1608 through 1612 removed outlier: 4.352A pdb=" N HIS B1612 " --> pdb=" O TYR B1609 " (cutoff:3.500A) Processing helix chain 'B' and resid 1650 through 1654 removed outlier: 3.612A pdb=" N GLU B1654 " --> pdb=" O LYS B1651 " (cutoff:3.500A) Processing helix chain 'B' and resid 1665 through 1669 removed outlier: 3.606A pdb=" N GLU B1668 " --> pdb=" O TRP B1665 " (cutoff:3.500A) removed outlier: 4.175A pdb=" N ASP B1669 " --> pdb=" O PHE B1666 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 1665 through 1669' Processing helix chain 'B' and resid 1685 through 1689 Processing helix chain 'B' and resid 1707 through 1713 removed outlier: 3.863A pdb=" N ILE B1712 " --> pdb=" O PRO B1708 " (cutoff:3.500A) removed outlier: 3.711A pdb=" N HIS B1713 " --> pdb=" O GLU B1709 " (cutoff:3.500A) Processing helix chain 'B' and resid 1768 through 1772 removed outlier: 4.038A pdb=" N LYS B1772 " --> pdb=" O SER B1769 " (cutoff:3.500A) Processing helix chain 'B' and resid 1860 through 1865 removed outlier: 3.833A pdb=" N ARG B1865 " --> pdb=" O GLY B1861 " (cutoff:3.500A) Processing helix chain 'B' and resid 1891 through 1893 No H-bonds generated for 'chain 'B' and resid 1891 through 1893' Processing helix chain 'B' and resid 1905 through 1909 Processing helix chain 'B' and resid 2050 through 2052 No H-bonds generated for 'chain 'B' and resid 2050 through 2052' Processing sheet with id=AA1, first strand: chain 'B' and resid 36 through 40 removed outlier: 3.670A pdb=" N GLN B 11 " --> pdb=" O VAL B 37 " (cutoff:3.500A) removed outlier: 7.235A pdb=" N ARG B 4 " --> pdb=" O LYS B 71 " (cutoff:3.500A) removed outlier: 7.440A pdb=" N HIS B 73 " --> pdb=" O ARG B 4 " (cutoff:3.500A) removed outlier: 7.143A pdb=" N PHE B 6 " --> pdb=" O HIS B 73 " (cutoff:3.500A) removed outlier: 3.544A pdb=" N ILE B 70 " --> pdb=" O TRP B 124 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'B' and resid 62 through 64 removed outlier: 6.768A pdb=" N LEU B 62 " --> pdb=" O LEU B 163 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'B' and resid 75 through 76 removed outlier: 3.519A pdb=" N ASN B 76 " --> pdb=" O ARG B 118 " (cutoff:3.500A) removed outlier: 4.017A pdb=" N ARG B 118 " --> pdb=" O ASN B 76 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA3 Processing sheet with id=AA4, first strand: chain 'B' and resid 185 through 189 removed outlier: 6.339A pdb=" N LEU B 185 " --> pdb=" O LEU B 230 " (cutoff:3.500A) removed outlier: 7.538A pdb=" N MET B 232 " --> pdb=" O LEU B 185 " (cutoff:3.500A) removed outlier: 6.650A pdb=" N PHE B 187 " --> pdb=" O MET B 232 " (cutoff:3.500A) removed outlier: 3.570A pdb=" N MET B 270 " --> pdb=" O ILE B 225 " (cutoff:3.500A) removed outlier: 3.562A pdb=" N ALA B 268 " --> pdb=" O TRP B 227 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'B' and resid 217 through 219 removed outlier: 3.521A pdb=" N ILE B 217 " --> pdb=" O TYR B 295 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'B' and resid 370 through 371 removed outlier: 4.167A pdb=" N PHE B 325 " --> pdb=" O TYR B 371 " (cutoff:3.500A) removed outlier: 6.282A pdb=" N TYR B 324 " --> pdb=" O VAL B 406 " (cutoff:3.500A) removed outlier: 7.819A pdb=" N LYS B 408 " --> pdb=" O TYR B 324 " (cutoff:3.500A) removed outlier: 7.453A pdb=" N ILE B 326 " --> pdb=" O LYS B 408 " (cutoff:3.500A) removed outlier: 3.665A pdb=" N PHE B 407 " --> pdb=" O TYR B 453 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 394 through 397 Processing sheet with id=AA8, first strand: chain 'B' and resid 515 through 520 removed outlier: 3.695A pdb=" N CYS B 585 " --> pdb=" O ALA B 518 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'B' and resid 522 through 524 Processing sheet with id=AB1, first strand: chain 'B' and resid 612 through 615 removed outlier: 3.815A pdb=" N LEU B 613 " --> pdb=" O ILE B 596 " (cutoff:3.500A) removed outlier: 3.640A pdb=" N LEU B 594 " --> pdb=" O LEU B 615 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'B' and resid 1550 through 1554 removed outlier: 5.491A pdb=" N ASN B1550 " --> pdb=" O VAL B1627 " (cutoff:3.500A) removed outlier: 7.847A pdb=" N GLN B1629 " --> pdb=" O ASN B1550 " (cutoff:3.500A) removed outlier: 6.199A pdb=" N ARG B1552 " --> pdb=" O GLN B1629 " (cutoff:3.500A) removed outlier: 7.845A pdb=" N ARG B1631 " --> pdb=" O ARG B1552 " (cutoff:3.500A) removed outlier: 7.424A pdb=" N TYR B1554 " --> pdb=" O ARG B1631 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'B' and resid 1556 through 1563 removed outlier: 3.907A pdb=" N ALA B1557 " --> pdb=" O ARG B1593 " (cutoff:3.500A) removed outlier: 3.767A pdb=" N ARG B1593 " --> pdb=" O ALA B1557 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'B' and resid 1697 through 1698 Processing sheet with id=AB5, first strand: chain 'B' and resid 1701 through 1703 Processing sheet with id=AB6, first strand: chain 'B' and resid 1776 through 1780 removed outlier: 3.513A pdb=" N MET B1746 " --> pdb=" O ALA B1779 " (cutoff:3.500A) Processing sheet with id=AB7, first strand: chain 'B' and resid 1791 through 1792 Processing sheet with id=AB8, first strand: chain 'B' and resid 1834 through 1836 Processing sheet with id=AB9, first strand: chain 'B' and resid 1881 through 1882 removed outlier: 6.390A pdb=" N GLN B2031 " --> pdb=" O THR B1946 " (cutoff:3.500A) removed outlier: 3.923A pdb=" N VAL B1943 " --> pdb=" O ALA B2010 " (cutoff:3.500A) removed outlier: 7.930A pdb=" N GLY B1947 " --> pdb=" O PRO B2006 " (cutoff:3.500A) removed outlier: 11.882A pdb=" N GLN B1949 " --> pdb=" O PHE B2004 " (cutoff:3.500A) removed outlier: 12.332A pdb=" N PHE B2004 " --> pdb=" O GLN B1949 " (cutoff:3.500A) Processing sheet with id=AC1, first strand: chain 'B' and resid 1895 through 1897 removed outlier: 6.746A pdb=" N GLU B1964 " --> pdb=" O THR B2018 " (cutoff:3.500A) removed outlier: 6.309A pdb=" N ALA B2020 " --> pdb=" O THR B1962 " (cutoff:3.500A) removed outlier: 6.213A pdb=" N THR B1962 " --> pdb=" O ALA B2020 " (cutoff:3.500A) Processing sheet with id=AC2, first strand: chain 'B' and resid 2040 through 2041 removed outlier: 3.788A pdb=" N THR B2040 " --> pdb=" O GLY B2192 " (cutoff:3.500A) removed outlier: 3.531A pdb=" N GLY B2192 " --> pdb=" O THR B2040 " (cutoff:3.500A) removed outlier: 6.912A pdb=" N GLU B2189 " --> pdb=" O ILE B2108 " (cutoff:3.500A) removed outlier: 4.637A pdb=" N ILE B2108 " --> pdb=" O GLU B2189 " (cutoff:3.500A) removed outlier: 6.620A pdb=" N PHE B2191 " --> pdb=" O THR B2106 " (cutoff:3.500A) removed outlier: 6.999A pdb=" N ARG B2171 " --> pdb=" O LEU B2099 " (cutoff:3.500A) removed outlier: 10.342A pdb=" N LYS B2101 " --> pdb=" O ILE B2169 " (cutoff:3.500A) removed outlier: 9.333A pdb=" N ILE B2169 " --> pdb=" O LYS B2101 " (cutoff:3.500A) removed outlier: 9.659A pdb=" N LYS B2103 " --> pdb=" O PRO B2167 " (cutoff:3.500A) removed outlier: 10.470A pdb=" N ILE B2105 " --> pdb=" O ASN B2165 " (cutoff:3.500A) removed outlier: 7.759A pdb=" N ASN B2165 " --> pdb=" O ILE B2105 " (cutoff:3.500A) removed outlier: 7.277A pdb=" N ALA B2107 " --> pdb=" O PHE B2163 " (cutoff:3.500A) removed outlier: 7.174A pdb=" N PHE B2163 " --> pdb=" O ALA B2107 " (cutoff:3.500A) removed outlier: 7.033A pdb=" N ILE B2109 " --> pdb=" O LYS B2161 " (cutoff:3.500A) removed outlier: 7.176A pdb=" N LYS B2161 " --> pdb=" O ILE B2109 " (cutoff:3.500A) removed outlier: 6.930A pdb=" N SER B2124 " --> pdb=" O LYS B2178 " (cutoff:3.500A) removed outlier: 6.488A pdb=" N TRP B2180 " --> pdb=" O VAL B2122 " (cutoff:3.500A) removed outlier: 6.649A pdb=" N VAL B2122 " --> pdb=" O TRP B2180 " (cutoff:3.500A) removed outlier: 3.785A pdb=" N MET B2120 " --> pdb=" O CYS B2113 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N CYS B2113 " --> pdb=" O MET B2120 " (cutoff:3.500A) Processing sheet with id=AC3, first strand: chain 'B' and resid 2054 through 2056 removed outlier: 6.999A pdb=" N ARG B2171 " --> pdb=" O LEU B2099 " (cutoff:3.500A) removed outlier: 10.342A pdb=" N LYS B2101 " --> pdb=" O ILE B2169 " (cutoff:3.500A) removed outlier: 9.333A pdb=" N ILE B2169 " --> pdb=" O LYS B2101 " (cutoff:3.500A) removed outlier: 9.659A pdb=" N LYS B2103 " --> pdb=" O PRO B2167 " (cutoff:3.500A) removed outlier: 10.470A pdb=" N ILE B2105 " --> pdb=" O ASN B2165 " (cutoff:3.500A) removed outlier: 7.759A pdb=" N ASN B2165 " --> pdb=" O ILE B2105 " (cutoff:3.500A) removed outlier: 7.277A pdb=" N ALA B2107 " --> pdb=" O PHE B2163 " (cutoff:3.500A) removed outlier: 7.174A pdb=" N PHE B2163 " --> pdb=" O ALA B2107 " (cutoff:3.500A) removed outlier: 7.033A pdb=" N ILE B2109 " --> pdb=" O LYS B2161 " (cutoff:3.500A) removed outlier: 7.176A pdb=" N LYS B2161 " --> pdb=" O ILE B2109 " (cutoff:3.500A) removed outlier: 6.930A pdb=" N SER B2124 " --> pdb=" O LYS B2178 " (cutoff:3.500A) removed outlier: 6.488A pdb=" N TRP B2180 " --> pdb=" O VAL B2122 " (cutoff:3.500A) removed outlier: 6.649A pdb=" N VAL B2122 " --> pdb=" O TRP B2180 " (cutoff:3.500A) 160 hydrogen bonds defined for protein. 384 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 0 stacking parallelities Total time for adding SS restraints: 4.29 Time building geometry restraints manager: 4.53 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.15 - 1.29: 1844 1.29 - 1.42: 3119 1.42 - 1.56: 6380 1.56 - 1.70: 1 1.70 - 1.83: 92 Bond restraints: 11436 Sorted by residual: bond pdb=" C VAL B 654 " pdb=" O VAL B 654 " ideal model delta sigma weight residual 1.231 1.153 0.078 1.29e-02 6.01e+03 3.63e+01 bond pdb=" C CYS B 656 " pdb=" O CYS B 656 " ideal model delta sigma weight residual 1.235 1.188 0.047 1.26e-02 6.30e+03 1.38e+01 bond pdb=" C ASP B 653 " pdb=" N VAL B 654 " ideal model delta sigma weight residual 1.331 1.290 0.041 1.20e-02 6.94e+03 1.19e+01 bond pdb=" C ASN B 541 " pdb=" N PRO B 542 " ideal model delta sigma weight residual 1.331 1.359 -0.027 7.90e-03 1.60e+04 1.18e+01 bond pdb=" CA VAL B 654 " pdb=" C VAL B 654 " ideal model delta sigma weight residual 1.526 1.477 0.049 1.43e-02 4.89e+03 1.16e+01 ... (remaining 11431 not shown) Histogram of bond angle deviations from ideal: 99.92 - 107.30: 500 107.30 - 114.69: 6526 114.69 - 122.07: 5929 122.07 - 129.45: 2423 129.45 - 136.84: 112 Bond angle restraints: 15490 Sorted by residual: angle pdb=" N VAL B 654 " pdb=" CA VAL B 654 " pdb=" CB VAL B 654 " ideal model delta sigma weight residual 112.60 100.51 12.09 1.29e+00 6.01e-01 8.78e+01 angle pdb=" N PHE B2059 " pdb=" CA PHE B2059 " pdb=" C PHE B2059 " ideal model delta sigma weight residual 112.72 103.01 9.71 1.36e+00 5.41e-01 5.10e+01 angle pdb=" N VAL B1706 " pdb=" CA VAL B1706 " pdb=" C VAL B1706 " ideal model delta sigma weight residual 112.12 107.30 4.82 8.40e-01 1.42e+00 3.29e+01 angle pdb=" N ASP B 401 " pdb=" CA ASP B 401 " pdb=" C ASP B 401 " ideal model delta sigma weight residual 112.57 106.64 5.93 1.13e+00 7.83e-01 2.75e+01 angle pdb=" C ASN B2036 " pdb=" N GLY B2037 " pdb=" CA GLY B2037 " ideal model delta sigma weight residual 121.51 129.34 -7.83 1.52e+00 4.33e-01 2.65e+01 ... (remaining 15485 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.94: 6070 17.94 - 35.88: 666 35.88 - 53.82: 92 53.82 - 71.76: 20 71.76 - 89.70: 8 Dihedral angle restraints: 6856 sinusoidal: 2828 harmonic: 4028 Sorted by residual: dihedral pdb=" CB CYS B 139 " pdb=" SG CYS B 139 " pdb=" SG CYS B 165 " pdb=" CB CYS B 165 " ideal model delta sinusoidal sigma weight residual 93.00 20.33 72.67 1 1.00e+01 1.00e-02 6.74e+01 dihedral pdb=" N VAL B 654 " pdb=" C VAL B 654 " pdb=" CA VAL B 654 " pdb=" CB VAL B 654 " ideal model delta harmonic sigma weight residual 123.40 108.17 15.23 0 2.50e+00 1.60e-01 3.71e+01 dihedral pdb=" CA ARG B 90 " pdb=" C ARG B 90 " pdb=" N TYR B 91 " pdb=" CA TYR B 91 " ideal model delta harmonic sigma weight residual -180.00 -155.50 -24.50 0 5.00e+00 4.00e-02 2.40e+01 ... (remaining 6853 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.100: 1519 0.100 - 0.199: 101 0.199 - 0.299: 3 0.299 - 0.398: 3 0.398 - 0.498: 1 Chirality restraints: 1627 Sorted by residual: chirality pdb=" CB VAL B 654 " pdb=" CA VAL B 654 " pdb=" CG1 VAL B 654 " pdb=" CG2 VAL B 654 " both_signs ideal model delta sigma weight residual False -2.63 -2.13 -0.50 2.00e-01 2.50e+01 6.19e+00 chirality pdb=" CA ASP B 660 " pdb=" N ASP B 660 " pdb=" C ASP B 660 " pdb=" CB ASP B 660 " both_signs ideal model delta sigma weight residual False 2.51 2.15 0.36 2.00e-01 2.50e+01 3.28e+00 chirality pdb=" CA VAL B 654 " pdb=" N VAL B 654 " pdb=" C VAL B 654 " pdb=" CB VAL B 654 " both_signs ideal model delta sigma weight residual False 2.44 2.79 -0.35 2.00e-01 2.50e+01 2.99e+00 ... (remaining 1624 not shown) Planarity restraints: 1999 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA CYS B2033 " -0.021 2.00e-02 2.50e+03 4.13e-02 1.71e+01 pdb=" C CYS B2033 " 0.071 2.00e-02 2.50e+03 pdb=" O CYS B2033 " -0.026 2.00e-02 2.50e+03 pdb=" N GLU B2034 " -0.025 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C GLU B 669 " 0.060 5.00e-02 4.00e+02 9.13e-02 1.33e+01 pdb=" N PRO B 670 " -0.158 5.00e-02 4.00e+02 pdb=" CA PRO B 670 " 0.048 5.00e-02 4.00e+02 pdb=" CD PRO B 670 " 0.050 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C ARG B2023 " -0.033 5.00e-02 4.00e+02 4.94e-02 3.91e+00 pdb=" N PRO B2024 " 0.085 5.00e-02 4.00e+02 pdb=" CA PRO B2024 " -0.025 5.00e-02 4.00e+02 pdb=" CD PRO B2024 " -0.028 5.00e-02 4.00e+02 ... (remaining 1996 not shown) Histogram of nonbonded interaction distances: 2.20 - 2.74: 1223 2.74 - 3.28: 11124 3.28 - 3.82: 17231 3.82 - 4.36: 20716 4.36 - 4.90: 36116 Nonbonded interactions: 86410 Sorted by model distance: nonbonded pdb=" OD1 ASP B 191 " pdb=" OG SER B 193 " model vdw 2.204 2.440 nonbonded pdb=" OH TYR B1609 " pdb=" OD2 ASP B2005 " model vdw 2.204 2.440 nonbonded pdb=" O SER B 349 " pdb=" NE2 GLN B 358 " model vdw 2.204 2.520 nonbonded pdb=" OH TYR B1554 " pdb=" O ARG B1621 " model vdw 2.212 2.440 nonbonded pdb=" OH TYR B 605 " pdb=" OE2 GLU B 621 " model vdw 2.214 2.440 ... (remaining 86405 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: 4 Type Number sf(0) Gaussians S 55 5.16 5 C 7074 2.51 5 N 1905 2.21 5 O 2108 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 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: 0.950 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.000 Extract box with map and model: 6.370 Check model and map are aligned: 0.170 Convert atoms to be neutral: 0.110 Process input model: 32.050 Find NCS groups from input model: 0.340 Set up NCS constraints: 0.050 Set refine NCS operators: 0.000 Set scattering table: 0.010 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.010 Load rotamer database and sin/cos tables:3.030 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 43.090 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7266 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.005 0.078 11436 Z= 0.319 Angle : 0.817 12.087 15490 Z= 0.482 Chirality : 0.053 0.498 1627 Planarity : 0.005 0.091 1999 Dihedral : 14.673 89.697 4250 Min Nonbonded Distance : 2.204 Molprobity Statistics. All-atom Clashscore : 17.50 Ramachandran Plot: Outliers : 0.36 % Allowed : 18.25 % Favored : 81.39 % Rotamer Outliers : 0.57 % Cbeta Deviations : 0.15 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % 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: -5.88 (0.17), residues: 1370 helix: -4.25 (0.41), residues: 45 sheet: -3.39 (0.23), residues: 373 loop : -4.44 (0.14), residues: 952 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 160 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 153 time to evaluate : 1.377 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 7 outliers final: 1 residues processed: 158 average time/residue: 0.2616 time to fit residues: 57.5418 Evaluate side-chains 115 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 114 time to evaluate : 1.297 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1059 time to fit residues: 1.9870 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 137 random chunks: chunk 115 optimal weight: 3.9990 chunk 103 optimal weight: 4.9990 chunk 57 optimal weight: 0.9990 chunk 35 optimal weight: 8.9990 chunk 70 optimal weight: 9.9990 chunk 55 optimal weight: 7.9990 chunk 107 optimal weight: 6.9990 chunk 41 optimal weight: 0.6980 chunk 65 optimal weight: 4.9990 chunk 79 optimal weight: 4.9990 chunk 124 optimal weight: 0.9990 overall best weight: 2.3388 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 11 GLN B 134 HIS B 252 HIS B 307 ASN B 315 GLN ** B 350 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 398 GLN B 562 ASN ** B 583 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** B 597 HIS B 609 HIS B 639 ASN B 700 ASN ** B1570 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B1643 HIS B1645 HIS B1670 ASN B1713 HIS B1778 HIS ** B1781 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B1819 GLN B1936 GLN B2002 ASN B2003 GLN B2093 GLN B2128 HIS B2162 ASN Total number of N/Q/H flips: 23 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7308 moved from start: 0.1326 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.048 11436 Z= 0.351 Angle : 0.745 9.830 15490 Z= 0.390 Chirality : 0.049 0.168 1627 Planarity : 0.006 0.096 1999 Dihedral : 6.679 36.613 1508 Min Nonbonded Distance : 2.124 Molprobity Statistics. All-atom Clashscore : 15.86 Ramachandran Plot: Outliers : 0.00 % Allowed : 17.59 % Favored : 82.41 % Rotamer Outliers : 2.69 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % Twisted Proline : 0.00 % Twisted General : 0.08 % 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: -5.39 (0.18), residues: 1370 helix: -3.57 (0.51), residues: 50 sheet: -2.99 (0.25), residues: 385 loop : -4.14 (0.15), residues: 935 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 147 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 33 poor density : 114 time to evaluate : 1.345 Fit side-chains revert: symmetry clash outliers start: 33 outliers final: 20 residues processed: 141 average time/residue: 0.2178 time to fit residues: 45.5635 Evaluate side-chains 122 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 20 poor density : 102 time to evaluate : 1.215 Switching outliers to nearest non-outliers outliers start: 20 outliers final: 0 residues processed: 20 average time/residue: 0.1298 time to fit residues: 6.0577 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 137 random chunks: chunk 69 optimal weight: 0.0000 chunk 38 optimal weight: 0.7980 chunk 103 optimal weight: 4.9990 chunk 84 optimal weight: 0.4980 chunk 34 optimal weight: 0.3980 chunk 124 optimal weight: 1.9990 chunk 134 optimal weight: 1.9990 chunk 110 optimal weight: 2.9990 chunk 123 optimal weight: 2.9990 chunk 42 optimal weight: 9.9990 chunk 99 optimal weight: 0.1980 overall best weight: 0.3784 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 244 ASN B 289 GLN ** B 350 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 379 HIS B 398 GLN B 583 HIS ** B1570 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B1670 ASN ** B1781 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B1824 ASN B1828 GLN B1830 GLN B1936 GLN B1973 GLN Total number of N/Q/H flips: 11 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7179 moved from start: 0.1850 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.040 11436 Z= 0.156 Angle : 0.625 12.103 15490 Z= 0.318 Chirality : 0.046 0.151 1627 Planarity : 0.005 0.088 1999 Dihedral : 5.745 30.209 1508 Min Nonbonded Distance : 2.137 Molprobity Statistics. All-atom Clashscore : 11.53 Ramachandran Plot: Outliers : 0.00 % Allowed : 12.41 % Favored : 87.59 % Rotamer Outliers : 1.31 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % 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.83 (0.19), residues: 1370 helix: -3.20 (0.56), residues: 51 sheet: -2.50 (0.26), residues: 367 loop : -3.78 (0.16), residues: 952 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 142 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 16 poor density : 126 time to evaluate : 1.273 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 16 outliers final: 5 residues processed: 137 average time/residue: 0.2168 time to fit residues: 44.4509 Evaluate side-chains 118 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 113 time to evaluate : 1.347 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.1100 time to fit residues: 2.8082 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 137 random chunks: chunk 123 optimal weight: 4.9990 chunk 93 optimal weight: 0.7980 chunk 64 optimal weight: 10.0000 chunk 13 optimal weight: 0.9980 chunk 59 optimal weight: 4.9990 chunk 83 optimal weight: 5.9990 chunk 125 optimal weight: 9.9990 chunk 132 optimal weight: 0.9990 chunk 65 optimal weight: 0.8980 chunk 118 optimal weight: 3.9990 chunk 35 optimal weight: 8.9990 overall best weight: 1.5384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 244 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 350 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1570 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B1670 ASN ** B1781 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7257 moved from start: 0.1883 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.049 11436 Z= 0.259 Angle : 0.661 10.083 15490 Z= 0.341 Chirality : 0.047 0.184 1627 Planarity : 0.005 0.093 1999 Dihedral : 6.021 32.483 1508 Min Nonbonded Distance : 2.162 Molprobity Statistics. All-atom Clashscore : 14.13 Ramachandran Plot: Outliers : 0.00 % Allowed : 16.13 % Favored : 83.87 % Rotamer Outliers : 2.04 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % 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.72 (0.19), residues: 1370 helix: -3.28 (0.54), residues: 52 sheet: -2.34 (0.27), residues: 369 loop : -3.74 (0.16), residues: 949 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 141 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 25 poor density : 116 time to evaluate : 1.471 Fit side-chains revert: symmetry clash outliers start: 25 outliers final: 14 residues processed: 134 average time/residue: 0.2153 time to fit residues: 43.7222 Evaluate side-chains 122 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 108 time to evaluate : 1.344 Switching outliers to nearest non-outliers outliers start: 14 outliers final: 0 residues processed: 14 average time/residue: 0.1066 time to fit residues: 4.7411 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 137 random chunks: chunk 110 optimal weight: 4.9990 chunk 75 optimal weight: 3.9990 chunk 1 optimal weight: 3.9990 chunk 98 optimal weight: 9.9990 chunk 54 optimal weight: 5.9990 chunk 113 optimal weight: 0.7980 chunk 91 optimal weight: 6.9990 chunk 0 optimal weight: 4.9990 chunk 67 optimal weight: 0.0770 chunk 118 optimal weight: 0.0870 chunk 33 optimal weight: 0.7980 overall best weight: 1.1518 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 350 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1570 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1781 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B1869 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7233 moved from start: 0.2067 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.038 11436 Z= 0.215 Angle : 0.640 10.164 15490 Z= 0.328 Chirality : 0.046 0.154 1627 Planarity : 0.005 0.092 1999 Dihedral : 5.875 31.503 1508 Min Nonbonded Distance : 2.157 Molprobity Statistics. All-atom Clashscore : 13.31 Ramachandran Plot: Outliers : 0.00 % Allowed : 14.38 % Favored : 85.62 % Rotamer Outliers : 1.88 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % 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.55 (0.20), residues: 1370 helix: -3.14 (0.57), residues: 52 sheet: -2.21 (0.26), residues: 373 loop : -3.63 (0.17), residues: 945 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 143 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 23 poor density : 120 time to evaluate : 1.432 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 23 outliers final: 11 residues processed: 140 average time/residue: 0.2051 time to fit residues: 44.7436 Evaluate side-chains 121 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 110 time to evaluate : 1.286 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 0 residues processed: 11 average time/residue: 0.1135 time to fit residues: 4.3078 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 137 random chunks: chunk 44 optimal weight: 10.0000 chunk 119 optimal weight: 0.8980 chunk 26 optimal weight: 3.9990 chunk 77 optimal weight: 0.9980 chunk 32 optimal weight: 0.9980 chunk 132 optimal weight: 3.9990 chunk 110 optimal weight: 2.9990 chunk 61 optimal weight: 1.9990 chunk 11 optimal weight: 4.9990 chunk 43 optimal weight: 1.9990 chunk 69 optimal weight: 9.9990 overall best weight: 1.3784 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 350 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1570 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1781 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7252 moved from start: 0.2167 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.046 11436 Z= 0.241 Angle : 0.655 10.129 15490 Z= 0.337 Chirality : 0.046 0.156 1627 Planarity : 0.005 0.092 1999 Dihedral : 5.952 32.187 1508 Min Nonbonded Distance : 2.116 Molprobity Statistics. All-atom Clashscore : 13.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 15.62 % Favored : 84.38 % Rotamer Outliers : 1.88 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % 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.47 (0.20), residues: 1370 helix: -3.11 (0.58), residues: 52 sheet: -2.13 (0.27), residues: 372 loop : -3.58 (0.17), residues: 946 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 134 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 23 poor density : 111 time to evaluate : 1.383 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 23 outliers final: 13 residues processed: 129 average time/residue: 0.2084 time to fit residues: 41.3976 Evaluate side-chains 119 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 106 time to evaluate : 1.322 Switching outliers to nearest non-outliers outliers start: 13 outliers final: 0 residues processed: 13 average time/residue: 0.1114 time to fit residues: 4.5088 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 137 random chunks: chunk 127 optimal weight: 5.9990 chunk 14 optimal weight: 0.8980 chunk 75 optimal weight: 7.9990 chunk 96 optimal weight: 4.9990 chunk 74 optimal weight: 3.9990 chunk 111 optimal weight: 0.7980 chunk 132 optimal weight: 0.9990 chunk 82 optimal weight: 5.9990 chunk 80 optimal weight: 1.9990 chunk 60 optimal weight: 1.9990 chunk 81 optimal weight: 0.0020 overall best weight: 0.9392 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 350 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1570 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1781 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7215 moved from start: 0.2341 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.040 11436 Z= 0.196 Angle : 0.634 10.056 15490 Z= 0.324 Chirality : 0.046 0.209 1627 Planarity : 0.005 0.090 1999 Dihedral : 5.702 30.090 1508 Min Nonbonded Distance : 2.110 Molprobity Statistics. All-atom Clashscore : 13.08 Ramachandran Plot: Outliers : 0.00 % Allowed : 14.45 % Favored : 85.55 % Rotamer Outliers : 1.80 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % 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.34 (0.20), residues: 1370 helix: -2.63 (0.62), residues: 58 sheet: -2.05 (0.27), residues: 367 loop : -3.50 (0.17), residues: 945 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 136 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 22 poor density : 114 time to evaluate : 1.383 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 22 outliers final: 5 residues processed: 133 average time/residue: 0.2301 time to fit residues: 46.9442 Evaluate side-chains 115 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 110 time to evaluate : 1.322 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.1356 time to fit residues: 3.0996 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 137 random chunks: chunk 52 optimal weight: 0.5980 chunk 78 optimal weight: 5.9990 chunk 39 optimal weight: 1.9990 chunk 25 optimal weight: 6.9990 chunk 83 optimal weight: 6.9990 chunk 89 optimal weight: 0.9980 chunk 65 optimal weight: 2.9990 chunk 12 optimal weight: 0.9980 chunk 103 optimal weight: 4.9990 chunk 120 optimal weight: 1.9990 chunk 126 optimal weight: 4.9990 overall best weight: 1.3184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 350 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1570 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1781 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7238 moved from start: 0.2356 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.043 11436 Z= 0.235 Angle : 0.658 10.067 15490 Z= 0.337 Chirality : 0.047 0.186 1627 Planarity : 0.005 0.091 1999 Dihedral : 5.845 30.993 1508 Min Nonbonded Distance : 2.126 Molprobity Statistics. All-atom Clashscore : 13.95 Ramachandran Plot: Outliers : 0.00 % Allowed : 15.55 % Favored : 84.45 % Rotamer Outliers : 1.39 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % 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.30 (0.20), residues: 1370 helix: -2.66 (0.61), residues: 58 sheet: -1.95 (0.27), residues: 363 loop : -3.50 (0.17), residues: 949 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 125 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 17 poor density : 108 time to evaluate : 1.314 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 17 outliers final: 12 residues processed: 121 average time/residue: 0.2113 time to fit residues: 39.3071 Evaluate side-chains 116 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 104 time to evaluate : 1.431 Switching outliers to nearest non-outliers outliers start: 12 outliers final: 0 residues processed: 12 average time/residue: 0.1190 time to fit residues: 4.4798 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 137 random chunks: chunk 115 optimal weight: 4.9990 chunk 123 optimal weight: 0.8980 chunk 126 optimal weight: 5.9990 chunk 74 optimal weight: 5.9990 chunk 53 optimal weight: 4.9990 chunk 96 optimal weight: 2.9990 chunk 37 optimal weight: 0.8980 chunk 111 optimal weight: 4.9990 chunk 116 optimal weight: 2.9990 chunk 122 optimal weight: 1.9990 chunk 80 optimal weight: 0.7980 overall best weight: 1.5184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 350 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1570 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1781 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B1824 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7253 moved from start: 0.2408 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.045 11436 Z= 0.257 Angle : 0.678 9.947 15490 Z= 0.347 Chirality : 0.047 0.167 1627 Planarity : 0.005 0.093 1999 Dihedral : 5.976 32.742 1508 Min Nonbonded Distance : 2.132 Molprobity Statistics. All-atom Clashscore : 14.36 Ramachandran Plot: Outliers : 0.00 % Allowed : 15.77 % Favored : 84.23 % Rotamer Outliers : 0.49 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % 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.26 (0.20), residues: 1370 helix: -2.62 (0.61), residues: 58 sheet: -1.91 (0.27), residues: 385 loop : -3.50 (0.17), residues: 927 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 113 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 107 time to evaluate : 1.289 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 6 outliers final: 2 residues processed: 112 average time/residue: 0.2102 time to fit residues: 36.4542 Evaluate side-chains 101 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 99 time to evaluate : 1.364 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.2648 time to fit residues: 2.5364 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 137 random chunks: chunk 130 optimal weight: 4.9990 chunk 79 optimal weight: 2.9990 chunk 61 optimal weight: 4.9990 chunk 90 optimal weight: 3.9990 chunk 136 optimal weight: 4.9990 chunk 125 optimal weight: 0.2980 chunk 108 optimal weight: 1.9990 chunk 11 optimal weight: 0.4980 chunk 83 optimal weight: 6.9990 chunk 66 optimal weight: 1.9990 chunk 86 optimal weight: 3.9990 overall best weight: 1.5586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 350 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1570 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1781 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7260 moved from start: 0.2430 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.044 11436 Z= 0.262 Angle : 0.685 10.007 15490 Z= 0.351 Chirality : 0.047 0.154 1627 Planarity : 0.005 0.093 1999 Dihedral : 6.047 33.501 1508 Min Nonbonded Distance : 2.136 Molprobity Statistics. All-atom Clashscore : 14.81 Ramachandran Plot: Outliers : 0.00 % Allowed : 15.84 % Favored : 84.16 % Rotamer Outliers : 0.41 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % 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.25 (0.20), residues: 1370 helix: -2.68 (0.60), residues: 58 sheet: -1.89 (0.27), residues: 385 loop : -3.50 (0.17), residues: 927 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2740 Ramachandran restraints generated. 1370 Oldfield, 0 Emsley, 1370 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 107 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 102 time to evaluate : 1.329 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 5 outliers final: 3 residues processed: 106 average time/residue: 0.2210 time to fit residues: 36.1039 Evaluate side-chains 102 residues out of total 1225 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 99 time to evaluate : 1.333 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1162 time to fit residues: 2.5593 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 137 random chunks: chunk 115 optimal weight: 4.9990 chunk 33 optimal weight: 0.9990 chunk 100 optimal weight: 0.9990 chunk 16 optimal weight: 0.9990 chunk 30 optimal weight: 3.9990 chunk 108 optimal weight: 0.5980 chunk 45 optimal weight: 0.9980 chunk 111 optimal weight: 4.9990 chunk 13 optimal weight: 0.5980 chunk 20 optimal weight: 3.9990 chunk 95 optimal weight: 7.9990 overall best weight: 0.8384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 350 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1570 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B1781 ASN ** 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.4158 r_free = 0.4158 target = 0.140836 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 37)----------------| | r_work = 0.3889 r_free = 0.3889 target = 0.121380 restraints weight = 25828.139| |-----------------------------------------------------------------------------| r_work (start): 0.3880 rms_B_bonded: 3.61 r_work: 0.3773 rms_B_bonded: 4.26 restraints_weight: 0.5000 r_work (final): 0.3773 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7348 moved from start: 0.2659 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.041 11436 Z= 0.191 Angle : 0.646 10.159 15490 Z= 0.329 Chirality : 0.046 0.156 1627 Planarity : 0.005 0.090 1999 Dihedral : 5.771 30.301 1508 Min Nonbonded Distance : 2.145 Molprobity Statistics. All-atom Clashscore : 13.58 Ramachandran Plot: Outliers : 0.00 % Allowed : 14.09 % Favored : 85.91 % Rotamer Outliers : 0.24 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.52 % Cis-general : 0.08 % 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.12 (0.20), residues: 1370 helix: -2.64 (0.61), residues: 58 sheet: -1.77 (0.27), residues: 374 loop : -3.41 (0.17), residues: 938 =============================================================================== Job complete usr+sys time: 2245.07 seconds wall clock time: 41 minutes 49.66 seconds (2509.66 seconds total)