Starting phenix.real_space_refine on Wed Jan 15 13:31:22 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/8vbe_43122/01_2025/8vbe_43122.cif Found real_map, /net/cci-nas-00/data/ceres_data/8vbe_43122/01_2025/8vbe_43122.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.6 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/8vbe_43122/01_2025/8vbe_43122.map" default_real_map = "/net/cci-nas-00/data/ceres_data/8vbe_43122/01_2025/8vbe_43122.map" model { file = "/net/cci-nas-00/data/ceres_data/8vbe_43122/01_2025/8vbe_43122.cif" } default_model = "/net/cci-nas-00/data/ceres_data/8vbe_43122/01_2025/8vbe_43122.cif" } resolution = 2.6 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= -0.000 sd= 0.009 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: 6 Type Number sf(0) Gaussians P 42 5.49 5 Mg 2 5.21 5 S 14 5.16 5 C 5312 2.51 5 N 1407 2.21 5 O 1644 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 27 residue(s): 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/chem_data/mon_lib" Total number of atoms: 8421 Number of models: 1 Model: "" Number of chains: 4 Chain: "A" Number of atoms: 4384 Number of conformers: 1 Conformer: "" Number of residues, atoms: 538, 4384 Classifications: {'peptide': 538} Incomplete info: {'truncation_to_alanine': 2} Link IDs: {'PCIS': 2, 'PTRANS': 35, 'TRANS': 500} Unresolved non-hydrogen bonds: 7 Unresolved non-hydrogen angles: 9 Unresolved non-hydrogen dihedrals: 5 Planarities with less than four sites: {'GLU:plan': 1, 'ASN:plan1': 1} Unresolved non-hydrogen planarities: 7 Chain: "B" Number of atoms: 3203 Number of conformers: 1 Conformer: "" Number of residues, atoms: 393, 3203 Classifications: {'peptide': 393} Incomplete info: {'truncation_to_alanine': 9} Link IDs: {'PTRANS': 28, 'TRANS': 364} Chain breaks: 2 Unresolved non-hydrogen bonds: 43 Unresolved non-hydrogen angles: 53 Unresolved non-hydrogen dihedrals: 40 Unresolved non-hydrogen chiralities: 2 Planarities with less than four sites: {'GLN:plan1': 1, 'ARG:plan': 2, 'TRP:plan': 1} Unresolved non-hydrogen planarities: 23 Chain: "F" Number of atoms: 802 Number of conformers: 1 Conformer: "" Number of residues, atoms: 39, 802 Classifications: {'DNA': 37, 'RNA': 2} Modifications used: {'rna3p': 2} Link IDs: {'rna3p': 38} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 6 Chain: "A" Number of atoms: 32 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 32 Unusual residues: {' MG': 2, 'DTP': 1} Classifications: {'undetermined': 3} Link IDs: {None: 2} Time building chain proxies: 6.02, per 1000 atoms: 0.71 Number of scatterers: 8421 At special positions: 0 Unit cell: (80.025, 95.7, 122.925, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 6 Type Number sf(0) S 14 16.00 P 42 15.00 Mg 2 11.99 O 1644 8.00 N 1407 7.00 C 5312 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.00 Amino acid : True - 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: 2.24 Conformation dependent library (CDL) restraints added in 904.5 milliseconds 1846 Ramachandran restraints generated. 923 Oldfield, 0 Emsley, 923 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1760 Finding SS restraints... Secondary structure from input PDB file: 32 helices and 12 sheets defined 37.8% alpha, 15.8% beta 13 base pairs and 27 stacking pairs defined. Time for finding SS restraints: 3.23 Creating SS restraints... Processing helix chain 'A' and resid 27 through 45 removed outlier: 3.598A pdb=" N ALA A 33 " --> pdb=" O GLU A 29 " (cutoff:3.500A) Processing helix chain 'A' and resid 77 through 84 Processing helix chain 'A' and resid 96 through 100 Processing helix chain 'A' and resid 114 through 118 Processing helix chain 'A' and resid 121 through 129 removed outlier: 4.471A pdb=" N LYS A 126 " --> pdb=" O ASP A 123 " (cutoff:3.500A) removed outlier: 3.511A pdb=" N TYR A 127 " --> pdb=" O PHE A 124 " (cutoff:3.500A) removed outlier: 3.860A pdb=" N ALA A 129 " --> pdb=" O LYS A 126 " (cutoff:3.500A) Processing helix chain 'A' and resid 155 through 175 Proline residue: A 170 - end of helix removed outlier: 3.596A pdb=" N GLN A 174 " --> pdb=" O PRO A 170 " (cutoff:3.500A) Processing helix chain 'A' and resid 194 through 212 removed outlier: 3.644A pdb=" N THR A 200 " --> pdb=" O GLY A 196 " (cutoff:3.500A) removed outlier: 3.552A pdb=" N LYS A 201 " --> pdb=" O GLN A 197 " (cutoff:3.500A) removed outlier: 3.651A pdb=" N TRP A 212 " --> pdb=" O HIS A 208 " (cutoff:3.500A) Processing helix chain 'A' and resid 217 through 221 Processing helix chain 'A' and resid 253 through 268 Processing helix chain 'A' and resid 276 through 281 Processing helix chain 'A' and resid 282 through 284 No H-bonds generated for 'chain 'A' and resid 282 through 284' Processing helix chain 'A' and resid 296 through 310 removed outlier: 3.557A pdb=" N GLU A 300 " --> pdb=" O THR A 296 " (cutoff:3.500A) removed outlier: 3.575A pdb=" N GLU A 302 " --> pdb=" O GLU A 298 " (cutoff:3.500A) Processing helix chain 'A' and resid 363 through 384 removed outlier: 3.718A pdb=" N GLN A 367 " --> pdb=" O ASN A 363 " (cutoff:3.500A) Processing helix chain 'A' and resid 394 through 405 removed outlier: 4.057A pdb=" N GLU A 404 " --> pdb=" O THR A 400 " (cutoff:3.500A) Processing helix chain 'A' and resid 473 through 488 removed outlier: 3.870A pdb=" N THR A 477 " --> pdb=" O THR A 473 " (cutoff:3.500A) removed outlier: 4.139A pdb=" N GLU A 478 " --> pdb=" O ASN A 474 " (cutoff:3.500A) removed outlier: 4.818A pdb=" N LEU A 479 " --> pdb=" O GLN A 475 " (cutoff:3.500A) Processing helix chain 'A' and resid 499 through 507 removed outlier: 3.661A pdb=" N LEU A 503 " --> pdb=" O SER A 499 " (cutoff:3.500A) Processing helix chain 'A' and resid 515 through 528 Processing helix chain 'B' and resid 27 through 44 removed outlier: 3.509A pdb=" N GLU B 40 " --> pdb=" O GLU B 36 " (cutoff:3.500A) Processing helix chain 'B' and resid 77 through 83 removed outlier: 3.767A pdb=" N ASN B 81 " --> pdb=" O PHE B 77 " (cutoff:3.500A) removed outlier: 3.737A pdb=" N LYS B 82 " --> pdb=" O ARG B 78 " (cutoff:3.500A) removed outlier: 3.660A pdb=" N ARG B 83 " --> pdb=" O GLU B 79 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 77 through 83' Processing helix chain 'B' and resid 112 through 118 removed outlier: 3.700A pdb=" N PHE B 116 " --> pdb=" O ASP B 113 " (cutoff:3.500A) removed outlier: 3.504A pdb=" N SER B 117 " --> pdb=" O ALA B 114 " (cutoff:3.500A) Processing helix chain 'B' and resid 124 through 127 Processing helix chain 'B' and resid 134 through 138 Processing helix chain 'B' and resid 156 through 165 removed outlier: 5.469A pdb=" N SER B 162 " --> pdb=" O ALA B 158 " (cutoff:3.500A) removed outlier: 4.799A pdb=" N SER B 163 " --> pdb=" O ILE B 159 " (cutoff:3.500A) Processing helix chain 'B' and resid 168 through 173 Processing helix chain 'B' and resid 197 through 212 Processing helix chain 'B' and resid 235 through 239 Processing helix chain 'B' and resid 253 through 271 removed outlier: 3.613A pdb=" N LYS B 263 " --> pdb=" O LYS B 259 " (cutoff:3.500A) removed outlier: 3.679A pdb=" N TRP B 266 " --> pdb=" O GLY B 262 " (cutoff:3.500A) removed outlier: 4.031A pdb=" N GLN B 269 " --> pdb=" O ASN B 265 " (cutoff:3.500A) removed outlier: 4.118A pdb=" N ILE B 270 " --> pdb=" O TRP B 266 " (cutoff:3.500A) Processing helix chain 'B' and resid 280 through 284 removed outlier: 3.776A pdb=" N ARG B 284 " --> pdb=" O LYS B 281 " (cutoff:3.500A) Processing helix chain 'B' and resid 296 through 309 Processing helix chain 'B' and resid 363 through 382 removed outlier: 3.909A pdb=" N GLN B 367 " --> pdb=" O ASN B 363 " (cutoff:3.500A) Processing helix chain 'B' and resid 394 through 402 Processing helix chain 'B' and resid 403 through 405 No H-bonds generated for 'chain 'B' and resid 403 through 405' Processing sheet with id=AA1, first strand: chain 'A' and resid 47 through 49 removed outlier: 3.556A pdb=" N ILE A 142 " --> pdb=" O ILE A 132 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 60 through 64 Processing sheet with id=AA3, first strand: chain 'A' and resid 105 through 110 Processing sheet with id=AA4, first strand: chain 'A' and resid 227 through 229 Processing sheet with id=AA5, first strand: chain 'A' and resid 347 through 354 removed outlier: 8.621A pdb=" N ASN A 348 " --> pdb=" O GLN A 343 " (cutoff:3.500A) removed outlier: 6.757A pdb=" N GLN A 343 " --> pdb=" O ASN A 348 " (cutoff:3.500A) removed outlier: 6.584A pdb=" N PHE A 389 " --> pdb=" O GLU A 415 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'A' and resid 361 through 362 Processing sheet with id=AA7, first strand: chain 'A' and resid 441 through 447 removed outlier: 3.786A pdb=" N LYS A 465 " --> pdb=" O TYR A 457 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'A' and resid 492 through 497 Processing sheet with id=AA9, first strand: chain 'B' and resid 47 through 49 Processing sheet with id=AB1, first strand: chain 'B' and resid 60 through 64 Processing sheet with id=AB2, first strand: chain 'B' and resid 105 through 110 Processing sheet with id=AB3, first strand: chain 'B' and resid 347 through 355 removed outlier: 8.581A pdb=" N ASN B 348 " --> pdb=" O GLN B 343 " (cutoff:3.500A) removed outlier: 7.023A pdb=" N GLN B 343 " --> pdb=" O ASN B 348 " (cutoff:3.500A) removed outlier: 3.925A pdb=" N ILE B 329 " --> pdb=" O LYS B 390 " (cutoff:3.500A) removed outlier: 5.868A pdb=" N PHE B 389 " --> pdb=" O GLU B 415 " (cutoff:3.500A) 300 hydrogen bonds defined for protein. 801 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 35 hydrogen bonds 70 hydrogen bond angles 0 basepair planarities 13 basepair parallelities 27 stacking parallelities Total time for adding SS restraints: 2.86 Time building geometry restraints manager: 2.56 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.20 - 1.32: 1392 1.32 - 1.45: 2403 1.45 - 1.57: 4819 1.57 - 1.69: 82 1.69 - 1.81: 24 Bond restraints: 8720 Sorted by residual: bond pdb=" N ILE A 132 " pdb=" CA ILE A 132 " ideal model delta sigma weight residual 1.460 1.490 -0.030 7.60e-03 1.73e+04 1.56e+01 bond pdb=" C ASP A 110 " pdb=" N VAL A 111 " ideal model delta sigma weight residual 1.329 1.271 0.058 1.69e-02 3.50e+03 1.19e+01 bond pdb=" N VAL B 381 " pdb=" CA VAL B 381 " ideal model delta sigma weight residual 1.461 1.502 -0.040 1.17e-02 7.31e+03 1.18e+01 bond pdb=" N VAL A 108 " pdb=" CA VAL A 108 " ideal model delta sigma weight residual 1.458 1.496 -0.038 1.17e-02 7.31e+03 1.04e+01 bond pdb=" N ILE A 380 " pdb=" CA ILE A 380 " ideal model delta sigma weight residual 1.460 1.495 -0.036 1.21e-02 6.83e+03 8.64e+00 ... (remaining 8715 not shown) Histogram of bond angle deviations from ideal: 0.00 - 1.89: 11698 1.89 - 3.78: 273 3.78 - 5.67: 41 5.67 - 7.56: 6 7.56 - 9.45: 3 Bond angle restraints: 12021 Sorted by residual: angle pdb=" CA ILE B 382 " pdb=" C ILE B 382 " pdb=" O ILE B 382 " ideal model delta sigma weight residual 121.68 116.33 5.35 7.90e-01 1.60e+00 4.58e+01 angle pdb=" O3' DG F 19 " pdb=" C3' DG F 19 " pdb=" C2' DG F 19 " ideal model delta sigma weight residual 111.50 102.05 9.45 1.50e+00 4.44e-01 3.97e+01 angle pdb=" O3' DG F 14 " pdb=" C3' DG F 14 " pdb=" C2' DG F 14 " ideal model delta sigma weight residual 111.50 103.51 7.99 1.50e+00 4.44e-01 2.84e+01 angle pdb=" C ILE A 382 " pdb=" N TRP A 383 " pdb=" CA TRP A 383 " ideal model delta sigma weight residual 122.23 114.01 8.22 1.57e+00 4.06e-01 2.74e+01 angle pdb=" PB DTP A 601 " pdb=" O3B DTP A 601 " pdb=" PG DTP A 601 " ideal model delta sigma weight residual 139.87 135.45 4.42 1.00e+00 1.00e+00 1.96e+01 ... (remaining 12016 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 27.96: 4617 27.96 - 55.93: 466 55.93 - 83.89: 48 83.89 - 111.86: 4 111.86 - 139.82: 1 Dihedral angle restraints: 5136 sinusoidal: 2411 harmonic: 2725 Sorted by residual: dihedral pdb=" O1B DTP A 601 " pdb=" O3B DTP A 601 " pdb=" PB DTP A 601 " pdb=" PG DTP A 601 " ideal model delta sinusoidal sigma weight residual 299.93 160.11 139.82 1 3.00e+01 1.11e-03 1.88e+01 dihedral pdb=" CA VAL B 10 " pdb=" C VAL B 10 " pdb=" N LYS B 11 " pdb=" CA LYS B 11 " ideal model delta harmonic sigma weight residual 180.00 161.07 18.93 0 5.00e+00 4.00e-02 1.43e+01 dihedral pdb=" CB GLU B 203 " pdb=" CG GLU B 203 " pdb=" CD GLU B 203 " pdb=" OE1 GLU B 203 " ideal model delta sinusoidal sigma weight residual 0.00 89.60 -89.60 1 3.00e+01 1.11e-03 1.06e+01 ... (remaining 5133 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.053: 1055 0.053 - 0.106: 212 0.106 - 0.158: 43 0.158 - 0.211: 3 0.211 - 0.264: 2 Chirality restraints: 1315 Sorted by residual: chirality pdb=" CA ILE B 382 " pdb=" N ILE B 382 " pdb=" C ILE B 382 " pdb=" CB ILE B 382 " both_signs ideal model delta sigma weight residual False 2.43 2.70 -0.26 2.00e-01 2.50e+01 1.74e+00 chirality pdb=" C3' DG F 19 " pdb=" C4' DG F 19 " pdb=" O3' DG F 19 " pdb=" C2' DG F 19 " both_signs ideal model delta sigma weight residual False -2.66 -2.91 0.25 2.00e-01 2.50e+01 1.52e+00 chirality pdb=" CA ILE A 132 " pdb=" N ILE A 132 " pdb=" C ILE A 132 " pdb=" CB ILE A 132 " both_signs ideal model delta sigma weight residual False 2.43 2.64 -0.21 2.00e-01 2.50e+01 1.06e+00 ... (remaining 1312 not shown) Planarity restraints: 1372 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA SER A 379 " -0.013 2.00e-02 2.50e+03 2.63e-02 6.90e+00 pdb=" C SER A 379 " 0.045 2.00e-02 2.50e+03 pdb=" O SER A 379 " -0.017 2.00e-02 2.50e+03 pdb=" N ILE A 380 " -0.015 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA THR A 107 " 0.011 2.00e-02 2.50e+03 2.20e-02 4.82e+00 pdb=" C THR A 107 " -0.038 2.00e-02 2.50e+03 pdb=" O THR A 107 " 0.014 2.00e-02 2.50e+03 pdb=" N VAL A 108 " 0.013 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C1' DA F 34 " -0.026 2.00e-02 2.50e+03 1.16e-02 3.69e+00 pdb=" N9 DA F 34 " 0.026 2.00e-02 2.50e+03 pdb=" C8 DA F 34 " 0.003 2.00e-02 2.50e+03 pdb=" N7 DA F 34 " -0.002 2.00e-02 2.50e+03 pdb=" C5 DA F 34 " 0.000 2.00e-02 2.50e+03 pdb=" C6 DA F 34 " -0.003 2.00e-02 2.50e+03 pdb=" N6 DA F 34 " -0.006 2.00e-02 2.50e+03 pdb=" N1 DA F 34 " -0.004 2.00e-02 2.50e+03 pdb=" C2 DA F 34 " 0.003 2.00e-02 2.50e+03 pdb=" N3 DA F 34 " 0.004 2.00e-02 2.50e+03 pdb=" C4 DA F 34 " 0.006 2.00e-02 2.50e+03 ... (remaining 1369 not shown) Histogram of nonbonded interaction distances: 2.05 - 2.62: 89 2.62 - 3.19: 6581 3.19 - 3.76: 12891 3.76 - 4.33: 18087 4.33 - 4.90: 29498 Nonbonded interactions: 67146 Sorted by model distance: nonbonded pdb=" O VAL A 111 " pdb="MG MG A 602 " model vdw 2.052 2.170 nonbonded pdb=" O1A DTP A 601 " pdb="MG MG A 602 " model vdw 2.128 2.170 nonbonded pdb=" O ARG B 356 " pdb=" NH1 ARG B 356 " model vdw 2.129 3.120 nonbonded pdb=" NZ LYS B 172 " pdb=" O ILE B 178 " model vdw 2.200 3.120 nonbonded pdb=" OD2 ASP A 185 " pdb="MG MG A 602 " model vdw 2.287 2.170 ... (remaining 67141 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.00 ========== 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: 3.640 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: 0.350 Check model and map are aligned: 0.060 Set scattering table: 0.080 Process input model: 25.220 Find NCS groups from input model: 0.090 Set up NCS constraints: 0.030 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.010 Load rotamer database and sin/cos tables:13.830 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 43.320 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8393 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.058 8720 Z= 0.235 Angle : 0.682 9.451 12021 Z= 0.415 Chirality : 0.046 0.264 1315 Planarity : 0.005 0.047 1372 Dihedral : 19.817 139.825 3376 Min Nonbonded Distance : 2.052 Molprobity Statistics. All-atom Clashscore : 4.72 Ramachandran Plot: Outliers : 0.54 % Allowed : 4.01 % Favored : 95.45 % Rotamer: Outliers : 4.46 % Allowed : 20.99 % Favored : 74.55 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.08 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.35 (0.28), residues: 923 helix: 1.03 (0.32), residues: 309 sheet: 0.06 (0.38), residues: 188 loop : -1.33 (0.28), residues: 426 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP B 410 HIS 0.003 0.001 HIS B 198 PHE 0.008 0.001 PHE B 61 TYR 0.010 0.001 TYR A 501 ARG 0.015 0.001 ARG B 277 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1846 Ramachandran restraints generated. 923 Oldfield, 0 Emsley, 923 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1846 Ramachandran restraints generated. 923 Oldfield, 0 Emsley, 923 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 88 residues out of total 840 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 37 poor density : 51 time to evaluate : 0.877 Fit side-chains REVERT: B 113 ASP cc_start: 0.8031 (OUTLIER) cc_final: 0.7711 (m-30) REVERT: B 166 LYS cc_start: 0.6225 (OUTLIER) cc_final: 0.5560 (pttm) REVERT: B 277 ARG cc_start: 0.7553 (ttp80) cc_final: 0.6950 (ttp80) REVERT: B 281 LYS cc_start: 0.7452 (ttmt) cc_final: 0.7029 (mtmm) outliers start: 37 outliers final: 27 residues processed: 86 average time/residue: 1.2928 time to fit residues: 118.7338 Evaluate side-chains 77 residues out of total 840 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 29 poor density : 48 time to evaluate : 1.226 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 27 THR Chi-restraints excluded: chain A residue 40 GLU Chi-restraints excluded: chain A residue 50 ILE Chi-restraints excluded: chain A residue 53 GLU Chi-restraints excluded: chain A residue 220 LYS Chi-restraints excluded: chain A residue 245 VAL Chi-restraints excluded: chain A residue 253 THR Chi-restraints excluded: chain A residue 301 LEU Chi-restraints excluded: chain A residue 372 VAL Chi-restraints excluded: chain A residue 459 THR Chi-restraints excluded: chain A residue 475 GLN Chi-restraints excluded: chain A residue 503 LEU Chi-restraints excluded: chain B residue 80 LEU Chi-restraints excluded: chain B residue 102 LYS Chi-restraints excluded: chain B residue 113 ASP Chi-restraints excluded: chain B residue 166 LYS Chi-restraints excluded: chain B residue 168 LEU Chi-restraints excluded: chain B residue 169 GLU Chi-restraints excluded: chain B residue 191 SER Chi-restraints excluded: chain B residue 195 ILE Chi-restraints excluded: chain B residue 245 VAL Chi-restraints excluded: chain B residue 250 ASP Chi-restraints excluded: chain B residue 296 THR Chi-restraints excluded: chain B residue 315 HIS Chi-restraints excluded: chain B residue 362 THR Chi-restraints excluded: chain B residue 386 THR Chi-restraints excluded: chain B residue 397 THR Chi-restraints excluded: chain B residue 417 VAL Chi-restraints excluded: chain B residue 419 THR Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 95 random chunks: chunk 80 optimal weight: 4.9990 chunk 72 optimal weight: 7.9990 chunk 39 optimal weight: 2.9990 chunk 24 optimal weight: 2.9990 chunk 48 optimal weight: 1.9990 chunk 38 optimal weight: 1.9990 chunk 74 optimal weight: 5.9990 chunk 28 optimal weight: 2.9990 chunk 45 optimal weight: 0.8980 chunk 55 optimal weight: 0.4980 chunk 86 optimal weight: 0.6980 overall best weight: 1.2184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 373 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3689 r_free = 0.3689 target = 0.126621 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 42)----------------| | r_work = 0.3307 r_free = 0.3307 target = 0.099157 restraints weight = 12601.321| |-----------------------------------------------------------------------------| r_work (start): 0.3301 rms_B_bonded: 2.32 r_work: 0.3178 rms_B_bonded: 2.82 restraints_weight: 0.5000 r_work: 0.3059 rms_B_bonded: 4.66 restraints_weight: 0.2500 r_work (final): 0.3059 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8454 moved from start: 0.0772 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.028 8720 Z= 0.295 Angle : 0.588 6.948 12021 Z= 0.316 Chirality : 0.044 0.168 1315 Planarity : 0.005 0.041 1372 Dihedral : 18.537 153.337 1475 Min Nonbonded Distance : 2.069 Molprobity Statistics. All-atom Clashscore : 3.99 Ramachandran Plot: Outliers : 0.43 % Allowed : 4.55 % Favored : 95.02 % Rotamer: Outliers : 5.19 % Allowed : 19.66 % Favored : 75.15 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.08 % 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.38 (0.28), residues: 923 helix: 1.10 (0.31), residues: 311 sheet: -0.08 (0.37), residues: 186 loop : -1.36 (0.28), residues: 426 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A 88 HIS 0.005 0.001 HIS A 198 PHE 0.010 0.001 PHE B 61 TYR 0.012 0.001 TYR A 457 ARG 0.006 0.001 ARG A 358 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1846 Ramachandran restraints generated. 923 Oldfield, 0 Emsley, 923 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1846 Ramachandran restraints generated. 923 Oldfield, 0 Emsley, 923 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 98 residues out of total 840 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 43 poor density : 55 time to evaluate : 0.877 Fit side-chains REVERT: A 50 ILE cc_start: 0.8512 (OUTLIER) cc_final: 0.8086 (mp) REVERT: A 413 GLU cc_start: 0.8811 (OUTLIER) cc_final: 0.8127 (tt0) REVERT: A 500 GLN cc_start: 0.8425 (OUTLIER) cc_final: 0.8146 (mt0) REVERT: B 374 LYS cc_start: 0.7937 (OUTLIER) cc_final: 0.7629 (mmmt) outliers start: 43 outliers final: 25 residues processed: 93 average time/residue: 1.3715 time to fit residues: 135.4645 Evaluate side-chains 79 residues out of total 840 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 29 poor density : 50 time to evaluate : 0.901 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 27 THR Chi-restraints excluded: chain A residue 40 GLU Chi-restraints excluded: chain A residue 50 ILE Chi-restraints excluded: chain A residue 53 GLU Chi-restraints excluded: chain A residue 178 ILE Chi-restraints excluded: chain A residue 245 VAL Chi-restraints excluded: chain A residue 253 THR Chi-restraints excluded: chain A residue 301 LEU Chi-restraints excluded: chain A residue 362 THR Chi-restraints excluded: chain A residue 413 GLU Chi-restraints excluded: chain A residue 475 GLN Chi-restraints excluded: chain A residue 500 GLN Chi-restraints excluded: chain A residue 503 LEU Chi-restraints excluded: chain B residue 12 LEU Chi-restraints excluded: chain B residue 80 LEU Chi-restraints excluded: chain B residue 102 LYS Chi-restraints excluded: chain B residue 117 SER Chi-restraints excluded: chain B residue 189 VAL Chi-restraints excluded: chain B residue 191 SER Chi-restraints excluded: chain B residue 296 THR Chi-restraints excluded: chain B residue 315 HIS Chi-restraints excluded: chain B residue 324 ASP Chi-restraints excluded: chain B residue 362 THR Chi-restraints excluded: chain B residue 374 LYS Chi-restraints excluded: chain B residue 382 ILE Chi-restraints excluded: chain B residue 386 THR Chi-restraints excluded: chain B residue 397 THR Chi-restraints excluded: chain B residue 417 VAL Chi-restraints excluded: chain B residue 419 THR Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 95 random chunks: chunk 22 optimal weight: 0.8980 chunk 54 optimal weight: 2.9990 chunk 59 optimal weight: 2.9990 chunk 75 optimal weight: 1.9990 chunk 62 optimal weight: 3.9990 chunk 72 optimal weight: 5.9990 chunk 46 optimal weight: 3.9990 chunk 2 optimal weight: 0.7980 chunk 31 optimal weight: 1.9990 chunk 45 optimal weight: 0.7980 chunk 87 optimal weight: 1.9990 overall best weight: 1.2984 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 373 GLN B 147 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3680 r_free = 0.3680 target = 0.126041 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 44)----------------| | r_work = 0.3296 r_free = 0.3296 target = 0.098576 restraints weight = 12686.003| |-----------------------------------------------------------------------------| r_work (start): 0.3289 rms_B_bonded: 2.34 r_work: 0.3165 rms_B_bonded: 2.83 restraints_weight: 0.5000 r_work: 0.3047 rms_B_bonded: 4.66 restraints_weight: 0.2500 r_work (final): 0.3047 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8464 moved from start: 0.0987 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.034 8720 Z= 0.308 Angle : 0.581 6.887 12021 Z= 0.311 Chirality : 0.044 0.162 1315 Planarity : 0.005 0.042 1372 Dihedral : 18.072 160.804 1453 Min Nonbonded Distance : 2.056 Molprobity Statistics. All-atom Clashscore : 4.12 Ramachandran Plot: Outliers : 0.43 % Allowed : 4.55 % Favored : 95.02 % Rotamer: Outliers : 5.07 % Allowed : 20.27 % Favored : 74.67 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.08 % 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.41 (0.27), residues: 923 helix: 1.13 (0.31), residues: 311 sheet: -0.15 (0.37), residues: 186 loop : -1.40 (0.28), residues: 426 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP A 402 HIS 0.005 0.001 HIS A 198 PHE 0.010 0.001 PHE B 61 TYR 0.012 0.001 TYR A 144 ARG 0.007 0.001 ARG A 358 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1846 Ramachandran restraints generated. 923 Oldfield, 0 Emsley, 923 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1846 Ramachandran restraints generated. 923 Oldfield, 0 Emsley, 923 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 95 residues out of total 840 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 42 poor density : 53 time to evaluate : 0.972 Fit side-chains REVERT: A 50 ILE cc_start: 0.8511 (OUTLIER) cc_final: 0.8068 (mp) REVERT: A 413 GLU cc_start: 0.8837 (OUTLIER) cc_final: 0.8162 (tt0) REVERT: B 277 ARG cc_start: 0.7409 (ttp80) cc_final: 0.7194 (ttp80) REVERT: B 374 LYS cc_start: 0.7903 (OUTLIER) cc_final: 0.7606 (mmmt) outliers start: 42 outliers final: 21 residues processed: 89 average time/residue: 1.4264 time to fit residues: 134.9645 Evaluate side-chains 75 residues out of total 840 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 24 poor density : 51 time to evaluate : 0.989 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 27 THR Chi-restraints excluded: chain A residue 40 GLU Chi-restraints excluded: chain A residue 50 ILE Chi-restraints excluded: chain A residue 53 GLU Chi-restraints excluded: chain A residue 178 ILE Chi-restraints excluded: chain A residue 245 VAL Chi-restraints excluded: chain A residue 253 THR Chi-restraints excluded: chain A residue 301 LEU Chi-restraints excluded: chain A residue 362 THR Chi-restraints excluded: chain A residue 413 GLU Chi-restraints excluded: chain A residue 475 GLN Chi-restraints excluded: chain A residue 536 VAL Chi-restraints excluded: chain B residue 68 SER Chi-restraints excluded: chain B residue 102 LYS Chi-restraints excluded: chain B residue 117 SER Chi-restraints excluded: chain B residue 189 VAL Chi-restraints excluded: chain B residue 275 LYS Chi-restraints excluded: chain B residue 317 VAL Chi-restraints excluded: chain B residue 362 THR Chi-restraints excluded: chain B residue 374 LYS Chi-restraints excluded: chain B residue 382 ILE Chi-restraints excluded: chain B residue 386 THR Chi-restraints excluded: chain B residue 397 THR Chi-restraints excluded: chain B residue 419 THR Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 95 random chunks: chunk 58 optimal weight: 1.9990 chunk 94 optimal weight: 3.9990 chunk 67 optimal weight: 0.7980 chunk 66 optimal weight: 1.9990 chunk 73 optimal weight: 2.9990 chunk 89 optimal weight: 2.9990 chunk 32 optimal weight: 0.9990 chunk 15 optimal weight: 1.9990 chunk 9 optimal weight: 1.9990 chunk 88 optimal weight: 0.2980 chunk 42 optimal weight: 0.3980 overall best weight: 0.8984 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 373 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3695 r_free = 0.3695 target = 0.127267 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 43)----------------| | r_work = 0.3315 r_free = 0.3315 target = 0.099829 restraints weight = 12583.796| |-----------------------------------------------------------------------------| r_work (start): 0.3311 rms_B_bonded: 2.33 r_work: 0.3188 rms_B_bonded: 2.83 restraints_weight: 0.5000 r_work: 0.3069 rms_B_bonded: 4.66 restraints_weight: 0.2500 r_work (final): 0.3069 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8445 moved from start: 0.1091 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.026 8720 Z= 0.224 Angle : 0.536 6.914 12021 Z= 0.287 Chirality : 0.042 0.155 1315 Planarity : 0.004 0.042 1372 Dihedral : 17.685 165.131 1444 Min Nonbonded Distance : 2.056 Molprobity Statistics. All-atom Clashscore : 4.30 Ramachandran Plot: Outliers : 0.43 % Allowed : 4.66 % Favored : 94.91 % Rotamer: Outliers : 5.07 % Allowed : 20.27 % Favored : 74.67 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 3.08 % 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.09 (0.28), residues: 923 helix: 1.65 (0.32), residues: 297 sheet: -0.03 (0.37), residues: 186 loop : -1.30 (0.28), residues: 440 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP A 406 HIS 0.005 0.001 HIS B 198 PHE 0.010 0.001 PHE B 61 TYR 0.011 0.001 TYR A 501 ARG 0.007 0.001 ARG B 284 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1846 Ramachandran restraints generated. 923 Oldfield, 0 Emsley, 923 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1846 Ramachandran restraints generated. 923 Oldfield, 0 Emsley, 923 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 99 residues out of total 840 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 42 poor density : 57 time to evaluate : 0.940 Fit side-chains REVERT: A 50 ILE cc_start: 0.8476 (OUTLIER) cc_final: 0.8079 (mp) REVERT: A 413 GLU cc_start: 0.8822 (OUTLIER) cc_final: 0.8155 (tt0) REVERT: B 206 ARG cc_start: 0.4485 (OUTLIER) cc_final: 0.3879 (ttp80) REVERT: B 278 GLN cc_start: 0.7573 (OUTLIER) cc_final: 0.7229 (tt0) REVERT: B 356 ARG cc_start: 0.7633 (pmt-80) cc_final: 0.7430 (pmt-80) outliers start: 42 outliers final: 20 residues processed: 90 average time/residue: 1.3771 time to fit residues: 131.6927 Evaluate side-chains 76 residues out of total 840 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 24 poor density : 52 time to evaluate : 0.981 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 27 THR Chi-restraints excluded: chain A residue 40 GLU Chi-restraints excluded: chain A residue 50 ILE Chi-restraints excluded: chain A residue 53 GLU Chi-restraints excluded: chain A residue 58 THR Chi-restraints excluded: chain A residue 178 ILE Chi-restraints excluded: chain A residue 253 THR Chi-restraints excluded: chain A residue 362 THR Chi-restraints excluded: chain A residue 413 GLU Chi-restraints excluded: chain A residue 475 GLN Chi-restraints excluded: chain A residue 503 LEU Chi-restraints excluded: chain B residue 102 LYS Chi-restraints excluded: chain B residue 169 GLU Chi-restraints excluded: chain B residue 189 VAL Chi-restraints excluded: chain B residue 206 ARG Chi-restraints excluded: chain B residue 275 LYS Chi-restraints excluded: chain B residue 278 GLN Chi-restraints excluded: chain B residue 296 THR Chi-restraints excluded: chain B residue 315 HIS Chi-restraints excluded: chain B residue 317 VAL Chi-restraints excluded: chain B residue 324 ASP Chi-restraints excluded: chain B residue 362 THR Chi-restraints excluded: chain B residue 386 THR Chi-restraints excluded: chain B residue 419 THR Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/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-5571/modules/cctbx_project/iotbx/cli_parser.py", line 946, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/phenix/phenix/programs/real_space_refine.py", line 210, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 306, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 702, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 767, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/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-5571/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-5571/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-5571/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-5571/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-5571/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-5571/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/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-5571/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-5571/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1525, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1427, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1308, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 69.1429 > 50: distance: 4 - 9: 29.629 distance: 9 - 10: 33.288 distance: 10 - 11: 48.678 distance: 10 - 13: 59.099 distance: 11 - 12: 46.519 distance: 11 - 17: 57.149 distance: 13 - 14: 21.648 distance: 13 - 15: 45.429 distance: 14 - 16: 43.078 distance: 17 - 18: 44.932 distance: 18 - 19: 38.974 distance: 19 - 20: 47.416 distance: 19 - 24: 58.142 distance: 21 - 22: 18.879 distance: 21 - 23: 40.538 distance: 24 - 25: 31.301 distance: 25 - 26: 45.215 distance: 25 - 28: 21.997 distance: 26 - 27: 39.080 distance: 26 - 32: 52.040 distance: 28 - 29: 19.574 distance: 29 - 30: 56.888 distance: 29 - 31: 43.678 distance: 32 - 33: 66.455 distance: 32 - 38: 51.258 distance: 33 - 34: 39.964 distance: 33 - 36: 32.266 distance: 34 - 35: 56.797 distance: 34 - 39: 50.539 distance: 36 - 37: 30.182 distance: 37 - 38: 39.604 distance: 39 - 40: 56.669 distance: 40 - 41: 49.304 distance: 40 - 43: 8.785 distance: 41 - 48: 58.457 distance: 43 - 44: 56.392 distance: 44 - 45: 45.073 distance: 45 - 46: 47.093 distance: 45 - 47: 47.099 distance: 48 - 49: 14.541 distance: 49 - 50: 29.537 distance: 49 - 52: 39.045 distance: 50 - 51: 32.320 distance: 50 - 57: 39.670 distance: 52 - 53: 14.484 distance: 53 - 54: 47.878 distance: 54 - 55: 39.498 distance: 55 - 56: 38.650 distance: 57 - 58: 40.399 distance: 58 - 59: 15.610 distance: 59 - 60: 22.508 distance: 59 - 65: 16.309 distance: 61 - 62: 16.338 distance: 62 - 63: 21.572 distance: 62 - 64: 44.631 distance: 66 - 69: 63.193 distance: 67 - 68: 56.321 distance: 67 - 71: 34.777 distance: 69 - 70: 43.390 distance: 71 - 72: 45.433 distance: 72 - 73: 44.611 distance: 72 - 75: 52.173 distance: 73 - 74: 23.229 distance: 73 - 85: 47.925 distance: 75 - 76: 8.976 distance: 76 - 77: 32.192 distance: 76 - 78: 48.372 distance: 77 - 79: 46.344 distance: 78 - 80: 66.172 distance: 79 - 80: 57.960 distance: 80 - 82: 46.744 distance: 82 - 84: 14.265 distance: 83 - 84: 39.476 distance: 85 - 86: 22.268 distance: 86 - 87: 44.912 distance: 86 - 89: 6.523 distance: 87 - 88: 29.547 distance: 87 - 92: 28.985 distance: 88 - 114: 34.395 distance: 89 - 90: 27.292 distance: 89 - 91: 34.932