Starting phenix.real_space_refine on Fri Dec 8 16:24:35 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wyb_32890/12_2023/7wyb_32890_neut.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wyb_32890/12_2023/7wyb_32890.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.97 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wyb_32890/12_2023/7wyb_32890.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wyb_32890/12_2023/7wyb_32890.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wyb_32890/12_2023/7wyb_32890_neut.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wyb_32890/12_2023/7wyb_32890_neut.pdb" } resolution = 2.97 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.001 sd= 0.019 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: 4 Type Number sf(0) Gaussians S 66 5.16 5 C 5587 2.51 5 N 1476 2.21 5 O 1628 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "E GLU 248": "OE1" <-> "OE2" Residue "R PHE 1027": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/chem_data/mon_lib" Total number of atoms: 8757 Number of models: 1 Model: "" Number of chains: 5 Chain: "B" Number of atoms: 1812 Number of conformers: 1 Conformer: "" Number of residues, atoms: 226, 1812 Classifications: {'peptide': 226} Link IDs: {'PTRANS': 2, 'TRANS': 223} Chain breaks: 1 Chain: "C" Number of atoms: 2607 Number of conformers: 1 Conformer: "" Number of residues, atoms: 339, 2607 Classifications: {'peptide': 339} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 5, 'TRANS': 333} Chain: "D" Number of atoms: 444 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 444 Classifications: {'peptide': 58} Link IDs: {'PTRANS': 4, 'TRANS': 53} Chain: "E" Number of atoms: 1771 Number of conformers: 1 Conformer: "" Number of residues, atoms: 230, 1771 Classifications: {'peptide': 230} Link IDs: {'PCIS': 1, 'PTRANS': 9, 'TRANS': 219} Chain breaks: 1 Chain: "R" Number of atoms: 2123 Number of conformers: 1 Conformer: "" Number of residues, atoms: 262, 2123 Classifications: {'peptide': 262} Link IDs: {'PTRANS': 3, 'TRANS': 258} Chain breaks: 1 Time building chain proxies: 4.78, per 1000 atoms: 0.55 Number of scatterers: 8757 At special positions: 0 Unit cell: (92.4, 122.1, 130.9, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 4 Type Number sf(0) S 66 16.00 O 1628 8.00 N 1476 7.00 C 5587 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=2, symmetry=0 Simple disulfide: pdb=" SG CYS E 161 " - pdb=" SG CYS E 231 " distance=2.04 Simple disulfide: pdb=" SG CYS R1008 " - pdb=" SG CYS R1080 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 3.26 Conformation dependent library (CDL) restraints added in 1.7 seconds 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2080 Finding SS restraints... Secondary structure from input PDB file: 32 helices and 12 sheets defined 35.3% alpha, 22.4% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.97 Creating SS restraints... Processing helix chain 'B' and resid 8 through 30 removed outlier: 3.551A pdb=" N ALA B 12 " --> pdb=" O GLU B 8 " (cutoff:3.500A) Processing helix chain 'B' and resid 45 through 54 Processing helix chain 'B' and resid 207 through 212 removed outlier: 4.407A pdb=" N TRP B 211 " --> pdb=" O GLU B 207 " (cutoff:3.500A) removed outlier: 3.928A pdb=" N ILE B 212 " --> pdb=" O ARG B 208 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 207 through 212' Processing helix chain 'B' and resid 213 through 216 removed outlier: 3.736A pdb=" N GLU B 216 " --> pdb=" O HIS B 213 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 213 through 216' Processing helix chain 'B' and resid 241 through 255 removed outlier: 3.842A pdb=" N ALA B 245 " --> pdb=" O ASN B 241 " (cutoff:3.500A) Processing helix chain 'B' and resid 256 through 260 removed outlier: 3.777A pdb=" N PHE B 259 " --> pdb=" O ASN B 256 " (cutoff:3.500A) removed outlier: 3.622A pdb=" N THR B 260 " --> pdb=" O LYS B 257 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 256 through 260' Processing helix chain 'B' and resid 271 through 279 Processing helix chain 'B' and resid 295 through 309 removed outlier: 3.654A pdb=" N ALA B 299 " --> pdb=" O THR B 295 " (cutoff:3.500A) Processing helix chain 'B' and resid 328 through 331 Processing helix chain 'B' and resid 332 through 348 Processing helix chain 'C' and resid 3 through 26 removed outlier: 3.774A pdb=" N ALA C 26 " --> pdb=" O ARG C 22 " (cutoff:3.500A) Processing helix chain 'C' and resid 29 through 33 Processing helix chain 'D' and resid 6 through 24 Processing helix chain 'D' and resid 29 through 43 Processing helix chain 'D' and resid 44 through 48 removed outlier: 3.521A pdb=" N GLU D 47 " --> pdb=" O HIS D 44 " (cutoff:3.500A) Processing helix chain 'E' and resid 28 through 32 removed outlier: 3.957A pdb=" N SER E 31 " --> pdb=" O ALA E 28 " (cutoff:3.500A) removed outlier: 3.618A pdb=" N PHE E 32 " --> pdb=" O PHE E 29 " (cutoff:3.500A) No H-bonds generated for 'chain 'E' and resid 28 through 32' Processing helix chain 'E' and resid 87 through 91 removed outlier: 4.034A pdb=" N THR E 91 " --> pdb=" O SER E 88 " (cutoff:3.500A) Processing helix chain 'R' and resid 924 through 929 Processing helix chain 'R' and resid 940 through 956 Processing helix chain 'R' and resid 958 through 970 removed outlier: 3.821A pdb=" N ILE R 962 " --> pdb=" O VAL R 958 " (cutoff:3.500A) removed outlier: 3.822A pdb=" N PHE R 970 " --> pdb=" O THR R 966 " (cutoff:3.500A) Processing helix chain 'R' and resid 971 through 974 removed outlier: 3.656A pdb=" N GLN R 974 " --> pdb=" O ARG R 971 " (cutoff:3.500A) No H-bonds generated for 'chain 'R' and resid 971 through 974' Processing helix chain 'R' and resid 975 through 999 removed outlier: 3.597A pdb=" N ILE R 999 " --> pdb=" O PHE R 995 " (cutoff:3.500A) Processing helix chain 'R' and resid 1007 through 1038 removed outlier: 3.930A pdb=" N PHE R1011 " --> pdb=" O ALA R1007 " (cutoff:3.500A) Processing helix chain 'R' and resid 1047 through 1070 removed outlier: 3.762A pdb=" N PHE R1051 " --> pdb=" O ARG R1047 " (cutoff:3.500A) Proline residue: R1059 - end of helix Processing helix chain 'R' and resid 1088 through 1091 Processing helix chain 'R' and resid 1092 through 1099 removed outlier: 3.658A pdb=" N ALA R1096 " --> pdb=" O PHE R1092 " (cutoff:3.500A) Processing helix chain 'R' and resid 1099 through 1117 removed outlier: 4.131A pdb=" N ILE R1105 " --> pdb=" O MET R1101 " (cutoff:3.500A) Processing helix chain 'R' and resid 1144 through 1154 removed outlier: 3.652A pdb=" N CYS R1152 " --> pdb=" O ILE R1148 " (cutoff:3.500A) removed outlier: 3.909A pdb=" N LEU R1153 " --> pdb=" O ALA R1149 " (cutoff:3.500A) removed outlier: 3.679A pdb=" N LEU R1154 " --> pdb=" O LEU R1150 " (cutoff:3.500A) Processing helix chain 'R' and resid 1155 through 1162 Processing helix chain 'R' and resid 1169 through 1181 removed outlier: 3.801A pdb=" N ALA R1173 " --> pdb=" O THR R1169 " (cutoff:3.500A) Processing helix chain 'R' and resid 1182 through 1193 removed outlier: 4.024A pdb=" N VAL R1193 " --> pdb=" O ILE R1189 " (cutoff:3.500A) Processing helix chain 'R' and resid 1195 through 1206 Processing sheet with id=AA1, first strand: chain 'B' and resid 185 through 191 removed outlier: 6.355A pdb=" N VAL B 34 " --> pdb=" O LYS B 197 " (cutoff:3.500A) removed outlier: 7.586A pdb=" N PHE B 199 " --> pdb=" O VAL B 34 " (cutoff:3.500A) removed outlier: 7.002A pdb=" N LEU B 36 " --> pdb=" O PHE B 199 " (cutoff:3.500A) removed outlier: 4.143A pdb=" N ALA B 220 " --> pdb=" O LYS B 35 " (cutoff:3.500A) removed outlier: 7.146A pdb=" N ILE B 264 " --> pdb=" O TYR B 320 " (cutoff:3.500A) removed outlier: 8.483A pdb=" N HIS B 322 " --> pdb=" O ILE B 264 " (cutoff:3.500A) removed outlier: 6.839A pdb=" N LEU B 266 " --> pdb=" O HIS B 322 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'C' and resid 47 through 51 removed outlier: 7.155A pdb=" N ILE C 338 " --> pdb=" O ARG C 48 " (cutoff:3.500A) removed outlier: 4.827A pdb=" N THR C 50 " --> pdb=" O LEU C 336 " (cutoff:3.500A) removed outlier: 6.670A pdb=" N LEU C 336 " --> pdb=" O THR C 50 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'C' and resid 58 through 63 removed outlier: 6.808A pdb=" N ALA C 73 " --> pdb=" O TYR C 59 " (cutoff:3.500A) removed outlier: 4.294A pdb=" N MET C 61 " --> pdb=" O VAL C 71 " (cutoff:3.500A) removed outlier: 6.817A pdb=" N VAL C 71 " --> pdb=" O MET C 61 " (cutoff:3.500A) removed outlier: 5.261A pdb=" N TRP C 63 " --> pdb=" O LEU C 69 " (cutoff:3.500A) removed outlier: 7.280A pdb=" N LEU C 69 " --> pdb=" O TRP C 63 " (cutoff:3.500A) removed outlier: 3.871A pdb=" N LYS C 78 " --> pdb=" O SER C 74 " (cutoff:3.500A) removed outlier: 6.770A pdb=" N LEU C 79 " --> pdb=" O ALA C 92 " (cutoff:3.500A) removed outlier: 4.805A pdb=" N ALA C 92 " --> pdb=" O LEU C 79 " (cutoff:3.500A) removed outlier: 6.639A pdb=" N ILE C 81 " --> pdb=" O VAL C 90 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'C' and resid 100 through 105 removed outlier: 7.284A pdb=" N GLY C 115 " --> pdb=" O MET C 101 " (cutoff:3.500A) removed outlier: 4.785A pdb=" N CYS C 103 " --> pdb=" O ALA C 113 " (cutoff:3.500A) removed outlier: 6.465A pdb=" N ALA C 113 " --> pdb=" O CYS C 103 " (cutoff:3.500A) removed outlier: 4.546A pdb=" N TYR C 105 " --> pdb=" O TYR C 111 " (cutoff:3.500A) removed outlier: 6.976A pdb=" N TYR C 111 " --> pdb=" O TYR C 105 " (cutoff:3.500A) removed outlier: 3.709A pdb=" N ARG C 137 " --> pdb=" O ILE C 123 " (cutoff:3.500A) removed outlier: 6.582A pdb=" N ASN C 125 " --> pdb=" O VAL C 135 " (cutoff:3.500A) removed outlier: 6.155A pdb=" N VAL C 135 " --> pdb=" O ASN C 125 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'C' and resid 146 through 153 removed outlier: 3.807A pdb=" N CYS C 148 " --> pdb=" O SER C 160 " (cutoff:3.500A) removed outlier: 6.573A pdb=" N GLN C 156 " --> pdb=" O LEU C 152 " (cutoff:3.500A) removed outlier: 3.741A pdb=" N SER C 161 " --> pdb=" O THR C 165 " (cutoff:3.500A) removed outlier: 5.121A pdb=" N THR C 165 " --> pdb=" O SER C 161 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'C' and resid 187 through 191 removed outlier: 6.829A pdb=" N GLY C 202 " --> pdb=" O MET C 188 " (cutoff:3.500A) removed outlier: 4.325A pdb=" N LEU C 190 " --> pdb=" O VAL C 200 " (cutoff:3.500A) removed outlier: 6.642A pdb=" N VAL C 200 " --> pdb=" O LEU C 190 " (cutoff:3.500A) removed outlier: 3.983A pdb=" N SER C 207 " --> pdb=" O ALA C 203 " (cutoff:3.500A) removed outlier: 6.628A pdb=" N ASP C 212 " --> pdb=" O CYS C 218 " (cutoff:3.500A) removed outlier: 6.060A pdb=" N CYS C 218 " --> pdb=" O ASP C 212 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'C' and resid 229 through 234 Processing sheet with id=AA8, first strand: chain 'C' and resid 273 through 277 removed outlier: 3.691A pdb=" N SER C 275 " --> pdb=" O GLY C 288 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'E' and resid 3 through 7 removed outlier: 6.545A pdb=" N GLY E 16 " --> pdb=" O SER E 85 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'E' and resid 10 through 12 removed outlier: 5.739A pdb=" N ARG E 38 " --> pdb=" O TRP E 47 " (cutoff:3.500A) removed outlier: 5.335A pdb=" N TRP E 47 " --> pdb=" O ARG E 38 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'E' and resid 142 through 144 Processing sheet with id=AB3, first strand: chain 'E' and resid 149 through 150 removed outlier: 3.898A pdb=" N GLU E 248 " --> pdb=" O VAL E 149 " (cutoff:3.500A) removed outlier: 6.126A pdb=" N LEU E 180 " --> pdb=" O LEU E 189 " (cutoff:3.500A) removed outlier: 5.769A pdb=" N LEU E 189 " --> pdb=" O LEU E 180 " (cutoff:3.500A) removed outlier: 3.513A pdb=" N TYR E 192 " --> pdb=" O ASN E 196 " (cutoff:3.500A) removed outlier: 3.625A pdb=" N ASN E 196 " --> pdb=" O TYR E 192 " (cutoff:3.500A) 402 hydrogen bonds defined for protein. 1104 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.63 Time building geometry restraints manager: 3.55 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.35: 2808 1.35 - 1.47: 2259 1.47 - 1.59: 3778 1.59 - 1.71: 0 1.71 - 1.83: 96 Bond restraints: 8941 Sorted by residual: bond pdb=" N THR E 212 " pdb=" CA THR E 212 " ideal model delta sigma weight residual 1.459 1.495 -0.036 1.32e-02 5.74e+03 7.61e+00 bond pdb=" N GLU E 236 " pdb=" CA GLU E 236 " ideal model delta sigma weight residual 1.455 1.488 -0.032 1.18e-02 7.18e+03 7.46e+00 bond pdb=" CA ALA C 299 " pdb=" C ALA C 299 " ideal model delta sigma weight residual 1.523 1.479 0.044 1.80e-02 3.09e+03 6.01e+00 bond pdb=" N TYR E 237 " pdb=" CA TYR E 237 " ideal model delta sigma weight residual 1.452 1.487 -0.035 1.41e-02 5.03e+03 5.99e+00 bond pdb=" C GLN C 44 " pdb=" N MET C 45 " ideal model delta sigma weight residual 1.330 1.298 0.032 1.31e-02 5.83e+03 5.87e+00 ... (remaining 8936 not shown) Histogram of bond angle deviations from ideal: 99.17 - 106.15: 150 106.15 - 113.13: 4805 113.13 - 120.11: 3197 120.11 - 127.09: 3855 127.09 - 134.07: 92 Bond angle restraints: 12099 Sorted by residual: angle pdb=" N SER R1168 " pdb=" CA SER R1168 " pdb=" C SER R1168 " ideal model delta sigma weight residual 113.23 107.67 5.56 1.22e+00 6.72e-01 2.08e+01 angle pdb=" N LYS B 349 " pdb=" CA LYS B 349 " pdb=" C LYS B 349 " ideal model delta sigma weight residual 114.04 108.94 5.10 1.24e+00 6.50e-01 1.69e+01 angle pdb=" C PRO E 14 " pdb=" CA PRO E 14 " pdb=" CB PRO E 14 " ideal model delta sigma weight residual 111.39 116.45 -5.06 1.28e+00 6.10e-01 1.57e+01 angle pdb=" N PHE R 970 " pdb=" CA PHE R 970 " pdb=" C PHE R 970 " ideal model delta sigma weight residual 109.15 114.77 -5.62 1.44e+00 4.82e-01 1.52e+01 angle pdb=" C LEU C 55 " pdb=" CA LEU C 55 " pdb=" CB LEU C 55 " ideal model delta sigma weight residual 109.65 103.09 6.56 1.68e+00 3.54e-01 1.52e+01 ... (remaining 12094 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.55: 4853 17.55 - 35.10: 360 35.10 - 52.65: 52 52.65 - 70.20: 3 70.20 - 87.75: 6 Dihedral angle restraints: 5274 sinusoidal: 2047 harmonic: 3227 Sorted by residual: dihedral pdb=" CB CYS E 161 " pdb=" SG CYS E 161 " pdb=" SG CYS E 231 " pdb=" CB CYS E 231 " ideal model delta sinusoidal sigma weight residual 93.00 135.64 -42.64 1 1.00e+01 1.00e-02 2.54e+01 dihedral pdb=" CB CYS R1008 " pdb=" SG CYS R1008 " pdb=" SG CYS R1080 " pdb=" CB CYS R1080 " ideal model delta sinusoidal sigma weight residual -86.00 -47.47 -38.53 1 1.00e+01 1.00e-02 2.09e+01 dihedral pdb=" CA ARG E 193 " pdb=" C ARG E 193 " pdb=" N MET E 194 " pdb=" CA MET E 194 " ideal model delta harmonic sigma weight residual 180.00 160.55 19.45 0 5.00e+00 4.00e-02 1.51e+01 ... (remaining 5271 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.051: 1023 0.051 - 0.102: 282 0.102 - 0.153: 56 0.153 - 0.204: 6 0.204 - 0.255: 2 Chirality restraints: 1369 Sorted by residual: chirality pdb=" CA LEU R 928 " pdb=" N LEU R 928 " pdb=" C LEU R 928 " pdb=" CB LEU R 928 " both_signs ideal model delta sigma weight residual False 2.51 2.26 0.26 2.00e-01 2.50e+01 1.63e+00 chirality pdb=" CA GLU B 216 " pdb=" N GLU B 216 " pdb=" C GLU B 216 " pdb=" CB GLU B 216 " both_signs ideal model delta sigma weight residual False 2.51 2.26 0.25 2.00e-01 2.50e+01 1.57e+00 chirality pdb=" CA ALA C 56 " pdb=" N ALA C 56 " pdb=" C ALA C 56 " pdb=" CB ALA C 56 " both_signs ideal model delta sigma weight residual False 2.48 2.68 -0.20 2.00e-01 2.50e+01 9.81e-01 ... (remaining 1366 not shown) Planarity restraints: 1525 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C ARG E 182 " -0.050 5.00e-02 4.00e+02 7.57e-02 9.18e+00 pdb=" N PRO E 183 " 0.131 5.00e-02 4.00e+02 pdb=" CA PRO E 183 " -0.039 5.00e-02 4.00e+02 pdb=" CD PRO E 183 " -0.041 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CA SER E 210 " -0.009 2.00e-02 2.50e+03 1.88e-02 3.53e+00 pdb=" C SER E 210 " 0.032 2.00e-02 2.50e+03 pdb=" O SER E 210 " -0.012 2.00e-02 2.50e+03 pdb=" N GLY E 211 " -0.011 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE C 222 " -0.007 2.00e-02 2.50e+03 1.37e-02 3.29e+00 pdb=" CG PHE C 222 " 0.028 2.00e-02 2.50e+03 pdb=" CD1 PHE C 222 " -0.019 2.00e-02 2.50e+03 pdb=" CD2 PHE C 222 " -0.005 2.00e-02 2.50e+03 pdb=" CE1 PHE C 222 " 0.008 2.00e-02 2.50e+03 pdb=" CE2 PHE C 222 " -0.006 2.00e-02 2.50e+03 pdb=" CZ PHE C 222 " 0.002 2.00e-02 2.50e+03 ... (remaining 1522 not shown) Histogram of nonbonded interaction distances: 2.18 - 2.72: 447 2.72 - 3.27: 8475 3.27 - 3.81: 13324 3.81 - 4.36: 16117 4.36 - 4.90: 28133 Nonbonded interactions: 66496 Sorted by model distance: nonbonded pdb=" O GLY E 211 " pdb=" OG1 THR E 212 " model vdw 2.176 2.440 nonbonded pdb=" OG1 THR C 274 " pdb=" O VAL C 315 " model vdw 2.258 2.440 nonbonded pdb=" O ALA B 41 " pdb=" OG SER B 44 " model vdw 2.259 2.440 nonbonded pdb=" O GLU R1167 " pdb=" OG SER R1168 " model vdw 2.276 2.440 nonbonded pdb=" O ARG C 314 " pdb=" OG SER C 331 " model vdw 2.293 2.440 ... (remaining 66491 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.810 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: 3.860 Check model and map are aligned: 0.100 Set scattering table: 0.080 Process input model: 26.100 Find NCS groups from input model: 0.290 Set up NCS constraints: 0.040 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.000 Load rotamer database and sin/cos tables:12.760 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 44.050 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5858 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.053 8941 Z= 0.354 Angle : 0.803 9.493 12099 Z= 0.464 Chirality : 0.049 0.255 1369 Planarity : 0.005 0.076 1525 Dihedral : 12.351 87.746 3188 Min Nonbonded Distance : 2.176 Molprobity Statistics. All-atom Clashscore : 5.06 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.01 % Favored : 92.99 % Rotamer: Outliers : 0.00 % Allowed : 1.79 % Favored : 98.21 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.17 % 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.33 (0.25), residues: 1099 helix: 0.32 (0.29), residues: 348 sheet: -0.19 (0.32), residues: 276 loop : -2.27 (0.25), residues: 475 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.018 0.002 TRP C 82 HIS 0.008 0.001 HIS R1016 PHE 0.028 0.003 PHE C 222 TYR 0.031 0.002 TYR C 105 ARG 0.006 0.000 ARG C 304 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 117 residues out of total 952 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 117 time to evaluate : 1.018 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 117 average time/residue: 0.9794 time to fit residues: 124.3884 Evaluate side-chains 105 residues out of total 952 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 105 time to evaluate : 0.918 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.3919 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 109 random chunks: chunk 92 optimal weight: 0.9980 chunk 82 optimal weight: 0.5980 chunk 45 optimal weight: 0.5980 chunk 28 optimal weight: 0.7980 chunk 55 optimal weight: 0.0030 chunk 44 optimal weight: 0.9990 chunk 85 optimal weight: 0.7980 chunk 33 optimal weight: 1.9990 chunk 51 optimal weight: 1.9990 chunk 63 optimal weight: 0.9980 chunk 98 optimal weight: 0.9980 overall best weight: 0.5590 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 213 HIS B 255 ASN ** C 119 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** E 173 ASN E 233 GLN R1180 ASN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5835 moved from start: 0.0990 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.034 8941 Z= 0.222 Angle : 0.613 9.682 12099 Z= 0.325 Chirality : 0.043 0.193 1369 Planarity : 0.004 0.056 1525 Dihedral : 5.194 47.235 1201 Min Nonbonded Distance : 2.107 Molprobity Statistics. All-atom Clashscore : 10.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.46 % Favored : 93.54 % Rotamer: Outliers : 1.58 % Allowed : 8.51 % Favored : 89.92 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.17 % 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.83 (0.25), residues: 1099 helix: 0.81 (0.29), residues: 348 sheet: 0.12 (0.32), residues: 271 loop : -2.03 (0.25), residues: 480 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.001 TRP C 211 HIS 0.005 0.001 HIS R1016 PHE 0.018 0.002 PHE C 253 TYR 0.020 0.002 TYR C 105 ARG 0.004 0.000 ARG C 304 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 125 residues out of total 952 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 110 time to evaluate : 1.026 Fit side-chains revert: symmetry clash outliers start: 15 outliers final: 6 residues processed: 115 average time/residue: 1.0745 time to fit residues: 133.4685 Evaluate side-chains 115 residues out of total 952 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 109 time to evaluate : 0.905 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 2 residues processed: 4 average time/residue: 0.9384 time to fit residues: 5.4329 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 109 random chunks: chunk 55 optimal weight: 9.9990 chunk 30 optimal weight: 0.9990 chunk 82 optimal weight: 0.9980 chunk 67 optimal weight: 0.8980 chunk 27 optimal weight: 0.9980 chunk 99 optimal weight: 0.0040 chunk 107 optimal weight: 0.9980 chunk 88 optimal weight: 0.4980 chunk 98 optimal weight: 0.7980 chunk 33 optimal weight: 1.9990 chunk 79 optimal weight: 0.8980 overall best weight: 0.6192 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 213 HIS B 255 ASN ** C 119 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5838 moved from start: 0.1341 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.031 8941 Z= 0.228 Angle : 0.597 10.207 12099 Z= 0.317 Chirality : 0.042 0.219 1369 Planarity : 0.004 0.049 1525 Dihedral : 5.013 46.195 1201 Min Nonbonded Distance : 2.126 Molprobity Statistics. All-atom Clashscore : 10.23 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.55 % Favored : 93.45 % Rotamer: Outliers : 1.89 % Allowed : 11.24 % Favored : 86.87 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.17 % 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.65 (0.25), residues: 1099 helix: 0.93 (0.29), residues: 354 sheet: 0.20 (0.32), residues: 283 loop : -1.97 (0.26), residues: 462 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP R1143 HIS 0.006 0.001 HIS R1016 PHE 0.018 0.002 PHE C 253 TYR 0.022 0.002 TYR C 105 ARG 0.003 0.000 ARG C 219 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 127 residues out of total 952 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 109 time to evaluate : 1.031 Fit side-chains revert: symmetry clash outliers start: 18 outliers final: 8 residues processed: 116 average time/residue: 1.1041 time to fit residues: 138.0070 Evaluate side-chains 112 residues out of total 952 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 104 time to evaluate : 0.997 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 3 residues processed: 5 average time/residue: 0.3215 time to fit residues: 3.4434 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 109 random chunks: chunk 97 optimal weight: 6.9990 chunk 74 optimal weight: 0.6980 chunk 51 optimal weight: 0.8980 chunk 10 optimal weight: 0.4980 chunk 47 optimal weight: 0.9990 chunk 66 optimal weight: 1.9990 chunk 99 optimal weight: 4.9990 chunk 105 optimal weight: 0.0010 chunk 94 optimal weight: 0.6980 chunk 28 optimal weight: 0.5980 chunk 87 optimal weight: 0.8980 overall best weight: 0.4986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 213 HIS B 255 ASN ** C 119 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** C 268 ASN ** C 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5815 moved from start: 0.1610 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.031 8941 Z= 0.193 Angle : 0.573 10.450 12099 Z= 0.301 Chirality : 0.041 0.153 1369 Planarity : 0.004 0.045 1525 Dihedral : 4.762 44.827 1201 Min Nonbonded Distance : 2.146 Molprobity Statistics. All-atom Clashscore : 9.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.55 % Favored : 94.45 % Rotamer: Outliers : 1.79 % Allowed : 13.66 % Favored : 84.56 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.17 % 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.48 (0.26), residues: 1099 helix: 1.11 (0.29), residues: 353 sheet: 0.17 (0.31), residues: 294 loop : -1.83 (0.27), residues: 452 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP C 211 HIS 0.007 0.001 HIS B 244 PHE 0.017 0.002 PHE C 253 TYR 0.018 0.002 TYR C 105 ARG 0.007 0.000 ARG C 52 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 124 residues out of total 952 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 17 poor density : 107 time to evaluate : 0.988 Fit side-chains revert: symmetry clash outliers start: 17 outliers final: 12 residues processed: 117 average time/residue: 0.9949 time to fit residues: 126.7508 Evaluate side-chains 116 residues out of total 952 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 104 time to evaluate : 0.956 Switching outliers to nearest non-outliers outliers start: 12 outliers final: 7 residues processed: 5 average time/residue: 0.5010 time to fit residues: 4.2856 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 109 random chunks: chunk 59 optimal weight: 0.0470 chunk 1 optimal weight: 1.9990 chunk 78 optimal weight: 0.9990 chunk 43 optimal weight: 2.9990 chunk 89 optimal weight: 0.9980 chunk 72 optimal weight: 0.9980 chunk 0 optimal weight: 3.9990 chunk 53 optimal weight: 0.0570 chunk 94 optimal weight: 0.9980 chunk 26 optimal weight: 0.6980 chunk 35 optimal weight: 0.8980 overall best weight: 0.5396 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 213 HIS B 255 ASN ** C 119 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 340 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** E 173 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5819 moved from start: 0.1760 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.030 8941 Z= 0.201 Angle : 0.569 10.176 12099 Z= 0.300 Chirality : 0.041 0.148 1369 Planarity : 0.004 0.043 1525 Dihedral : 4.698 45.389 1201 Min Nonbonded Distance : 2.143 Molprobity Statistics. All-atom Clashscore : 9.48 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.64 % Favored : 94.36 % Rotamer: Outliers : 1.89 % Allowed : 14.29 % Favored : 83.82 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.17 % 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.40 (0.26), residues: 1099 helix: 1.19 (0.29), residues: 352 sheet: 0.14 (0.31), residues: 294 loop : -1.75 (0.27), residues: 453 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.009 0.001 TRP C 211 HIS 0.005 0.001 HIS R1016 PHE 0.017 0.002 PHE C 241 TYR 0.019 0.001 TYR C 105 ARG 0.005 0.000 ARG C 52 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2198 Ramachandran restraints generated. 1099 Oldfield, 0 Emsley, 1099 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 124 residues out of total 952 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 106 time to evaluate : 1.016 Fit side-chains revert: symmetry clash outliers start: 18 outliers final: 11 residues processed: 117 average time/residue: 1.0715 time to fit residues: 135.2675 Evaluate side-chains 115 residues out of total 952 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 104 time to evaluate : 0.910 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 7 residues processed: 4 average time/residue: 0.3412 time to fit residues: 2.9803 ------------------------------------------------------------------------------- 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-1.21rc1-5172/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-1.21rc1-5172/modules/cctbx_project/iotbx/cli_parser.py", line 931, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/phenix/phenix/programs/real_space_refine.py", line 191, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 312, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 737, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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: 63.4474 > 50: distance: 53 - 69: 22.362 distance: 57 - 78: 20.131 distance: 61 - 86: 26.655 distance: 64 - 69: 25.609 distance: 65 - 98: 28.755 distance: 69 - 70: 12.540 distance: 70 - 71: 55.451 distance: 71 - 72: 45.213 distance: 71 - 78: 44.901 distance: 72 - 106: 40.603 distance: 73 - 74: 12.323 distance: 74 - 75: 35.656 distance: 75 - 76: 31.503 distance: 75 - 77: 32.763 distance: 78 - 79: 33.989 distance: 79 - 80: 14.025 distance: 79 - 82: 18.870 distance: 80 - 81: 16.343 distance: 80 - 86: 12.535 distance: 81 - 114: 24.779 distance: 82 - 83: 19.154 distance: 83 - 84: 40.710 distance: 83 - 85: 53.823 distance: 86 - 87: 23.877 distance: 87 - 88: 29.392 distance: 87 - 90: 18.284 distance: 88 - 89: 31.507 distance: 88 - 98: 43.512 distance: 89 - 122: 44.763 distance: 90 - 91: 34.255 distance: 91 - 92: 50.638 distance: 91 - 93: 33.996 distance: 92 - 94: 35.313 distance: 93 - 95: 14.475 distance: 94 - 96: 5.966 distance: 95 - 96: 23.317 distance: 96 - 97: 6.417 distance: 98 - 99: 41.478 distance: 99 - 100: 41.928 distance: 99 - 102: 20.056 distance: 100 - 101: 17.818 distance: 100 - 106: 25.010 distance: 102 - 103: 25.586 distance: 102 - 104: 16.135 distance: 103 - 105: 32.239 distance: 106 - 107: 37.155 distance: 107 - 108: 37.845 distance: 107 - 110: 24.791 distance: 108 - 109: 32.720 distance: 108 - 114: 21.527 distance: 110 - 111: 42.114 distance: 111 - 112: 27.140 distance: 112 - 113: 17.391 distance: 114 - 115: 19.495 distance: 115 - 116: 38.719 distance: 115 - 118: 20.580 distance: 116 - 117: 11.061 distance: 116 - 122: 48.692 distance: 118 - 119: 33.559 distance: 119 - 120: 19.595 distance: 119 - 121: 38.660 distance: 122 - 123: 33.953 distance: 123 - 124: 21.239 distance: 123 - 126: 12.099 distance: 124 - 125: 11.423 distance: 124 - 129: 40.767 distance: 126 - 127: 23.159 distance: 126 - 128: 16.025 distance: 129 - 130: 7.224 distance: 130 - 131: 35.392 distance: 130 - 133: 23.106 distance: 131 - 132: 57.227 distance: 131 - 138: 18.201 distance: 133 - 134: 32.823 distance: 134 - 135: 6.652 distance: 135 - 136: 8.308 distance: 135 - 137: 31.500 distance: 138 - 139: 24.925 distance: 139 - 140: 53.685 distance: 139 - 142: 41.252 distance: 140 - 141: 23.768 distance: 140 - 145: 40.612 distance: 142 - 143: 43.677 distance: 142 - 144: 27.190 distance: 145 - 146: 18.691 distance: 146 - 147: 25.210 distance: 146 - 149: 35.113 distance: 147 - 148: 32.417 distance: 147 - 156: 33.880 distance: 149 - 150: 13.952 distance: 150 - 151: 31.103 distance: 150 - 152: 18.244 distance: 151 - 153: 32.440 distance: 152 - 154: 15.985 distance: 153 - 155: 22.319 distance: 154 - 155: 26.451