Starting phenix.real_space_refine (version: dev) on Fri Feb 17 05:14:11 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7u16_26279/02_2023/7u16_26279.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7u16_26279/02_2023/7u16_26279.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.7 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7u16_26279/02_2023/7u16_26279.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7u16_26279/02_2023/7u16_26279.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7u16_26279/02_2023/7u16_26279.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7u16_26279/02_2023/7u16_26279.pdb" } resolution = 2.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 None Time to flip residues: 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/chem_data/mon_lib" Total number of atoms: 3456 Number of models: 1 Model: "" Number of chains: 6 Chain: "A" Number of atoms: 1096 Number of conformers: 2 Conformer: "A" Number of residues, atoms: 135, 1087 Classifications: {'peptide': 135} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 2, 'TRANS': 132} Conformer: "B" Number of residues, atoms: 135, 1087 Classifications: {'peptide': 135} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 2, 'TRANS': 132} bond proxies already assigned to first conformer: 1094 Chain: "B" Number of atoms: 1096 Number of conformers: 2 Conformer: "A" Number of residues, atoms: 135, 1087 Classifications: {'peptide': 135} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 2, 'TRANS': 132} Conformer: "B" Number of residues, atoms: 135, 1087 Classifications: {'peptide': 135} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 2, 'TRANS': 132} bond proxies already assigned to first conformer: 1094 Chain: "C" Number of atoms: 1096 Number of conformers: 2 Conformer: "A" Number of residues, atoms: 135, 1087 Classifications: {'peptide': 135} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 2, 'TRANS': 132} Conformer: "B" Number of residues, atoms: 135, 1087 Classifications: {'peptide': 135} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 2, 'TRANS': 132} bond proxies already assigned to first conformer: 1094 Chain: "A" Number of atoms: 56 Number of conformers: 1 Conformer: "" Number of residues, atoms: 4, 56 Unusual residues: {'NAG': 4} Classifications: {'undetermined': 4} Link IDs: {None: 3} Unresolved non-hydrogen bonds: 4 Unresolved non-hydrogen angles: 8 Unresolved non-hydrogen chiralities: 4 Chain: "B" Number of atoms: 56 Number of conformers: 1 Conformer: "" Number of residues, atoms: 4, 56 Unusual residues: {'NAG': 4} Classifications: {'undetermined': 4} Link IDs: {None: 3} Unresolved non-hydrogen bonds: 4 Unresolved non-hydrogen angles: 8 Unresolved non-hydrogen chiralities: 4 Chain: "C" Number of atoms: 56 Number of conformers: 1 Conformer: "" Number of residues, atoms: 4, 56 Unusual residues: {'NAG': 4} Classifications: {'undetermined': 4} Link IDs: {None: 3} Unresolved non-hydrogen bonds: 4 Unresolved non-hydrogen angles: 8 Unresolved non-hydrogen chiralities: 4 Time building chain proxies: 3.76, per 1000 atoms: 1.09 Number of scatterers: 3456 At special positions: 0 Unit cell: (105.252, 99.882, 32.22, 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 21 16.00 O 702 8.00 N 531 7.00 C 2202 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=3, symmetry=0 Simple disulfide: pdb=" SG CYS A 214 " - pdb=" SG CYS A 253 " distance=2.03 Simple disulfide: pdb=" SG CYS B 214 " - pdb=" SG CYS B 253 " distance=2.03 Simple disulfide: pdb=" SG CYS C 214 " - pdb=" SG CYS C 253 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Links applied NAG-ASN " NAG A 301 " - " ASN A 183 " " NAG A 302 " - " ASN A 151 " " NAG A 303 " - " ASN A 145 " " NAG A 304 " - " ASN A 164 " " NAG B 301 " - " ASN B 183 " " NAG B 302 " - " ASN B 151 " " NAG B 303 " - " ASN B 145 " " NAG B 304 " - " ASN B 164 " " NAG C 301 " - " ASN C 183 " " NAG C 302 " - " ASN C 151 " " NAG C 303 " - " ASN C 145 " " NAG C 304 " - " ASN C 164 " Time building additional restraints: 1.19 Conformation dependent library (CDL) restraints added in 833.9 milliseconds 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 786 Finding SS restraints... Secondary structure from input PDB file: 0 helices and 11 sheets defined 0.0% alpha, 67.4% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.20 Creating SS restraints... Processing sheet with id=AA1, first strand: chain 'A' and resid 122 through 126 Processing sheet with id=AA2, first strand: chain 'A' and resid 129 through 137 removed outlier: 6.371A pdb=" N LYS A 129 " --> pdb=" O SER B 130 " (cutoff:3.500A) removed outlier: 8.229A pdb=" N TYR B 132 " --> pdb=" O LYS A 129 " (cutoff:3.500A) removed outlier: 9.088A pdb=" N ALA A 131 " --> pdb=" O TYR B 132 " (cutoff:3.500A) removed outlier: 6.498A pdb=" N SER B 134 " --> pdb=" O ALA A 131 " (cutoff:3.500A) removed outlier: 9.018A pdb=" N VAL A 133 " --> pdb=" O SER B 134 " (cutoff:3.500A) removed outlier: 9.432A pdb=" N ASP B 136 " --> pdb=" O VAL A 133 " (cutoff:3.500A) removed outlier: 8.667A pdb=" N TYR A 135 " --> pdb=" O ASP B 136 " (cutoff:3.500A) removed outlier: 6.392A pdb=" N LYS B 129 " --> pdb=" O SER C 130 " (cutoff:3.500A) removed outlier: 8.254A pdb=" N TYR C 132 " --> pdb=" O LYS B 129 " (cutoff:3.500A) removed outlier: 9.092A pdb=" N ALA B 131 " --> pdb=" O TYR C 132 " (cutoff:3.500A) removed outlier: 6.505A pdb=" N SER C 134 " --> pdb=" O ALA B 131 " (cutoff:3.500A) removed outlier: 8.999A pdb=" N VAL B 133 " --> pdb=" O SER C 134 " (cutoff:3.500A) removed outlier: 9.419A pdb=" N ASP C 136 " --> pdb=" O VAL B 133 " (cutoff:3.500A) removed outlier: 8.650A pdb=" N TYR B 135 " --> pdb=" O ASP C 136 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA2 Processing sheet with id=AA3, first strand: chain 'A' and resid 142 through 165 removed outlier: 6.342A pdb=" N ILE A 142 " --> pdb=" O TYR B 143 " (cutoff:3.500A) removed outlier: 7.798A pdb=" N ASN B 145 " --> pdb=" O ILE A 142 " (cutoff:3.500A) removed outlier: 8.827A pdb=" N LEU A 144 " --> pdb=" O ASN B 145 " (cutoff:3.500A) removed outlier: 8.022A pdb=" N THR B 147 " --> pdb=" O LEU A 144 " (cutoff:3.500A) removed outlier: 9.195A pdb=" N ILE A 146 " --> pdb=" O THR B 147 " (cutoff:3.500A) removed outlier: 10.073A pdb=" N THR B 149 " --> pdb=" O ILE A 146 " (cutoff:3.500A) removed outlier: 9.033A pdb=" N ASN A 148 " --> pdb=" O THR B 149 " (cutoff:3.500A) removed outlier: 10.294A pdb=" N ASN B 151 " --> pdb=" O ASN A 148 " (cutoff:3.500A) removed outlier: 8.986A pdb=" N LEU A 150 " --> pdb=" O ASN B 151 " (cutoff:3.500A) removed outlier: 10.178A pdb=" N THR B 153 " --> pdb=" O LEU A 150 " (cutoff:3.500A) removed outlier: 9.199A pdb=" N ILE A 152 " --> pdb=" O THR B 153 " (cutoff:3.500A) removed outlier: 9.734A pdb=" N ASN B 155 " --> pdb=" O ILE A 152 " (cutoff:3.500A) removed outlier: 8.872A pdb=" N ASN A 154 " --> pdb=" O ASN B 155 " (cutoff:3.500A) removed outlier: 10.389A pdb=" N TYR B 157 " --> pdb=" O ASN A 154 " (cutoff:3.500A) removed outlier: 9.051A pdb=" N ASN A 156 " --> pdb=" O TYR B 157 " (cutoff:3.500A) removed outlier: 10.322A pdb=" N SER B 159 " --> pdb=" O ASN A 156 " (cutoff:3.500A) removed outlier: 9.519A pdb=" N TYR A 158 " --> pdb=" O SER B 159 " (cutoff:3.500A) removed outlier: 10.246A pdb=" N GLU B 161 " --> pdb=" O TYR A 158 " (cutoff:3.500A) removed outlier: 6.018A pdb=" N VAL A 160 " --> pdb=" O GLU B 161 " (cutoff:3.500A) removed outlier: 8.788A pdb=" N GLU B 163 " --> pdb=" O VAL A 160 " (cutoff:3.500A) removed outlier: 6.595A pdb=" N VAL A 162 " --> pdb=" O GLU B 163 " (cutoff:3.500A) removed outlier: 7.846A pdb=" N ILE B 165 " --> pdb=" O VAL A 162 " (cutoff:3.500A) removed outlier: 5.974A pdb=" N ASN A 164 " --> pdb=" O ILE B 165 " (cutoff:3.500A) removed outlier: 6.379A pdb=" N ILE B 142 " --> pdb=" O TYR C 143 " (cutoff:3.500A) removed outlier: 7.842A pdb=" N ASN C 145 " --> pdb=" O ILE B 142 " (cutoff:3.500A) removed outlier: 8.843A pdb=" N LEU B 144 " --> pdb=" O ASN C 145 " (cutoff:3.500A) removed outlier: 8.030A pdb=" N THR C 147 " --> pdb=" O LEU B 144 " (cutoff:3.500A) removed outlier: 9.182A pdb=" N ILE B 146 " --> pdb=" O THR C 147 " (cutoff:3.500A) removed outlier: 10.050A pdb=" N THR C 149 " --> pdb=" O ILE B 146 " (cutoff:3.500A) removed outlier: 9.001A pdb=" N ASN B 148 " --> pdb=" O THR C 149 " (cutoff:3.500A) removed outlier: 10.256A pdb=" N ASN C 151 " --> pdb=" O ASN B 148 " (cutoff:3.500A) removed outlier: 8.938A pdb=" N LEU B 150 " --> pdb=" O ASN C 151 " (cutoff:3.500A) removed outlier: 10.126A pdb=" N THR C 153 " --> pdb=" O LEU B 150 " (cutoff:3.500A) removed outlier: 9.141A pdb=" N ILE B 152 " --> pdb=" O THR C 153 " (cutoff:3.500A) removed outlier: 9.671A pdb=" N ASN C 155 " --> pdb=" O ILE B 152 " (cutoff:3.500A) removed outlier: 8.802A pdb=" N ASN B 154 " --> pdb=" O ASN C 155 " (cutoff:3.500A) removed outlier: 10.322A pdb=" N TYR C 157 " --> pdb=" O ASN B 154 " (cutoff:3.500A) removed outlier: 8.976A pdb=" N ASN B 156 " --> pdb=" O TYR C 157 " (cutoff:3.500A) removed outlier: 10.247A pdb=" N SER C 159 " --> pdb=" O ASN B 156 " (cutoff:3.500A) removed outlier: 9.437A pdb=" N TYR B 158 " --> pdb=" O SER C 159 " (cutoff:3.500A) removed outlier: 10.166A pdb=" N GLU C 161 " --> pdb=" O TYR B 158 " (cutoff:3.500A) removed outlier: 5.955A pdb=" N VAL B 160 " --> pdb=" O GLU C 161 " (cutoff:3.500A) removed outlier: 8.720A pdb=" N GLU C 163 " --> pdb=" O VAL B 160 " (cutoff:3.500A) removed outlier: 6.541A pdb=" N VAL B 162 " --> pdb=" O GLU C 163 " (cutoff:3.500A) removed outlier: 7.797A pdb=" N ILE C 165 " --> pdb=" O VAL B 162 " (cutoff:3.500A) removed outlier: 5.912A pdb=" N ASN B 164 " --> pdb=" O ILE C 165 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA3 Processing sheet with id=AA4, first strand: chain 'A' and resid 168 through 177 Processing sheet with id=AA5, first strand: chain 'A' and resid 180 through 181 removed outlier: 6.572A pdb=" N ARG A 180 " --> pdb=" O LEU B 181 " (cutoff:3.500A) removed outlier: 6.560A pdb=" N ARG B 180 " --> pdb=" O LEU C 181 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA5 Processing sheet with id=AA6, first strand: chain 'A' and resid 184 through 187 Processing sheet with id=AA7, first strand: chain 'A' and resid 191 through 198 removed outlier: 8.944A pdb=" N MET A 192 " --> pdb=" O LYS B 193 " (cutoff:3.500A) removed outlier: 9.058A pdb=" N ILE B 195 " --> pdb=" O MET A 192 " (cutoff:3.500A) removed outlier: 6.798A pdb=" N GLN A 194 " --> pdb=" O ILE B 195 " (cutoff:3.500A) removed outlier: 7.984A pdb=" N TYR B 197 " --> pdb=" O GLN A 194 " (cutoff:3.500A) removed outlier: 6.438A pdb=" N ASP A 196 " --> pdb=" O TYR B 197 " (cutoff:3.500A) removed outlier: 8.933A pdb=" N MET B 192 " --> pdb=" O LYS C 193 " (cutoff:3.500A) removed outlier: 9.066A pdb=" N ILE C 195 " --> pdb=" O MET B 192 " (cutoff:3.500A) removed outlier: 6.830A pdb=" N GLN B 194 " --> pdb=" O ILE C 195 " (cutoff:3.500A) removed outlier: 8.015A pdb=" N TYR C 197 " --> pdb=" O GLN B 194 " (cutoff:3.500A) removed outlier: 6.470A pdb=" N ASP B 196 " --> pdb=" O TYR C 197 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA7 Processing sheet with id=AA8, first strand: chain 'A' and resid 201 through 204 Processing sheet with id=AA9, first strand: chain 'A' and resid 209 through 228 removed outlier: 6.946A pdb=" N MET A 210 " --> pdb=" O TYR B 209 " (cutoff:3.500A) removed outlier: 6.954A pdb=" N TYR B 211 " --> pdb=" O MET A 210 " (cutoff:3.500A) removed outlier: 6.974A pdb=" N ASP A 212 " --> pdb=" O TYR B 211 " (cutoff:3.500A) removed outlier: 6.941A pdb=" N PHE B 213 " --> pdb=" O ASP A 212 " (cutoff:3.500A) removed outlier: 6.744A pdb=" N VAL B 221 " --> pdb=" O LYS A 220 " (cutoff:3.500A) removed outlier: 7.012A pdb=" N HIS A 222 " --> pdb=" O VAL B 221 " (cutoff:3.500A) removed outlier: 6.873A pdb=" N ASN B 223 " --> pdb=" O HIS A 222 " (cutoff:3.500A) removed outlier: 6.102A pdb=" N TYR B 209 " --> pdb=" O MET C 210 " (cutoff:3.500A) removed outlier: 7.336A pdb=" N ASP C 212 " --> pdb=" O TYR B 209 " (cutoff:3.500A) removed outlier: 6.357A pdb=" N TYR B 211 " --> pdb=" O ASP C 212 " (cutoff:3.500A) removed outlier: 7.255A pdb=" N CYS C 214 " --> pdb=" O TYR B 211 " (cutoff:3.500A) removed outlier: 8.371A pdb=" N PHE B 213 " --> pdb=" O CYS C 214 " (cutoff:3.500A) removed outlier: 10.442A pdb=" N LEU C 216 " --> pdb=" O PHE B 213 " (cutoff:3.500A) removed outlier: 8.876A pdb=" N THR B 215 " --> pdb=" O LEU C 216 " (cutoff:3.500A) removed outlier: 9.799A pdb=" N SER C 218 " --> pdb=" O THR B 215 " (cutoff:3.500A) removed outlier: 8.795A pdb=" N ILE B 217 " --> pdb=" O SER C 218 " (cutoff:3.500A) removed outlier: 9.946A pdb=" N LYS C 220 " --> pdb=" O ILE B 217 " (cutoff:3.500A) removed outlier: 6.071A pdb=" N ILE B 219 " --> pdb=" O LYS C 220 " (cutoff:3.500A) removed outlier: 8.299A pdb=" N HIS C 222 " --> pdb=" O ILE B 219 " (cutoff:3.500A) removed outlier: 5.887A pdb=" N VAL B 221 " --> pdb=" O HIS C 222 " (cutoff:3.500A) removed outlier: 7.500A pdb=" N ILE C 224 " --> pdb=" O VAL B 221 " (cutoff:3.500A) removed outlier: 7.857A pdb=" N ASN B 223 " --> pdb=" O ILE C 224 " (cutoff:3.500A) removed outlier: 9.821A pdb=" N LEU C 226 " --> pdb=" O ASN B 223 " (cutoff:3.500A) removed outlier: 8.747A pdb=" N VAL B 225 " --> pdb=" O LEU C 226 " (cutoff:3.500A) removed outlier: 10.088A pdb=" N MET C 228 " --> pdb=" O VAL B 225 " (cutoff:3.500A) removed outlier: 9.225A pdb=" N MET B 227 " --> pdb=" O MET C 228 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'A' and resid 233 through 245 Processing sheet with id=AB2, first strand: chain 'A' and resid 251 through 253 78 hydrogen bonds defined for protein. 234 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: 0.82 Time building geometry restraints manager: 1.30 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.23 - 1.34: 1053 1.34 - 1.46: 791 1.46 - 1.58: 1636 1.58 - 1.69: 0 1.69 - 1.81: 36 Bond restraints: 3516 Sorted by residual: bond pdb=" C1 NAG A 304 " pdb=" O5 NAG A 304 " ideal model delta sigma weight residual 1.426 1.478 -0.052 3.20e-02 9.77e+02 2.67e+00 bond pdb=" C1 NAG C 304 " pdb=" O5 NAG C 304 " ideal model delta sigma weight residual 1.426 1.478 -0.052 3.20e-02 9.77e+02 2.64e+00 bond pdb=" C1 NAG B 304 " pdb=" O5 NAG B 304 " ideal model delta sigma weight residual 1.426 1.478 -0.052 3.20e-02 9.77e+02 2.64e+00 bond pdb=" CA TYR B 236 " pdb=" CB TYR B 236 " ideal model delta sigma weight residual 1.534 1.551 -0.017 1.76e-02 3.23e+03 9.66e-01 bond pdb=" CB THR B 234 " pdb=" CG2 THR B 234 " ideal model delta sigma weight residual 1.521 1.489 0.032 3.30e-02 9.18e+02 9.44e-01 ... (remaining 3511 not shown) Histogram of bond angle deviations from ideal: 100.54 - 106.79: 75 106.79 - 113.04: 1984 113.04 - 119.29: 878 119.29 - 125.54: 1839 125.54 - 131.80: 12 Bond angle restraints: 4788 Sorted by residual: angle pdb=" C1 NAG B 304 " pdb=" O5 NAG B 304 " pdb=" C5 NAG B 304 " ideal model delta sigma weight residual 113.34 117.83 -4.49 1.74e+00 3.30e-01 6.64e+00 angle pdb=" C1 NAG A 304 " pdb=" O5 NAG A 304 " pdb=" C5 NAG A 304 " ideal model delta sigma weight residual 113.34 117.78 -4.44 1.74e+00 3.30e-01 6.50e+00 angle pdb=" C1 NAG C 304 " pdb=" O5 NAG C 304 " pdb=" C5 NAG C 304 " ideal model delta sigma weight residual 113.34 117.77 -4.43 1.74e+00 3.30e-01 6.49e+00 angle pdb=" CA ARG C 247 " pdb=" CB ARG C 247 " pdb=" CG ARG C 247 " ideal model delta sigma weight residual 114.10 118.77 -4.67 2.00e+00 2.50e-01 5.45e+00 angle pdb=" CA ARG B 247 " pdb=" CB ARG B 247 " pdb=" CG ARG B 247 " ideal model delta sigma weight residual 114.10 118.75 -4.65 2.00e+00 2.50e-01 5.40e+00 ... (remaining 4783 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 21.88: 1957 21.88 - 43.76: 113 43.76 - 65.64: 12 65.64 - 87.51: 6 87.51 - 109.39: 3 Dihedral angle restraints: 2091 sinusoidal: 858 harmonic: 1233 Sorted by residual: dihedral pdb=" C2 NAG C 304 " pdb=" C1 NAG C 304 " pdb=" O5 NAG C 304 " pdb=" C5 NAG C 304 " ideal model delta sinusoidal sigma weight residual -50.00 59.39 -109.39 1 2.00e+01 2.50e-03 3.20e+01 dihedral pdb=" C2 NAG B 304 " pdb=" C1 NAG B 304 " pdb=" O5 NAG B 304 " pdb=" C5 NAG B 304 " ideal model delta sinusoidal sigma weight residual -50.00 59.34 -109.34 1 2.00e+01 2.50e-03 3.19e+01 dihedral pdb=" C2 NAG A 304 " pdb=" C1 NAG A 304 " pdb=" O5 NAG A 304 " pdb=" C5 NAG A 304 " ideal model delta sinusoidal sigma weight residual -50.00 59.31 -109.31 1 2.00e+01 2.50e-03 3.19e+01 ... (remaining 2088 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.038: 381 0.038 - 0.075: 114 0.075 - 0.113: 77 0.113 - 0.151: 28 0.151 - 0.188: 6 Chirality restraints: 606 Sorted by residual: chirality pdb=" C1 NAG A 301 " pdb=" ND2 ASN A 183 " pdb=" C2 NAG A 301 " pdb=" O5 NAG A 301 " both_signs ideal model delta sigma weight residual False -2.40 -2.21 -0.19 2.00e-01 2.50e+01 8.87e-01 chirality pdb=" C1 NAG C 301 " pdb=" ND2 ASN C 183 " pdb=" C2 NAG C 301 " pdb=" O5 NAG C 301 " both_signs ideal model delta sigma weight residual False -2.40 -2.22 -0.18 2.00e-01 2.50e+01 8.40e-01 chirality pdb=" CA TYR C 236 " pdb=" N TYR C 236 " pdb=" C TYR C 236 " pdb=" CB TYR C 236 " both_signs ideal model delta sigma weight residual False 2.51 2.33 0.18 2.00e-01 2.50e+01 8.27e-01 ... (remaining 603 not shown) Planarity restraints: 591 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB TYR B 236 " -0.005 2.00e-02 2.50e+03 2.80e-02 1.57e+01 pdb=" CG TYR B 236 " 0.056 2.00e-02 2.50e+03 pdb=" CD1 TYR B 236 " -0.046 2.00e-02 2.50e+03 pdb=" CD2 TYR B 236 " -0.023 2.00e-02 2.50e+03 pdb=" CE1 TYR B 236 " 0.020 2.00e-02 2.50e+03 pdb=" CE2 TYR B 236 " -0.002 2.00e-02 2.50e+03 pdb=" CZ TYR B 236 " -0.003 2.00e-02 2.50e+03 pdb=" OH TYR B 236 " 0.003 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TYR C 236 " -0.005 2.00e-02 2.50e+03 2.79e-02 1.56e+01 pdb=" CG TYR C 236 " 0.056 2.00e-02 2.50e+03 pdb=" CD1 TYR C 236 " -0.046 2.00e-02 2.50e+03 pdb=" CD2 TYR C 236 " -0.024 2.00e-02 2.50e+03 pdb=" CE1 TYR C 236 " 0.020 2.00e-02 2.50e+03 pdb=" CE2 TYR C 236 " -0.002 2.00e-02 2.50e+03 pdb=" CZ TYR C 236 " -0.003 2.00e-02 2.50e+03 pdb=" OH TYR C 236 " 0.003 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TYR A 236 " -0.005 2.00e-02 2.50e+03 2.79e-02 1.56e+01 pdb=" CG TYR A 236 " 0.056 2.00e-02 2.50e+03 pdb=" CD1 TYR A 236 " -0.046 2.00e-02 2.50e+03 pdb=" CD2 TYR A 236 " -0.023 2.00e-02 2.50e+03 pdb=" CE1 TYR A 236 " 0.020 2.00e-02 2.50e+03 pdb=" CE2 TYR A 236 " -0.002 2.00e-02 2.50e+03 pdb=" CZ TYR A 236 " -0.003 2.00e-02 2.50e+03 pdb=" OH TYR A 236 " 0.003 2.00e-02 2.50e+03 ... (remaining 588 not shown) Histogram of nonbonded interaction distances: 2.31 - 2.83: 883 2.83 - 3.35: 2777 3.35 - 3.86: 5897 3.86 - 4.38: 6080 4.38 - 4.90: 11940 Nonbonded interactions: 27577 Sorted by model distance: nonbonded pdb=" NZ LYS A 129 " pdb=" OD2 ASP A 136 " model vdw 2.310 2.520 nonbonded pdb=" NZ LYS B 129 " pdb=" OD2 ASP B 136 " model vdw 2.311 2.520 nonbonded pdb=" NZ LYS C 129 " pdb=" OD2 ASP C 136 " model vdw 2.311 2.520 nonbonded pdb=" OD1 ASN C 145 " pdb=" N2 NAG C 303 " model vdw 2.361 2.520 nonbonded pdb=" OD1 ASN B 145 " pdb=" N2 NAG B 303 " model vdw 2.362 2.520 ... (remaining 27572 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = (chain 'A' and (resid 120 through 210 or resid 212 through 254 or resid 301 thro \ ugh 304)) selection = (chain 'B' and (resid 120 through 210 or resid 212 through 254 or resid 301 thro \ ugh 304)) selection = (chain 'C' and (resid 120 through 210 or resid 212 through 254 or resid 301 thro \ ugh 304)) } 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 21 5.16 5 C 2202 2.51 5 N 531 2.21 5 O 702 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=0.50 max=1.00 mean=0.99 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 2.040 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: 9.360 Check model and map are aligned: 0.050 Process input model: 15.820 Find NCS groups from input model: 0.220 Set up NCS constraints: 0.020 Set refine NCS operators: 0.000 Set scattering table: 0.040 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:11.560 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 39.110 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6834 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.052 3516 Z= 0.208 Angle : 0.550 4.668 4788 Z= 0.311 Chirality : 0.055 0.188 606 Planarity : 0.004 0.028 579 Dihedral : 14.720 109.392 1296 Min Nonbonded Distance : 2.310 Molprobity Statistics. All-atom Clashscore : 8.58 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.77 % Favored : 93.23 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.67 (0.39), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.03 (0.30), residues: 402 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 60 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 60 time to evaluate : 0.442 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 60 average time/residue: 1.1538 time to fit residues: 71.5000 Evaluate side-chains 47 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 47 time to evaluate : 0.397 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5609 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 39 random chunks: chunk 32 optimal weight: 1.9990 chunk 29 optimal weight: 3.9990 chunk 16 optimal weight: 0.8980 chunk 10 optimal weight: 7.9990 chunk 19 optimal weight: 0.9990 chunk 15 optimal weight: 6.9990 chunk 30 optimal weight: 3.9990 chunk 11 optimal weight: 4.9990 chunk 18 optimal weight: 3.9990 chunk 22 optimal weight: 9.9990 chunk 35 optimal weight: 0.9980 overall best weight: 1.7786 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7016 moved from start: 0.2250 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.052 3516 Z= 0.295 Angle : 0.549 5.396 4788 Z= 0.309 Chirality : 0.057 0.417 606 Planarity : 0.003 0.022 579 Dihedral : 12.881 109.102 531 Min Nonbonded Distance : 2.297 Molprobity Statistics. All-atom Clashscore : 7.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.27 % Favored : 93.73 % Rotamer Outliers : 2.89 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.57 (0.39), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.95 (0.30), residues: 402 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 55 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 44 time to evaluate : 0.368 Fit side-chains outliers start: 11 outliers final: 4 residues processed: 47 average time/residue: 1.2098 time to fit residues: 58.6890 Evaluate side-chains 45 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 41 time to evaluate : 0.388 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 1 residues processed: 3 average time/residue: 0.0337 time to fit residues: 0.7242 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 39 random chunks: chunk 19 optimal weight: 2.9990 chunk 10 optimal weight: 5.9990 chunk 29 optimal weight: 6.9990 chunk 24 optimal weight: 0.7980 chunk 9 optimal weight: 4.9990 chunk 35 optimal weight: 7.9990 chunk 38 optimal weight: 2.9990 chunk 31 optimal weight: 0.9980 chunk 12 optimal weight: 3.9990 chunk 28 optimal weight: 6.9990 chunk 26 optimal weight: 2.9990 overall best weight: 2.1586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7090 moved from start: 0.3116 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.055 3516 Z= 0.321 Angle : 0.558 5.019 4788 Z= 0.314 Chirality : 0.061 0.610 606 Planarity : 0.003 0.022 579 Dihedral : 12.861 108.709 531 Min Nonbonded Distance : 2.234 Molprobity Statistics. All-atom Clashscore : 7.12 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.02 % Favored : 91.98 % Rotamer Outliers : 3.67 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.58 (0.38), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.96 (0.29), residues: 402 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 61 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 47 time to evaluate : 0.374 Fit side-chains revert: symmetry clash outliers start: 14 outliers final: 6 residues processed: 53 average time/residue: 1.2446 time to fit residues: 68.0642 Evaluate side-chains 48 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 42 time to evaluate : 0.388 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 4 residues processed: 2 average time/residue: 0.0346 time to fit residues: 0.6637 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 39 random chunks: chunk 18 optimal weight: 3.9990 chunk 3 optimal weight: 2.9990 chunk 16 optimal weight: 3.9990 chunk 23 optimal weight: 5.9990 chunk 35 optimal weight: 8.9990 chunk 37 optimal weight: 2.9990 chunk 33 optimal weight: 7.9990 chunk 10 optimal weight: 2.9990 chunk 31 optimal weight: 3.9990 chunk 21 optimal weight: 4.9990 chunk 0 optimal weight: 5.9990 overall best weight: 3.3990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7158 moved from start: 0.3764 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.007 0.060 3516 Z= 0.465 Angle : 0.630 5.990 4788 Z= 0.353 Chirality : 0.069 0.711 606 Planarity : 0.003 0.026 579 Dihedral : 12.864 107.258 531 Min Nonbonded Distance : 2.183 Molprobity Statistics. All-atom Clashscore : 8.00 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.02 % Favored : 90.98 % Rotamer Outliers : 3.41 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.63 (0.38), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.00 (0.29), residues: 402 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 58 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 45 time to evaluate : 0.414 Fit side-chains revert: symmetry clash outliers start: 13 outliers final: 9 residues processed: 52 average time/residue: 1.6155 time to fit residues: 86.0963 Evaluate side-chains 52 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 43 time to evaluate : 0.414 Switching outliers to nearest non-outliers outliers start: 9 outliers final: 6 residues processed: 3 average time/residue: 0.5703 time to fit residues: 2.3468 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 39 random chunks: chunk 28 optimal weight: 2.9990 chunk 15 optimal weight: 3.9990 chunk 32 optimal weight: 3.9990 chunk 26 optimal weight: 3.9990 chunk 0 optimal weight: 4.9990 chunk 19 optimal weight: 5.9990 chunk 33 optimal weight: 5.9990 chunk 9 optimal weight: 7.9990 chunk 12 optimal weight: 0.7980 chunk 7 optimal weight: 1.9990 chunk 22 optimal weight: 5.9990 overall best weight: 2.7588 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 229 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7152 moved from start: 0.3890 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.006 0.059 3516 Z= 0.383 Angle : 0.576 5.451 4788 Z= 0.322 Chirality : 0.062 0.543 606 Planarity : 0.003 0.027 579 Dihedral : 12.794 107.202 531 Min Nonbonded Distance : 2.209 Molprobity Statistics. All-atom Clashscore : 8.58 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.27 % Favored : 90.73 % Rotamer Outliers : 2.62 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.66 (0.38), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.02 (0.29), residues: 402 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 53 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 43 time to evaluate : 0.351 Fit side-chains revert: symmetry clash outliers start: 10 outliers final: 8 residues processed: 49 average time/residue: 1.5905 time to fit residues: 79.9429 Evaluate side-chains 50 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 42 time to evaluate : 0.387 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 6 residues processed: 2 average time/residue: 0.0337 time to fit residues: 0.6534 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 39 random chunks: chunk 9 optimal weight: 4.9990 chunk 37 optimal weight: 0.9990 chunk 31 optimal weight: 5.9990 chunk 17 optimal weight: 0.9980 chunk 3 optimal weight: 3.9990 chunk 12 optimal weight: 3.9990 chunk 19 optimal weight: 5.9990 chunk 36 optimal weight: 1.9990 chunk 4 optimal weight: 0.1980 chunk 21 optimal weight: 0.9980 chunk 27 optimal weight: 7.9990 overall best weight: 1.0384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7101 moved from start: 0.3838 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.060 3516 Z= 0.201 Angle : 0.493 5.076 4788 Z= 0.276 Chirality : 0.055 0.342 606 Planarity : 0.003 0.036 579 Dihedral : 12.713 108.603 531 Min Nonbonded Distance : 2.286 Molprobity Statistics. All-atom Clashscore : 7.70 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.26 % Favored : 94.74 % Rotamer Outliers : 2.62 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.58 (0.38), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.96 (0.29), residues: 402 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 50 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 40 time to evaluate : 0.343 Fit side-chains revert: symmetry clash outliers start: 10 outliers final: 7 residues processed: 46 average time/residue: 1.4689 time to fit residues: 69.4150 Evaluate side-chains 45 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 38 time to evaluate : 0.321 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 5 residues processed: 2 average time/residue: 0.0325 time to fit residues: 0.6598 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 39 random chunks: chunk 21 optimal weight: 9.9990 chunk 31 optimal weight: 3.9990 chunk 37 optimal weight: 1.9990 chunk 23 optimal weight: 7.9990 chunk 22 optimal weight: 1.9990 chunk 17 optimal weight: 4.9990 chunk 15 optimal weight: 5.9990 chunk 11 optimal weight: 5.9990 chunk 7 optimal weight: 1.9990 chunk 25 optimal weight: 3.9990 chunk 18 optimal weight: 3.9990 overall best weight: 2.7990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7152 moved from start: 0.4105 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.006 0.059 3516 Z= 0.390 Angle : 0.583 5.351 4788 Z= 0.327 Chirality : 0.063 0.582 606 Planarity : 0.003 0.031 579 Dihedral : 12.782 107.374 531 Min Nonbonded Distance : 2.209 Molprobity Statistics. All-atom Clashscore : 8.72 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.02 % Favored : 90.98 % Rotamer Outliers : 2.62 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.65 (0.38), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.01 (0.29), residues: 402 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 51 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 41 time to evaluate : 0.407 Fit side-chains outliers start: 10 outliers final: 8 residues processed: 48 average time/residue: 1.5645 time to fit residues: 77.0564 Evaluate side-chains 48 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 40 time to evaluate : 0.388 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 5 residues processed: 3 average time/residue: 0.4529 time to fit residues: 1.9738 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 39 random chunks: chunk 3 optimal weight: 3.9990 chunk 29 optimal weight: 3.9990 chunk 34 optimal weight: 0.9980 chunk 36 optimal weight: 0.6980 chunk 32 optimal weight: 2.9990 chunk 35 optimal weight: 4.9990 chunk 21 optimal weight: 3.9990 chunk 15 optimal weight: 1.9990 chunk 27 optimal weight: 3.9990 chunk 10 optimal weight: 2.9990 chunk 31 optimal weight: 5.9990 overall best weight: 1.9386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7136 moved from start: 0.4095 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.061 3516 Z= 0.290 Angle : 0.538 5.336 4788 Z= 0.300 Chirality : 0.058 0.447 606 Planarity : 0.003 0.032 579 Dihedral : 12.729 107.853 531 Min Nonbonded Distance : 2.252 Molprobity Statistics. All-atom Clashscore : 8.43 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.27 % Favored : 93.73 % Rotamer Outliers : 1.84 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.66 (0.38), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.02 (0.29), residues: 402 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 44 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 37 time to evaluate : 0.413 Fit side-chains revert: symmetry clash outliers start: 7 outliers final: 6 residues processed: 43 average time/residue: 1.5730 time to fit residues: 69.4768 Evaluate side-chains 43 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 37 time to evaluate : 0.390 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 5 residues processed: 1 average time/residue: 0.0345 time to fit residues: 0.6072 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 39 random chunks: chunk 33 optimal weight: 5.9990 chunk 34 optimal weight: 0.7980 chunk 23 optimal weight: 1.9990 chunk 37 optimal weight: 0.8980 chunk 22 optimal weight: 7.9990 chunk 17 optimal weight: 1.9990 chunk 25 optimal weight: 4.9990 chunk 38 optimal weight: 2.9990 chunk 35 optimal weight: 0.9980 chunk 30 optimal weight: 2.9990 chunk 3 optimal weight: 3.9990 overall best weight: 1.3384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7117 moved from start: 0.4048 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.061 3516 Z= 0.230 Angle : 0.509 5.226 4788 Z= 0.284 Chirality : 0.056 0.386 606 Planarity : 0.003 0.033 579 Dihedral : 12.697 108.397 531 Min Nonbonded Distance : 2.291 Molprobity Statistics. All-atom Clashscore : 7.27 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.52 % Favored : 93.48 % Rotamer Outliers : 1.31 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.59 (0.38), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.96 (0.29), residues: 402 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 44 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 39 time to evaluate : 0.413 Fit side-chains revert: symmetry clash outliers start: 5 outliers final: 5 residues processed: 44 average time/residue: 1.5225 time to fit residues: 69.0304 Evaluate side-chains 43 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 38 time to evaluate : 0.491 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 5 residues processed: 0 time to fit residues: 0.5888 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 39 random chunks: chunk 23 optimal weight: 6.9990 chunk 18 optimal weight: 4.9990 chunk 24 optimal weight: 0.8980 chunk 32 optimal weight: 3.9990 chunk 9 optimal weight: 4.9990 chunk 28 optimal weight: 0.0970 chunk 4 optimal weight: 4.9990 chunk 8 optimal weight: 0.3980 chunk 30 optimal weight: 7.9990 chunk 12 optimal weight: 4.9990 chunk 31 optimal weight: 2.9990 overall best weight: 1.6782 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7129 moved from start: 0.4118 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.060 3516 Z= 0.262 Angle : 0.525 5.263 4788 Z= 0.293 Chirality : 0.057 0.413 606 Planarity : 0.003 0.032 579 Dihedral : 12.702 108.234 531 Min Nonbonded Distance : 2.269 Molprobity Statistics. All-atom Clashscore : 7.70 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.27 % Favored : 92.73 % Rotamer Outliers : 1.84 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.61 (0.38), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.98 (0.29), residues: 402 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 804 Ramachandran restraints generated. 402 Oldfield, 0 Emsley, 402 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 46 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 39 time to evaluate : 0.375 Fit side-chains revert: symmetry clash outliers start: 7 outliers final: 7 residues processed: 44 average time/residue: 1.4189 time to fit residues: 64.2793 Evaluate side-chains 46 residues out of total 375 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 39 time to evaluate : 0.357 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 5 residues processed: 2 average time/residue: 0.0352 time to fit residues: 0.6691 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 39 random chunks: chunk 3 optimal weight: 4.9990 chunk 5 optimal weight: 10.0000 chunk 27 optimal weight: 7.9990 chunk 1 optimal weight: 0.7980 chunk 22 optimal weight: 0.9980 chunk 35 optimal weight: 0.5980 chunk 20 optimal weight: 2.9990 chunk 26 optimal weight: 1.9990 chunk 24 optimal weight: 0.9990 chunk 23 optimal weight: 0.9990 chunk 15 optimal weight: 1.9990 overall best weight: 0.8784 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4364 r_free = 0.4364 target = 0.132144 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 40)----------------| | r_work = 0.4104 r_free = 0.4104 target = 0.111656 restraints weight = 10837.804| |-----------------------------------------------------------------------------| r_work (start): 0.4107 rms_B_bonded: 4.20 r_work: 0.3926 rms_B_bonded: 4.08 restraints_weight: 0.5000 r_work (final): 0.3926 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.4364 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.4364 r_free = 0.4364 target_work(ls_wunit_k1) = 0.132 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.4364 r_free = 0.4364 target_work(ls_wunit_k1) = 0.132 | | occupancies: max = 1.00 min = 0.49 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| r_final: 0.4364 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7958 moved from start: 0.4073 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.061 3516 Z= 0.189 Angle : 0.489 5.179 4788 Z= 0.274 Chirality : 0.054 0.334 606 Planarity : 0.003 0.033 579 Dihedral : 12.674 108.970 531 Min Nonbonded Distance : 2.142 Molprobity Statistics. All-atom Clashscore : 7.27 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.76 % Favored : 95.24 % Rotamer Outliers : 1.57 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.57 (0.38), residues: 402 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.95 (0.29), residues: 402 =============================================================================== Job complete usr+sys time: 1818.74 seconds wall clock time: 33 minutes 10.34 seconds (1990.34 seconds total)