Starting phenix.real_space_refine on Thu Dec 7 16:35:22 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xo3_33335/12_2023/7xo3_33335.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xo3_33335/12_2023/7xo3_33335.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/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xo3_33335/12_2023/7xo3_33335.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xo3_33335/12_2023/7xo3_33335.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xo3_33335/12_2023/7xo3_33335.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xo3_33335/12_2023/7xo3_33335.pdb" } 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.001 sd= 0.004 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: 3 Type Number sf(0) Gaussians C 1626 2.51 5 N 450 2.21 5 O 522 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "B PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "E PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "F PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/chem_data/mon_lib" Total number of atoms: 2598 Number of models: 1 Model: "" Number of chains: 6 Chain: "B" Number of atoms: 433 Number of conformers: 1 Conformer: "" Number of residues, atoms: 63, 433 Classifications: {'peptide': 63} Link IDs: {'TRANS': 62} Chain: "A" Number of atoms: 433 Number of conformers: 1 Conformer: "" Number of residues, atoms: 63, 433 Classifications: {'peptide': 63} Link IDs: {'TRANS': 62} Chain: "E" Number of atoms: 433 Number of conformers: 1 Conformer: "" Number of residues, atoms: 63, 433 Classifications: {'peptide': 63} Link IDs: {'TRANS': 62} Chain: "C" Number of atoms: 433 Number of conformers: 1 Conformer: "" Number of residues, atoms: 63, 433 Classifications: {'peptide': 63} Link IDs: {'TRANS': 62} Chain: "F" Number of atoms: 433 Number of conformers: 1 Conformer: "" Number of residues, atoms: 63, 433 Classifications: {'peptide': 63} Link IDs: {'TRANS': 62} Chain: "D" Number of atoms: 433 Number of conformers: 1 Conformer: "" Number of residues, atoms: 63, 433 Classifications: {'peptide': 63} Link IDs: {'TRANS': 62} Time building chain proxies: 1.79, per 1000 atoms: 0.69 Number of scatterers: 2598 At special positions: 0 Unit cell: (79.5, 117.66, 32.86, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 3 Type Number sf(0) O 522 8.00 N 450 7.00 C 1626 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 0.98 Conformation dependent library (CDL) restraints added in 524.9 milliseconds 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 636 Finding SS restraints... Secondary structure from input PDB file: 0 helices and 14 sheets defined 0.0% alpha, 42.9% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.38 Creating SS restraints... Processing sheet with id=AA1, first strand: chain 'E' and resid 37 through 42 removed outlier: 6.386A pdb=" N LEU B 38 " --> pdb=" O TYR E 39 " (cutoff:3.500A) removed outlier: 7.745A pdb=" N GLY E 41 " --> pdb=" O LEU B 38 " (cutoff:3.500A) removed outlier: 6.479A pdb=" N VAL B 40 " --> pdb=" O GLY E 41 " (cutoff:3.500A) removed outlier: 6.375A pdb=" N VAL B 37 " --> pdb=" O LEU F 38 " (cutoff:3.500A) removed outlier: 7.690A pdb=" N VAL F 40 " --> pdb=" O VAL B 37 " (cutoff:3.500A) removed outlier: 6.403A pdb=" N TYR B 39 " --> pdb=" O VAL F 40 " (cutoff:3.500A) removed outlier: 7.790A pdb=" N SER F 42 " --> pdb=" O TYR B 39 " (cutoff:3.500A) removed outlier: 6.477A pdb=" N GLY B 41 " --> pdb=" O SER F 42 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA1 Processing sheet with id=AA2, first strand: chain 'E' and resid 48 through 57 removed outlier: 6.293A pdb=" N VAL B 49 " --> pdb=" O HIS E 50 " (cutoff:3.500A) removed outlier: 8.236A pdb=" N VAL E 52 " --> pdb=" O VAL B 49 " (cutoff:3.500A) removed outlier: 7.702A pdb=" N GLY B 51 " --> pdb=" O VAL E 52 " (cutoff:3.500A) removed outlier: 9.531A pdb=" N THR E 54 " --> pdb=" O GLY B 51 " (cutoff:3.500A) removed outlier: 8.686A pdb=" N ALA B 53 " --> pdb=" O THR E 54 " (cutoff:3.500A) removed outlier: 10.093A pdb=" N ALA E 56 " --> pdb=" O ALA B 53 " (cutoff:3.500A) removed outlier: 8.846A pdb=" N VAL B 55 " --> pdb=" O ALA E 56 " (cutoff:3.500A) removed outlier: 6.989A pdb=" N VAL B 49 " --> pdb=" O VAL F 48 " (cutoff:3.500A) removed outlier: 7.012A pdb=" N HIS F 50 " --> pdb=" O VAL B 49 " (cutoff:3.500A) removed outlier: 6.908A pdb=" N GLY B 51 " --> pdb=" O HIS F 50 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'E' and resid 61 through 67 removed outlier: 6.057A pdb=" N GLU B 61 " --> pdb=" O GLN F 62 " (cutoff:3.500A) removed outlier: 7.646A pdb=" N THR F 64 " --> pdb=" O GLU B 61 " (cutoff:3.500A) removed outlier: 6.581A pdb=" N VAL B 63 " --> pdb=" O THR F 64 " (cutoff:3.500A) removed outlier: 7.924A pdb=" N VAL F 66 " --> pdb=" O VAL B 63 " (cutoff:3.500A) removed outlier: 6.289A pdb=" N ASN B 65 " --> pdb=" O VAL F 66 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'E' and resid 71 through 72 removed outlier: 6.466A pdb=" N VAL B 71 " --> pdb=" O THR F 72 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'E' and resid 75 through 79 removed outlier: 5.994A pdb=" N THR B 75 " --> pdb=" O ALA E 76 " (cutoff:3.500A) removed outlier: 7.473A pdb=" N ALA E 78 " --> pdb=" O THR B 75 " (cutoff:3.500A) removed outlier: 6.274A pdb=" N VAL B 77 " --> pdb=" O ALA E 78 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'E' and resid 82 through 83 removed outlier: 6.362A pdb=" N VAL B 82 " --> pdb=" O GLU F 83 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'E' and resid 87 through 88 removed outlier: 6.898A pdb=" N SER B 87 " --> pdb=" O ILE E 88 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'C' and resid 37 through 42 removed outlier: 6.433A pdb=" N LEU A 38 " --> pdb=" O TYR C 39 " (cutoff:3.500A) removed outlier: 7.786A pdb=" N GLY C 41 " --> pdb=" O LEU A 38 " (cutoff:3.500A) removed outlier: 6.532A pdb=" N VAL A 40 " --> pdb=" O GLY C 41 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'C' and resid 48 through 57 removed outlier: 8.988A pdb=" N VAL A 48 " --> pdb=" O VAL C 49 " (cutoff:3.500A) removed outlier: 10.297A pdb=" N GLY C 51 " --> pdb=" O VAL A 48 " (cutoff:3.500A) removed outlier: 6.664A pdb=" N HIS A 50 " --> pdb=" O GLY C 51 " (cutoff:3.500A) removed outlier: 8.808A pdb=" N ALA C 53 " --> pdb=" O HIS A 50 " (cutoff:3.500A) removed outlier: 6.372A pdb=" N VAL A 52 " --> pdb=" O ALA C 53 " (cutoff:3.500A) removed outlier: 7.952A pdb=" N VAL C 55 " --> pdb=" O VAL A 52 " (cutoff:3.500A) removed outlier: 6.399A pdb=" N THR A 54 " --> pdb=" O VAL C 55 " (cutoff:3.500A) removed outlier: 8.449A pdb=" N GLU C 57 " --> pdb=" O THR A 54 " (cutoff:3.500A) removed outlier: 6.821A pdb=" N ALA A 56 " --> pdb=" O GLU C 57 " (cutoff:3.500A) removed outlier: 6.960A pdb=" N VAL A 49 " --> pdb=" O VAL D 48 " (cutoff:3.500A) removed outlier: 6.995A pdb=" N HIS D 50 " --> pdb=" O VAL A 49 " (cutoff:3.500A) removed outlier: 6.889A pdb=" N GLY A 51 " --> pdb=" O HIS D 50 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'C' and resid 61 through 67 removed outlier: 6.050A pdb=" N GLU A 61 " --> pdb=" O GLN D 62 " (cutoff:3.500A) removed outlier: 7.635A pdb=" N THR D 64 " --> pdb=" O GLU A 61 " (cutoff:3.500A) removed outlier: 6.570A pdb=" N VAL A 63 " --> pdb=" O THR D 64 " (cutoff:3.500A) removed outlier: 7.913A pdb=" N VAL D 66 " --> pdb=" O VAL A 63 " (cutoff:3.500A) removed outlier: 6.277A pdb=" N ASN A 65 " --> pdb=" O VAL D 66 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'C' and resid 71 through 72 removed outlier: 6.445A pdb=" N VAL A 71 " --> pdb=" O THR D 72 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'C' and resid 75 through 79 removed outlier: 6.033A pdb=" N THR A 75 " --> pdb=" O ALA C 76 " (cutoff:3.500A) removed outlier: 7.508A pdb=" N ALA C 78 " --> pdb=" O THR A 75 " (cutoff:3.500A) removed outlier: 6.300A pdb=" N VAL A 77 " --> pdb=" O ALA C 78 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'C' and resid 82 through 83 removed outlier: 6.374A pdb=" N VAL A 82 " --> pdb=" O GLU D 83 " (cutoff:3.500A) Processing sheet with id=AB5, first strand: chain 'C' and resid 87 through 88 removed outlier: 6.889A pdb=" N SER A 87 " --> pdb=" O ILE C 88 " (cutoff:3.500A) 43 hydrogen bonds defined for protein. 129 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.61 Time building geometry restraints manager: 1.08 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.23 - 1.29: 456 1.29 - 1.35: 407 1.35 - 1.42: 103 1.42 - 1.48: 426 1.48 - 1.55: 1218 Bond restraints: 2610 Sorted by residual: bond pdb=" C GLY B 67 " pdb=" O GLY B 67 " ideal model delta sigma weight residual 1.233 1.227 0.006 6.90e-03 2.10e+04 7.30e-01 bond pdb=" C GLY E 67 " pdb=" O GLY E 67 " ideal model delta sigma weight residual 1.233 1.227 0.006 6.90e-03 2.10e+04 6.58e-01 bond pdb=" C GLY C 67 " pdb=" O GLY C 67 " ideal model delta sigma weight residual 1.233 1.227 0.005 6.90e-03 2.10e+04 6.26e-01 bond pdb=" C GLY A 67 " pdb=" O GLY A 67 " ideal model delta sigma weight residual 1.233 1.227 0.005 6.90e-03 2.10e+04 6.24e-01 bond pdb=" C GLY D 67 " pdb=" O GLY D 67 " ideal model delta sigma weight residual 1.233 1.227 0.005 6.90e-03 2.10e+04 5.82e-01 ... (remaining 2605 not shown) Histogram of bond angle deviations from ideal: 106.10 - 111.13: 1174 111.13 - 116.17: 710 116.17 - 121.20: 824 121.20 - 126.24: 814 126.24 - 131.28: 6 Bond angle restraints: 3528 Sorted by residual: angle pdb=" CA GLY C 84 " pdb=" C GLY C 84 " pdb=" N ALA C 85 " ideal model delta sigma weight residual 115.25 116.89 -1.64 1.43e+00 4.89e-01 1.31e+00 angle pdb=" CA GLY B 84 " pdb=" C GLY B 84 " pdb=" N ALA B 85 " ideal model delta sigma weight residual 115.25 116.86 -1.61 1.43e+00 4.89e-01 1.27e+00 angle pdb=" CA GLY F 84 " pdb=" C GLY F 84 " pdb=" N ALA F 85 " ideal model delta sigma weight residual 115.25 116.85 -1.60 1.43e+00 4.89e-01 1.25e+00 angle pdb=" CA GLY D 84 " pdb=" C GLY D 84 " pdb=" N ALA D 85 " ideal model delta sigma weight residual 115.25 116.84 -1.59 1.43e+00 4.89e-01 1.23e+00 angle pdb=" CA GLY A 84 " pdb=" C GLY A 84 " pdb=" N ALA A 85 " ideal model delta sigma weight residual 115.25 116.83 -1.58 1.43e+00 4.89e-01 1.22e+00 ... (remaining 3523 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 10.58: 1314 10.58 - 21.16: 128 21.16 - 31.74: 46 31.74 - 42.32: 12 42.32 - 52.90: 18 Dihedral angle restraints: 1518 sinusoidal: 504 harmonic: 1014 Sorted by residual: dihedral pdb=" CA LEU B 38 " pdb=" CB LEU B 38 " pdb=" CG LEU B 38 " pdb=" CD1 LEU B 38 " ideal model delta sinusoidal sigma weight residual 180.00 139.66 40.34 3 1.50e+01 4.44e-03 7.18e+00 dihedral pdb=" CA LEU A 38 " pdb=" CB LEU A 38 " pdb=" CG LEU A 38 " pdb=" CD1 LEU A 38 " ideal model delta sinusoidal sigma weight residual 180.00 139.68 40.32 3 1.50e+01 4.44e-03 7.18e+00 dihedral pdb=" CA LEU D 38 " pdb=" CB LEU D 38 " pdb=" CG LEU D 38 " pdb=" CD1 LEU D 38 " ideal model delta sinusoidal sigma weight residual 180.00 139.68 40.32 3 1.50e+01 4.44e-03 7.18e+00 ... (remaining 1515 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.026: 224 0.026 - 0.051: 118 0.051 - 0.077: 48 0.077 - 0.103: 20 0.103 - 0.129: 52 Chirality restraints: 462 Sorted by residual: chirality pdb=" CA ILE A 88 " pdb=" N ILE A 88 " pdb=" C ILE A 88 " pdb=" CB ILE A 88 " both_signs ideal model delta sigma weight residual False 2.43 2.56 -0.13 2.00e-01 2.50e+01 4.13e-01 chirality pdb=" CA ILE D 88 " pdb=" N ILE D 88 " pdb=" C ILE D 88 " pdb=" CB ILE D 88 " both_signs ideal model delta sigma weight residual False 2.43 2.56 -0.13 2.00e-01 2.50e+01 4.09e-01 chirality pdb=" CA ILE B 88 " pdb=" N ILE B 88 " pdb=" C ILE B 88 " pdb=" CB ILE B 88 " both_signs ideal model delta sigma weight residual False 2.43 2.56 -0.13 2.00e-01 2.50e+01 4.06e-01 ... (remaining 459 not shown) Planarity restraints: 438 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA GLY A 68 " 0.002 2.00e-02 2.50e+03 4.77e-03 2.28e-01 pdb=" C GLY A 68 " -0.008 2.00e-02 2.50e+03 pdb=" O GLY A 68 " 0.003 2.00e-02 2.50e+03 pdb=" N ALA A 69 " 0.003 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA GLY D 68 " 0.002 2.00e-02 2.50e+03 4.74e-03 2.25e-01 pdb=" C GLY D 68 " -0.008 2.00e-02 2.50e+03 pdb=" O GLY D 68 " 0.003 2.00e-02 2.50e+03 pdb=" N ALA D 69 " 0.003 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA GLY B 68 " -0.002 2.00e-02 2.50e+03 4.73e-03 2.24e-01 pdb=" C GLY B 68 " 0.008 2.00e-02 2.50e+03 pdb=" O GLY B 68 " -0.003 2.00e-02 2.50e+03 pdb=" N ALA B 69 " -0.003 2.00e-02 2.50e+03 ... (remaining 435 not shown) Histogram of nonbonded interaction distances: 2.22 - 2.76: 223 2.76 - 3.29: 2264 3.29 - 3.83: 3987 3.83 - 4.36: 4368 4.36 - 4.90: 8731 Nonbonded interactions: 19573 Sorted by model distance: nonbonded pdb=" NE2 GLN C 62 " pdb=" OG1 THR C 64 " model vdw 2.223 2.520 nonbonded pdb=" NE2 GLN D 62 " pdb=" OG1 THR D 64 " model vdw 2.223 2.520 nonbonded pdb=" NE2 GLN F 62 " pdb=" OG1 THR F 64 " model vdw 2.223 2.520 nonbonded pdb=" NE2 GLN E 62 " pdb=" OG1 THR E 64 " model vdw 2.223 2.520 nonbonded pdb=" NE2 GLN B 62 " pdb=" OG1 THR B 64 " model vdw 2.223 2.520 ... (remaining 19568 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = chain 'A' selection = chain 'B' selection = chain 'C' selection = chain 'D' selection = chain 'E' selection = chain 'F' } 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: 1.000 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: 4.910 Check model and map are aligned: 0.040 Set scattering table: 0.030 Process input model: 12.300 Find NCS groups from input model: 0.200 Set up NCS constraints: 0.020 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:12.260 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 30.760 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7881 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.002 0.020 2610 Z= 0.125 Angle : 0.408 1.796 3528 Z= 0.248 Chirality : 0.051 0.129 462 Planarity : 0.001 0.005 438 Dihedral : 12.373 52.896 882 Min Nonbonded Distance : 2.223 Molprobity Statistics. All-atom Clashscore : 4.70 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.20 % Favored : 91.80 % Rotamer: Outliers : 0.00 % Allowed : 4.55 % Favored : 95.45 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.31 (0.34), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.52 (0.26), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.001 0.001 HIS B 50 PHE 0.007 0.002 PHE F 94 TYR 0.001 0.000 TYR A 39 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 49 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 49 time to evaluate : 0.302 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 49 average time/residue: 1.5808 time to fit residues: 78.8976 Evaluate side-chains 28 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 28 time to evaluate : 0.301 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.3766 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 36 random chunks: chunk 30 optimal weight: 10.0000 chunk 27 optimal weight: 10.0000 chunk 15 optimal weight: 10.0000 chunk 9 optimal weight: 8.9990 chunk 18 optimal weight: 6.9990 chunk 14 optimal weight: 5.9990 chunk 28 optimal weight: 3.9990 chunk 10 optimal weight: 4.9990 chunk 17 optimal weight: 8.9990 chunk 21 optimal weight: 1.9990 chunk 32 optimal weight: 9.9990 overall best weight: 4.7990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 62 GLN D 62 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7977 moved from start: 0.1938 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.027 2610 Z= 0.390 Angle : 0.581 5.123 3528 Z= 0.334 Chirality : 0.055 0.141 462 Planarity : 0.003 0.010 438 Dihedral : 5.719 15.954 378 Min Nonbonded Distance : 2.289 Molprobity Statistics. All-atom Clashscore : 3.19 Ramachandran Plot: Outliers : 0.00 % Allowed : 15.03 % Favored : 84.97 % Rotamer: Outliers : 1.52 % Allowed : 18.56 % Favored : 79.92 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.51 (0.40), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.90 (0.31), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.001 HIS D 50 PHE 0.013 0.003 PHE A 94 TYR 0.007 0.002 TYR D 39 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 37 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 33 time to evaluate : 0.328 Fit side-chains outliers start: 4 outliers final: 1 residues processed: 35 average time/residue: 1.5723 time to fit residues: 56.2406 Evaluate side-chains 29 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 28 time to evaluate : 0.308 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 1 residues processed: 0 time to fit residues: 0.3857 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 36 random chunks: chunk 18 optimal weight: 6.9990 chunk 10 optimal weight: 7.9990 chunk 27 optimal weight: 7.9990 chunk 22 optimal weight: 0.8980 chunk 9 optimal weight: 7.9990 chunk 32 optimal weight: 8.9990 chunk 35 optimal weight: 5.9990 chunk 29 optimal weight: 2.9990 chunk 11 optimal weight: 5.9990 chunk 26 optimal weight: 7.9990 chunk 24 optimal weight: 8.9990 overall best weight: 4.5788 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 62 GLN C 62 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7966 moved from start: 0.2048 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.027 2610 Z= 0.375 Angle : 0.596 5.031 3528 Z= 0.344 Chirality : 0.054 0.138 462 Planarity : 0.003 0.010 438 Dihedral : 5.778 17.356 378 Min Nonbonded Distance : 2.249 Molprobity Statistics. All-atom Clashscore : 3.01 Ramachandran Plot: Outliers : 0.00 % Allowed : 13.39 % Favored : 86.61 % Rotamer: Outliers : 1.89 % Allowed : 21.21 % Favored : 76.89 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.30 (0.41), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.74 (0.31), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.001 HIS D 50 PHE 0.011 0.003 PHE A 94 TYR 0.007 0.002 TYR D 39 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 32 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 27 time to evaluate : 0.280 Fit side-chains outliers start: 5 outliers final: 4 residues processed: 32 average time/residue: 1.5272 time to fit residues: 49.9388 Evaluate side-chains 30 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 26 time to evaluate : 0.304 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 4 residues processed: 0 time to fit residues: 0.3838 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 36 random chunks: chunk 16 optimal weight: 4.9990 chunk 3 optimal weight: 5.9990 chunk 15 optimal weight: 2.9990 chunk 21 optimal weight: 0.7980 chunk 32 optimal weight: 4.9990 chunk 34 optimal weight: 0.8980 chunk 17 optimal weight: 9.9990 chunk 31 optimal weight: 10.0000 chunk 9 optimal weight: 3.9990 chunk 28 optimal weight: 2.9990 chunk 19 optimal weight: 10.0000 overall best weight: 2.3386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 62 GLN C 62 GLN D 62 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7926 moved from start: 0.2062 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.017 2610 Z= 0.209 Angle : 0.494 4.189 3528 Z= 0.283 Chirality : 0.052 0.139 462 Planarity : 0.002 0.015 438 Dihedral : 5.260 16.474 378 Min Nonbonded Distance : 2.300 Molprobity Statistics. All-atom Clashscore : 2.63 Ramachandran Plot: Outliers : 0.00 % Allowed : 10.66 % Favored : 89.34 % Rotamer: Outliers : 2.65 % Allowed : 21.97 % Favored : 75.38 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.12 (0.42), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.61 (0.32), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.002 0.001 HIS F 50 PHE 0.007 0.002 PHE D 94 TYR 0.004 0.001 TYR A 39 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 37 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 30 time to evaluate : 0.301 Fit side-chains outliers start: 7 outliers final: 5 residues processed: 35 average time/residue: 1.4945 time to fit residues: 53.4605 Evaluate side-chains 33 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 28 time to evaluate : 0.292 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 4 residues processed: 1 average time/residue: 1.3795 time to fit residues: 1.7731 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 36 random chunks: chunk 0 optimal weight: 10.0000 chunk 25 optimal weight: 5.9990 chunk 14 optimal weight: 7.9990 chunk 29 optimal weight: 9.9990 chunk 24 optimal weight: 4.9990 chunk 17 optimal weight: 10.0000 chunk 31 optimal weight: 10.0000 chunk 8 optimal weight: 9.9990 chunk 11 optimal weight: 5.9990 chunk 6 optimal weight: 10.0000 chunk 20 optimal weight: 1.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 62 GLN C 62 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7977 moved from start: 0.2358 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.031 2610 Z= 0.430 Angle : 0.636 5.072 3528 Z= 0.367 Chirality : 0.055 0.141 462 Planarity : 0.003 0.011 438 Dihedral : 6.046 18.413 378 Min Nonbonded Distance : 2.235 Molprobity Statistics. All-atom Clashscore : 3.01 Ramachandran Plot: Outliers : 0.00 % Allowed : 16.12 % Favored : 83.88 % Rotamer: Outliers : 3.03 % Allowed : 21.97 % Favored : 75.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: -2.19 (0.42), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.66 (0.32), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.001 HIS D 50 PHE 0.014 0.004 PHE D 94 TYR 0.007 0.003 TYR D 39 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 38 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 30 time to evaluate : 0.316 Fit side-chains outliers start: 8 outliers final: 7 residues processed: 37 average time/residue: 1.5704 time to fit residues: 59.4338 Evaluate side-chains 37 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 30 time to evaluate : 0.309 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 7 residues processed: 0 time to fit residues: 0.3882 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 36 random chunks: chunk 8 optimal weight: 5.9990 chunk 34 optimal weight: 0.9980 chunk 28 optimal weight: 3.9990 chunk 16 optimal weight: 6.9990 chunk 2 optimal weight: 9.9990 chunk 11 optimal weight: 0.7980 chunk 18 optimal weight: 5.9990 chunk 33 optimal weight: 2.9990 chunk 3 optimal weight: 0.6980 chunk 19 optimal weight: 6.9990 chunk 25 optimal weight: 7.9990 overall best weight: 1.8984 ------------------------------------------------------------------------------- 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.7908 moved from start: 0.2172 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.016 2610 Z= 0.179 Angle : 0.485 3.855 3528 Z= 0.279 Chirality : 0.051 0.140 462 Planarity : 0.002 0.005 438 Dihedral : 5.172 16.405 378 Min Nonbonded Distance : 2.303 Molprobity Statistics. All-atom Clashscore : 2.63 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.29 % Favored : 90.71 % Rotamer: Outliers : 3.03 % Allowed : 21.97 % Favored : 75.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: -1.99 (0.42), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.50 (0.32), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.002 0.001 HIS F 50 PHE 0.007 0.002 PHE D 94 TYR 0.004 0.001 TYR C 39 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 35 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 27 time to evaluate : 0.312 Fit side-chains outliers start: 8 outliers final: 7 residues processed: 33 average time/residue: 1.2607 time to fit residues: 42.7193 Evaluate side-chains 33 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 26 time to evaluate : 0.278 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 6 residues processed: 1 average time/residue: 1.4345 time to fit residues: 1.8262 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 36 random chunks: chunk 19 optimal weight: 9.9990 chunk 29 optimal weight: 9.9990 chunk 34 optimal weight: 0.0270 chunk 21 optimal weight: 6.9990 chunk 16 optimal weight: 8.9990 chunk 13 optimal weight: 1.9990 chunk 20 optimal weight: 9.9990 chunk 10 optimal weight: 1.9990 chunk 6 optimal weight: 9.9990 chunk 22 optimal weight: 1.9990 chunk 23 optimal weight: 0.7980 overall best weight: 1.3644 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 62 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7886 moved from start: 0.2254 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.016 2610 Z= 0.140 Angle : 0.438 3.509 3528 Z= 0.252 Chirality : 0.051 0.140 462 Planarity : 0.001 0.004 438 Dihedral : 4.767 14.985 378 Min Nonbonded Distance : 2.330 Molprobity Statistics. All-atom Clashscore : 2.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 10.38 % Favored : 89.62 % Rotamer: Outliers : 3.03 % Allowed : 21.97 % Favored : 75.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: -1.89 (0.43), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.43 (0.33), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.001 0.001 HIS D 50 PHE 0.005 0.001 PHE F 94 TYR 0.005 0.001 TYR C 39 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 34 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 26 time to evaluate : 0.300 Fit side-chains outliers start: 8 outliers final: 7 residues processed: 33 average time/residue: 1.3714 time to fit residues: 46.3956 Evaluate side-chains 32 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 25 time to evaluate : 0.310 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 6 residues processed: 1 average time/residue: 1.3941 time to fit residues: 1.8003 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 36 random chunks: chunk 17 optimal weight: 0.7980 chunk 3 optimal weight: 10.0000 chunk 27 optimal weight: 0.4980 chunk 31 optimal weight: 10.0000 chunk 33 optimal weight: 7.9990 chunk 30 optimal weight: 8.9990 chunk 32 optimal weight: 9.9990 chunk 19 optimal weight: 10.0000 chunk 14 optimal weight: 0.0020 chunk 25 optimal weight: 4.9990 chunk 9 optimal weight: 2.9990 overall best weight: 1.8592 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 62 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7900 moved from start: 0.2254 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.014 2610 Z= 0.174 Angle : 0.460 3.461 3528 Z= 0.266 Chirality : 0.051 0.140 462 Planarity : 0.002 0.005 438 Dihedral : 4.931 16.160 378 Min Nonbonded Distance : 2.298 Molprobity Statistics. All-atom Clashscore : 2.44 Ramachandran Plot: Outliers : 0.00 % Allowed : 10.38 % Favored : 89.62 % Rotamer: Outliers : 2.65 % Allowed : 21.97 % Favored : 75.38 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.80 (0.43), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.37 (0.33), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.002 0.001 HIS D 50 PHE 0.006 0.002 PHE D 94 TYR 0.003 0.001 TYR A 39 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 32 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 25 time to evaluate : 0.304 Fit side-chains outliers start: 7 outliers final: 6 residues processed: 31 average time/residue: 1.5321 time to fit residues: 48.6338 Evaluate side-chains 31 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 25 time to evaluate : 0.299 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 6 residues processed: 0 time to fit residues: 0.3776 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 36 random chunks: chunk 29 optimal weight: 8.9990 chunk 30 optimal weight: 6.9990 chunk 32 optimal weight: 0.0270 chunk 21 optimal weight: 0.0970 chunk 34 optimal weight: 4.9990 chunk 20 optimal weight: 0.9980 chunk 16 optimal weight: 6.9990 chunk 23 optimal weight: 8.9990 chunk 35 optimal weight: 5.9990 chunk 33 optimal weight: 1.9990 chunk 28 optimal weight: 1.9990 overall best weight: 1.0240 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 62 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7857 moved from start: 0.2308 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.015 2610 Z= 0.118 Angle : 0.419 3.271 3528 Z= 0.242 Chirality : 0.051 0.141 462 Planarity : 0.001 0.004 438 Dihedral : 4.486 14.201 378 Min Nonbonded Distance : 2.331 Molprobity Statistics. All-atom Clashscore : 2.07 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.74 % Favored : 91.26 % Rotamer: Outliers : 3.03 % Allowed : 21.21 % Favored : 75.76 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.75 (0.43), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.32 (0.33), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.001 0.001 HIS D 50 PHE 0.004 0.001 PHE B 94 TYR 0.003 0.001 TYR A 39 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 33 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 25 time to evaluate : 0.274 Fit side-chains outliers start: 8 outliers final: 6 residues processed: 32 average time/residue: 1.4055 time to fit residues: 46.0713 Evaluate side-chains 30 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 24 time to evaluate : 0.363 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 6 residues processed: 0 time to fit residues: 0.3849 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 36 random chunks: chunk 2 optimal weight: 9.9990 chunk 22 optimal weight: 7.9990 chunk 17 optimal weight: 5.9990 chunk 30 optimal weight: 4.9990 chunk 8 optimal weight: 7.9990 chunk 26 optimal weight: 0.9990 chunk 4 optimal weight: 10.0000 chunk 7 optimal weight: 9.9990 chunk 28 optimal weight: 6.9990 chunk 11 optimal weight: 4.9990 chunk 29 optimal weight: 0.6980 overall best weight: 3.5388 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 62 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7945 moved from start: 0.2358 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.025 2610 Z= 0.298 Angle : 0.548 5.515 3528 Z= 0.316 Chirality : 0.053 0.144 462 Planarity : 0.002 0.007 438 Dihedral : 5.475 17.465 378 Min Nonbonded Distance : 2.258 Molprobity Statistics. All-atom Clashscore : 3.01 Ramachandran Plot: Outliers : 0.00 % Allowed : 15.57 % Favored : 84.43 % Rotamer: Outliers : 3.03 % Allowed : 21.21 % Favored : 75.76 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.74 (0.44), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.32 (0.33), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.002 0.001 HIS D 50 PHE 0.010 0.003 PHE D 94 TYR 0.006 0.002 TYR A 39 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 732 Ramachandran restraints generated. 366 Oldfield, 0 Emsley, 366 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 34 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 26 time to evaluate : 0.316 Fit side-chains outliers start: 8 outliers final: 8 residues processed: 34 average time/residue: 1.3507 time to fit residues: 47.1162 Evaluate side-chains 34 residues out of total 264 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 26 time to evaluate : 0.310 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 8 residues processed: 0 time to fit residues: 0.3902 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 36 random chunks: chunk 3 optimal weight: 9.9990 chunk 5 optimal weight: 10.0000 chunk 25 optimal weight: 9.9990 chunk 1 optimal weight: 2.9990 chunk 20 optimal weight: 0.9990 chunk 32 optimal weight: 0.8980 chunk 19 optimal weight: 10.0000 chunk 24 optimal weight: 6.9990 chunk 0 optimal weight: 8.9990 chunk 22 optimal weight: 6.9990 chunk 21 optimal weight: 0.0980 overall best weight: 2.3986 ------------------------------------------------------------------------------- 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.4108 r_free = 0.4108 target = 0.137595 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 50)----------------| | r_work = 0.3783 r_free = 0.3783 target = 0.115767 restraints weight = 3459.130| |-----------------------------------------------------------------------------| r_work (start): 0.3794 rms_B_bonded: 2.00 r_work: 0.3680 rms_B_bonded: 2.35 restraints_weight: 0.5000 r_work: 0.3573 rms_B_bonded: 3.91 restraints_weight: 0.2500 r_work (final): 0.3573 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8247 moved from start: 0.2281 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.020 2610 Z= 0.216 Angle : 0.504 3.784 3528 Z= 0.291 Chirality : 0.051 0.143 462 Planarity : 0.002 0.007 438 Dihedral : 5.205 16.596 378 Min Nonbonded Distance : 2.293 Molprobity Statistics. All-atom Clashscore : 2.44 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.84 % Favored : 90.16 % Rotamer: Outliers : 3.03 % Allowed : 21.21 % Favored : 75.76 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.76 (0.43), residues: 366 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -1.33 (0.33), residues: 366 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.002 0.001 HIS D 50 PHE 0.007 0.002 PHE F 94 TYR 0.004 0.001 TYR D 39 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 1488.86 seconds wall clock time: 27 minutes 0.48 seconds (1620.48 seconds total)