Starting phenix.real_space_refine (version: dev) on Sun Dec 18 12:31:59 2022 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ch7_30370/12_2022/7ch7_30370.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ch7_30370/12_2022/7ch7_30370.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.9 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ch7_30370/12_2022/7ch7_30370.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ch7_30370/12_2022/7ch7_30370.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ch7_30370/12_2022/7ch7_30370.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ch7_30370/12_2022/7ch7_30370.pdb" } resolution = 3.9 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.002 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Convert atoms to be neutral Process input model Symmetric amino acids flipped Residue "A ARG 22": "NH1" <-> "NH2" Residue "A ARG 97": "NH1" <-> "NH2" Residue "A ARG 111": "NH1" <-> "NH2" Residue "A ARG 147": "NH1" <-> "NH2" Residue "A PHE 184": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 237": "NH1" <-> "NH2" Residue "A PHE 258": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 22": "NH1" <-> "NH2" Residue "B ARG 97": "NH1" <-> "NH2" Residue "B ARG 111": "NH1" <-> "NH2" Residue "B ARG 147": "NH1" <-> "NH2" Residue "B PHE 184": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 237": "NH1" <-> "NH2" Residue "B PHE 258": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C ARG 18": "NH1" <-> "NH2" Residue "C ARG 21": "NH1" <-> "NH2" Residue "C ASP 25": "OD1" <-> "OD2" Residue "C ARG 33": "NH1" <-> "NH2" Residue "C PHE 67": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C ARG 77": "NH1" <-> "NH2" Residue "C ARG 84": "NH1" <-> "NH2" Residue "C ARG 86": "NH1" <-> "NH2" Residue "C ASP 220": "OD1" <-> "OD2" Residue "C ARG 241": "NH1" <-> "NH2" Residue "C ARG 255": "NH1" <-> "NH2" Residue "C TYR 256": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C ASP 260": "OD1" <-> "OD2" Residue "D ARG 18": "NH1" <-> "NH2" Residue "D ARG 21": "NH1" <-> "NH2" Residue "D ASP 25": "OD1" <-> "OD2" Residue "D ARG 33": "NH1" <-> "NH2" Residue "D PHE 67": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D ARG 77": "NH1" <-> "NH2" Residue "D ARG 84": "NH1" <-> "NH2" Residue "D ARG 86": "NH1" <-> "NH2" Residue "D ASP 220": "OD1" <-> "OD2" Residue "D ARG 241": "NH1" <-> "NH2" Residue "D ARG 255": "NH1" <-> "NH2" Residue "D ASP 260": "OD1" <-> "OD2" Residue "E ARG 32": "NH1" <-> "NH2" Residue "E ARG 49": "NH1" <-> "NH2" Residue "E TYR 82": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "F ARG 32": "NH1" <-> "NH2" Residue "F ARG 49": "NH1" <-> "NH2" Residue "F TYR 82": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.05s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/chem_data/mon_lib" Total number of atoms: 9198 Number of models: 1 Model: "" Number of chains: 6 Chain: "A" Number of atoms: 1871 Number of conformers: 1 Conformer: "" Number of residues, atoms: 248, 1871 Classifications: {'peptide': 248} Link IDs: {'PTRANS': 7, 'TRANS': 240} Chain: "B" Number of atoms: 1871 Number of conformers: 1 Conformer: "" Number of residues, atoms: 248, 1871 Classifications: {'peptide': 248} Link IDs: {'PTRANS': 7, 'TRANS': 240} Chain: "C" Number of atoms: 2006 Number of conformers: 1 Conformer: "" Number of residues, atoms: 265, 2006 Classifications: {'peptide': 265} Link IDs: {'PTRANS': 18, 'TRANS': 246} Chain: "D" Number of atoms: 2006 Number of conformers: 1 Conformer: "" Number of residues, atoms: 265, 2006 Classifications: {'peptide': 265} Link IDs: {'PTRANS': 18, 'TRANS': 246} Chain: "E" Number of atoms: 722 Number of conformers: 1 Conformer: "" Number of residues, atoms: 94, 722 Classifications: {'peptide': 94} Link IDs: {'PTRANS': 3, 'TRANS': 90} Chain: "F" Number of atoms: 722 Number of conformers: 1 Conformer: "" Number of residues, atoms: 94, 722 Classifications: {'peptide': 94} Link IDs: {'PTRANS': 3, 'TRANS': 90} Time building chain proxies: 5.76, per 1000 atoms: 0.63 Number of scatterers: 9198 At special positions: 0 Unit cell: (99, 111, 110, 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 58 16.00 O 1654 8.00 N 1572 7.00 C 5914 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 3.87 Conformation dependent library (CDL) restraints added in 1.5 seconds 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2212 Finding SS restraints... Secondary structure from input PDB file: 50 helices and 6 sheets defined 57.8% alpha, 7.2% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.44 Creating SS restraints... Processing helix chain 'A' and resid 17 through 29 removed outlier: 3.787A pdb=" N LEU A 27 " --> pdb=" O ALA A 23 " (cutoff:3.500A) Processing helix chain 'A' and resid 39 through 52 removed outlier: 4.387A pdb=" N LEU A 43 " --> pdb=" O LYS A 39 " (cutoff:3.500A) Processing helix chain 'A' and resid 55 through 79 removed outlier: 3.552A pdb=" N ILE A 59 " --> pdb=" O SER A 55 " (cutoff:3.500A) removed outlier: 4.240A pdb=" N SER A 62 " --> pdb=" O ILE A 58 " (cutoff:3.500A) removed outlier: 3.785A pdb=" N ILE A 66 " --> pdb=" O SER A 62 " (cutoff:3.500A) removed outlier: 3.621A pdb=" N VAL A 77 " --> pdb=" O GLN A 73 " (cutoff:3.500A) Processing helix chain 'A' and resid 86 through 98 removed outlier: 4.154A pdb=" N LEU A 93 " --> pdb=" O MET A 89 " (cutoff:3.500A) removed outlier: 3.630A pdb=" N SER A 94 " --> pdb=" O LEU A 90 " (cutoff:3.500A) Processing helix chain 'A' and resid 99 through 112 removed outlier: 4.169A pdb=" N ALA A 105 " --> pdb=" O PRO A 101 " (cutoff:3.500A) Processing helix chain 'A' and resid 116 through 127 removed outlier: 4.560A pdb=" N LEU A 122 " --> pdb=" O ALA A 118 " (cutoff:3.500A) Processing helix chain 'A' and resid 127 through 135 Processing helix chain 'A' and resid 138 through 175 removed outlier: 4.742A pdb=" N SER A 145 " --> pdb=" O ARG A 141 " (cutoff:3.500A) Proline residue: A 146 - end of helix removed outlier: 3.541A pdb=" N ILE A 153 " --> pdb=" O TRP A 149 " (cutoff:3.500A) Proline residue: A 156 - end of helix Processing helix chain 'A' and resid 189 through 193 Processing helix chain 'A' and resid 195 through 224 removed outlier: 3.715A pdb=" N LEU A 199 " --> pdb=" O TRP A 195 " (cutoff:3.500A) removed outlier: 4.109A pdb=" N LEU A 218 " --> pdb=" O THR A 214 " (cutoff:3.500A) removed outlier: 3.577A pdb=" N ASN A 220 " --> pdb=" O ILE A 216 " (cutoff:3.500A) removed outlier: 4.707A pdb=" N GLY A 221 " --> pdb=" O SER A 217 " (cutoff:3.500A) removed outlier: 3.728A pdb=" N TYR A 222 " --> pdb=" O LEU A 218 " (cutoff:3.500A) Processing helix chain 'A' and resid 230 through 259 removed outlier: 3.917A pdb=" N SER A 243 " --> pdb=" O VAL A 239 " (cutoff:3.500A) removed outlier: 4.724A pdb=" N VAL A 246 " --> pdb=" O SER A 242 " (cutoff:3.500A) Processing helix chain 'B' and resid 17 through 29 removed outlier: 3.787A pdb=" N LEU B 27 " --> pdb=" O ALA B 23 " (cutoff:3.500A) Processing helix chain 'B' and resid 39 through 52 removed outlier: 4.388A pdb=" N LEU B 43 " --> pdb=" O LYS B 39 " (cutoff:3.500A) Processing helix chain 'B' and resid 55 through 79 removed outlier: 3.551A pdb=" N ILE B 59 " --> pdb=" O SER B 55 " (cutoff:3.500A) removed outlier: 4.240A pdb=" N SER B 62 " --> pdb=" O ILE B 58 " (cutoff:3.500A) removed outlier: 3.785A pdb=" N ILE B 66 " --> pdb=" O SER B 62 " (cutoff:3.500A) removed outlier: 3.622A pdb=" N VAL B 77 " --> pdb=" O GLN B 73 " (cutoff:3.500A) Processing helix chain 'B' and resid 86 through 98 removed outlier: 4.152A pdb=" N LEU B 93 " --> pdb=" O MET B 89 " (cutoff:3.500A) removed outlier: 3.630A pdb=" N SER B 94 " --> pdb=" O LEU B 90 " (cutoff:3.500A) Processing helix chain 'B' and resid 99 through 112 removed outlier: 4.170A pdb=" N ALA B 105 " --> pdb=" O PRO B 101 " (cutoff:3.500A) Processing helix chain 'B' and resid 116 through 127 removed outlier: 4.561A pdb=" N LEU B 122 " --> pdb=" O ALA B 118 " (cutoff:3.500A) Processing helix chain 'B' and resid 127 through 135 Processing helix chain 'B' and resid 138 through 175 removed outlier: 4.742A pdb=" N SER B 145 " --> pdb=" O ARG B 141 " (cutoff:3.500A) Proline residue: B 146 - end of helix removed outlier: 3.541A pdb=" N ILE B 153 " --> pdb=" O TRP B 149 " (cutoff:3.500A) Proline residue: B 156 - end of helix Processing helix chain 'B' and resid 189 through 193 Processing helix chain 'B' and resid 195 through 224 removed outlier: 3.714A pdb=" N LEU B 199 " --> pdb=" O TRP B 195 " (cutoff:3.500A) removed outlier: 4.109A pdb=" N LEU B 218 " --> pdb=" O THR B 214 " (cutoff:3.500A) removed outlier: 3.577A pdb=" N ASN B 220 " --> pdb=" O ILE B 216 " (cutoff:3.500A) removed outlier: 4.708A pdb=" N GLY B 221 " --> pdb=" O SER B 217 " (cutoff:3.500A) removed outlier: 3.729A pdb=" N TYR B 222 " --> pdb=" O LEU B 218 " (cutoff:3.500A) Processing helix chain 'B' and resid 230 through 259 removed outlier: 3.916A pdb=" N SER B 243 " --> pdb=" O VAL B 239 " (cutoff:3.500A) removed outlier: 4.724A pdb=" N VAL B 246 " --> pdb=" O SER B 242 " (cutoff:3.500A) Processing helix chain 'C' and resid 46 through 55 Processing helix chain 'C' and resid 78 through 85 Processing helix chain 'C' and resid 101 through 114 Proline residue: C 109 - end of helix Processing helix chain 'C' and resid 117 through 133 removed outlier: 3.710A pdb=" N GLU C 130 " --> pdb=" O MET C 126 " (cutoff:3.500A) Processing helix chain 'C' and resid 134 through 138 Processing helix chain 'C' and resid 146 through 160 removed outlier: 3.726A pdb=" N ARG C 151 " --> pdb=" O GLY C 147 " (cutoff:3.500A) removed outlier: 3.686A pdb=" N ALA C 160 " --> pdb=" O ALA C 156 " (cutoff:3.500A) Processing helix chain 'C' and resid 176 through 192 removed outlier: 3.725A pdb=" N SER C 188 " --> pdb=" O VAL C 184 " (cutoff:3.500A) removed outlier: 3.895A pdb=" N ASN C 191 " --> pdb=" O ILE C 187 " (cutoff:3.500A) Processing helix chain 'C' and resid 204 through 209 Processing helix chain 'C' and resid 228 through 234 Processing helix chain 'C' and resid 237 through 246 removed outlier: 3.683A pdb=" N PHE C 243 " --> pdb=" O ARG C 239 " (cutoff:3.500A) removed outlier: 3.862A pdb=" N ASP C 245 " --> pdb=" O ARG C 241 " (cutoff:3.500A) Processing helix chain 'D' and resid 46 through 55 Processing helix chain 'D' and resid 78 through 85 Processing helix chain 'D' and resid 101 through 114 Proline residue: D 109 - end of helix Processing helix chain 'D' and resid 117 through 133 removed outlier: 3.709A pdb=" N GLU D 130 " --> pdb=" O MET D 126 " (cutoff:3.500A) Processing helix chain 'D' and resid 134 through 138 Processing helix chain 'D' and resid 146 through 160 removed outlier: 3.726A pdb=" N ARG D 151 " --> pdb=" O GLY D 147 " (cutoff:3.500A) removed outlier: 3.686A pdb=" N ALA D 160 " --> pdb=" O ALA D 156 " (cutoff:3.500A) Processing helix chain 'D' and resid 176 through 192 removed outlier: 3.725A pdb=" N SER D 188 " --> pdb=" O VAL D 184 " (cutoff:3.500A) removed outlier: 3.896A pdb=" N ASN D 191 " --> pdb=" O ILE D 187 " (cutoff:3.500A) Processing helix chain 'D' and resid 204 through 209 Processing helix chain 'D' and resid 228 through 234 Processing helix chain 'D' and resid 237 through 246 removed outlier: 3.682A pdb=" N PHE D 243 " --> pdb=" O ARG D 239 " (cutoff:3.500A) removed outlier: 3.862A pdb=" N ASP D 245 " --> pdb=" O ARG D 241 " (cutoff:3.500A) Processing helix chain 'E' and resid 25 through 30 removed outlier: 3.696A pdb=" N GLU E 30 " --> pdb=" O LEU E 26 " (cutoff:3.500A) Processing helix chain 'E' and resid 31 through 36 Processing helix chain 'E' and resid 51 through 66 Processing helix chain 'E' and resid 78 through 88 Processing helix chain 'F' and resid 25 through 30 removed outlier: 3.697A pdb=" N GLU F 30 " --> pdb=" O LEU F 26 " (cutoff:3.500A) Processing helix chain 'F' and resid 31 through 36 Processing helix chain 'F' and resid 51 through 66 Processing helix chain 'F' and resid 78 through 88 Processing sheet with id=AA1, first strand: chain 'C' and resid 27 through 32 removed outlier: 4.512A pdb=" N LEU C 29 " --> pdb=" O ARG C 12 " (cutoff:3.500A) removed outlier: 4.580A pdb=" N ARG C 12 " --> pdb=" O LEU C 29 " (cutoff:3.500A) removed outlier: 7.260A pdb=" N VAL C 31 " --> pdb=" O ASP C 10 " (cutoff:3.500A) removed outlier: 5.272A pdb=" N ASP C 10 " --> pdb=" O VAL C 31 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'C' and resid 90 through 91 removed outlier: 5.894A pdb=" N PHE C 168 " --> pdb=" O VAL C 201 " (cutoff:3.500A) removed outlier: 6.752A pdb=" N THR C 37 " --> pdb=" O VAL C 200 " (cutoff:3.500A) removed outlier: 7.776A pdb=" N SER C 202 " --> pdb=" O THR C 37 " (cutoff:3.500A) removed outlier: 5.714A pdb=" N ILE C 39 " --> pdb=" O SER C 202 " (cutoff:3.500A) removed outlier: 5.612A pdb=" N ALA C 219 " --> pdb=" O ILE C 223 " (cutoff:3.500A) removed outlier: 5.091A pdb=" N ILE C 223 " --> pdb=" O ALA C 219 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'D' and resid 27 through 32 removed outlier: 4.511A pdb=" N LEU D 29 " --> pdb=" O ARG D 12 " (cutoff:3.500A) removed outlier: 4.580A pdb=" N ARG D 12 " --> pdb=" O LEU D 29 " (cutoff:3.500A) removed outlier: 7.260A pdb=" N VAL D 31 " --> pdb=" O ASP D 10 " (cutoff:3.500A) removed outlier: 5.271A pdb=" N ASP D 10 " --> pdb=" O VAL D 31 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'D' and resid 90 through 91 removed outlier: 5.894A pdb=" N PHE D 168 " --> pdb=" O VAL D 201 " (cutoff:3.500A) removed outlier: 6.752A pdb=" N THR D 37 " --> pdb=" O VAL D 200 " (cutoff:3.500A) removed outlier: 7.776A pdb=" N SER D 202 " --> pdb=" O THR D 37 " (cutoff:3.500A) removed outlier: 5.714A pdb=" N ILE D 39 " --> pdb=" O SER D 202 " (cutoff:3.500A) removed outlier: 5.612A pdb=" N ALA D 219 " --> pdb=" O ILE D 223 " (cutoff:3.500A) removed outlier: 5.091A pdb=" N ILE D 223 " --> pdb=" O ALA D 219 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'E' and resid 5 through 9 removed outlier: 3.893A pdb=" N LEU E 16 " --> pdb=" O ASP E 43 " (cutoff:3.500A) removed outlier: 6.600A pdb=" N ILE E 42 " --> pdb=" O GLN E 75 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'F' and resid 5 through 9 removed outlier: 3.893A pdb=" N LEU F 16 " --> pdb=" O ASP F 43 " (cutoff:3.500A) removed outlier: 6.600A pdb=" N ILE F 42 " --> pdb=" O GLN F 75 " (cutoff:3.500A) 472 hydrogen bonds defined for protein. 1374 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: 2.59 Time building geometry restraints manager: 4.19 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.34: 2996 1.34 - 1.46: 1730 1.46 - 1.57: 4536 1.57 - 1.69: 0 1.69 - 1.81: 108 Bond restraints: 9370 Sorted by residual: bond pdb=" CA GLU C 162 " pdb=" C GLU C 162 " ideal model delta sigma weight residual 1.531 1.517 0.014 7.40e-03 1.83e+04 3.68e+00 bond pdb=" CA GLU D 162 " pdb=" C GLU D 162 " ideal model delta sigma weight residual 1.531 1.518 0.013 7.40e-03 1.83e+04 3.06e+00 bond pdb=" CA GLU D 170 " pdb=" C GLU D 170 " ideal model delta sigma weight residual 1.531 1.520 0.012 7.40e-03 1.83e+04 2.43e+00 bond pdb=" CA SER B 114 " pdb=" C SER B 114 " ideal model delta sigma weight residual 1.525 1.510 0.015 1.02e-02 9.61e+03 2.27e+00 bond pdb=" CA GLU C 170 " pdb=" C GLU C 170 " ideal model delta sigma weight residual 1.531 1.520 0.011 7.40e-03 1.83e+04 2.24e+00 ... (remaining 9365 not shown) Histogram of bond angle deviations from ideal: 99.06 - 106.04: 258 106.04 - 113.02: 5394 113.02 - 120.01: 3230 120.01 - 126.99: 3731 126.99 - 133.98: 117 Bond angle restraints: 12730 Sorted by residual: angle pdb=" C PRO C 251 " pdb=" N VAL C 252 " pdb=" CA VAL C 252 " ideal model delta sigma weight residual 122.59 119.80 2.79 7.20e-01 1.93e+00 1.50e+01 angle pdb=" C PRO D 251 " pdb=" N VAL D 252 " pdb=" CA VAL D 252 " ideal model delta sigma weight residual 122.59 119.80 2.79 7.20e-01 1.93e+00 1.50e+01 angle pdb=" N ILE C 211 " pdb=" CA ILE C 211 " pdb=" C ILE C 211 " ideal model delta sigma weight residual 113.71 110.19 3.52 9.50e-01 1.11e+00 1.37e+01 angle pdb=" N GLY A 113 " pdb=" CA GLY A 113 " pdb=" C GLY A 113 " ideal model delta sigma weight residual 114.66 110.08 4.58 1.24e+00 6.50e-01 1.36e+01 angle pdb=" N ILE D 211 " pdb=" CA ILE D 211 " pdb=" C ILE D 211 " ideal model delta sigma weight residual 113.71 110.22 3.49 9.50e-01 1.11e+00 1.35e+01 ... (remaining 12725 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 15.81: 4813 15.81 - 31.61: 534 31.61 - 47.42: 174 47.42 - 63.23: 79 63.23 - 79.04: 4 Dihedral angle restraints: 5604 sinusoidal: 2164 harmonic: 3440 Sorted by residual: dihedral pdb=" CA LEU F 95 " pdb=" C LEU F 95 " pdb=" N PRO F 96 " pdb=" CA PRO F 96 " ideal model delta harmonic sigma weight residual -180.00 -150.99 -29.01 0 5.00e+00 4.00e-02 3.37e+01 dihedral pdb=" CA LEU E 95 " pdb=" C LEU E 95 " pdb=" N PRO E 96 " pdb=" CA PRO E 96 " ideal model delta harmonic sigma weight residual -180.00 -151.00 -29.00 0 5.00e+00 4.00e-02 3.36e+01 dihedral pdb=" CA GLU D 162 " pdb=" C GLU D 162 " pdb=" N PRO D 163 " pdb=" CA PRO D 163 " ideal model delta harmonic sigma weight residual -180.00 -159.17 -20.83 0 5.00e+00 4.00e-02 1.74e+01 ... (remaining 5601 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.030: 717 0.030 - 0.060: 536 0.060 - 0.089: 179 0.089 - 0.119: 80 0.119 - 0.149: 16 Chirality restraints: 1528 Sorted by residual: chirality pdb=" CA PRO D 163 " pdb=" N PRO D 163 " pdb=" C PRO D 163 " pdb=" CB PRO D 163 " both_signs ideal model delta sigma weight residual False 2.72 2.57 0.15 2.00e-01 2.50e+01 5.56e-01 chirality pdb=" CA PRO C 163 " pdb=" N PRO C 163 " pdb=" C PRO C 163 " pdb=" CB PRO C 163 " both_signs ideal model delta sigma weight residual False 2.72 2.57 0.15 2.00e-01 2.50e+01 5.37e-01 chirality pdb=" CA PRO E 96 " pdb=" N PRO E 96 " pdb=" C PRO E 96 " pdb=" CB PRO E 96 " both_signs ideal model delta sigma weight residual False 2.72 2.58 0.14 2.00e-01 2.50e+01 4.89e-01 ... (remaining 1525 not shown) Planarity restraints: 1604 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C LEU F 95 " 0.074 5.00e-02 4.00e+02 1.12e-01 1.99e+01 pdb=" N PRO F 96 " -0.193 5.00e-02 4.00e+02 pdb=" CA PRO F 96 " 0.061 5.00e-02 4.00e+02 pdb=" CD PRO F 96 " 0.058 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C LEU E 95 " -0.074 5.00e-02 4.00e+02 1.11e-01 1.99e+01 pdb=" N PRO E 96 " 0.193 5.00e-02 4.00e+02 pdb=" CA PRO E 96 " -0.061 5.00e-02 4.00e+02 pdb=" CD PRO E 96 " -0.058 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CA LEU B 140 " -0.007 2.00e-02 2.50e+03 1.42e-02 2.03e+00 pdb=" C LEU B 140 " 0.025 2.00e-02 2.50e+03 pdb=" O LEU B 140 " -0.009 2.00e-02 2.50e+03 pdb=" N ARG B 141 " -0.008 2.00e-02 2.50e+03 ... (remaining 1601 not shown) Histogram of nonbonded interaction distances: 2.10 - 2.66: 301 2.66 - 3.22: 10345 3.22 - 3.78: 15437 3.78 - 4.34: 20279 4.34 - 4.90: 32211 Nonbonded interactions: 78573 Sorted by model distance: nonbonded pdb=" O TYR B 75 " pdb=" OG1 THR B 79 " model vdw 2.096 2.440 nonbonded pdb=" O TYR A 75 " pdb=" OG1 THR A 79 " model vdw 2.096 2.440 nonbonded pdb=" OG SER A 228 " pdb=" OG SER B 228 " model vdw 2.142 2.440 nonbonded pdb=" O VAL A 51 " pdb=" OG SER A 55 " model vdw 2.170 2.440 nonbonded pdb=" O TRP B 168 " pdb=" OG SER B 171 " model vdw 2.170 2.440 ... (remaining 78568 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = chain 'A' selection = chain 'B' } ncs_group { reference = chain 'C' selection = chain 'D' } ncs_group { reference = chain 'E' selection = chain 'F' } 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 58 5.16 5 C 5914 2.51 5 N 1572 2.21 5 O 1654 1.98 5 sf(0) = scattering factor at diffraction angle 0. Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.920 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.010 Extract box with map and model: 4.390 Check model and map are aligned: 0.130 Convert atoms to be neutral: 0.080 Process input model: 27.780 Find NCS groups from input model: 0.390 Set up NCS constraints: 0.040 Set refine NCS operators: 0.000 Set scattering table: 0.010 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:2.160 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 35.910 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8114 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.007 0.046 9370 Z= 0.491 Angle : 0.939 13.756 12730 Z= 0.514 Chirality : 0.047 0.149 1528 Planarity : 0.006 0.112 1604 Dihedral : 16.386 79.037 3392 Min Nonbonded Distance : 2.096 Molprobity Statistics. All-atom Clashscore : 26.29 Ramachandran Plot: Outliers : 0.00 % Allowed : 13.56 % Favored : 86.44 % Rotamer Outliers : 18.11 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.81 (0.20), residues: 1202 helix: -3.04 (0.15), residues: 646 sheet: -2.97 (0.57), residues: 66 loop : -3.26 (0.26), residues: 490 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 482 residues out of total 994 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 180 poor density : 302 time to evaluate : 0.809 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 revert: symmetry clash outliers start: 180 outliers final: 83 residues processed: 420 average time/residue: 0.1427 time to fit residues: 83.6409 Evaluate side-chains 308 residues out of total 994 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 83 poor density : 225 time to evaluate : 0.717 Switching outliers to nearest non-outliers outliers start: 83 outliers final: 0 residues processed: 83 average time/residue: 0.0643 time to fit residues: 10.0979 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 99 optimal weight: 3.9990 chunk 89 optimal weight: 0.9990 chunk 49 optimal weight: 5.9990 chunk 30 optimal weight: 2.9990 chunk 60 optimal weight: 1.9990 chunk 47 optimal weight: 0.9980 chunk 92 optimal weight: 1.9990 chunk 35 optimal weight: 1.9990 chunk 56 optimal weight: 2.9990 chunk 68 optimal weight: 2.9990 chunk 107 optimal weight: 8.9990 overall best weight: 1.5988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 128 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 190 GLN A 220 ASN ** B 73 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 128 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 190 GLN ** C 101 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 105 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** C 175 GLN ** C 226 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 101 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 175 GLN ** D 226 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** E 60 HIS F 60 HIS Total number of N/Q/H flips: 7 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8013 moved from start: 0.2771 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.051 9370 Z= 0.293 Angle : 0.830 10.514 12730 Z= 0.394 Chirality : 0.047 0.257 1528 Planarity : 0.006 0.102 1604 Dihedral : 5.515 28.892 1294 Min Nonbonded Distance : 2.026 Molprobity Statistics. All-atom Clashscore : 23.40 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.90 % Favored : 92.10 % Rotamer Outliers : 6.34 % 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: -3.06 (0.22), residues: 1202 helix: -1.46 (0.19), residues: 636 sheet: -1.67 (0.60), residues: 62 loop : -2.78 (0.25), residues: 504 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 344 residues out of total 994 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 63 poor density : 281 time to evaluate : 1.121 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 revert: symmetry clash outliers start: 63 outliers final: 27 residues processed: 320 average time/residue: 0.1815 time to fit residues: 82.5816 Evaluate side-chains 252 residues out of total 994 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 27 poor density : 225 time to evaluate : 1.100 Switching outliers to nearest non-outliers outliers start: 27 outliers final: 0 residues processed: 27 average time/residue: 0.0847 time to fit residues: 5.7170 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 59 optimal weight: 0.8980 chunk 33 optimal weight: 0.7980 chunk 89 optimal weight: 4.9990 chunk 72 optimal weight: 0.0570 chunk 29 optimal weight: 1.9990 chunk 107 optimal weight: 8.9990 chunk 115 optimal weight: 0.2980 chunk 95 optimal weight: 7.9990 chunk 106 optimal weight: 5.9990 chunk 36 optimal weight: 0.9990 chunk 86 optimal weight: 0.6980 overall best weight: 0.5498 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 128 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 220 ASN ** B 128 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 220 ASN ** C 105 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** C 175 GLN D 226 HIS Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7961 moved from start: 0.3502 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.046 9370 Z= 0.228 Angle : 0.769 9.120 12730 Z= 0.363 Chirality : 0.048 0.350 1528 Planarity : 0.005 0.096 1604 Dihedral : 5.241 28.007 1294 Min Nonbonded Distance : 2.084 Molprobity Statistics. All-atom Clashscore : 19.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.15 % Favored : 92.85 % Rotamer Outliers : 4.73 % 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.40 (0.23), residues: 1202 helix: -0.87 (0.19), residues: 650 sheet: -2.07 (0.53), residues: 80 loop : -2.44 (0.27), residues: 472 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 306 residues out of total 994 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 47 poor density : 259 time to evaluate : 1.003 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 revert: symmetry clash outliers start: 47 outliers final: 20 residues processed: 289 average time/residue: 0.1867 time to fit residues: 76.0818 Evaluate side-chains 249 residues out of total 994 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 20 poor density : 229 time to evaluate : 1.109 Switching outliers to nearest non-outliers outliers start: 20 outliers final: 0 residues processed: 20 average time/residue: 0.0862 time to fit residues: 4.6268 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 106 optimal weight: 7.9990 chunk 80 optimal weight: 0.9980 chunk 55 optimal weight: 6.9990 chunk 11 optimal weight: 7.9990 chunk 51 optimal weight: 3.9990 chunk 72 optimal weight: 2.9990 chunk 107 optimal weight: 6.9990 chunk 114 optimal weight: 2.9990 chunk 56 optimal weight: 3.9990 chunk 102 optimal weight: 0.9990 chunk 30 optimal weight: 6.9990 overall best weight: 2.3988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 128 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 220 ASN ** B 128 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 105 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8023 moved from start: 0.3656 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.044 9370 Z= 0.309 Angle : 0.779 11.150 12730 Z= 0.371 Chirality : 0.048 0.344 1528 Planarity : 0.006 0.097 1604 Dihedral : 5.169 27.037 1294 Min Nonbonded Distance : 2.085 Molprobity Statistics. All-atom Clashscore : 22.07 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.65 % Favored : 92.35 % Rotamer Outliers : 5.53 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.26 (0.23), residues: 1202 helix: -0.70 (0.20), residues: 652 sheet: -2.51 (0.47), residues: 90 loop : -2.35 (0.28), residues: 460 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 291 residues out of total 994 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 55 poor density : 236 time to evaluate : 1.129 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 revert: symmetry clash outliers start: 55 outliers final: 30 residues processed: 267 average time/residue: 0.1917 time to fit residues: 71.6698 Evaluate side-chains 248 residues out of total 994 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 30 poor density : 218 time to evaluate : 1.116 Switching outliers to nearest non-outliers outliers start: 30 outliers final: 0 residues processed: 30 average time/residue: 0.0847 time to fit residues: 6.1194 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 94 optimal weight: 3.9990 chunk 64 optimal weight: 1.9990 chunk 1 optimal weight: 0.9990 chunk 84 optimal weight: 0.6980 chunk 47 optimal weight: 0.8980 chunk 97 optimal weight: 1.9990 chunk 78 optimal weight: 0.9990 chunk 0 optimal weight: 5.9990 chunk 58 optimal weight: 0.8980 chunk 102 optimal weight: 2.9990 chunk 28 optimal weight: 0.3980 overall best weight: 0.7782 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 128 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 220 ASN ** B 128 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 105 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 191 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7950 moved from start: 0.4159 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.046 9370 Z= 0.209 Angle : 0.730 9.097 12730 Z= 0.346 Chirality : 0.045 0.336 1528 Planarity : 0.005 0.090 1604 Dihedral : 5.007 26.831 1294 Min Nonbonded Distance : 2.094 Molprobity Statistics. All-atom Clashscore : 18.54 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.40 % Favored : 92.60 % Rotamer Outliers : 3.52 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.94 (0.23), residues: 1202 helix: -0.48 (0.19), residues: 664 sheet: -2.17 (0.49), residues: 90 loop : -2.24 (0.28), residues: 448 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2404 Ramachandran restraints generated. 1202 Oldfield, 0 Emsley, 1202 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 283 residues out of total 994 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 35 poor density : 248 time to evaluate : 1.087 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: 35 outliers final: 17 residues processed: 268 average time/residue: 0.1902 time to fit residues: 71.4364 Evaluate side-chains 247 residues out of total 994 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 17 poor density : 230 time to evaluate : 1.062 Switching outliers to nearest non-outliers outliers start: 17 outliers final: 0 residues processed: 17 average time/residue: 0.0861 time to fit residues: 4.1339 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 8, in run_program(real_space_refine.Program) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/iotbx/cli_parser.py", line 864, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/phenix/phenix/programs/real_space_refine.py", line 191, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 297, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 312, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 673, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 736, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 51, in __init__ gradients_method = gradients_method) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/mmtbx/refinement/real_space/weight.py", line 87, in __init__ rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 368, in refine rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 112, in refine states_collector = self.states_accumulator) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 174, in __init__ exception_handling_params=lbfgs_exception_handling_params) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 131, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 245, in compute_functional_and_gradients compute_gradients=True) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1479, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1381, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1262, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 61.3906 > 50: distance: 58 - 64: 14.255 distance: 64 - 65: 3.852 distance: 65 - 66: 3.172 distance: 66 - 67: 5.745 distance: 68 - 69: 10.520 distance: 69 - 70: 7.653 distance: 70 - 71: 13.532 distance: 70 - 72: 4.021 distance: 73 - 74: 15.290 distance: 74 - 75: 8.496 distance: 74 - 77: 16.557 distance: 75 - 76: 15.457 distance: 75 - 81: 14.678 distance: 77 - 78: 17.168 distance: 78 - 79: 6.173 distance: 81 - 82: 11.020 distance: 82 - 83: 15.794 distance: 82 - 85: 21.666 distance: 83 - 84: 14.086 distance: 83 - 88: 26.155 distance: 85 - 86: 12.795 distance: 85 - 87: 34.470 distance: 88 - 89: 10.891 distance: 89 - 90: 48.932 distance: 89 - 92: 28.158 distance: 90 - 91: 12.556 distance: 90 - 96: 17.681 distance: 91 - 119: 17.817 distance: 92 - 93: 13.592 distance: 93 - 94: 25.541 distance: 93 - 95: 32.194 distance: 96 - 97: 23.345 distance: 97 - 98: 50.533 distance: 97 - 100: 56.436 distance: 98 - 99: 37.414 distance: 98 - 104: 22.152 distance: 100 - 101: 14.006 distance: 101 - 102: 13.113 distance: 101 - 103: 15.197 distance: 104 - 105: 35.118 distance: 104 - 110: 22.166 distance: 105 - 106: 16.888 distance: 105 - 108: 3.242 distance: 106 - 111: 10.860 distance: 108 - 109: 21.919 distance: 109 - 110: 16.149 distance: 111 - 112: 24.130 distance: 112 - 113: 6.650 distance: 112 - 115: 19.185 distance: 113 - 114: 31.494 distance: 113 - 119: 21.770 distance: 115 - 116: 27.892 distance: 116 - 117: 12.932 distance: 116 - 118: 10.790 distance: 119 - 120: 19.771 distance: 120 - 121: 20.135 distance: 120 - 123: 15.034 distance: 121 - 122: 20.344 distance: 121 - 133: 26.897 distance: 123 - 124: 15.146 distance: 124 - 125: 8.712 distance: 124 - 126: 7.224 distance: 125 - 127: 7.647 distance: 126 - 128: 21.973 distance: 126 - 129: 13.443 distance: 127 - 128: 19.853 distance: 128 - 130: 8.856 distance: 129 - 131: 6.580 distance: 130 - 132: 11.270 distance: 131 - 132: 11.825 distance: 133 - 134: 40.074 distance: 134 - 135: 37.462 distance: 134 - 137: 17.855 distance: 135 - 142: 43.217 distance: 137 - 138: 17.138 distance: 138 - 139: 26.326 distance: 139 - 140: 12.927 distance: 139 - 141: 5.902 distance: 142 - 143: 8.850 distance: 143 - 144: 10.748 distance: 143 - 146: 13.962 distance: 144 - 145: 12.224 distance: 144 - 150: 22.673 distance: 145 - 179: 22.930 distance: 146 - 147: 18.486 distance: 147 - 148: 11.177 distance: 148 - 149: 37.453