Starting phenix.real_space_refine on Thu Jul 18 17:10:43 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m65_23689/07_2024/7m65_23689.cif Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m65_23689/07_2024/7m65_23689.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=4.1 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m65_23689/07_2024/7m65_23689.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m65_23689/07_2024/7m65_23689.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m65_23689/07_2024/7m65_23689.cif" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m65_23689/07_2024/7m65_23689.cif" } resolution = 4.1 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 ------------------------------------------------------------------------------- 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 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 5 5.16 5 C 1235 2.51 5 N 370 2.21 5 O 415 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped None Time to flip residues: 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/modules/chem_data/mon_lib" Total number of atoms: 2025 Number of models: 1 Model: "" Number of chains: 15 Chain: "A" Number of atoms: 237 Number of conformers: 1 Conformer: "" Number of residues, atoms: 32, 237 Classifications: {'peptide': 32} Link IDs: {'TRANS': 31} Chain: "A" Number of atoms: 1 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 1 Unusual residues: {'NH2': 1} Classifications: {'undetermined': 1} Chain: "B" Number of atoms: 167 Number of conformers: 1 Conformer: "" Number of residues, atoms: 23, 167 Classifications: {'peptide': 23} Link IDs: {'TRANS': 22} Chain: "C" Number of atoms: 237 Number of conformers: 1 Conformer: "" Number of residues, atoms: 32, 237 Classifications: {'peptide': 32} Link IDs: {'TRANS': 31} Chain: "C" Number of atoms: 1 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 1 Unusual residues: {'NH2': 1} Classifications: {'undetermined': 1} Chain: "D" Number of atoms: 167 Number of conformers: 1 Conformer: "" Number of residues, atoms: 23, 167 Classifications: {'peptide': 23} Link IDs: {'TRANS': 22} Chain: "E" Number of atoms: 237 Number of conformers: 1 Conformer: "" Number of residues, atoms: 32, 237 Classifications: {'peptide': 32} Link IDs: {'TRANS': 31} Chain: "E" Number of atoms: 1 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 1 Unusual residues: {'NH2': 1} Classifications: {'undetermined': 1} Chain: "F" Number of atoms: 167 Number of conformers: 1 Conformer: "" Number of residues, atoms: 23, 167 Classifications: {'peptide': 23} Link IDs: {'TRANS': 22} Chain: "G" Number of atoms: 237 Number of conformers: 1 Conformer: "" Number of residues, atoms: 32, 237 Classifications: {'peptide': 32} Link IDs: {'TRANS': 31} Chain: "G" Number of atoms: 1 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 1 Unusual residues: {'NH2': 1} Classifications: {'undetermined': 1} Chain: "H" Number of atoms: 167 Number of conformers: 1 Conformer: "" Number of residues, atoms: 23, 167 Classifications: {'peptide': 23} Link IDs: {'TRANS': 22} Chain: "I" Number of atoms: 237 Number of conformers: 1 Conformer: "" Number of residues, atoms: 32, 237 Classifications: {'peptide': 32} Link IDs: {'TRANS': 31} Chain: "I" Number of atoms: 1 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 1 Unusual residues: {'NH2': 1} Classifications: {'undetermined': 1} Chain: "J" Number of atoms: 167 Number of conformers: 1 Conformer: "" Number of residues, atoms: 23, 167 Classifications: {'peptide': 23} Link IDs: {'TRANS': 22} Time building chain proxies: 1.46, per 1000 atoms: 0.72 Number of scatterers: 2025 At special positions: 0 Unit cell: (62.245, 68.575, 45.365, 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 5 16.00 O 415 8.00 N 370 7.00 C 1235 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=5, symmetry=0 Number of additional bonds: simple=5, symmetry=0 Coordination: Other bonds: Time building additional restraints: 2.05 Conformation dependent library (CDL) restraints added in 368.6 milliseconds 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 510 Finding SS restraints... Secondary structure from input PDB file: 0 helices and 1 sheets defined 0.0% alpha, 5.5% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.14 Creating SS restraints... Processing sheet with id=AA1, first strand: chain 'C' and resid 7 through 10 removed outlier: 3.604A pdb=" N THR E 9 " --> pdb=" O ALA C 8 " (cutoff:3.500A) removed outlier: 3.604A pdb=" N THR G 9 " --> pdb=" O ALA E 8 " (cutoff:3.500A) removed outlier: 5.879A pdb=" N CYS A 7 " --> pdb=" O ALA G 8 " (cutoff:3.500A) removed outlier: 8.182A pdb=" N GLN G 10 " --> pdb=" O CYS A 7 " (cutoff:3.500A) removed outlier: 6.768A pdb=" N THR A 9 " --> pdb=" O GLN G 10 " (cutoff:3.500A) removed outlier: 3.604A pdb=" N THR I 9 " --> pdb=" O ALA A 8 " (cutoff:3.500A) 6 hydrogen bonds defined for protein. 18 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.27 Time building geometry restraints manager: 1.05 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.19 - 1.31: 337 1.31 - 1.44: 613 1.44 - 1.56: 1090 1.56 - 1.68: 0 1.68 - 1.81: 5 Bond restraints: 2045 Sorted by residual: bond pdb=" CG ASN C 35 " pdb=" OD1 ASN C 35 " ideal model delta sigma weight residual 1.231 1.189 0.042 1.90e-02 2.77e+03 4.93e+00 bond pdb=" CG ASN E 35 " pdb=" OD1 ASN E 35 " ideal model delta sigma weight residual 1.231 1.189 0.042 1.90e-02 2.77e+03 4.78e+00 bond pdb=" CG ASN I 35 " pdb=" OD1 ASN I 35 " ideal model delta sigma weight residual 1.231 1.190 0.041 1.90e-02 2.77e+03 4.74e+00 bond pdb=" CG ASN G 35 " pdb=" OD1 ASN G 35 " ideal model delta sigma weight residual 1.231 1.190 0.041 1.90e-02 2.77e+03 4.59e+00 bond pdb=" CG ASN A 35 " pdb=" OD1 ASN A 35 " ideal model delta sigma weight residual 1.231 1.190 0.041 1.90e-02 2.77e+03 4.55e+00 ... (remaining 2040 not shown) Histogram of bond angle deviations from ideal: 103.62 - 109.16: 224 109.16 - 114.71: 1001 114.71 - 120.26: 585 120.26 - 125.81: 955 125.81 - 131.36: 10 Bond angle restraints: 2775 Sorted by residual: angle pdb=" N PHE F 15 " pdb=" CA PHE F 15 " pdb=" C PHE F 15 " ideal model delta sigma weight residual 112.41 119.39 -6.98 1.30e+00 5.92e-01 2.88e+01 angle pdb=" N PHE J 15 " pdb=" CA PHE J 15 " pdb=" C PHE J 15 " ideal model delta sigma weight residual 112.41 119.38 -6.97 1.30e+00 5.92e-01 2.87e+01 angle pdb=" N PHE B 15 " pdb=" CA PHE B 15 " pdb=" C PHE B 15 " ideal model delta sigma weight residual 112.41 119.37 -6.96 1.30e+00 5.92e-01 2.87e+01 angle pdb=" N PHE D 15 " pdb=" CA PHE D 15 " pdb=" C PHE D 15 " ideal model delta sigma weight residual 112.41 119.36 -6.95 1.30e+00 5.92e-01 2.86e+01 angle pdb=" N PHE H 15 " pdb=" CA PHE H 15 " pdb=" C PHE H 15 " ideal model delta sigma weight residual 112.41 119.34 -6.93 1.30e+00 5.92e-01 2.84e+01 ... (remaining 2770 not shown) Histogram of dihedral angle deviations from ideal: 0.01 - 8.26: 870 8.26 - 16.50: 195 16.50 - 24.75: 55 24.75 - 33.00: 35 33.00 - 41.24: 10 Dihedral angle restraints: 1165 sinusoidal: 385 harmonic: 780 Sorted by residual: dihedral pdb=" CA PHE H 15 " pdb=" C PHE H 15 " pdb=" N LEU H 16 " pdb=" CA LEU H 16 " ideal model delta harmonic sigma weight residual 180.00 152.44 27.56 0 5.00e+00 4.00e-02 3.04e+01 dihedral pdb=" CA PHE J 15 " pdb=" C PHE J 15 " pdb=" N LEU J 16 " pdb=" CA LEU J 16 " ideal model delta harmonic sigma weight residual 180.00 152.47 27.53 0 5.00e+00 4.00e-02 3.03e+01 dihedral pdb=" CA PHE D 15 " pdb=" C PHE D 15 " pdb=" N LEU D 16 " pdb=" CA LEU D 16 " ideal model delta harmonic sigma weight residual 180.00 152.49 27.51 0 5.00e+00 4.00e-02 3.03e+01 ... (remaining 1162 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.132: 312 0.132 - 0.264: 23 0.264 - 0.396: 0 0.396 - 0.528: 0 0.528 - 0.660: 5 Chirality restraints: 340 Sorted by residual: chirality pdb=" CB ILE B 26 " pdb=" CA ILE B 26 " pdb=" CG1 ILE B 26 " pdb=" CG2 ILE B 26 " both_signs ideal model delta sigma weight residual False 2.64 1.99 0.66 2.00e-01 2.50e+01 1.09e+01 chirality pdb=" CB ILE J 26 " pdb=" CA ILE J 26 " pdb=" CG1 ILE J 26 " pdb=" CG2 ILE J 26 " both_signs ideal model delta sigma weight residual False 2.64 1.99 0.66 2.00e-01 2.50e+01 1.09e+01 chirality pdb=" CB ILE F 26 " pdb=" CA ILE F 26 " pdb=" CG1 ILE F 26 " pdb=" CG2 ILE F 26 " both_signs ideal model delta sigma weight residual False 2.64 1.99 0.66 2.00e-01 2.50e+01 1.08e+01 ... (remaining 337 not shown) Planarity restraints: 360 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB PHE E 15 " -0.015 2.00e-02 2.50e+03 1.24e-02 2.70e+00 pdb=" CG PHE E 15 " 0.028 2.00e-02 2.50e+03 pdb=" CD1 PHE E 15 " 0.001 2.00e-02 2.50e+03 pdb=" CD2 PHE E 15 " -0.005 2.00e-02 2.50e+03 pdb=" CE1 PHE E 15 " -0.006 2.00e-02 2.50e+03 pdb=" CE2 PHE E 15 " 0.000 2.00e-02 2.50e+03 pdb=" CZ PHE E 15 " -0.003 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE I 15 " -0.015 2.00e-02 2.50e+03 1.24e-02 2.69e+00 pdb=" CG PHE I 15 " 0.028 2.00e-02 2.50e+03 pdb=" CD1 PHE I 15 " 0.001 2.00e-02 2.50e+03 pdb=" CD2 PHE I 15 " -0.005 2.00e-02 2.50e+03 pdb=" CE1 PHE I 15 " -0.006 2.00e-02 2.50e+03 pdb=" CE2 PHE I 15 " 0.000 2.00e-02 2.50e+03 pdb=" CZ PHE I 15 " -0.003 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE A 15 " -0.015 2.00e-02 2.50e+03 1.24e-02 2.68e+00 pdb=" CG PHE A 15 " 0.028 2.00e-02 2.50e+03 pdb=" CD1 PHE A 15 " 0.001 2.00e-02 2.50e+03 pdb=" CD2 PHE A 15 " -0.005 2.00e-02 2.50e+03 pdb=" CE1 PHE A 15 " -0.006 2.00e-02 2.50e+03 pdb=" CE2 PHE A 15 " 0.000 2.00e-02 2.50e+03 pdb=" CZ PHE A 15 " -0.003 2.00e-02 2.50e+03 ... (remaining 357 not shown) Histogram of nonbonded interaction distances: 2.39 - 2.89: 936 2.89 - 3.40: 2346 3.40 - 3.90: 3188 3.90 - 4.40: 3720 4.40 - 4.90: 6916 Nonbonded interactions: 17106 Sorted by model distance: nonbonded pdb=" O LEU J 16 " pdb=" CD2 LEU J 16 " model vdw 2.392 3.460 nonbonded pdb=" O LEU D 16 " pdb=" CD2 LEU D 16 " model vdw 2.392 3.460 nonbonded pdb=" O LEU B 16 " pdb=" CD2 LEU B 16 " model vdw 2.392 3.460 nonbonded pdb=" O LEU F 16 " pdb=" CD2 LEU F 16 " model vdw 2.392 3.460 nonbonded pdb=" O LEU H 16 " pdb=" CD2 LEU H 16 " model vdw 2.393 3.460 ... (remaining 17101 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = chain 'A' selection = chain 'C' selection = chain 'E' selection = chain 'G' selection = chain 'I' } ncs_group { reference = chain 'B' selection = chain 'D' selection = chain 'F' selection = chain 'H' selection = chain 'J' } 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 group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.130 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.000 Extract box with map and model: 0.100 Check model and map are aligned: 0.010 Set scattering table: 0.030 Process input model: 10.590 Find NCS groups from input model: 0.070 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.370 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 13.310 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7384 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.009 0.063 2045 Z= 0.595 Angle : 1.401 11.663 2775 Z= 0.744 Chirality : 0.102 0.660 340 Planarity : 0.004 0.012 360 Dihedral : 12.046 41.242 655 Min Nonbonded Distance : 2.392 Molprobity Statistics. All-atom Clashscore : 33.63 Ramachandran Plot: Outliers : 0.00 % Allowed : 17.65 % Favored : 82.35 % Rotamer: Outliers : 0.00 % Allowed : 4.26 % Favored : 95.74 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -6.09 (0.25), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.63 (0.19), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.006 0.004 HIS H 18 PHE 0.028 0.005 PHE E 15 TYR 0.008 0.002 TYR E 37 ARG 0.001 0.000 ARG G 11 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 54 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 54 time to evaluate : 0.623 Fit side-chains REVERT: C 15 PHE cc_start: 0.7309 (p90) cc_final: 0.5719 (t80) REVERT: C 28 SER cc_start: 0.8833 (t) cc_final: 0.8519 (p) REVERT: E 28 SER cc_start: 0.8926 (t) cc_final: 0.8483 (p) REVERT: G 28 SER cc_start: 0.9166 (t) cc_final: 0.8770 (p) outliers start: 0 outliers final: 0 residues processed: 54 average time/residue: 0.1427 time to fit residues: 8.9681 Evaluate side-chains 34 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 34 time to evaluate : 0.260 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 25 random chunks: chunk 21 optimal weight: 8.9990 chunk 18 optimal weight: 20.0000 chunk 10 optimal weight: 10.0000 chunk 6 optimal weight: 10.0000 chunk 12 optimal weight: 10.0000 chunk 19 optimal weight: 1.9990 chunk 7 optimal weight: 4.9990 chunk 11 optimal weight: 7.9990 chunk 14 optimal weight: 3.9990 chunk 22 optimal weight: 10.0000 chunk 15 optimal weight: 7.9990 overall best weight: 5.3990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 21 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 21 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 21 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** G 21 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7324 moved from start: 0.3154 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 2045 Z= 0.223 Angle : 0.835 12.079 2775 Z= 0.387 Chirality : 0.044 0.134 340 Planarity : 0.003 0.015 360 Dihedral : 7.480 23.495 275 Min Nonbonded Distance : 2.255 Molprobity Statistics. All-atom Clashscore : 23.69 Ramachandran Plot: Outliers : 0.00 % Allowed : 17.25 % Favored : 82.75 % Rotamer: Outliers : 0.00 % Allowed : 9.36 % Favored : 90.64 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -5.83 (0.28), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.44 (0.22), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.004 0.002 HIS I 18 PHE 0.012 0.002 PHE J 15 TYR 0.006 0.001 TYR G 37 ARG 0.001 0.001 ARG E 11 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 39 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 39 time to evaluate : 0.243 Fit side-chains REVERT: C 28 SER cc_start: 0.8726 (t) cc_final: 0.8494 (p) REVERT: E 27 LEU cc_start: 0.8867 (mp) cc_final: 0.8613 (mp) REVERT: G 28 SER cc_start: 0.9103 (t) cc_final: 0.8700 (p) REVERT: J 22 ASN cc_start: 0.7821 (p0) cc_final: 0.7182 (t0) outliers start: 0 outliers final: 0 residues processed: 39 average time/residue: 0.0826 time to fit residues: 4.2216 Evaluate side-chains 30 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 30 time to evaluate : 0.256 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 25 random chunks: chunk 6 optimal weight: 10.0000 chunk 22 optimal weight: 20.0000 chunk 24 optimal weight: 30.0000 chunk 20 optimal weight: 10.0000 chunk 7 optimal weight: 5.9990 chunk 18 optimal weight: 40.0000 chunk 17 optimal weight: 20.0000 chunk 11 optimal weight: 20.0000 chunk 2 optimal weight: 8.9990 chunk 10 optimal weight: 10.0000 chunk 15 optimal weight: 4.9990 overall best weight: 7.9994 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 21 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** C 35 ASN I 35 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7322 moved from start: 0.3856 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.027 2045 Z= 0.251 Angle : 0.752 9.739 2775 Z= 0.371 Chirality : 0.043 0.139 340 Planarity : 0.003 0.021 360 Dihedral : 7.286 22.752 275 Min Nonbonded Distance : 2.225 Molprobity Statistics. All-atom Clashscore : 21.66 Ramachandran Plot: Outliers : 0.00 % Allowed : 18.04 % Favored : 81.96 % Rotamer: Outliers : 0.00 % Allowed : 7.23 % Favored : 92.77 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -5.91 (0.29), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.50 (0.22), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.005 0.002 HIS B 18 PHE 0.013 0.003 PHE J 23 TYR 0.005 0.001 TYR E 37 ARG 0.003 0.000 ARG I 11 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 37 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 37 time to evaluate : 0.243 Fit side-chains REVERT: G 28 SER cc_start: 0.9138 (t) cc_final: 0.8796 (p) REVERT: J 22 ASN cc_start: 0.7897 (p0) cc_final: 0.7112 (t0) outliers start: 0 outliers final: 0 residues processed: 37 average time/residue: 0.0770 time to fit residues: 3.7238 Evaluate side-chains 25 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 25 time to evaluate : 0.186 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 25 random chunks: chunk 22 optimal weight: 20.0000 chunk 24 optimal weight: 20.0000 chunk 11 optimal weight: 8.9990 chunk 21 optimal weight: 8.9990 chunk 6 optimal weight: 9.9990 chunk 20 optimal weight: 6.9990 chunk 13 optimal weight: 50.0000 chunk 0 optimal weight: 30.0000 chunk 17 optimal weight: 9.9990 chunk 9 optimal weight: 2.9990 chunk 16 optimal weight: 7.9990 overall best weight: 7.1990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 21 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7298 moved from start: 0.4338 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.026 2045 Z= 0.223 Angle : 0.779 10.377 2775 Z= 0.363 Chirality : 0.042 0.140 340 Planarity : 0.003 0.018 360 Dihedral : 7.084 22.068 275 Min Nonbonded Distance : 2.190 Molprobity Statistics. All-atom Clashscore : 20.13 Ramachandran Plot: Outliers : 0.00 % Allowed : 17.65 % Favored : 82.35 % Rotamer: Outliers : 0.00 % Allowed : 3.40 % Favored : 96.60 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -6.02 (0.28), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.58 (0.22), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.004 0.002 HIS B 18 PHE 0.011 0.003 PHE G 15 TYR 0.003 0.001 TYR G 37 ARG 0.002 0.000 ARG I 11 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 35 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 35 time to evaluate : 0.203 Fit side-chains REVERT: E 10 GLN cc_start: 0.7932 (mp10) cc_final: 0.7545 (mp10) REVERT: J 22 ASN cc_start: 0.7862 (p0) cc_final: 0.6951 (t0) outliers start: 0 outliers final: 0 residues processed: 35 average time/residue: 0.0716 time to fit residues: 3.2876 Evaluate side-chains 26 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 26 time to evaluate : 0.200 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 25 random chunks: chunk 0 optimal weight: 50.0000 chunk 12 optimal weight: 40.0000 chunk 21 optimal weight: 10.0000 chunk 6 optimal weight: 7.9990 chunk 8 optimal weight: 1.9990 chunk 4 optimal weight: 3.9990 chunk 14 optimal weight: 40.0000 chunk 5 optimal weight: 20.0000 chunk 24 optimal weight: 0.9990 chunk 20 optimal weight: 5.9990 chunk 11 optimal weight: 7.9990 overall best weight: 4.1990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 18 HIS A 35 ASN I 18 HIS Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7215 moved from start: 0.5027 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.022 2045 Z= 0.162 Angle : 0.682 10.724 2775 Z= 0.323 Chirality : 0.041 0.138 340 Planarity : 0.002 0.013 360 Dihedral : 6.267 20.942 275 Min Nonbonded Distance : 2.178 Molprobity Statistics. All-atom Clashscore : 18.34 Ramachandran Plot: Outliers : 0.00 % Allowed : 14.51 % Favored : 85.49 % Rotamer: Outliers : 0.00 % Allowed : 5.96 % Favored : 94.04 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -5.93 (0.28), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.51 (0.22), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.001 HIS B 18 PHE 0.011 0.002 PHE E 15 TYR 0.003 0.001 TYR E 37 ARG 0.001 0.000 ARG I 11 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 41 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 41 time to evaluate : 0.251 Fit side-chains REVERT: E 10 GLN cc_start: 0.8034 (mp10) cc_final: 0.7698 (mp10) REVERT: J 22 ASN cc_start: 0.7499 (p0) cc_final: 0.6925 (t0) outliers start: 0 outliers final: 0 residues processed: 41 average time/residue: 0.0851 time to fit residues: 4.5021 Evaluate side-chains 26 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 26 time to evaluate : 0.554 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 25 random chunks: chunk 2 optimal weight: 50.0000 chunk 8 optimal weight: 0.8980 chunk 12 optimal weight: 7.9990 chunk 23 optimal weight: 40.0000 chunk 13 optimal weight: 8.9990 chunk 17 optimal weight: 8.9990 chunk 20 optimal weight: 5.9990 chunk 24 optimal weight: 6.9990 chunk 15 optimal weight: 5.9990 chunk 14 optimal weight: 8.9990 chunk 11 optimal weight: 1.9990 overall best weight: 4.3788 ------------------------------------------------------------------------------- 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.7222 moved from start: 0.5356 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.024 2045 Z= 0.173 Angle : 0.743 10.541 2775 Z= 0.339 Chirality : 0.043 0.140 340 Planarity : 0.002 0.011 360 Dihedral : 6.208 20.930 275 Min Nonbonded Distance : 2.191 Molprobity Statistics. All-atom Clashscore : 16.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 15.29 % Favored : 84.71 % Rotamer: Outliers : 0.00 % Allowed : 3.83 % Favored : 96.17 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -5.93 (0.28), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.51 (0.21), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.004 0.002 HIS B 18 PHE 0.012 0.002 PHE G 15 TYR 0.003 0.001 TYR E 37 ARG 0.002 0.000 ARG I 11 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 35 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 35 time to evaluate : 0.194 Fit side-chains REVERT: E 10 GLN cc_start: 0.7952 (mp10) cc_final: 0.7740 (mp10) REVERT: J 22 ASN cc_start: 0.7494 (p0) cc_final: 0.6699 (t0) outliers start: 0 outliers final: 0 residues processed: 35 average time/residue: 0.0537 time to fit residues: 2.6229 Evaluate side-chains 25 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 25 time to evaluate : 0.187 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 25 random chunks: chunk 14 optimal weight: 30.0000 chunk 9 optimal weight: 20.0000 chunk 7 optimal weight: 6.9990 chunk 4 optimal weight: 30.0000 chunk 15 optimal weight: 40.0000 chunk 16 optimal weight: 30.0000 chunk 11 optimal weight: 20.0000 chunk 2 optimal weight: 30.0000 chunk 18 optimal weight: 30.0000 chunk 21 optimal weight: 3.9990 chunk 23 optimal weight: 30.0000 overall best weight: 16.1996 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 22 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7406 moved from start: 0.5173 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.037 2045 Z= 0.434 Angle : 0.910 10.558 2775 Z= 0.450 Chirality : 0.048 0.192 340 Planarity : 0.004 0.036 360 Dihedral : 8.121 22.051 275 Min Nonbonded Distance : 2.139 Molprobity Statistics. All-atom Clashscore : 30.06 Ramachandran Plot: Outliers : 0.00 % Allowed : 21.96 % Favored : 78.04 % Rotamer: Outliers : 0.00 % Allowed : 2.13 % Favored : 97.87 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -6.33 (0.27), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.82 (0.20), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.006 0.003 HIS B 18 PHE 0.019 0.004 PHE F 15 TYR 0.004 0.001 TYR E 37 ARG 0.007 0.001 ARG I 11 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 27 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 27 time to evaluate : 0.393 Fit side-chains REVERT: E 10 GLN cc_start: 0.8316 (mp10) cc_final: 0.7884 (mp10) outliers start: 0 outliers final: 0 residues processed: 27 average time/residue: 0.0521 time to fit residues: 2.2317 Evaluate side-chains 22 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 22 time to evaluate : 0.220 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 25 random chunks: chunk 21 optimal weight: 5.9990 chunk 22 optimal weight: 10.0000 chunk 23 optimal weight: 20.0000 chunk 13 optimal weight: 40.0000 chunk 9 optimal weight: 4.9990 chunk 17 optimal weight: 20.0000 chunk 6 optimal weight: 7.9990 chunk 20 optimal weight: 10.0000 chunk 14 optimal weight: 7.9990 chunk 11 optimal weight: 7.9990 chunk 16 optimal weight: 9.9990 overall best weight: 6.9990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: E 18 HIS G 18 HIS Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7289 moved from start: 0.5570 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.025 2045 Z= 0.224 Angle : 0.761 10.753 2775 Z= 0.358 Chirality : 0.043 0.145 340 Planarity : 0.003 0.017 360 Dihedral : 7.139 22.492 275 Min Nonbonded Distance : 2.137 Molprobity Statistics. All-atom Clashscore : 19.11 Ramachandran Plot: Outliers : 0.00 % Allowed : 13.73 % Favored : 86.27 % Rotamer: Outliers : 0.00 % Allowed : 1.70 % Favored : 98.30 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -6.18 (0.27), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.71 (0.21), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.004 0.002 HIS H 18 PHE 0.012 0.002 PHE H 15 TYR 0.004 0.001 TYR E 37 ARG 0.004 0.001 ARG I 11 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 31 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 31 time to evaluate : 0.260 Fit side-chains REVERT: E 10 GLN cc_start: 0.8206 (mp10) cc_final: 0.7805 (mp10) outliers start: 0 outliers final: 0 residues processed: 31 average time/residue: 0.0665 time to fit residues: 2.9792 Evaluate side-chains 25 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 25 time to evaluate : 0.247 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 25 random chunks: chunk 24 optimal weight: 30.0000 chunk 22 optimal weight: 30.0000 chunk 19 optimal weight: 20.0000 chunk 2 optimal weight: 3.9990 chunk 15 optimal weight: 9.9990 chunk 12 optimal weight: 7.9990 chunk 21 optimal weight: 10.0000 chunk 6 optimal weight: 0.0060 chunk 18 optimal weight: 30.0000 chunk 5 optimal weight: 9.9990 chunk 8 optimal weight: 0.0050 overall best weight: 4.4016 ------------------------------------------------------------------------------- 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.7197 moved from start: 0.5996 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.021 2045 Z= 0.174 Angle : 0.739 11.584 2775 Z= 0.337 Chirality : 0.043 0.140 340 Planarity : 0.002 0.012 360 Dihedral : 6.377 20.848 275 Min Nonbonded Distance : 2.231 Molprobity Statistics. All-atom Clashscore : 15.80 Ramachandran Plot: Outliers : 0.00 % Allowed : 15.29 % Favored : 84.71 % Rotamer: Outliers : 0.00 % Allowed : 0.85 % Favored : 99.15 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -6.09 (0.27), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.63 (0.21), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.004 0.002 HIS G 18 PHE 0.013 0.002 PHE E 15 TYR 0.003 0.000 TYR E 37 ARG 0.002 0.000 ARG I 11 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 34 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 34 time to evaluate : 0.247 Fit side-chains REVERT: E 10 GLN cc_start: 0.8064 (mp10) cc_final: 0.7755 (mp10) REVERT: E 32 VAL cc_start: 0.8376 (p) cc_final: 0.8163 (p) outliers start: 0 outliers final: 0 residues processed: 34 average time/residue: 0.0565 time to fit residues: 2.7373 Evaluate side-chains 27 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 27 time to evaluate : 0.198 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 25 random chunks: chunk 20 optimal weight: 6.9990 chunk 2 optimal weight: 30.0000 chunk 3 optimal weight: 0.1980 chunk 17 optimal weight: 20.0000 chunk 1 optimal weight: 7.9990 chunk 14 optimal weight: 1.9990 chunk 22 optimal weight: 8.9990 chunk 13 optimal weight: 9.9990 chunk 0 optimal weight: 5.9990 chunk 15 optimal weight: 0.6980 chunk 9 optimal weight: 9.9990 overall best weight: 3.1786 ------------------------------------------------------------------------------- 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.7149 moved from start: 0.6273 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.018 2045 Z= 0.151 Angle : 0.718 10.701 2775 Z= 0.326 Chirality : 0.043 0.133 340 Planarity : 0.002 0.011 360 Dihedral : 6.026 20.731 275 Min Nonbonded Distance : 2.212 Molprobity Statistics. All-atom Clashscore : 15.03 Ramachandran Plot: Outliers : 0.00 % Allowed : 11.37 % Favored : 88.63 % Rotamer: Outliers : 0.00 % Allowed : 0.85 % Favored : 99.15 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -5.94 (0.27), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.52 (0.21), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.001 HIS B 18 PHE 0.016 0.002 PHE G 15 TYR 0.002 0.000 TYR E 37 ARG 0.002 0.000 ARG I 11 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 510 Ramachandran restraints generated. 255 Oldfield, 0 Emsley, 255 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 34 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 34 time to evaluate : 0.198 Fit side-chains REVERT: E 10 GLN cc_start: 0.8073 (mp10) cc_final: 0.7866 (mp10) outliers start: 0 outliers final: 0 residues processed: 34 average time/residue: 0.0647 time to fit residues: 3.0569 Evaluate side-chains 27 residues out of total 235 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 27 time to evaluate : 0.248 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 25 random chunks: chunk 9 optimal weight: 0.0470 chunk 24 optimal weight: 30.0000 chunk 0 optimal weight: 50.0000 chunk 4 optimal weight: 9.9990 chunk 3 optimal weight: 4.9990 chunk 5 optimal weight: 40.0000 chunk 20 optimal weight: 5.9990 chunk 23 optimal weight: 6.9990 chunk 10 optimal weight: 10.0000 chunk 2 optimal weight: 7.9990 chunk 6 optimal weight: 20.0000 overall best weight: 5.2086 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3920 r_free = 0.3920 target = 0.098460 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 34)----------------| | r_work = 0.3839 r_free = 0.3839 target = 0.090199 restraints weight = 9882.376| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 29)----------------| | r_work = 0.3876 r_free = 0.3876 target = 0.091969 restraints weight = 5106.860| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 33)----------------| | r_work = 0.3902 r_free = 0.3902 target = 0.093126 restraints weight = 3089.370| |-----------------------------------------------------------------------------| r_work (final): 0.3859 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7261 moved from start: 0.6275 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.020 2045 Z= 0.183 Angle : 0.736 10.003 2775 Z= 0.336 Chirality : 0.043 0.144 340 Planarity : 0.002 0.011 360 Dihedral : 6.159 20.087 275 Min Nonbonded Distance : 2.270 Molprobity Statistics. All-atom Clashscore : 17.83 Ramachandran Plot: Outliers : 0.00 % Allowed : 14.51 % Favored : 85.49 % Rotamer: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -6.00 (0.27), residues: 255 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -4.56 (0.20), residues: 255 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.004 0.001 HIS B 18 PHE 0.013 0.002 PHE G 15 TYR 0.002 0.001 TYR E 37 ARG 0.002 0.000 ARG I 11 =============================================================================== Job complete usr+sys time: 671.02 seconds wall clock time: 19 minutes 8.91 seconds (1148.91 seconds total)