Starting phenix.real_space_refine on Thu Dec 7 15:49:01 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7us9_26729/12_2023/7us9_26729.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7us9_26729/12_2023/7us9_26729.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.8 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7us9_26729/12_2023/7us9_26729.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7us9_26729/12_2023/7us9_26729.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7us9_26729/12_2023/7us9_26729.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7us9_26729/12_2023/7us9_26729.pdb" } resolution = 3.8 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= 0.000 sd= 0.019 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 4 5.16 5 C 744 2.51 5 N 185 2.21 5 O 194 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A PHE 39": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 43": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 58": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 59": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 112": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 159": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 180": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 197": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 202": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 236": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.00s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/chem_data/mon_lib" Total number of atoms: 1127 Number of models: 1 Model: "" Number of chains: 3 Chain: "A" Number of atoms: 1071 Number of conformers: 1 Conformer: "" Number of residues, atoms: 145, 1071 Classifications: {'peptide': 145} Incomplete info: {'truncation_to_alanine': 52} Link IDs: {'PCIS': 1, 'PTRANS': 6, 'TRANS': 137} Chain breaks: 7 Unresolved non-hydrogen bonds: 134 Unresolved non-hydrogen angles: 176 Unresolved non-hydrogen dihedrals: 94 Unresolved non-hydrogen chiralities: 27 Planarities with less than four sites: {'GLN:plan1': 1, 'GLU:plan': 3, 'HIS:plan': 1, 'ASN:plan1': 5, 'ASP:plan': 4} Unresolved non-hydrogen planarities: 48 Chain: "B" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Chain: "A" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Time building chain proxies: 1.32, per 1000 atoms: 1.17 Number of scatterers: 1127 At special positions: 0 Unit cell: (50.58, 48.894, 59.853, 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 4 16.00 O 194 8.00 N 185 7.00 C 744 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=1, symmetry=0 Simple disulfide: pdb=" SG CYS A 133 " - pdb=" SG CYS A 153 " distance=2.09 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 Links applied BETA1-4 " NAG B 1 " - " NAG B 2 " NAG-ASN " NAG A 301 " - " ASN A 42 " " NAG A 302 " - " ASN A 242 " " NAG B 1 " - " ASN A 249 " Time building additional restraints: 0.42 Conformation dependent library (CDL) restraints added in 343.3 milliseconds 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 274 Finding SS restraints... Secondary structure from input PDB file: 1 helices and 3 sheets defined 9.7% alpha, 43.4% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.16 Creating SS restraints... Processing helix chain 'A' and resid 32 through 45 removed outlier: 4.150A pdb=" N ILE A 36 " --> pdb=" O HIS A 32 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 52 through 57 Processing sheet with id=AA2, first strand: chain 'A' and resid 78 through 81 removed outlier: 7.641A pdb=" N PHE A 202 " --> pdb=" O LYS A 110 " (cutoff:3.500A) removed outlier: 7.859A pdb=" N LYS A 110 " --> pdb=" O PHE A 202 " (cutoff:3.500A) removed outlier: 5.324A pdb=" N THR A 204 " --> pdb=" O LYS A 108 " (cutoff:3.500A) removed outlier: 3.596A pdb=" N GLY A 109 " --> pdb=" O TYR A 112 " (cutoff:3.500A) removed outlier: 4.225A pdb=" N TYR A 112 " --> pdb=" O GLY A 109 " (cutoff:3.500A) removed outlier: 6.788A pdb=" N ILE A 130 " --> pdb=" O ASN A 156 " (cutoff:3.500A) removed outlier: 4.145A pdb=" N ASN A 156 " --> pdb=" O ILE A 130 " (cutoff:3.500A) removed outlier: 6.369A pdb=" N ILE A 132 " --> pdb=" O ARG A 154 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 191 through 196 removed outlier: 3.514A pdb=" N ALA A 85 " --> pdb=" O VAL A 230 " (cutoff:3.500A) 54 hydrogen bonds defined for protein. 141 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.30 Time building geometry restraints manager: 0.47 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.20 - 1.33: 251 1.33 - 1.46: 358 1.46 - 1.58: 546 1.58 - 1.71: 1 1.71 - 1.84: 6 Bond restraints: 1162 Sorted by residual: bond pdb=" N PRO A 111 " pdb=" CD PRO A 111 " ideal model delta sigma weight residual 1.474 1.381 0.093 1.40e-02 5.10e+03 4.37e+01 bond pdb=" CZ2 TRP A 225 " pdb=" CH2 TRP A 225 " ideal model delta sigma weight residual 1.368 1.455 -0.087 1.90e-02 2.77e+03 2.10e+01 bond pdb=" CD2 TRP A 225 " pdb=" CE3 TRP A 225 " ideal model delta sigma weight residual 1.398 1.468 -0.070 1.60e-02 3.91e+03 1.90e+01 bond pdb=" CB LEU A 76 " pdb=" CG LEU A 76 " ideal model delta sigma weight residual 1.530 1.615 -0.085 2.00e-02 2.50e+03 1.83e+01 bond pdb=" C LEU A 47 " pdb=" N PRO A 48 " ideal model delta sigma weight residual 1.329 1.371 -0.042 1.18e-02 7.18e+03 1.24e+01 ... (remaining 1157 not shown) Histogram of bond angle deviations from ideal: 98.41 - 105.51: 32 105.51 - 112.61: 561 112.61 - 119.71: 360 119.71 - 126.81: 608 126.81 - 133.91: 28 Bond angle restraints: 1589 Sorted by residual: angle pdb=" N ILE A 195 " pdb=" CA ILE A 195 " pdb=" C ILE A 195 " ideal model delta sigma weight residual 108.11 117.48 -9.37 1.40e+00 5.10e-01 4.47e+01 angle pdb=" C PHE A 159 " pdb=" N PRO A 160 " pdb=" CA PRO A 160 " ideal model delta sigma weight residual 120.31 126.85 -6.54 9.80e-01 1.04e+00 4.45e+01 angle pdb=" C PRO A 48 " pdb=" N PRO A 49 " pdb=" CA PRO A 49 " ideal model delta sigma weight residual 119.84 127.79 -7.95 1.25e+00 6.40e-01 4.05e+01 angle pdb=" C LEU A 47 " pdb=" N PRO A 48 " pdb=" CA PRO A 48 " ideal model delta sigma weight residual 120.38 126.66 -6.28 1.03e+00 9.43e-01 3.72e+01 angle pdb=" N VAL A 78 " pdb=" CA VAL A 78 " pdb=" CB VAL A 78 " ideal model delta sigma weight residual 110.72 104.39 6.33 1.09e+00 8.42e-01 3.37e+01 ... (remaining 1584 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 8.01: 623 8.01 - 16.02: 44 16.02 - 24.03: 10 24.03 - 32.04: 7 32.04 - 40.05: 4 Dihedral angle restraints: 688 sinusoidal: 266 harmonic: 422 Sorted by residual: dihedral pdb=" C GLN A 63 " pdb=" N GLN A 63 " pdb=" CA GLN A 63 " pdb=" CB GLN A 63 " ideal model delta harmonic sigma weight residual -122.60 -130.32 7.72 0 2.50e+00 1.60e-01 9.54e+00 dihedral pdb=" CA ILE A 195 " pdb=" CB ILE A 195 " pdb=" CG1 ILE A 195 " pdb=" CD1 ILE A 195 " ideal model delta sinusoidal sigma weight residual 60.00 100.05 -40.05 3 1.50e+01 4.44e-03 7.12e+00 dihedral pdb=" C LYS A 110 " pdb=" N LYS A 110 " pdb=" CA LYS A 110 " pdb=" CB LYS A 110 " ideal model delta harmonic sigma weight residual -122.60 -116.09 -6.51 0 2.50e+00 1.60e-01 6.77e+00 ... (remaining 685 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.078: 122 0.078 - 0.157: 30 0.157 - 0.235: 13 0.235 - 0.313: 9 0.313 - 0.392: 4 Chirality restraints: 178 Sorted by residual: chirality pdb=" C1 NAG B 2 " pdb=" O4 NAG B 1 " pdb=" C2 NAG B 2 " pdb=" O5 NAG B 2 " both_signs ideal model delta sigma weight residual False -2.40 -2.57 0.17 2.00e-02 2.50e+03 7.42e+01 chirality pdb=" C2 NAG B 1 " pdb=" C1 NAG B 1 " pdb=" C3 NAG B 1 " pdb=" N2 NAG B 1 " both_signs ideal model delta sigma weight residual False -2.49 -2.88 0.39 2.00e-01 2.50e+01 3.84e+00 chirality pdb=" CG LEU A 76 " pdb=" CB LEU A 76 " pdb=" CD1 LEU A 76 " pdb=" CD2 LEU A 76 " both_signs ideal model delta sigma weight residual False -2.59 -2.97 0.38 2.00e-01 2.50e+01 3.59e+00 ... (remaining 175 not shown) Planarity restraints: 196 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB ASN A 42 " -0.046 2.00e-02 2.50e+03 4.80e-02 2.88e+01 pdb=" CG ASN A 42 " 0.020 2.00e-02 2.50e+03 pdb=" OD1 ASN A 42 " 0.009 2.00e-02 2.50e+03 pdb=" ND2 ASN A 42 " 0.075 2.00e-02 2.50e+03 pdb=" C1 NAG A 301 " -0.058 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB ASN A 242 " -0.039 2.00e-02 2.50e+03 4.07e-02 2.07e+01 pdb=" CG ASN A 242 " 0.017 2.00e-02 2.50e+03 pdb=" OD1 ASN A 242 " 0.008 2.00e-02 2.50e+03 pdb=" ND2 ASN A 242 " 0.063 2.00e-02 2.50e+03 pdb=" C1 NAG A 302 " -0.049 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA PRO A 64 " -0.008 2.00e-02 2.50e+03 1.70e-02 2.90e+00 pdb=" C PRO A 64 " 0.029 2.00e-02 2.50e+03 pdb=" O PRO A 64 " -0.011 2.00e-02 2.50e+03 pdb=" N TRP A 65 " -0.010 2.00e-02 2.50e+03 ... (remaining 193 not shown) Histogram of nonbonded interaction distances: 2.49 - 2.97: 519 2.97 - 3.45: 846 3.45 - 3.94: 1656 3.94 - 4.42: 1735 4.42 - 4.90: 3029 Nonbonded interactions: 7785 Sorted by model distance: nonbonded pdb=" O LYS A 110 " pdb=" CD PRO A 111 " model vdw 2.489 2.752 nonbonded pdb=" O GLN A 178 " pdb=" CA TRP A 179 " model vdw 2.627 2.776 nonbonded pdb=" N GLY A 109 " pdb=" O GLY A 109 " model vdw 2.636 2.496 nonbonded pdb=" N MET A 80 " pdb=" O MET A 80 " model vdw 2.641 2.496 nonbonded pdb=" O HIS A 189 " pdb=" ND1 HIS A 189 " model vdw 2.643 2.520 ... (remaining 7780 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 group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 5.420 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: 25.530 Check model and map are aligned: 0.010 Set scattering table: 0.020 Process input model: 8.910 Find NCS groups from input model: 0.030 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.870 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 42.800 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7302 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.015 0.093 1162 Z= 0.922 Angle : 1.582 11.720 1589 Z= 1.016 Chirality : 0.112 0.392 178 Planarity : 0.005 0.017 193 Dihedral : 7.489 40.047 411 Min Nonbonded Distance : 2.489 Molprobity Statistics. All-atom Clashscore : 3.45 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.33 % Favored : 97.67 % Rotamer: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 14.29 % 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.07 (0.73), residues: 129 helix: 3.26 (1.42), residues: 12 sheet: 1.80 (0.75), residues: 47 loop : -0.57 (0.71), residues: 70 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.015 0.004 TRP A 65 HIS 0.004 0.002 HIS A 189 PHE 0.011 0.003 PHE A 227 TYR 0.015 0.004 TYR A 43 ARG 0.002 0.000 ARG A 154 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 33 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 33 time to evaluate : 0.124 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 33 average time/residue: 0.1203 time to fit residues: 4.4836 Evaluate side-chains 21 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 21 time to evaluate : 0.127 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.1682 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 14 random chunks: chunk 11 optimal weight: 0.0570 chunk 10 optimal weight: 0.5980 chunk 5 optimal weight: 0.7980 chunk 3 optimal weight: 0.0970 chunk 7 optimal weight: 0.9980 chunk 4 optimal weight: 0.6980 chunk 6 optimal weight: 0.1980 chunk 8 optimal weight: 0.3980 chunk 12 optimal weight: 0.0570 chunk 13 optimal weight: 0.1980 chunk 9 optimal weight: 0.1980 overall best weight: 0.1214 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 199 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7362 moved from start: 0.2864 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 1162 Z= 0.204 Angle : 0.649 6.087 1589 Z= 0.344 Chirality : 0.052 0.195 178 Planarity : 0.004 0.029 193 Dihedral : 7.665 53.577 235 Min Nonbonded Distance : 2.421 Molprobity Statistics. All-atom Clashscore : 5.92 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.55 % Favored : 98.45 % Rotamer: Outliers : 2.47 % Allowed : 8.64 % Favored : 88.89 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 14.29 % 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.43 (0.79), residues: 129 helix: 3.83 (1.25), residues: 12 sheet: 1.48 (0.59), residues: 75 loop : -1.00 (1.12), residues: 42 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.001 TRP A 179 HIS 0.000 0.000 HIS A 189 PHE 0.008 0.001 PHE A 180 TYR 0.008 0.001 TYR A 43 ARG 0.001 0.000 ARG A 154 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 22 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 20 time to evaluate : 0.127 Fit side-chains outliers start: 2 outliers final: 1 residues processed: 22 average time/residue: 0.0671 time to fit residues: 1.9043 Evaluate side-chains 17 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 16 time to evaluate : 0.113 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0163 time to fit residues: 0.1935 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 14 random chunks: chunk 6 optimal weight: 0.3980 chunk 1 optimal weight: 0.3980 chunk 8 optimal weight: 0.9990 chunk 12 optimal weight: 0.9990 chunk 13 optimal weight: 0.1980 chunk 3 optimal weight: 0.7980 chunk 11 optimal weight: 0.9980 chunk 7 optimal weight: 0.6980 chunk 0 optimal weight: 0.4980 chunk 10 optimal weight: 0.9990 chunk 5 optimal weight: 0.8980 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 199 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7560 moved from start: 0.3777 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.034 1162 Z= 0.398 Angle : 0.809 8.811 1589 Z= 0.422 Chirality : 0.055 0.185 178 Planarity : 0.007 0.066 193 Dihedral : 7.806 42.420 235 Min Nonbonded Distance : 2.354 Molprobity Statistics. All-atom Clashscore : 5.92 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.88 % Favored : 96.12 % Rotamer: Outliers : 3.70 % Allowed : 9.88 % Favored : 86.42 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 14.29 % 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.07 (0.73), residues: 129 helix: 1.61 (1.36), residues: 12 sheet: 0.79 (0.59), residues: 77 loop : -1.98 (0.86), residues: 40 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.002 TRP A 179 HIS 0.003 0.001 HIS A 189 PHE 0.016 0.003 PHE A 197 TYR 0.015 0.002 TYR A 235 ARG 0.001 0.000 ARG A 194 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 19 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 16 time to evaluate : 0.138 Fit side-chains outliers start: 3 outliers final: 3 residues processed: 17 average time/residue: 0.1026 time to fit residues: 2.1173 Evaluate side-chains 16 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 13 time to evaluate : 0.129 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0167 time to fit residues: 0.2469 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 14 random chunks: chunk 11 optimal weight: 0.9980 chunk 9 optimal weight: 0.6980 chunk 0 optimal weight: 0.5980 chunk 6 optimal weight: 0.2980 chunk 12 optimal weight: 0.9980 chunk 3 optimal weight: 0.0030 chunk 4 optimal weight: 0.0980 chunk 2 optimal weight: 0.0870 chunk 7 optimal weight: 0.0980 chunk 13 optimal weight: 0.0980 chunk 1 optimal weight: 0.5980 overall best weight: 0.0768 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 199 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7471 moved from start: 0.3946 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.025 1162 Z= 0.171 Angle : 0.633 5.095 1589 Z= 0.322 Chirality : 0.051 0.240 178 Planarity : 0.005 0.045 193 Dihedral : 6.987 39.337 235 Min Nonbonded Distance : 2.430 Molprobity Statistics. All-atom Clashscore : 5.42 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.10 % Favored : 96.90 % Rotamer: Outliers : 0.00 % Allowed : 12.35 % Favored : 87.65 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 14.29 % 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.51 (0.74), residues: 129 helix: 2.27 (1.36), residues: 12 sheet: 0.88 (0.62), residues: 73 loop : -1.14 (0.87), residues: 44 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP A 179 HIS 0.001 0.000 HIS A 189 PHE 0.024 0.002 PHE A 39 TYR 0.005 0.001 TYR A 43 ARG 0.000 0.000 ARG A 194 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 19 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 19 time to evaluate : 0.133 Fit side-chains revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 19 average time/residue: 0.0695 time to fit residues: 1.7138 Evaluate side-chains 16 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 16 time to evaluate : 0.127 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.1772 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 14 random chunks: chunk 4 optimal weight: 0.1980 chunk 7 optimal weight: 1.9990 chunk 13 optimal weight: 0.5980 chunk 1 optimal weight: 0.1980 chunk 9 optimal weight: 0.1980 chunk 11 optimal weight: 0.1980 chunk 8 optimal weight: 0.7980 chunk 6 optimal weight: 0.7980 chunk 5 optimal weight: 0.0980 chunk 2 optimal weight: 0.0980 chunk 10 optimal weight: 0.6980 overall best weight: 0.1580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 199 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7499 moved from start: 0.4477 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.026 1162 Z= 0.190 Angle : 0.575 8.009 1589 Z= 0.293 Chirality : 0.049 0.164 178 Planarity : 0.005 0.038 193 Dihedral : 6.784 44.178 235 Min Nonbonded Distance : 2.462 Molprobity Statistics. All-atom Clashscore : 5.92 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.65 % Favored : 95.35 % Rotamer: Outliers : 0.00 % Allowed : 14.81 % Favored : 85.19 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 14.29 % 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.54 (0.75), residues: 129 helix: 2.41 (1.32), residues: 12 sheet: 0.88 (0.61), residues: 78 loop : -1.40 (0.95), residues: 39 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.004 0.001 TRP A 179 HIS 0.001 0.000 HIS A 157 PHE 0.020 0.002 PHE A 39 TYR 0.006 0.001 TYR A 43 ARG 0.001 0.000 ARG A 154 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 17 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 17 time to evaluate : 0.126 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 17 average time/residue: 0.1300 time to fit residues: 2.5793 Evaluate side-chains 16 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 16 time to evaluate : 0.130 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.2132 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 14 random chunks: chunk 12 optimal weight: 0.9990 chunk 11 optimal weight: 0.5980 chunk 7 optimal weight: 3.9990 chunk 5 optimal weight: 0.6980 chunk 9 optimal weight: 0.8980 chunk 3 optimal weight: 0.5980 chunk 8 optimal weight: 0.1980 chunk 13 optimal weight: 0.4980 chunk 6 optimal weight: 0.9990 chunk 1 optimal weight: 0.3980 chunk 10 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 199 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7626 moved from start: 0.5057 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.034 1162 Z= 0.403 Angle : 0.710 5.443 1589 Z= 0.369 Chirality : 0.058 0.293 178 Planarity : 0.007 0.077 193 Dihedral : 8.599 52.702 235 Min Nonbonded Distance : 2.375 Molprobity Statistics. All-atom Clashscore : 7.40 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.65 % Favored : 95.35 % Rotamer: Outliers : 1.23 % Allowed : 17.28 % Favored : 81.48 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 14.29 % 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.55 (0.71), residues: 129 helix: 0.56 (1.30), residues: 12 sheet: 0.49 (0.59), residues: 77 loop : -2.39 (0.81), residues: 40 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.002 TRP A 179 HIS 0.003 0.001 HIS A 189 PHE 0.011 0.002 PHE A 59 TYR 0.015 0.002 TYR A 43 ARG 0.002 0.001 ARG A 237 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 20 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 19 time to evaluate : 0.128 Fit side-chains outliers start: 1 outliers final: 1 residues processed: 20 average time/residue: 0.1418 time to fit residues: 3.2228 Evaluate side-chains 18 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 17 time to evaluate : 0.124 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0155 time to fit residues: 0.1881 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 14 random chunks: chunk 1 optimal weight: 0.1980 chunk 3 optimal weight: 0.0470 chunk 11 optimal weight: 0.9980 chunk 4 optimal weight: 0.0980 chunk 2 optimal weight: 0.3980 chunk 9 optimal weight: 0.5980 chunk 0 optimal weight: 0.0980 chunk 8 optimal weight: 0.0970 chunk 12 optimal weight: 0.3980 chunk 7 optimal weight: 1.9990 chunk 5 optimal weight: 0.2980 overall best weight: 0.1076 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 199 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7521 moved from start: 0.5345 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.027 1162 Z= 0.186 Angle : 0.655 8.311 1589 Z= 0.353 Chirality : 0.053 0.231 178 Planarity : 0.007 0.058 193 Dihedral : 6.388 36.550 235 Min Nonbonded Distance : 2.477 Molprobity Statistics. All-atom Clashscore : 4.93 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.33 % Favored : 97.67 % Rotamer: Outliers : 0.00 % Allowed : 20.99 % Favored : 79.01 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 14.29 % 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.03 (0.72), residues: 129 helix: 2.83 (1.30), residues: 12 sheet: 0.67 (0.59), residues: 78 loop : -2.47 (0.78), residues: 39 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A 65 HIS 0.001 0.000 HIS A 157 PHE 0.007 0.001 PHE A 39 TYR 0.017 0.001 TYR A 236 ARG 0.001 0.000 ARG A 237 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 17 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 17 time to evaluate : 0.290 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 17 average time/residue: 0.2985 time to fit residues: 5.8583 Evaluate side-chains 15 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 15 time to evaluate : 0.285 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.3558 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 14 random chunks: chunk 5 optimal weight: 0.0980 chunk 13 optimal weight: 0.3980 chunk 0 optimal weight: 0.5980 chunk 2 optimal weight: 0.0010 chunk 1 optimal weight: 0.0970 chunk 11 optimal weight: 0.0980 chunk 3 optimal weight: 0.5980 chunk 9 optimal weight: 0.6980 chunk 4 optimal weight: 0.1980 chunk 7 optimal weight: 2.9990 chunk 10 optimal weight: 0.7980 overall best weight: 0.0984 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 199 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7499 moved from start: 0.5393 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.018 1162 Z= 0.146 Angle : 0.558 4.556 1589 Z= 0.292 Chirality : 0.052 0.223 178 Planarity : 0.006 0.046 193 Dihedral : 7.036 45.948 235 Min Nonbonded Distance : 2.492 Molprobity Statistics. All-atom Clashscore : 5.92 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.10 % Favored : 96.90 % Rotamer: Outliers : 0.00 % Allowed : 22.22 % Favored : 77.78 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 14.29 % 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.29 (0.74), residues: 129 helix: 3.50 (1.21), residues: 12 sheet: 0.81 (0.60), residues: 78 loop : -2.27 (0.81), residues: 39 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP A 65 HIS 0.001 0.000 HIS A 157 PHE 0.006 0.001 PHE A 39 TYR 0.015 0.001 TYR A 236 ARG 0.000 0.000 ARG A 237 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 17 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 17 time to evaluate : 0.158 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 17 average time/residue: 0.1602 time to fit residues: 3.1016 Evaluate side-chains 16 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 16 time to evaluate : 0.123 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.1771 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 14 random chunks: chunk 11 optimal weight: 0.9980 chunk 12 optimal weight: 0.9990 chunk 2 optimal weight: 0.0980 chunk 9 optimal weight: 0.8980 chunk 13 optimal weight: 0.0980 chunk 0 optimal weight: 0.1980 chunk 4 optimal weight: 0.9980 chunk 3 optimal weight: 0.4980 chunk 8 optimal weight: 0.1980 chunk 6 optimal weight: 0.9980 chunk 5 optimal weight: 0.4980 overall best weight: 0.2180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 199 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 240 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7525 moved from start: 0.5579 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.033 1162 Z= 0.266 Angle : 0.940 15.722 1589 Z= 0.505 Chirality : 0.064 0.307 178 Planarity : 0.008 0.070 193 Dihedral : 7.496 41.540 235 Min Nonbonded Distance : 2.225 Molprobity Statistics. All-atom Clashscore : 8.38 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.10 % Favored : 96.90 % Rotamer: Outliers : 1.23 % Allowed : 22.22 % Favored : 76.54 % Cbeta Deviations : 0.73 % Peptide Plane: Cis-proline : 14.29 % 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.72), residues: 129 helix: 3.17 (1.20), residues: 12 sheet: 0.79 (0.59), residues: 78 loop : -2.85 (0.73), residues: 39 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP A 65 HIS 0.002 0.001 HIS A 189 PHE 0.007 0.001 PHE A 202 TYR 0.012 0.001 TYR A 236 ARG 0.000 0.000 ARG A 194 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 17 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 16 time to evaluate : 0.134 Fit side-chains revert: symmetry clash outliers start: 1 outliers final: 1 residues processed: 17 average time/residue: 0.1511 time to fit residues: 2.9440 Evaluate side-chains 17 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 16 time to evaluate : 0.129 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0197 time to fit residues: 0.1892 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 14 random chunks: chunk 7 optimal weight: 1.9990 chunk 4 optimal weight: 0.1980 chunk 5 optimal weight: 0.0870 chunk 11 optimal weight: 0.1980 chunk 3 optimal weight: 0.4980 chunk 0 optimal weight: 0.5980 chunk 10 optimal weight: 0.5980 chunk 13 optimal weight: 0.0570 chunk 8 optimal weight: 0.1980 chunk 9 optimal weight: 0.1980 chunk 2 optimal weight: 0.1980 overall best weight: 0.1476 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 199 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 240 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7461 moved from start: 0.5798 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.025 1162 Z= 0.193 Angle : 0.717 9.636 1589 Z= 0.369 Chirality : 0.054 0.196 178 Planarity : 0.007 0.060 193 Dihedral : 7.187 47.031 235 Min Nonbonded Distance : 2.491 Molprobity Statistics. All-atom Clashscore : 7.89 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.88 % Favored : 96.12 % Rotamer: Outliers : 0.00 % Allowed : 20.99 % Favored : 79.01 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 14.29 % 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.05 (0.73), residues: 129 helix: 3.53 (1.23), residues: 12 sheet: 0.80 (0.60), residues: 77 loop : -2.76 (0.72), residues: 40 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A 65 HIS 0.001 0.000 HIS A 189 PHE 0.003 0.001 PHE A 180 TYR 0.014 0.001 TYR A 236 ARG 0.001 0.000 ARG A 237 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 258 Ramachandran restraints generated. 129 Oldfield, 0 Emsley, 129 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 17 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 17 time to evaluate : 0.124 Fit side-chains revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 17 average time/residue: 0.1485 time to fit residues: 2.8793 Evaluate side-chains 14 residues out of total 133 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 14 time to evaluate : 0.116 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.1535 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 14 random chunks: chunk 4 optimal weight: 0.8980 chunk 3 optimal weight: 0.4980 chunk 8 optimal weight: 0.1980 chunk 5 optimal weight: 0.9980 chunk 12 optimal weight: 0.5980 chunk 7 optimal weight: 0.4980 chunk 9 optimal weight: 0.6980 chunk 0 optimal weight: 0.4980 chunk 6 optimal weight: 0.0980 chunk 2 optimal weight: 0.0370 chunk 10 optimal weight: 0.0980 overall best weight: 0.1858 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 199 ASN ** 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.3828 r_free = 0.3828 target = 0.162551 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 40)----------------| | r_work = 0.3560 r_free = 0.3560 target = 0.139466 restraints weight = 1813.908| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 26)----------------| | r_work = 0.3582 r_free = 0.3582 target = 0.141557 restraints weight = 1249.068| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 42)----------------| | r_work = 0.3607 r_free = 0.3607 target = 0.143781 restraints weight = 981.647| |-----------------------------------------------------------------------------| r_work (final): 0.3595 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7683 moved from start: 0.5863 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.026 1162 Z= 0.225 Angle : 0.753 12.267 1589 Z= 0.385 Chirality : 0.054 0.191 178 Planarity : 0.007 0.055 193 Dihedral : 6.577 38.014 235 Min Nonbonded Distance : 2.421 Molprobity Statistics. All-atom Clashscore : 8.38 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.88 % Favored : 96.12 % Rotamer: Outliers : 0.00 % Allowed : 24.69 % Favored : 75.31 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 14.29 % 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.72), residues: 129 helix: 3.14 (1.21), residues: 12 sheet: 0.78 (0.59), residues: 78 loop : -2.82 (0.74), residues: 39 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.004 0.001 TRP A 179 HIS 0.001 0.000 HIS A 189 PHE 0.006 0.001 PHE A 202 TYR 0.012 0.001 TYR A 236 ARG 0.008 0.001 ARG A 194 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 1196.72 seconds wall clock time: 26 minutes 33.33 seconds (1593.33 seconds total)