Starting phenix.real_space_refine (version: 1.21rc1) on Tue Oct 10 12:39:07 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8d45_27175/10_2023/8d45_27175.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8d45_27175/10_2023/8d45_27175.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.62 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8d45_27175/10_2023/8d45_27175.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8d45_27175/10_2023/8d45_27175.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8d45_27175/10_2023/8d45_27175.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8d45_27175/10_2023/8d45_27175.pdb" } resolution = 2.62 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= -0.000 sd= 0.052 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 5 Type Number sf(0) Gaussians S 56 5.16 5 C 6716 2.51 5 N 1846 2.21 5 O 1944 1.98 5 H 10370 0.53 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A GLU 21": "OE1" <-> "OE2" Residue "A GLU 47": "OE1" <-> "OE2" Residue "A GLU 159": "OE1" <-> "OE2" Residue "A GLU 317": "OE1" <-> "OE2" Residue "B GLU 303": "OE1" <-> "OE2" Residue "B GLU 317": "OE1" <-> "OE2" Residue "C GLU 21": "OE1" <-> "OE2" Residue "C GLU 47": "OE1" <-> "OE2" Residue "D GLU 21": "OE1" <-> "OE2" Residue "D GLU 303": "OE1" <-> "OE2" Residue "D GLU 317": "OE1" <-> "OE2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5097/modules/chem_data/mon_lib" Total number of atoms: 20932 Number of models: 1 Model: "" Number of chains: 4 Chain: "A" Number of atoms: 5225 Number of conformers: 2 Conformer: "A" Number of residues, atoms: 339, 5201 Classifications: {'peptide': 339} Link IDs: {'PCIS': 1, 'PTRANS': 21, 'TRANS': 316} Chain breaks: 2 Conformer: "B" Number of residues, atoms: 339, 5201 Classifications: {'peptide': 339} Link IDs: {'PCIS': 1, 'PTRANS': 21, 'TRANS': 316} Chain breaks: 2 bond proxies already assigned to first conformer: 5239 Chain: "B" Number of atoms: 5225 Number of conformers: 2 Conformer: "A" Number of residues, atoms: 339, 5201 Classifications: {'peptide': 339} Link IDs: {'PCIS': 1, 'PTRANS': 21, 'TRANS': 316} Chain breaks: 2 Conformer: "B" Number of residues, atoms: 339, 5201 Classifications: {'peptide': 339} Link IDs: {'PCIS': 1, 'PTRANS': 21, 'TRANS': 316} Chain breaks: 2 bond proxies already assigned to first conformer: 5239 Chain: "C" Number of atoms: 5225 Number of conformers: 2 Conformer: "A" Number of residues, atoms: 339, 5201 Classifications: {'peptide': 339} Link IDs: {'PCIS': 1, 'PTRANS': 21, 'TRANS': 316} Chain breaks: 2 Conformer: "B" Number of residues, atoms: 339, 5201 Classifications: {'peptide': 339} Link IDs: {'PCIS': 1, 'PTRANS': 21, 'TRANS': 316} Chain breaks: 2 bond proxies already assigned to first conformer: 5239 Chain: "D" Number of atoms: 5257 Number of conformers: 2 Conformer: "A" Number of residues, atoms: 341, 5233 Classifications: {'peptide': 341} Link IDs: {'PCIS': 1, 'PTRANS': 21, 'TRANS': 318} Chain breaks: 2 Conformer: "B" Number of residues, atoms: 341, 5233 Classifications: {'peptide': 341} Link IDs: {'PCIS': 1, 'PTRANS': 21, 'TRANS': 318} Chain breaks: 2 bond proxies already assigned to first conformer: 5272 Time building chain proxies: 16.22, per 1000 atoms: 0.77 Number of scatterers: 20932 At special positions: 0 Unit cell: (103.79, 63.13, 110.21, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) S 56 16.00 O 1944 8.00 N 1846 7.00 C 6716 6.00 H 10370 1.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 20.24 Conformation dependent library (CDL) restraints added in 3.6 seconds 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2476 Finding SS restraints... Secondary structure from input PDB file: 62 helices and 4 sheets defined 54.0% alpha, 6.7% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.64 Creating SS restraints... Processing helix chain 'A' and resid 12 through 20 Processing helix chain 'A' and resid 28 through 36 Processing helix chain 'A' and resid 47 through 52 Processing helix chain 'A' and resid 52 through 67 Processing helix chain 'A' and resid 96 through 113 removed outlier: 3.752A pdb=" N ASN A 100 " --> pdb=" O GLY A 96 " (cutoff:3.500A) removed outlier: 3.669A pdb=" N GLU A 113 " --> pdb=" O GLN A 109 " (cutoff:3.500A) Processing helix chain 'A' and resid 134 through 152 removed outlier: 3.906A pdb=" N VAL A 147 " --> pdb=" O GLN A 143 " (cutoff:3.500A) removed outlier: 3.629A pdb=" N PHE A 148 " --> pdb=" O GLN A 144 " (cutoff:3.500A) Processing helix chain 'A' and resid 163 through 176 removed outlier: 3.700A pdb=" N ALA A 176 " --> pdb=" O GLU A 172 " (cutoff:3.500A) Processing helix chain 'A' and resid 197 through 209 Processing helix chain 'A' and resid 220 through 239 Processing helix chain 'A' and resid 260 through 264 removed outlier: 3.853A pdb=" N ASP A 263 " --> pdb=" O GLY A 260 " (cutoff:3.500A) Processing helix chain 'A' and resid 277 through 292 removed outlier: 3.739A pdb=" N ARG A 286 " --> pdb=" O GLN A 282 " (cutoff:3.500A) Processing helix chain 'A' and resid 303 through 319 Proline residue: A 316 - end of helix Processing helix chain 'A' and resid 322 through 328 Processing helix chain 'A' and resid 332 through 338 removed outlier: 4.006A pdb=" N MET A 337 " --> pdb=" O GLY A 334 " (cutoff:3.500A) Processing helix chain 'A' and resid 340 through 346 Processing helix chain 'A' and resid 348 through 354 Processing helix chain 'B' and resid 12 through 20 Processing helix chain 'B' and resid 28 through 33 Processing helix chain 'B' and resid 47 through 52 Processing helix chain 'B' and resid 52 through 67 Processing helix chain 'B' and resid 96 through 112 Processing helix chain 'B' and resid 134 through 151 removed outlier: 3.871A pdb=" N VAL B 147 " --> pdb=" O GLN B 143 " (cutoff:3.500A) removed outlier: 3.606A pdb=" N PHE B 148 " --> pdb=" O GLN B 144 " (cutoff:3.500A) Processing helix chain 'B' and resid 163 through 176 removed outlier: 3.535A pdb=" N ALA B 167 " --> pdb=" O HIS B 163 " (cutoff:3.500A) removed outlier: 3.539A pdb=" N TRP B 169 " --> pdb=" O GLU B 165 " (cutoff:3.500A) removed outlier: 3.687A pdb=" N ALA B 176 " --> pdb=" O GLU B 172 " (cutoff:3.500A) Processing helix chain 'B' and resid 197 through 209 Processing helix chain 'B' and resid 222 through 239 Processing helix chain 'B' and resid 277 through 292 removed outlier: 4.011A pdb=" N ARG B 286 " --> pdb=" O GLN B 282 " (cutoff:3.500A) Processing helix chain 'B' and resid 303 through 319 Proline residue: B 316 - end of helix Processing helix chain 'B' and resid 322 through 328 Processing helix chain 'B' and resid 332 through 338 removed outlier: 3.959A pdb=" N MET B 337 " --> pdb=" O GLY B 334 " (cutoff:3.500A) Processing helix chain 'B' and resid 340 through 346 Processing helix chain 'B' and resid 348 through 355 removed outlier: 3.515A pdb=" N TRP B 352 " --> pdb=" O ARG B 348 " (cutoff:3.500A) Processing helix chain 'C' and resid 12 through 19 Processing helix chain 'C' and resid 28 through 36 removed outlier: 3.708A pdb=" N GLU C 34 " --> pdb=" O VAL C 30 " (cutoff:3.500A) Processing helix chain 'C' and resid 47 through 52 Processing helix chain 'C' and resid 52 through 67 Processing helix chain 'C' and resid 96 through 113 removed outlier: 3.644A pdb=" N ASN C 100 " --> pdb=" O GLY C 96 " (cutoff:3.500A) removed outlier: 3.646A pdb=" N GLU C 113 " --> pdb=" O GLN C 109 " (cutoff:3.500A) Processing helix chain 'C' and resid 134 through 152 removed outlier: 3.902A pdb=" N VAL C 147 " --> pdb=" O GLN C 143 " (cutoff:3.500A) removed outlier: 3.608A pdb=" N PHE C 148 " --> pdb=" O GLN C 144 " (cutoff:3.500A) Processing helix chain 'C' and resid 163 through 176 removed outlier: 3.659A pdb=" N ALA C 167 " --> pdb=" O HIS C 163 " (cutoff:3.500A) removed outlier: 3.729A pdb=" N ALA C 176 " --> pdb=" O GLU C 172 " (cutoff:3.500A) Processing helix chain 'C' and resid 197 through 209 Processing helix chain 'C' and resid 222 through 238 Processing helix chain 'C' and resid 279 through 292 removed outlier: 4.000A pdb=" N ARG C 286 " --> pdb=" O GLN C 282 " (cutoff:3.500A) Processing helix chain 'C' and resid 303 through 319 Proline residue: C 316 - end of helix Processing helix chain 'C' and resid 322 through 328 Processing helix chain 'C' and resid 332 through 336 removed outlier: 3.536A pdb=" N ASP C 336 " --> pdb=" O SER C 333 " (cutoff:3.500A) Processing helix chain 'C' and resid 340 through 346 Processing helix chain 'C' and resid 348 through 355 removed outlier: 3.574A pdb=" N TRP C 352 " --> pdb=" O ARG C 348 " (cutoff:3.500A) Processing helix chain 'D' and resid 12 through 20 Processing helix chain 'D' and resid 28 through 33 Processing helix chain 'D' and resid 47 through 52 Processing helix chain 'D' and resid 52 through 67 Processing helix chain 'D' and resid 96 through 112 Processing helix chain 'D' and resid 134 through 152 removed outlier: 3.871A pdb=" N VAL D 147 " --> pdb=" O GLN D 143 " (cutoff:3.500A) removed outlier: 3.613A pdb=" N PHE D 148 " --> pdb=" O GLN D 144 " (cutoff:3.500A) Processing helix chain 'D' and resid 163 through 176 removed outlier: 3.557A pdb=" N ALA D 167 " --> pdb=" O HIS D 163 " (cutoff:3.500A) removed outlier: 3.731A pdb=" N ALA D 176 " --> pdb=" O GLU D 172 " (cutoff:3.500A) Processing helix chain 'D' and resid 197 through 209 Processing helix chain 'D' and resid 220 through 238 Processing helix chain 'D' and resid 259 through 264 removed outlier: 3.845A pdb=" N ILE D 262 " --> pdb=" O GLN D 259 " (cutoff:3.500A) removed outlier: 3.870A pdb=" N ASP D 263 " --> pdb=" O GLY D 260 " (cutoff:3.500A) Processing helix chain 'D' and resid 277 through 292 removed outlier: 4.009A pdb=" N ARG D 286 " --> pdb=" O GLN D 282 " (cutoff:3.500A) Processing helix chain 'D' and resid 303 through 319 Proline residue: D 316 - end of helix Processing helix chain 'D' and resid 322 through 328 Processing helix chain 'D' and resid 332 through 338 removed outlier: 3.922A pdb=" N MET D 337 " --> pdb=" O GLY D 334 " (cutoff:3.500A) Processing helix chain 'D' and resid 340 through 345 Processing helix chain 'D' and resid 348 through 355 removed outlier: 3.697A pdb=" N TRP D 352 " --> pdb=" O ARG D 348 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 122 through 123 removed outlier: 5.833A pdb=" N VAL A 122 " --> pdb=" O GLU A 159 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'B' and resid 70 through 72 removed outlier: 3.569A pdb=" N ALA B 119 " --> pdb=" O MET B 71 " (cutoff:3.500A) removed outlier: 7.157A pdb=" N LEU B 244 " --> pdb=" O TYR B 295 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'C' and resid 70 through 72 removed outlier: 6.619A pdb=" N LEU C 156 " --> pdb=" O ALA C 182 " (cutoff:3.500A) removed outlier: 7.934A pdb=" N THR C 184 " --> pdb=" O LEU C 156 " (cutoff:3.500A) removed outlier: 6.061A pdb=" N ALA C 158 " --> pdb=" O THR C 184 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'D' and resid 70 through 72 518 hydrogen bonds defined for protein. 1440 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: 7.84 Time building geometry restraints manager: 20.69 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 0.82 - 1.02: 10342 1.02 - 1.22: 33 1.22 - 1.42: 4492 1.42 - 1.61: 6246 1.61 - 1.81: 84 Bond restraints: 21197 Sorted by residual: bond pdb=" NE ARG B 108 " pdb=" HE ARG B 108 " ideal model delta sigma weight residual 0.860 0.971 -0.111 2.00e-02 2.50e+03 3.08e+01 bond pdb=" C TRP B 44 " pdb=" N THR B 45 " ideal model delta sigma weight residual 1.331 1.384 -0.053 2.07e-02 2.33e+03 6.62e+00 bond pdb=" C TRP C 44 " pdb=" N THR C 45 " ideal model delta sigma weight residual 1.331 1.378 -0.047 2.07e-02 2.33e+03 5.09e+00 bond pdb=" N THR B 222 " pdb=" CA THR B 222 " ideal model delta sigma weight residual 1.457 1.483 -0.026 1.25e-02 6.40e+03 4.23e+00 bond pdb=" N THR C 222 " pdb=" CA THR C 222 " ideal model delta sigma weight residual 1.458 1.481 -0.023 1.27e-02 6.20e+03 3.26e+00 ... (remaining 21192 not shown) Histogram of bond angle deviations from ideal: 97.73 - 104.98: 292 104.98 - 112.24: 24061 112.24 - 119.49: 5709 119.49 - 126.75: 8125 126.75 - 134.00: 205 Bond angle restraints: 38392 Sorted by residual: angle pdb=" CG LYS A 139 " pdb=" CD LYS A 139 " pdb=" CE LYS A 139 " ideal model delta sigma weight residual 111.30 118.65 -7.35 2.30e+00 1.89e-01 1.02e+01 angle pdb=" N ASN A 152 " pdb=" CA ASN A 152 " pdb=" C ASN A 152 " ideal model delta sigma weight residual 110.91 107.56 3.35 1.17e+00 7.31e-01 8.18e+00 angle pdb=" N ARG A 239 " pdb=" CA ARG A 239 " pdb=" C ARG A 239 " ideal model delta sigma weight residual 111.24 107.32 3.92 1.38e+00 5.25e-01 8.09e+00 angle pdb=" N ARG B 239 " pdb=" CA ARG B 239 " pdb=" C ARG B 239 " ideal model delta sigma weight residual 111.24 107.40 3.84 1.38e+00 5.25e-01 7.74e+00 angle pdb=" N ASN C 152 " pdb=" CA ASN C 152 " pdb=" C ASN C 152 " ideal model delta sigma weight residual 110.91 107.68 3.23 1.17e+00 7.31e-01 7.63e+00 ... (remaining 38387 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.62: 7373 17.62 - 35.24: 774 35.24 - 52.87: 206 52.87 - 70.49: 40 70.49 - 88.11: 27 Dihedral angle restraints: 8420 sinusoidal: 3853 harmonic: 4567 Sorted by residual: dihedral pdb=" CA ASP C 336 " pdb=" C ASP C 336 " pdb=" N MET C 337 " pdb=" CA MET C 337 " ideal model delta harmonic sigma weight residual 180.00 155.09 24.91 0 5.00e+00 4.00e-02 2.48e+01 dihedral pdb=" CA PHE C 74 " pdb=" C PHE C 74 " pdb=" N THR C 75 " pdb=" CA THR C 75 " ideal model delta harmonic sigma weight residual -180.00 -162.49 -17.51 0 5.00e+00 4.00e-02 1.23e+01 dihedral pdb=" CA PHE D 74 " pdb=" C PHE D 74 " pdb=" N THR D 75 " pdb=" CA THR D 75 " ideal model delta harmonic sigma weight residual 180.00 -162.95 -17.05 0 5.00e+00 4.00e-02 1.16e+01 ... (remaining 8417 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.029: 853 0.029 - 0.058: 422 0.058 - 0.087: 176 0.087 - 0.115: 80 0.115 - 0.144: 23 Chirality restraints: 1554 Sorted by residual: chirality pdb=" CA PRO A 46 " pdb=" N PRO A 46 " pdb=" C PRO A 46 " pdb=" CB PRO A 46 " both_signs ideal model delta sigma weight residual False 2.72 2.57 0.14 2.00e-01 2.50e+01 5.21e-01 chirality pdb=" CA PRO B 46 " pdb=" N PRO B 46 " pdb=" C PRO B 46 " pdb=" CB PRO B 46 " both_signs ideal model delta sigma weight residual False 2.72 2.58 0.14 2.00e-01 2.50e+01 5.09e-01 chirality pdb=" CA VAL A 196 " pdb=" N VAL A 196 " pdb=" C VAL A 196 " pdb=" CB VAL A 196 " both_signs ideal model delta sigma weight residual False 2.44 2.58 -0.14 2.00e-01 2.50e+01 5.02e-01 ... (remaining 1551 not shown) Planarity restraints: 3201 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CG GLN B 109 " -0.077 2.00e-02 2.50e+03 8.49e-02 1.08e+02 pdb=" CD GLN B 109 " 0.004 2.00e-02 2.50e+03 pdb=" OE1 GLN B 109 " 0.077 2.00e-02 2.50e+03 pdb=" NE2 GLN B 109 " -0.004 2.00e-02 2.50e+03 pdb="HE21 GLN B 109 " 0.125 2.00e-02 2.50e+03 pdb="HE22 GLN B 109 " -0.125 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C ALA D 315 " 0.048 5.00e-02 4.00e+02 7.23e-02 8.36e+00 pdb=" N PRO D 316 " -0.125 5.00e-02 4.00e+02 pdb=" CA PRO D 316 " 0.036 5.00e-02 4.00e+02 pdb=" CD PRO D 316 " 0.041 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C ALA C 315 " -0.048 5.00e-02 4.00e+02 7.16e-02 8.19e+00 pdb=" N PRO C 316 " 0.124 5.00e-02 4.00e+02 pdb=" CA PRO C 316 " -0.036 5.00e-02 4.00e+02 pdb=" CD PRO C 316 " -0.040 5.00e-02 4.00e+02 ... (remaining 3198 not shown) Histogram of nonbonded interaction distances: 1.61 - 2.21: 1588 2.21 - 2.81: 45468 2.81 - 3.41: 62571 3.41 - 4.00: 84710 4.00 - 4.60: 129854 Nonbonded interactions: 324191 Sorted by model distance: nonbonded pdb=" H ALA C 158 " pdb=" O ALA C 182 " model vdw 1.613 1.850 nonbonded pdb=" O ALA C 158 " pdb=" H THR C 184 " model vdw 1.626 1.850 nonbonded pdb=" O ILE D 24 " pdb=" H ASN D 69 " model vdw 1.686 1.850 nonbonded pdb=" O ILE A 357 " pdb="HD22 ASN B 364 " model vdw 1.686 1.850 nonbonded pdb=" O ILE B 24 " pdb=" H ASN B 69 " model vdw 1.700 1.850 ... (remaining 324186 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 9 through 14 or resid 16 through 325 or resid 327 through \ 372)) selection = (chain 'B' and (resid 9 through 14 or resid 16 through 325 or resid 327 through \ 372)) selection = (chain 'C' and (resid 9 through 14 or resid 16 through 325 or resid 327 through \ 372)) selection = (chain 'D' and (resid 9 through 14 or resid 16 through 127 or resid 134 through \ 325 or resid 327 through 372)) } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=0.40 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 1.320 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.320 Extract box with map and model: 7.060 Check model and map are aligned: 0.290 Set scattering table: 0.180 Process input model: 81.290 Find NCS groups from input model: 1.100 Set up NCS constraints: 0.100 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.010 Load rotamer database and sin/cos tables:11.450 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 103.120 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7063 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.004 0.053 10827 Z= 0.288 Angle : 0.788 7.348 14666 Z= 0.445 Chirality : 0.044 0.144 1554 Planarity : 0.007 0.072 1935 Dihedral : 17.508 88.110 4022 Min Nonbonded Distance : 2.203 Molprobity Statistics. All-atom Clashscore : 4.67 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.25 % Favored : 94.75 % Rotamer: Outliers : 0.92 % Allowed : 18.26 % Favored : 80.83 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.55 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.43 (0.22), residues: 1342 helix: -0.69 (0.19), residues: 676 sheet: -1.97 (0.32), residues: 132 loop : -0.76 (0.28), residues: 534 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 146 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 138 time to evaluate : 1.701 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 8 outliers final: 7 residues processed: 145 average time/residue: 2.0533 time to fit residues: 329.9112 Evaluate side-chains 138 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 131 time to evaluate : 1.706 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 7 residues processed: 0 time to fit residues: 2.0723 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 112 optimal weight: 2.9990 chunk 100 optimal weight: 1.9990 chunk 55 optimal weight: 3.9990 chunk 34 optimal weight: 1.9990 chunk 67 optimal weight: 1.9990 chunk 53 optimal weight: 3.9990 chunk 104 optimal weight: 0.2980 chunk 40 optimal weight: 3.9990 chunk 63 optimal weight: 7.9990 chunk 77 optimal weight: 2.9990 chunk 120 optimal weight: 1.9990 overall best weight: 1.6588 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 152 ASN A 290 ASN B 152 ASN B 290 ASN C 152 ASN C 354 ASN D 152 ASN D 259 GLN D 290 ASN Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7107 moved from start: 0.0874 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.066 10827 Z= 0.348 Angle : 0.689 5.294 14666 Z= 0.361 Chirality : 0.044 0.152 1554 Planarity : 0.007 0.059 1935 Dihedral : 4.827 20.509 1475 Min Nonbonded Distance : 2.094 Molprobity Statistics. All-atom Clashscore : 7.49 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.22 % Favored : 93.78 % Rotamer: Outliers : 3.03 % Allowed : 15.78 % Favored : 81.19 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.55 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.37 (0.22), residues: 1342 helix: -0.56 (0.19), residues: 652 sheet: -2.18 (0.32), residues: 132 loop : -0.78 (0.27), residues: 558 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 156 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 31 poor density : 125 time to evaluate : 1.678 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 31 outliers final: 12 residues processed: 144 average time/residue: 2.0855 time to fit residues: 331.5387 Evaluate side-chains 140 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 128 time to evaluate : 1.709 Switching outliers to nearest non-outliers outliers start: 12 outliers final: 8 residues processed: 4 average time/residue: 1.0660 time to fit residues: 7.1377 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 67 optimal weight: 1.9990 chunk 37 optimal weight: 0.9980 chunk 100 optimal weight: 1.9990 chunk 82 optimal weight: 0.0870 chunk 33 optimal weight: 2.9990 chunk 120 optimal weight: 0.9980 chunk 130 optimal weight: 2.9990 chunk 107 optimal weight: 1.9990 chunk 119 optimal weight: 2.9990 chunk 41 optimal weight: 0.9980 chunk 97 optimal weight: 0.6980 overall best weight: 0.7558 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 152 ASN C 152 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7060 moved from start: 0.1098 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.043 10827 Z= 0.220 Angle : 0.600 4.820 14666 Z= 0.308 Chirality : 0.041 0.138 1554 Planarity : 0.005 0.048 1935 Dihedral : 4.500 18.007 1475 Min Nonbonded Distance : 2.083 Molprobity Statistics. All-atom Clashscore : 6.68 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.02 % Favored : 94.98 % Rotamer: Outliers : 2.48 % Allowed : 15.96 % Favored : 81.56 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.55 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.72 (0.22), residues: 1342 helix: 0.05 (0.20), residues: 650 sheet: -2.24 (0.33), residues: 138 loop : -0.43 (0.28), residues: 554 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 159 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 25 poor density : 134 time to evaluate : 1.624 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 25 outliers final: 12 residues processed: 149 average time/residue: 1.9912 time to fit residues: 328.9841 Evaluate side-chains 137 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 125 time to evaluate : 1.968 Switching outliers to nearest non-outliers outliers start: 12 outliers final: 8 residues processed: 4 average time/residue: 0.6421 time to fit residues: 5.5273 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 119 optimal weight: 0.9980 chunk 90 optimal weight: 0.9980 chunk 62 optimal weight: 5.9990 chunk 13 optimal weight: 2.9990 chunk 57 optimal weight: 0.9990 chunk 81 optimal weight: 1.9990 chunk 121 optimal weight: 2.9990 chunk 128 optimal weight: 0.9990 chunk 63 optimal weight: 0.6980 chunk 115 optimal weight: 1.9990 chunk 34 optimal weight: 1.9990 overall best weight: 0.9384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 152 ASN C 152 ASN C 163 HIS Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7078 moved from start: 0.1383 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.042 10827 Z= 0.240 Angle : 0.598 4.749 14666 Z= 0.306 Chirality : 0.042 0.138 1554 Planarity : 0.005 0.047 1935 Dihedral : 4.401 18.446 1475 Min Nonbonded Distance : 2.052 Molprobity Statistics. All-atom Clashscore : 6.72 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.55 % Favored : 94.45 % Rotamer: Outliers : 2.39 % Allowed : 16.79 % Favored : 80.83 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.55 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.56 (0.23), residues: 1342 helix: 0.20 (0.20), residues: 656 sheet: -2.23 (0.33), residues: 138 loop : -0.37 (0.28), residues: 548 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 157 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 24 poor density : 133 time to evaluate : 1.816 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 24 outliers final: 11 residues processed: 149 average time/residue: 1.9648 time to fit residues: 328.0880 Evaluate side-chains 136 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 125 time to evaluate : 1.659 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 8 residues processed: 3 average time/residue: 0.8167 time to fit residues: 5.2345 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 107 optimal weight: 1.9990 chunk 72 optimal weight: 2.9990 chunk 1 optimal weight: 0.7980 chunk 95 optimal weight: 3.9990 chunk 53 optimal weight: 2.9990 chunk 109 optimal weight: 8.9990 chunk 88 optimal weight: 0.0020 chunk 0 optimal weight: 1.9990 chunk 65 optimal weight: 1.9990 chunk 115 optimal weight: 2.9990 chunk 32 optimal weight: 0.5980 overall best weight: 1.0792 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 152 ASN C 152 ASN C 163 HIS Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7080 moved from start: 0.1564 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.036 10827 Z= 0.256 Angle : 0.600 5.076 14666 Z= 0.308 Chirality : 0.042 0.138 1554 Planarity : 0.006 0.087 1935 Dihedral : 4.400 19.701 1475 Min Nonbonded Distance : 2.032 Molprobity Statistics. All-atom Clashscore : 7.15 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.32 % Favored : 94.68 % Rotamer: Outliers : 3.03 % Allowed : 16.06 % Favored : 80.92 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.55 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.37 (0.23), residues: 1342 helix: 0.36 (0.20), residues: 635 sheet: -2.23 (0.33), residues: 133 loop : -0.26 (0.28), residues: 574 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 162 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 31 poor density : 131 time to evaluate : 1.723 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 31 outliers final: 12 residues processed: 153 average time/residue: 2.0718 time to fit residues: 350.8894 Evaluate side-chains 140 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 128 time to evaluate : 1.774 Switching outliers to nearest non-outliers outliers start: 12 outliers final: 8 residues processed: 4 average time/residue: 0.6886 time to fit residues: 5.7900 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 43 optimal weight: 0.9990 chunk 115 optimal weight: 1.9990 chunk 25 optimal weight: 1.9990 chunk 75 optimal weight: 2.9990 chunk 31 optimal weight: 3.9990 chunk 128 optimal weight: 0.0980 chunk 106 optimal weight: 1.9990 chunk 59 optimal weight: 1.9990 chunk 10 optimal weight: 0.5980 chunk 42 optimal weight: 1.9990 chunk 67 optimal weight: 1.9990 overall best weight: 1.1386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 152 ASN C 152 ASN C 163 HIS D 354 ASN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7090 moved from start: 0.1736 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.042 10827 Z= 0.268 Angle : 0.604 5.288 14666 Z= 0.310 Chirality : 0.042 0.140 1554 Planarity : 0.006 0.088 1935 Dihedral : 4.395 20.170 1475 Min Nonbonded Distance : 2.022 Molprobity Statistics. All-atom Clashscore : 7.58 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.00 % Favored : 94.00 % Rotamer: Outliers : 2.75 % Allowed : 16.79 % Favored : 80.46 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.55 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.35 (0.23), residues: 1342 helix: 0.40 (0.20), residues: 635 sheet: -2.27 (0.33), residues: 133 loop : -0.25 (0.28), residues: 574 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 162 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 28 poor density : 134 time to evaluate : 1.868 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 28 outliers final: 13 residues processed: 151 average time/residue: 2.0112 time to fit residues: 336.2064 Evaluate side-chains 142 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 129 time to evaluate : 1.710 Switching outliers to nearest non-outliers outliers start: 13 outliers final: 8 residues processed: 5 average time/residue: 0.6103 time to fit residues: 6.1856 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 124 optimal weight: 0.9980 chunk 14 optimal weight: 0.2980 chunk 73 optimal weight: 2.9990 chunk 93 optimal weight: 1.9990 chunk 72 optimal weight: 1.9990 chunk 108 optimal weight: 2.9990 chunk 71 optimal weight: 0.8980 chunk 128 optimal weight: 0.8980 chunk 80 optimal weight: 0.9990 chunk 78 optimal weight: 0.9980 chunk 59 optimal weight: 1.9990 overall best weight: 0.8180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 152 ASN C 163 HIS Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7072 moved from start: 0.1842 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.036 10827 Z= 0.223 Angle : 0.583 5.748 14666 Z= 0.298 Chirality : 0.041 0.138 1554 Planarity : 0.005 0.085 1935 Dihedral : 4.271 19.952 1475 Min Nonbonded Distance : 2.019 Molprobity Statistics. All-atom Clashscore : 7.44 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.02 % Favored : 94.98 % Rotamer: Outliers : 1.93 % Allowed : 17.34 % Favored : 80.73 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.55 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.12 (0.23), residues: 1342 helix: 0.61 (0.20), residues: 635 sheet: -2.31 (0.33), residues: 133 loop : -0.12 (0.28), residues: 574 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 144 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 19 poor density : 125 time to evaluate : 2.098 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 19 outliers final: 7 residues processed: 138 average time/residue: 2.0440 time to fit residues: 312.6270 Evaluate side-chains 131 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 124 time to evaluate : 1.602 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 6 residues processed: 1 average time/residue: 0.3824 time to fit residues: 2.6314 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 79 optimal weight: 1.9990 chunk 51 optimal weight: 2.9990 chunk 76 optimal weight: 3.9990 chunk 38 optimal weight: 1.9990 chunk 25 optimal weight: 1.9990 chunk 24 optimal weight: 0.9990 chunk 81 optimal weight: 2.9990 chunk 87 optimal weight: 1.9990 chunk 63 optimal weight: 0.9990 chunk 11 optimal weight: 1.9990 chunk 100 optimal weight: 1.9990 overall best weight: 1.5990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 152 ASN C 163 HIS C 364 ASN D 259 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7112 moved from start: 0.1929 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.056 10827 Z= 0.339 Angle : 0.645 5.931 14666 Z= 0.334 Chirality : 0.044 0.150 1554 Planarity : 0.006 0.083 1935 Dihedral : 4.510 20.965 1475 Min Nonbonded Distance : 2.014 Molprobity Statistics. All-atom Clashscore : 8.44 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.15 % Favored : 93.85 % Rotamer: Outliers : 1.65 % Allowed : 17.61 % Favored : 80.73 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.55 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.54 (0.23), residues: 1342 helix: 0.21 (0.20), residues: 635 sheet: -2.29 (0.33), residues: 133 loop : -0.33 (0.28), residues: 574 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 143 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 16 poor density : 127 time to evaluate : 1.725 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 16 outliers final: 7 residues processed: 137 average time/residue: 2.0816 time to fit residues: 316.6515 Evaluate side-chains 131 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 124 time to evaluate : 1.613 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 6 residues processed: 1 average time/residue: 1.6437 time to fit residues: 3.8322 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 116 optimal weight: 2.9990 chunk 122 optimal weight: 0.5980 chunk 112 optimal weight: 0.9980 chunk 119 optimal weight: 0.9990 chunk 71 optimal weight: 1.9990 chunk 52 optimal weight: 0.9990 chunk 93 optimal weight: 0.4980 chunk 36 optimal weight: 0.2980 chunk 107 optimal weight: 1.9990 chunk 78 optimal weight: 1.9990 chunk 126 optimal weight: 0.9980 overall best weight: 0.6780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 152 ASN C 163 HIS D 259 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7066 moved from start: 0.2004 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.037 10827 Z= 0.211 Angle : 0.584 6.263 14666 Z= 0.298 Chirality : 0.041 0.139 1554 Planarity : 0.005 0.083 1935 Dihedral : 4.265 19.940 1475 Min Nonbonded Distance : 2.039 Molprobity Statistics. All-atom Clashscore : 8.20 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.80 % Favored : 95.20 % Rotamer: Outliers : 1.38 % Allowed : 18.35 % Favored : 80.28 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.55 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.06 (0.23), residues: 1342 helix: 0.67 (0.20), residues: 635 sheet: -2.33 (0.34), residues: 133 loop : -0.09 (0.28), residues: 574 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2684 Ramachandran restraints generated. 1342 Oldfield, 0 Emsley, 1342 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 137 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 124 time to evaluate : 1.671 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 13 outliers final: 5 residues processed: 136 average time/residue: 2.1460 time to fit residues: 321.9563 Evaluate side-chains 129 residues out of total 1074 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 124 time to evaluate : 1.700 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 5 residues processed: 0 time to fit residues: 2.1435 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 77 optimal weight: 0.4980 chunk 59 optimal weight: 1.9990 chunk 87 optimal weight: 1.9990 chunk 132 optimal weight: 0.7980 chunk 121 optimal weight: 0.7980 chunk 105 optimal weight: 2.9990 chunk 10 optimal weight: 3.9990 chunk 81 optimal weight: 2.9990 chunk 64 optimal weight: 0.8980 chunk 83 optimal weight: 1.9990 chunk 112 optimal weight: 1.9990 overall best weight: 0.9982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Sorry: Reduce crashed with command 'molprobity.reduce -quiet -trim -'. Dumping stdin to file 'reduce_failure.pdb'. Return code: -15 Dumping stderr: