Starting phenix.real_space_refine on Sun Mar 10 15:59:17 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8esa_28573/03_2024/8esa_28573.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8esa_28573/03_2024/8esa_28573.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.4 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8esa_28573/03_2024/8esa_28573.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8esa_28573/03_2024/8esa_28573.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8esa_28573/03_2024/8esa_28573.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8esa_28573/03_2024/8esa_28573.pdb" } resolution = 3.4 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= 1.000 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Anisotropic ADP refinement not supported. Converting 72 atoms to isotropic. 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 14 5.16 5 C 1970 2.51 5 N 563 2.21 5 O 593 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "B PHE 30": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 63": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/chem_data/mon_lib" Total number of atoms: 3140 Number of models: 1 Model: "" Number of chains: 3 Chain: "A" Number of atoms: 2223 Number of conformers: 1 Conformer: "" Number of residues, atoms: 273, 2223 Classifications: {'peptide': 273} Link IDs: {'PCIS': 1, 'PTRANS': 11, 'TRANS': 260} Chain breaks: 1 Chain: "B" Number of atoms: 828 Number of conformers: 1 Conformer: "" Number of residues, atoms: 99, 828 Classifications: {'peptide': 99} Link IDs: {'PCIS': 1, 'PTRANS': 4, 'TRANS': 93} Chain: "C" Number of atoms: 89 Number of conformers: 1 Conformer: "" Number of residues, atoms: 12, 89 Classifications: {'peptide': 12} Link IDs: {'TRANS': 11} Time building chain proxies: 2.21, per 1000 atoms: 0.70 Number of scatterers: 3140 At special positions: 0 Unit cell: (84.15, 65.45, 67.15, 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 14 16.00 O 593 8.00 N 563 7.00 C 1970 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=4, symmetry=0 Simple disulfide: pdb=" SG CYS A 84 " - pdb=" SG CYS C 12 " distance=2.03 Simple disulfide: pdb=" SG CYS A 101 " - pdb=" SG CYS A 164 " distance=2.03 Simple disulfide: pdb=" SG CYS A 203 " - pdb=" SG CYS A 259 " distance=2.03 Simple disulfide: pdb=" SG CYS B 25 " - pdb=" SG CYS B 80 " 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: 1.12 Conformation dependent library (CDL) restraints added in 560.2 milliseconds 752 Ramachandran restraints generated. 376 Oldfield, 0 Emsley, 376 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 716 Finding SS restraints... Secondary structure from input PDB file: 3 helices and 7 sheets defined 17.4% alpha, 33.3% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.36 Creating SS restraints... Processing helix chain 'A' and resid 56 through 85 removed outlier: 3.820A pdb=" N TRP A 60 " --> pdb=" O GLY A 56 " (cutoff:3.500A) removed outlier: 3.593A pdb=" N VAL A 76 " --> pdb=" O GLN A 72 " (cutoff:3.500A) Processing helix chain 'A' and resid 137 through 148 removed outlier: 3.699A pdb=" N HIS A 145 " --> pdb=" O GLN A 141 " (cutoff:3.500A) Processing helix chain 'A' and resid 151 through 175 removed outlier: 3.511A pdb=" N GLN A 155 " --> pdb=" O HIS A 151 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 31 through 34 removed outlier: 6.965A pdb=" N GLY A 26 " --> pdb=" O PHE A 33 " (cutoff:3.500A) removed outlier: 3.576A pdb=" N CYS A 101 " --> pdb=" O MET A 5 " (cutoff:3.500A) removed outlier: 3.546A pdb=" N TYR A 99 " --> pdb=" O TYR A 7 " (cutoff:3.500A) removed outlier: 6.875A pdb=" N THR A 94 " --> pdb=" O TYR A 118 " (cutoff:3.500A) removed outlier: 4.680A pdb=" N TYR A 118 " --> pdb=" O THR A 94 " (cutoff:3.500A) removed outlier: 7.018A pdb=" N GLN A 96 " --> pdb=" O TYR A 116 " (cutoff:3.500A) removed outlier: 5.057A pdb=" N TYR A 116 " --> pdb=" O GLN A 96 " (cutoff:3.500A) removed outlier: 6.835A pdb=" N MET A 98 " --> pdb=" O HIS A 114 " (cutoff:3.500A) removed outlier: 4.290A pdb=" N HIS A 114 " --> pdb=" O MET A 98 " (cutoff:3.500A) removed outlier: 6.924A pdb=" N GLY A 100 " --> pdb=" O GLY A 112 " (cutoff:3.500A) removed outlier: 4.871A pdb=" N GLY A 112 " --> pdb=" O GLY A 100 " (cutoff:3.500A) removed outlier: 6.761A pdb=" N ASP A 102 " --> pdb=" O LEU A 110 " (cutoff:3.500A) removed outlier: 7.153A pdb=" N HIS A 114 " --> pdb=" O ALA A 125 " (cutoff:3.500A) removed outlier: 4.644A pdb=" N ALA A 125 " --> pdb=" O HIS A 114 " (cutoff:3.500A) removed outlier: 7.094A pdb=" N TYR A 116 " --> pdb=" O TYR A 123 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 186 through 190 removed outlier: 4.072A pdb=" N ARG A 202 " --> pdb=" O THR A 190 " (cutoff:3.500A) removed outlier: 3.589A pdb=" N LEU A 201 " --> pdb=" O VAL A 247 " (cutoff:3.500A) removed outlier: 5.811A pdb=" N PHE A 241 " --> pdb=" O SER A 207 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 186 through 190 removed outlier: 4.072A pdb=" N ARG A 202 " --> pdb=" O THR A 190 " (cutoff:3.500A) removed outlier: 3.589A pdb=" N LEU A 201 " --> pdb=" O VAL A 247 " (cutoff:3.500A) removed outlier: 5.811A pdb=" N PHE A 241 " --> pdb=" O SER A 207 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'A' and resid 222 through 223 removed outlier: 3.801A pdb=" N GLU A 222 " --> pdb=" O ARG A 219 " (cutoff:3.500A) removed outlier: 3.568A pdb=" N THR A 214 " --> pdb=" O GLN A 262 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'B' and resid 6 through 11 Processing sheet with id=AA6, first strand: chain 'B' and resid 44 through 45 Processing sheet with id=AA7, first strand: chain 'B' and resid 55 through 56 removed outlier: 3.981A pdb=" N SER B 55 " --> pdb=" O TYR B 63 " (cutoff:3.500A) removed outlier: 3.557A pdb=" N TYR B 63 " --> pdb=" O SER B 55 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA7 134 hydrogen bonds defined for protein. 354 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 0 stacking parallelities Total time for adding SS restraints: 0.78 Time building geometry restraints manager: 1.24 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.23 - 1.34: 1049 1.34 - 1.46: 794 1.46 - 1.58: 1365 1.58 - 1.70: 0 1.70 - 1.81: 20 Bond restraints: 3228 Sorted by residual: bond pdb=" CB PRO B 32 " pdb=" CG PRO B 32 " ideal model delta sigma weight residual 1.506 1.468 0.038 3.90e-02 6.57e+02 9.41e-01 bond pdb=" CB PRO A 210 " pdb=" CG PRO A 210 " ideal model delta sigma weight residual 1.506 1.468 0.038 3.90e-02 6.57e+02 9.38e-01 bond pdb=" CB VAL C 2 " pdb=" CG2 VAL C 2 " ideal model delta sigma weight residual 1.521 1.490 0.031 3.30e-02 9.18e+02 8.64e-01 bond pdb=" CB ARG A 35 " pdb=" CG ARG A 35 " ideal model delta sigma weight residual 1.520 1.496 0.024 3.00e-02 1.11e+03 6.39e-01 bond pdb=" CG PRO A 210 " pdb=" CD PRO A 210 " ideal model delta sigma weight residual 1.512 1.491 0.021 2.70e-02 1.37e+03 5.95e-01 ... (remaining 3223 not shown) Histogram of bond angle deviations from ideal: 100.57 - 107.26: 133 107.26 - 113.95: 1722 113.95 - 120.63: 1219 120.63 - 127.32: 1243 127.32 - 134.01: 58 Bond angle restraints: 4375 Sorted by residual: angle pdb=" C CYS A 164 " pdb=" N VAL A 165 " pdb=" CA VAL A 165 " ideal model delta sigma weight residual 120.42 117.04 3.38 1.42e+00 4.96e-01 5.66e+00 angle pdb=" CA TYR A 159 " pdb=" CB TYR A 159 " pdb=" CG TYR A 159 " ideal model delta sigma weight residual 113.90 110.50 3.40 1.80e+00 3.09e-01 3.56e+00 angle pdb=" C GLN A 72 " pdb=" N THR A 73 " pdb=" CA THR A 73 " ideal model delta sigma weight residual 120.68 118.37 2.31 1.70e+00 3.46e-01 1.85e+00 angle pdb=" C THR A 163 " pdb=" N CYS A 164 " pdb=" CA CYS A 164 " ideal model delta sigma weight residual 120.29 118.39 1.90 1.42e+00 4.96e-01 1.79e+00 angle pdb=" CA CYS A 164 " pdb=" C CYS A 164 " pdb=" N VAL A 165 " ideal model delta sigma weight residual 117.39 118.94 -1.55 1.24e+00 6.50e-01 1.57e+00 ... (remaining 4370 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 15.45: 1642 15.45 - 30.89: 200 30.89 - 46.33: 45 46.33 - 61.78: 12 61.78 - 77.22: 5 Dihedral angle restraints: 1904 sinusoidal: 788 harmonic: 1116 Sorted by residual: dihedral pdb=" CA TYR A 27 " pdb=" C TYR A 27 " pdb=" N VAL A 28 " pdb=" CA VAL A 28 " ideal model delta harmonic sigma weight residual 180.00 162.65 17.35 0 5.00e+00 4.00e-02 1.20e+01 dihedral pdb=" CA ARG B 81 " pdb=" CB ARG B 81 " pdb=" CG ARG B 81 " pdb=" CD ARG B 81 " ideal model delta sinusoidal sigma weight residual -60.00 -111.49 51.49 3 1.50e+01 4.44e-03 9.02e+00 dihedral pdb=" CA ARG A 256 " pdb=" CB ARG A 256 " pdb=" CG ARG A 256 " pdb=" CD ARG A 256 " ideal model delta sinusoidal sigma weight residual -180.00 -129.33 -50.67 3 1.50e+01 4.44e-03 8.93e+00 ... (remaining 1901 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.028: 245 0.028 - 0.056: 129 0.056 - 0.084: 45 0.084 - 0.111: 17 0.111 - 0.139: 8 Chirality restraints: 444 Sorted by residual: chirality pdb=" CA VAL A 261 " pdb=" N VAL A 261 " pdb=" C VAL A 261 " pdb=" CB VAL A 261 " both_signs ideal model delta sigma weight residual False 2.44 2.58 -0.14 2.00e-01 2.50e+01 4.84e-01 chirality pdb=" CA ILE A 213 " pdb=" N ILE A 213 " pdb=" C ILE A 213 " pdb=" CB ILE A 213 " both_signs ideal model delta sigma weight residual False 2.43 2.57 -0.13 2.00e-01 2.50e+01 4.34e-01 chirality pdb=" CB VAL A 165 " pdb=" CA VAL A 165 " pdb=" CG1 VAL A 165 " pdb=" CG2 VAL A 165 " both_signs ideal model delta sigma weight residual False -2.63 -2.50 -0.13 2.00e-01 2.50e+01 4.32e-01 ... (remaining 441 not shown) Planarity restraints: 571 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA LEU A 160 " 0.014 2.00e-02 2.50e+03 2.85e-02 8.10e+00 pdb=" C LEU A 160 " -0.049 2.00e-02 2.50e+03 pdb=" O LEU A 160 " 0.018 2.00e-02 2.50e+03 pdb=" N GLU A 161 " 0.017 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA GLU A 161 " 0.011 2.00e-02 2.50e+03 2.15e-02 4.61e+00 pdb=" C GLU A 161 " -0.037 2.00e-02 2.50e+03 pdb=" O GLU A 161 " 0.014 2.00e-02 2.50e+03 pdb=" N GLY A 162 " 0.013 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C ALA A 49 " 0.016 5.00e-02 4.00e+02 2.42e-02 9.38e-01 pdb=" N PRO A 50 " -0.042 5.00e-02 4.00e+02 pdb=" CA PRO A 50 " 0.012 5.00e-02 4.00e+02 pdb=" CD PRO A 50 " 0.014 5.00e-02 4.00e+02 ... (remaining 568 not shown) Histogram of nonbonded interaction distances: 2.24 - 2.77: 588 2.77 - 3.31: 2685 3.31 - 3.84: 4624 3.84 - 4.37: 5574 4.37 - 4.90: 10145 Nonbonded interactions: 23616 Sorted by model distance: nonbonded pdb=" OH TYR A 59 " pdb=" OH TYR A 171 " model vdw 2.243 2.440 nonbonded pdb=" O GLN A 255 " pdb=" NH2 ARG A 273 " model vdw 2.351 2.520 nonbonded pdb=" O SER A 207 " pdb=" OG SER A 207 " model vdw 2.391 2.440 nonbonded pdb=" OE2 GLU B 36 " pdb=" NH1 ARG B 81 " model vdw 2.415 2.520 nonbonded pdb=" NH2 ARG A 14 " pdb=" O GLU A 19 " model vdw 2.459 2.520 ... (remaining 23611 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: 2.040 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.010 Construct map_model_manager: 0.010 Extract box with map and model: 9.230 Check model and map are aligned: 0.050 Set scattering table: 0.030 Process input model: 13.350 Find NCS groups from input model: 0.080 Set up NCS constraints: 0.010 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:2.780 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 27.590 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8135 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.038 3228 Z= 0.256 Angle : 0.493 4.865 4375 Z= 0.260 Chirality : 0.041 0.139 444 Planarity : 0.004 0.028 571 Dihedral : 14.826 77.222 1176 Min Nonbonded Distance : 2.243 Molprobity Statistics. All-atom Clashscore : 5.75 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.99 % Favored : 96.01 % Rotamer: Outliers : 0.30 % Allowed : 0.30 % Favored : 99.40 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 11.76 % 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.99 (0.43), residues: 376 helix: 1.18 (0.60), residues: 66 sheet: -0.23 (0.52), residues: 110 loop : -1.58 (0.42), residues: 200 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A 244 HIS 0.003 0.001 HIS B 84 PHE 0.014 0.002 PHE B 56 TYR 0.010 0.001 TYR A 159 ARG 0.005 0.000 ARG A 219 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 752 Ramachandran restraints generated. 376 Oldfield, 0 Emsley, 376 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 752 Ramachandran restraints generated. 376 Oldfield, 0 Emsley, 376 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 52 residues out of total 332 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 51 time to evaluate : 0.388 Fit side-chains REVERT: A 110 LEU cc_start: 0.7718 (tt) cc_final: 0.7508 (tt) REVERT: A 154 GLU cc_start: 0.7217 (mm-30) cc_final: 0.6901 (pt0) REVERT: B 45 ARG cc_start: 0.7243 (ttp80) cc_final: 0.5821 (mmt90) outliers start: 1 outliers final: 0 residues processed: 52 average time/residue: 0.2122 time to fit residues: 12.9866 Evaluate side-chains 47 residues out of total 332 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 47 time to evaluate : 0.337 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 37 random chunks: chunk 31 optimal weight: 1.9990 chunk 28 optimal weight: 1.9990 chunk 15 optimal weight: 0.5980 chunk 9 optimal weight: 0.5980 chunk 18 optimal weight: 0.9990 chunk 14 optimal weight: 1.9990 chunk 29 optimal weight: 0.7980 chunk 11 optimal weight: 0.9980 chunk 17 optimal weight: 0.4980 chunk 21 optimal weight: 1.9990 chunk 33 optimal weight: 0.5980 overall best weight: 0.6180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8119 moved from start: 0.0977 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.033 3228 Z= 0.232 Angle : 0.529 6.167 4375 Z= 0.280 Chirality : 0.042 0.152 444 Planarity : 0.004 0.026 571 Dihedral : 4.257 16.127 429 Min Nonbonded Distance : 2.252 Molprobity Statistics. All-atom Clashscore : 6.57 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.72 % Favored : 96.28 % Rotamer: Outliers : 0.60 % Allowed : 6.02 % Favored : 93.37 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 11.76 % 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.01 (0.43), residues: 376 helix: 0.87 (0.61), residues: 72 sheet: -0.73 (0.46), residues: 134 loop : -1.21 (0.45), residues: 170 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.015 0.001 TRP A 204 HIS 0.003 0.001 HIS B 84 PHE 0.014 0.002 PHE B 56 TYR 0.012 0.001 TYR A 159 ARG 0.004 0.000 ARG A 157 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 752 Ramachandran restraints generated. 376 Oldfield, 0 Emsley, 376 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 752 Ramachandran restraints generated. 376 Oldfield, 0 Emsley, 376 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 51 residues out of total 332 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 49 time to evaluate : 0.372 Fit side-chains REVERT: A 110 LEU cc_start: 0.7805 (tt) cc_final: 0.7513 (tt) REVERT: A 154 GLU cc_start: 0.7182 (mm-30) cc_final: 0.6974 (pt0) REVERT: A 157 ARG cc_start: 0.7693 (ttt-90) cc_final: 0.7447 (ptm-80) REVERT: A 159 TYR cc_start: 0.8178 (t80) cc_final: 0.7839 (t80) REVERT: B 45 ARG cc_start: 0.7416 (ttp80) cc_final: 0.5807 (mmt90) outliers start: 2 outliers final: 1 residues processed: 51 average time/residue: 0.1994 time to fit residues: 12.0883 Evaluate side-chains 48 residues out of total 332 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 47 time to evaluate : 0.369 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 34 VAL Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/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-5265/modules/cctbx_project/iotbx/cli_parser.py", line 944, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/phenix/phenix/programs/real_space_refine.py", line 189, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/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-5265/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 700, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 764, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/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-5265/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-5265/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-5265/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-5265/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-5265/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-5265/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/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-5265/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-5265/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1479, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1381, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/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: 67.6788 > 50: distance: 30 - 34: 17.886 distance: 34 - 35: 32.833 distance: 35 - 36: 22.749 distance: 35 - 38: 24.911 distance: 36 - 37: 48.855 distance: 36 - 43: 12.433 distance: 38 - 39: 25.245 distance: 39 - 40: 9.293 distance: 40 - 41: 7.815 distance: 40 - 42: 14.350 distance: 43 - 44: 30.917 distance: 44 - 45: 35.327 distance: 44 - 47: 34.847 distance: 45 - 46: 6.450 distance: 45 - 54: 50.030 distance: 47 - 48: 38.164 distance: 48 - 49: 18.627 distance: 49 - 50: 9.722 distance: 51 - 52: 5.011 distance: 51 - 53: 11.386 distance: 54 - 55: 36.161 distance: 55 - 56: 22.225 distance: 55 - 58: 18.085 distance: 56 - 57: 19.776 distance: 56 - 61: 24.920 distance: 58 - 59: 27.088 distance: 58 - 60: 24.628 distance: 61 - 62: 35.175 distance: 62 - 63: 54.505 distance: 62 - 65: 26.352 distance: 63 - 64: 67.110 distance: 63 - 69: 49.983 distance: 65 - 66: 28.838 distance: 66 - 67: 11.800 distance: 66 - 68: 20.088 distance: 69 - 70: 50.860 distance: 70 - 71: 46.956 distance: 70 - 73: 16.074 distance: 71 - 72: 38.574 distance: 71 - 74: 34.343 distance: 74 - 75: 19.591 distance: 74 - 80: 25.188 distance: 75 - 76: 14.786 distance: 75 - 78: 29.618 distance: 76 - 77: 16.642 distance: 76 - 81: 25.287 distance: 78 - 79: 33.900 distance: 79 - 80: 3.439 distance: 81 - 82: 17.024 distance: 81 - 151: 4.660 distance: 82 - 83: 26.005 distance: 82 - 85: 18.232 distance: 83 - 84: 10.846 distance: 83 - 90: 33.259 distance: 84 - 148: 22.094 distance: 85 - 86: 32.849 distance: 86 - 87: 10.889 distance: 87 - 88: 15.392 distance: 88 - 89: 12.435 distance: 90 - 91: 29.377 distance: 91 - 92: 51.830 distance: 91 - 94: 10.465 distance: 92 - 93: 47.491 distance: 92 - 97: 48.943 distance: 94 - 95: 40.548 distance: 94 - 96: 34.328 distance: 97 - 98: 22.554 distance: 97 - 137: 29.685 distance: 98 - 99: 26.707 distance: 98 - 101: 11.457 distance: 99 - 100: 24.756 distance: 99 - 107: 27.626 distance: 100 - 134: 24.089 distance: 101 - 102: 17.615 distance: 102 - 103: 23.716 distance: 102 - 104: 14.277 distance: 103 - 105: 13.272 distance: 104 - 106: 11.815 distance: 105 - 106: 12.927 distance: 107 - 108: 16.731 distance: 108 - 109: 33.697 distance: 108 - 111: 29.585 distance: 109 - 110: 24.543 distance: 109 - 115: 22.551 distance: 111 - 112: 12.848 distance: 112 - 113: 25.664 distance: 113 - 114: 6.796