Starting phenix.real_space_refine on Thu Dec 7 15:48:52 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8jmt_36426/12_2023/8jmt_36426_neut.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8jmt_36426/12_2023/8jmt_36426.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.3 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8jmt_36426/12_2023/8jmt_36426.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8jmt_36426/12_2023/8jmt_36426.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8jmt_36426/12_2023/8jmt_36426_neut.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8jmt_36426/12_2023/8jmt_36426_neut.pdb" } resolution = 3.3 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.054 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 17 5.16 5 C 1284 2.51 5 N 296 2.21 5 O 297 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A PHE 884": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 914": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 1176": "OE1" <-> "OE2" 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: 1894 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 1894 Number of conformers: 1 Conformer: "" Number of residues, atoms: 234, 1894 Classifications: {'peptide': 234} Link IDs: {'PTRANS': 3, 'TRANS': 230} Chain breaks: 3 Time building chain proxies: 1.63, per 1000 atoms: 0.86 Number of scatterers: 1894 At special positions: 0 Unit cell: (63.13, 67.41, 65.27, 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 17 16.00 O 297 8.00 N 296 7.00 C 1284 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=2, symmetry=0 Simple disulfide: pdb=" SG CYS A 882 " - pdb=" SG CYS A1201 " distance=2.03 Simple disulfide: pdb=" SG CYS A 927 " - pdb=" SG CYS A 999 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 0.90 Conformation dependent library (CDL) restraints added in 311.8 milliseconds 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 442 Finding SS restraints... Secondary structure from input PDB file: 8 helices and 0 sheets defined 81.2% alpha, 0.0% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.18 Creating SS restraints... Processing helix chain 'A' and resid 861 through 882 removed outlier: 3.957A pdb=" N ILE A 867 " --> pdb=" O LEU A 863 " (cutoff:3.500A) removed outlier: 3.625A pdb=" N VAL A 870 " --> pdb=" O VAL A 866 " (cutoff:3.500A) removed outlier: 3.953A pdb=" N CYS A 882 " --> pdb=" O CYS A 878 " (cutoff:3.500A) Processing helix chain 'A' and resid 894 through 917 removed outlier: 3.880A pdb=" N ILE A 899 " --> pdb=" O ASP A 895 " (cutoff:3.500A) removed outlier: 3.694A pdb=" N HIS A 900 " --> pdb=" O ARG A 896 " (cutoff:3.500A) removed outlier: 3.681A pdb=" N CYS A 904 " --> pdb=" O HIS A 900 " (cutoff:3.500A) Processing helix chain 'A' and resid 923 through 958 removed outlier: 3.763A pdb=" N LEU A 939 " --> pdb=" O HIS A 935 " (cutoff:3.500A) Processing helix chain 'A' and resid 966 through 988 removed outlier: 3.645A pdb=" N MET A 977 " --> pdb=" O VAL A 973 " (cutoff:3.500A) Proline residue: A 978 - end of helix Processing helix chain 'A' and resid 1007 through 1036 Proline residue: A1014 - end of helix Processing helix chain 'A' and resid 1148 through 1168 Processing helix chain 'A' and resid 1181 through 1200 removed outlier: 3.751A pdb=" N THR A1186 " --> pdb=" O ALA A1182 " (cutoff:3.500A) removed outlier: 3.886A pdb=" N ILE A1187 " --> pdb=" O TYR A1183 " (cutoff:3.500A) removed outlier: 3.662A pdb=" N HIS A1200 " --> pdb=" O ILE A1196 " (cutoff:3.500A) Processing helix chain 'A' and resid 1206 through 1219 141 hydrogen bonds defined for protein. 423 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.41 Time building geometry restraints manager: 0.96 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.34: 545 1.34 - 1.46: 473 1.46 - 1.58: 900 1.58 - 1.70: 0 1.70 - 1.82: 25 Bond restraints: 1943 Sorted by residual: bond pdb=" CG LEU A1021 " pdb=" CD2 LEU A1021 " ideal model delta sigma weight residual 1.521 1.478 0.043 3.30e-02 9.18e+02 1.73e+00 bond pdb=" CA MET A 955 " pdb=" C MET A 955 " ideal model delta sigma weight residual 1.522 1.504 0.018 1.72e-02 3.38e+03 1.08e+00 bond pdb=" CA LEU A1021 " pdb=" C LEU A1021 " ideal model delta sigma weight residual 1.523 1.533 -0.010 1.37e-02 5.33e+03 5.85e-01 bond pdb=" N PHE A1025 " pdb=" CA PHE A1025 " ideal model delta sigma weight residual 1.459 1.468 -0.009 1.21e-02 6.83e+03 5.72e-01 bond pdb=" CA LYS A 901 " pdb=" C LYS A 901 " ideal model delta sigma weight residual 1.525 1.516 0.009 1.29e-02 6.01e+03 4.64e-01 ... (remaining 1938 not shown) Histogram of bond angle deviations from ideal: 98.44 - 105.56: 18 105.56 - 112.69: 1037 112.69 - 119.81: 640 119.81 - 126.94: 913 126.94 - 134.06: 25 Bond angle restraints: 2633 Sorted by residual: angle pdb=" N VAL A1179 " pdb=" CA VAL A1179 " pdb=" C VAL A1179 " ideal model delta sigma weight residual 113.00 109.13 3.87 1.30e+00 5.92e-01 8.88e+00 angle pdb=" N LYS A1206 " pdb=" CA LYS A1206 " pdb=" C LYS A1206 " ideal model delta sigma weight residual 113.43 110.03 3.40 1.26e+00 6.30e-01 7.28e+00 angle pdb=" N CYS A 887 " pdb=" CA CYS A 887 " pdb=" C CYS A 887 " ideal model delta sigma weight residual 114.56 111.15 3.41 1.27e+00 6.20e-01 7.19e+00 angle pdb=" CA LEU A1021 " pdb=" CB LEU A1021 " pdb=" CG LEU A1021 " ideal model delta sigma weight residual 116.30 125.25 -8.95 3.50e+00 8.16e-02 6.54e+00 angle pdb=" CA SER A1177 " pdb=" C SER A1177 " pdb=" O SER A1177 " ideal model delta sigma weight residual 121.32 118.53 2.79 1.16e+00 7.43e-01 5.79e+00 ... (remaining 2628 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.92: 986 17.92 - 35.84: 100 35.84 - 53.76: 19 53.76 - 71.68: 1 71.68 - 89.60: 2 Dihedral angle restraints: 1108 sinusoidal: 427 harmonic: 681 Sorted by residual: dihedral pdb=" CB CYS A 927 " pdb=" SG CYS A 927 " pdb=" SG CYS A 999 " pdb=" CB CYS A 999 " ideal model delta sinusoidal sigma weight residual -86.00 -127.16 41.16 1 1.00e+01 1.00e-02 2.37e+01 dihedral pdb=" CA TRP A1167 " pdb=" C TRP A1167 " pdb=" N ALA A1168 " pdb=" CA ALA A1168 " ideal model delta harmonic sigma weight residual 180.00 159.73 20.27 0 5.00e+00 4.00e-02 1.64e+01 dihedral pdb=" CA CYS A 882 " pdb=" C CYS A 882 " pdb=" N ILE A 883 " pdb=" CA ILE A 883 " ideal model delta harmonic sigma weight residual 180.00 160.31 19.69 0 5.00e+00 4.00e-02 1.55e+01 ... (remaining 1105 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.024: 155 0.024 - 0.048: 77 0.048 - 0.072: 57 0.072 - 0.096: 15 0.096 - 0.119: 8 Chirality restraints: 312 Sorted by residual: chirality pdb=" CG LEU A1021 " pdb=" CB LEU A1021 " pdb=" CD1 LEU A1021 " pdb=" CD2 LEU A1021 " both_signs ideal model delta sigma weight residual False -2.59 -2.71 0.12 2.00e-01 2.50e+01 3.57e-01 chirality pdb=" CA MET A 945 " pdb=" N MET A 945 " pdb=" C MET A 945 " pdb=" CB MET A 945 " both_signs ideal model delta sigma weight residual False 2.51 2.39 0.12 2.00e-01 2.50e+01 3.52e-01 chirality pdb=" CA PHE A1007 " pdb=" N PHE A1007 " pdb=" C PHE A1007 " pdb=" CB PHE A1007 " both_signs ideal model delta sigma weight residual False 2.51 2.40 0.11 2.00e-01 2.50e+01 3.29e-01 ... (remaining 309 not shown) Planarity restraints: 308 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C MET A 977 " -0.022 5.00e-02 4.00e+02 3.37e-02 1.82e+00 pdb=" N PRO A 978 " 0.058 5.00e-02 4.00e+02 pdb=" CA PRO A 978 " -0.017 5.00e-02 4.00e+02 pdb=" CD PRO A 978 " -0.019 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CA GLY A1164 " 0.006 2.00e-02 2.50e+03 1.21e-02 1.46e+00 pdb=" C GLY A1164 " -0.021 2.00e-02 2.50e+03 pdb=" O GLY A1164 " 0.008 2.00e-02 2.50e+03 pdb=" N LEU A1165 " 0.007 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE A1025 " 0.009 2.00e-02 2.50e+03 8.72e-03 1.33e+00 pdb=" CG PHE A1025 " -0.020 2.00e-02 2.50e+03 pdb=" CD1 PHE A1025 " 0.005 2.00e-02 2.50e+03 pdb=" CD2 PHE A1025 " 0.002 2.00e-02 2.50e+03 pdb=" CE1 PHE A1025 " -0.000 2.00e-02 2.50e+03 pdb=" CE2 PHE A1025 " 0.003 2.00e-02 2.50e+03 pdb=" CZ PHE A1025 " 0.001 2.00e-02 2.50e+03 ... (remaining 305 not shown) Histogram of nonbonded interaction distances: 2.20 - 2.74: 148 2.74 - 3.28: 2002 3.28 - 3.82: 3149 3.82 - 4.36: 3605 4.36 - 4.90: 6170 Nonbonded interactions: 15074 Sorted by model distance: nonbonded pdb=" O PHE A1185 " pdb=" ND2 ASN A1189 " model vdw 2.204 2.520 nonbonded pdb=" O TYR A1031 " pdb=" ND1 HIS A1035 " model vdw 2.302 2.520 nonbonded pdb=" O PHE A 884 " pdb=" N PHE A 888 " model vdw 2.346 2.520 nonbonded pdb=" OG SER A 875 " pdb=" O ASN A1189 " model vdw 2.385 2.440 nonbonded pdb=" O ASP A1004 " pdb=" OG1 THR A1005 " model vdw 2.478 2.440 ... (remaining 15069 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.850 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.010 Extract box with map and model: 3.420 Check model and map are aligned: 0.030 Set scattering table: 0.030 Process input model: 10.730 Find NCS groups from input model: 0.050 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.790 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 17.920 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5677 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.003 0.043 1943 Z= 0.196 Angle : 0.745 8.952 2633 Z= 0.456 Chirality : 0.040 0.119 312 Planarity : 0.004 0.034 308 Dihedral : 14.381 89.604 660 Min Nonbonded Distance : 2.204 Molprobity Statistics. All-atom Clashscore : 8.56 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.54 % Favored : 96.46 % 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: -0.01 (0.55), residues: 226 helix: 0.64 (0.37), residues: 184 sheet: None (None), residues: 0 loop : -2.71 (0.92), residues: 42 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP A 944 HIS 0.003 0.001 HIS A 935 PHE 0.020 0.002 PHE A1025 TYR 0.005 0.001 TYR A 969 ARG 0.002 0.000 ARG A 920 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 18 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 18 time to evaluate : 0.207 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 18 average time/residue: 0.0943 time to fit residues: 2.3850 Evaluate side-chains 16 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 16 time to evaluate : 0.203 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.2937 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 23 random chunks: chunk 19 optimal weight: 0.8980 chunk 17 optimal weight: 0.7980 chunk 9 optimal weight: 0.5980 chunk 5 optimal weight: 0.7980 chunk 11 optimal weight: 0.0970 chunk 18 optimal weight: 0.9990 chunk 6 optimal weight: 0.7980 chunk 10 optimal weight: 0.4980 chunk 13 optimal weight: 0.8980 chunk 20 optimal weight: 0.0980 chunk 14 optimal weight: 0.9980 overall best weight: 0.4178 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 900 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1204 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5783 moved from start: 0.1417 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 1943 Z= 0.204 Angle : 0.572 5.757 2633 Z= 0.293 Chirality : 0.041 0.117 312 Planarity : 0.004 0.029 308 Dihedral : 4.590 20.911 247 Min Nonbonded Distance : 2.211 Molprobity Statistics. All-atom Clashscore : 8.56 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.98 % Favored : 96.02 % Rotamer: Outliers : 0.49 % Allowed : 4.93 % Favored : 94.58 % 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: 0.53 (0.55), residues: 226 helix: 1.04 (0.37), residues: 187 sheet: None (None), residues: 0 loop : -2.98 (0.90), residues: 39 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.001 TRP A 944 HIS 0.003 0.001 HIS A 900 PHE 0.012 0.002 PHE A1199 TYR 0.004 0.001 TYR A 953 ARG 0.001 0.000 ARG A 920 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 19 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 18 time to evaluate : 0.214 Fit side-chains outliers start: 1 outliers final: 0 residues processed: 18 average time/residue: 0.1021 time to fit residues: 2.4621 Evaluate side-chains 17 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 17 time to evaluate : 0.201 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.2923 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 23 random chunks: chunk 5 optimal weight: 0.7980 chunk 20 optimal weight: 0.9990 chunk 22 optimal weight: 0.8980 chunk 18 optimal weight: 1.9990 chunk 7 optimal weight: 0.5980 chunk 16 optimal weight: 0.3980 chunk 15 optimal weight: 0.9980 chunk 10 optimal weight: 0.0980 chunk 2 optimal weight: 0.5980 chunk 9 optimal weight: 0.8980 chunk 14 optimal weight: 0.5980 overall best weight: 0.4580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 897 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 900 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1192 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A1204 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5837 moved from start: 0.1956 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.032 1943 Z= 0.209 Angle : 0.543 5.077 2633 Z= 0.288 Chirality : 0.040 0.119 312 Planarity : 0.003 0.030 308 Dihedral : 4.472 21.131 247 Min Nonbonded Distance : 2.227 Molprobity Statistics. All-atom Clashscore : 9.60 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.54 % Favored : 96.46 % Rotamer: Outliers : 0.00 % Allowed : 12.32 % Favored : 87.68 % 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: 0.88 (0.56), residues: 226 helix: 1.33 (0.37), residues: 187 sheet: None (None), residues: 0 loop : -3.13 (0.88), residues: 39 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.001 TRP A 944 HIS 0.003 0.001 HIS A 900 PHE 0.011 0.002 PHE A1199 TYR 0.006 0.001 TYR A1031 ARG 0.001 0.000 ARG A 896 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 16 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 16 time to evaluate : 0.201 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 16 average time/residue: 0.1025 time to fit residues: 2.2366 Evaluate side-chains 15 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 15 time to evaluate : 0.202 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.2917 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 23 random chunks: chunk 20 optimal weight: 0.9990 chunk 22 optimal weight: 0.4980 chunk 10 optimal weight: 0.0020 chunk 19 optimal weight: 3.9990 chunk 5 optimal weight: 1.9990 chunk 18 optimal weight: 0.5980 chunk 12 optimal weight: 0.8980 chunk 0 optimal weight: 5.9990 chunk 16 optimal weight: 0.0970 chunk 9 optimal weight: 0.5980 chunk 15 optimal weight: 0.5980 overall best weight: 0.3586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 897 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5816 moved from start: 0.2162 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.030 1943 Z= 0.176 Angle : 0.506 4.919 2633 Z= 0.267 Chirality : 0.039 0.103 312 Planarity : 0.003 0.030 308 Dihedral : 4.284 20.002 247 Min Nonbonded Distance : 2.228 Molprobity Statistics. All-atom Clashscore : 9.08 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.98 % Favored : 96.02 % Rotamer: Outliers : 1.48 % Allowed : 12.81 % Favored : 85.71 % 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.27 (0.56), residues: 226 helix: 1.59 (0.37), residues: 187 sheet: None (None), residues: 0 loop : -2.94 (0.85), residues: 39 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.007 0.001 TRP A 944 HIS 0.002 0.001 HIS A1035 PHE 0.010 0.001 PHE A 937 TYR 0.005 0.001 TYR A 953 ARG 0.001 0.000 ARG A 967 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 19 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 16 time to evaluate : 0.190 Fit side-chains outliers start: 3 outliers final: 0 residues processed: 17 average time/residue: 0.0975 time to fit residues: 2.3009 Evaluate side-chains 15 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 15 time to evaluate : 0.187 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.2852 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 23 random chunks: chunk 0 optimal weight: 8.9990 chunk 11 optimal weight: 0.0870 chunk 19 optimal weight: 2.9990 chunk 5 optimal weight: 0.5980 chunk 7 optimal weight: 0.3980 chunk 20 optimal weight: 0.7980 chunk 4 optimal weight: 0.5980 chunk 13 optimal weight: 0.9990 chunk 22 optimal weight: 0.4980 chunk 18 optimal weight: 0.6980 chunk 10 optimal weight: 0.9990 overall best weight: 0.4358 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 897 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5861 moved from start: 0.2452 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.030 1943 Z= 0.198 Angle : 0.531 4.941 2633 Z= 0.277 Chirality : 0.040 0.114 312 Planarity : 0.003 0.031 308 Dihedral : 4.344 21.284 247 Min Nonbonded Distance : 2.231 Molprobity Statistics. All-atom Clashscore : 8.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.42 % Favored : 95.58 % Rotamer: Outliers : 1.48 % Allowed : 15.27 % Favored : 83.25 % 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.37 (0.55), residues: 226 helix: 1.67 (0.37), residues: 187 sheet: None (None), residues: 0 loop : -3.00 (0.85), residues: 39 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.007 0.001 TRP A 944 HIS 0.003 0.001 HIS A1035 PHE 0.011 0.001 PHE A1195 TYR 0.006 0.001 TYR A1031 ARG 0.002 0.000 ARG A 896 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 20 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 17 time to evaluate : 0.224 Fit side-chains outliers start: 3 outliers final: 2 residues processed: 19 average time/residue: 0.0988 time to fit residues: 2.5468 Evaluate side-chains 19 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 17 time to evaluate : 0.199 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0310 time to fit residues: 0.3837 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 23 random chunks: chunk 1 optimal weight: 0.6980 chunk 7 optimal weight: 0.7980 chunk 11 optimal weight: 0.5980 chunk 21 optimal weight: 0.4980 chunk 2 optimal weight: 0.1980 chunk 12 optimal weight: 0.5980 chunk 16 optimal weight: 0.3980 chunk 18 optimal weight: 0.9990 chunk 22 optimal weight: 0.9980 chunk 13 optimal weight: 0.1980 chunk 10 optimal weight: 0.9980 overall best weight: 0.3780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 897 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5841 moved from start: 0.2525 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.030 1943 Z= 0.185 Angle : 0.558 9.149 2633 Z= 0.277 Chirality : 0.039 0.120 312 Planarity : 0.003 0.031 308 Dihedral : 4.232 20.010 247 Min Nonbonded Distance : 2.237 Molprobity Statistics. All-atom Clashscore : 10.89 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.98 % Favored : 96.02 % Rotamer: Outliers : 0.49 % Allowed : 16.75 % Favored : 82.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.56 (0.56), residues: 226 helix: 1.82 (0.37), residues: 187 sheet: None (None), residues: 0 loop : -2.98 (0.83), residues: 39 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.007 0.001 TRP A 944 HIS 0.003 0.001 HIS A1035 PHE 0.011 0.001 PHE A 937 TYR 0.008 0.001 TYR A 953 ARG 0.001 0.000 ARG A 967 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 18 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 17 time to evaluate : 0.213 Fit side-chains outliers start: 1 outliers final: 1 residues processed: 18 average time/residue: 0.1022 time to fit residues: 2.4711 Evaluate side-chains 18 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 17 time to evaluate : 0.161 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0238 time to fit residues: 0.3497 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 23 random chunks: chunk 13 optimal weight: 0.5980 chunk 8 optimal weight: 0.5980 chunk 6 optimal weight: 2.9990 chunk 4 optimal weight: 0.3980 chunk 14 optimal weight: 0.6980 chunk 15 optimal weight: 0.5980 chunk 10 optimal weight: 0.4980 chunk 2 optimal weight: 0.8980 chunk 17 optimal weight: 0.2980 chunk 20 optimal weight: 0.5980 chunk 21 optimal weight: 0.7980 overall best weight: 0.4780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 897 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5887 moved from start: 0.2773 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.032 1943 Z= 0.216 Angle : 0.587 8.223 2633 Z= 0.293 Chirality : 0.040 0.122 312 Planarity : 0.003 0.030 308 Dihedral : 4.362 21.642 247 Min Nonbonded Distance : 2.236 Molprobity Statistics. All-atom Clashscore : 11.15 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.98 % Favored : 96.02 % Rotamer: Outliers : 0.99 % Allowed : 17.73 % Favored : 81.28 % 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.44 (0.55), residues: 226 helix: 1.75 (0.37), residues: 186 sheet: None (None), residues: 0 loop : -3.04 (0.81), residues: 40 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A 944 HIS 0.003 0.001 HIS A1036 PHE 0.011 0.001 PHE A 884 TYR 0.008 0.001 TYR A 953 ARG 0.002 0.000 ARG A 896 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 21 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 19 time to evaluate : 0.198 Fit side-chains outliers start: 2 outliers final: 1 residues processed: 21 average time/residue: 0.0945 time to fit residues: 2.6560 Evaluate side-chains 20 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 19 time to evaluate : 0.197 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0283 time to fit residues: 0.3333 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 23 random chunks: chunk 19 optimal weight: 0.4980 chunk 20 optimal weight: 0.5980 chunk 21 optimal weight: 0.7980 chunk 12 optimal weight: 0.9980 chunk 8 optimal weight: 0.5980 chunk 16 optimal weight: 0.0870 chunk 6 optimal weight: 3.9990 chunk 18 optimal weight: 0.6980 chunk 13 optimal weight: 0.5980 chunk 10 optimal weight: 0.6980 chunk 15 optimal weight: 0.8980 overall best weight: 0.4758 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 897 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1192 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5884 moved from start: 0.2798 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.032 1943 Z= 0.224 Angle : 0.626 8.628 2633 Z= 0.306 Chirality : 0.042 0.122 312 Planarity : 0.003 0.031 308 Dihedral : 4.376 20.887 247 Min Nonbonded Distance : 2.239 Molprobity Statistics. All-atom Clashscore : 12.19 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.42 % Favored : 95.58 % Rotamer: Outliers : 0.00 % Allowed : 18.72 % Favored : 81.28 % 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.48 (0.55), residues: 226 helix: 1.79 (0.37), residues: 185 sheet: None (None), residues: 0 loop : -2.96 (0.78), residues: 41 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A1009 HIS 0.002 0.001 HIS A1036 PHE 0.015 0.002 PHE A1007 TYR 0.011 0.001 TYR A 953 ARG 0.002 0.000 ARG A 896 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 18 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 18 time to evaluate : 0.198 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 18 average time/residue: 0.2144 time to fit residues: 5.0265 Evaluate side-chains 17 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 17 time to evaluate : 0.561 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.7663 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 23 random chunks: chunk 22 optimal weight: 0.6980 chunk 21 optimal weight: 0.2980 chunk 18 optimal weight: 0.7980 chunk 1 optimal weight: 0.6980 chunk 14 optimal weight: 0.9980 chunk 11 optimal weight: 0.2980 chunk 19 optimal weight: 1.9990 chunk 5 optimal weight: 0.6980 chunk 16 optimal weight: 0.4980 chunk 2 optimal weight: 0.0870 chunk 7 optimal weight: 0.3980 overall best weight: 0.3158 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A1192 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5834 moved from start: 0.2849 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.030 1943 Z= 0.180 Angle : 0.579 7.979 2633 Z= 0.283 Chirality : 0.040 0.128 312 Planarity : 0.003 0.030 308 Dihedral : 4.199 19.657 247 Min Nonbonded Distance : 2.239 Molprobity Statistics. All-atom Clashscore : 11.67 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.54 % Favored : 96.46 % Rotamer: Outliers : 0.00 % Allowed : 19.21 % Favored : 80.79 % 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.77 (0.56), residues: 226 helix: 2.00 (0.37), residues: 186 sheet: None (None), residues: 0 loop : -2.99 (0.79), residues: 40 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A1009 HIS 0.002 0.001 HIS A 935 PHE 0.015 0.001 PHE A1007 TYR 0.009 0.001 TYR A 953 ARG 0.001 0.000 ARG A 967 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 17 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 17 time to evaluate : 0.213 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 17 average time/residue: 0.1016 time to fit residues: 2.3726 Evaluate side-chains 16 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 16 time to evaluate : 0.210 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.3036 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 23 random chunks: chunk 18 optimal weight: 0.5980 chunk 2 optimal weight: 0.8980 chunk 3 optimal weight: 0.9990 chunk 16 optimal weight: 0.5980 chunk 1 optimal weight: 0.7980 chunk 13 optimal weight: 0.7980 chunk 20 optimal weight: 0.6980 chunk 12 optimal weight: 3.9990 chunk 15 optimal weight: 0.1980 chunk 0 optimal weight: 2.9990 chunk 14 optimal weight: 0.0980 overall best weight: 0.4380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 897 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1192 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5866 moved from start: 0.2941 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.031 1943 Z= 0.213 Angle : 0.607 7.875 2633 Z= 0.297 Chirality : 0.041 0.126 312 Planarity : 0.003 0.030 308 Dihedral : 4.346 20.846 247 Min Nonbonded Distance : 2.236 Molprobity Statistics. All-atom Clashscore : 12.19 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.98 % Favored : 96.02 % Rotamer: Outliers : 0.00 % Allowed : 19.21 % Favored : 80.79 % 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.60 (0.55), residues: 226 helix: 1.88 (0.37), residues: 186 sheet: None (None), residues: 0 loop : -3.02 (0.79), residues: 40 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A 944 HIS 0.002 0.001 HIS A 935 PHE 0.013 0.002 PHE A1007 TYR 0.011 0.001 TYR A 953 ARG 0.002 0.000 ARG A 896 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 452 Ramachandran restraints generated. 226 Oldfield, 0 Emsley, 226 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 16 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 16 time to evaluate : 0.215 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 16 average time/residue: 0.0961 time to fit residues: 2.1761 Evaluate side-chains 16 residues out of total 203 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 16 time to evaluate : 0.198 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.2864 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 23 random chunks: chunk 13 optimal weight: 0.4980 chunk 8 optimal weight: 0.7980 chunk 22 optimal weight: 0.7980 chunk 0 optimal weight: 8.9990 chunk 4 optimal weight: 0.5980 chunk 2 optimal weight: 0.7980 chunk 18 optimal weight: 0.5980 chunk 21 optimal weight: 0.9990 chunk 9 optimal weight: 0.7980 chunk 5 optimal weight: 0.0970 chunk 14 optimal weight: 0.8980 overall best weight: 0.5178 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A1192 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4512 r_free = 0.4512 target = 0.243339 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 49)----------------| | r_work = 0.4233 r_free = 0.4233 target = 0.204777 restraints weight = 2063.646| |-----------------------------------------------------------------------------| r_work (start): 0.4219 rms_B_bonded: 2.09 r_work: 0.4092 rms_B_bonded: 3.47 restraints_weight: 0.5000 r_work (final): 0.4092 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5789 moved from start: 0.3061 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.032 1943 Z= 0.235 Angle : 0.623 7.612 2633 Z= 0.306 Chirality : 0.042 0.128 312 Planarity : 0.003 0.030 308 Dihedral : 4.457 21.366 247 Min Nonbonded Distance : 2.237 Molprobity Statistics. All-atom Clashscore : 11.67 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.98 % Favored : 96.02 % Rotamer: Outliers : 0.00 % Allowed : 19.21 % Favored : 80.79 % 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.51 (0.55), residues: 226 helix: 1.82 (0.37), residues: 185 sheet: None (None), residues: 0 loop : -2.99 (0.77), residues: 41 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.007 0.001 TRP A 944 HIS 0.002 0.001 HIS A 935 PHE 0.014 0.002 PHE A1007 TYR 0.011 0.001 TYR A 953 ARG 0.003 0.000 ARG A 896 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 964.07 seconds wall clock time: 23 minutes 35.14 seconds (1415.14 seconds total)