Starting phenix.real_space_refine on Fri Nov 15 05:28:30 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8cr9_16805/11_2024/8cr9_16805.cif Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8cr9_16805/11_2024/8cr9_16805.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=4.2 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8cr9_16805/11_2024/8cr9_16805.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8cr9_16805/11_2024/8cr9_16805.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8cr9_16805/11_2024/8cr9_16805.cif" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8cr9_16805/11_2024/8cr9_16805.cif" } resolution = 4.2 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.002 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: 5 Type Number sf(0) Gaussians S 11 5.16 5 C 2730 2.51 5 N 736 2.21 5 O 872 1.98 5 H 4272 0.53 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 20 residue(s): 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5513/modules/chem_data/mon_lib" Total number of atoms: 8621 Number of models: 1 Model: "" Number of chains: 3 Chain: "A" Number of atoms: 3300 Number of conformers: 1 Conformer: "" Number of residues, atoms: 222, 3300 Classifications: {'peptide': 222} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PCIS': 2, 'PTRANS': 10, 'TRANS': 209} Chain: "B" Number of atoms: 3290 Number of conformers: 1 Conformer: "" Number of residues, atoms: 217, 3290 Classifications: {'peptide': 217} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PCIS': 2, 'PTRANS': 9, 'TRANS': 205} Chain: "C" Number of atoms: 2031 Number of conformers: 1 Conformer: "" Number of residues, atoms: 132, 2031 Classifications: {'peptide': 132} Link IDs: {'PCIS': 1, 'PTRANS': 6, 'TRANS': 124} Time building chain proxies: 4.50, per 1000 atoms: 0.52 Number of scatterers: 8621 At special positions: 0 Unit cell: (64.6, 82.45, 105.4, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) S 11 16.00 O 872 8.00 N 736 7.00 C 2730 6.00 H 4272 1.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.00 Amino acid : True - 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: 2.58 Conformation dependent library (CDL) restraints added in 878.6 milliseconds 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1054 Finding SS restraints... Secondary structure from input PDB file: 15 helices and 10 sheets defined 21.7% alpha, 22.8% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.71 Creating SS restraints... Processing helix chain 'A' and resid 45 through 49 removed outlier: 3.918A pdb=" N VAL A 48 " --> pdb=" O LEU A 45 " (cutoff:3.500A) removed outlier: 4.262A pdb=" N ALA A 49 " --> pdb=" O GLU A 46 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 45 through 49' Processing helix chain 'A' and resid 87 through 91 Processing helix chain 'A' and resid 136 through 140 Processing helix chain 'A' and resid 196 through 198 No H-bonds generated for 'chain 'A' and resid 196 through 198' Processing helix chain 'A' and resid 210 through 213 Processing helix chain 'B' and resid 124 through 130 Processing helix chain 'B' and resid 185 through 192 Processing helix chain 'B' and resid 215 through 217 No H-bonds generated for 'chain 'B' and resid 215 through 217' Processing helix chain 'C' and resid 130 through 130 No H-bonds generated for 'chain 'C' and resid 130 through 130' Processing helix chain 'C' and resid 131 through 158 Processing helix chain 'C' and resid 173 through 178 Processing helix chain 'C' and resid 184 through 189 removed outlier: 4.243A pdb=" N ASP C 188 " --> pdb=" O PRO C 184 " (cutoff:3.500A) Processing helix chain 'C' and resid 189 through 197 removed outlier: 3.797A pdb=" N LEU C 195 " --> pdb=" O GLU C 191 " (cutoff:3.500A) Processing helix chain 'C' and resid 206 through 212 Processing helix chain 'C' and resid 235 through 261 removed outlier: 3.587A pdb=" N ARG C 247 " --> pdb=" O THR C 243 " (cutoff:3.500A) Proline residue: C 250 - end of helix Processing sheet with id=AA1, first strand: chain 'A' and resid 11 through 12 Processing sheet with id=AA2, first strand: chain 'A' and resid 18 through 19 removed outlier: 3.717A pdb=" N LEU A 18 " --> pdb=" O LEU A 81 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 34 through 39 Processing sheet with id=AA4, first strand: chain 'A' and resid 129 through 133 removed outlier: 4.841A pdb=" N ALA A 146 " --> pdb=" O VAL A 193 " (cutoff:3.500A) removed outlier: 6.450A pdb=" N VAL A 193 " --> pdb=" O ALA A 146 " (cutoff:3.500A) removed outlier: 5.740A pdb=" N GLY A 148 " --> pdb=" O VAL A 191 " (cutoff:3.500A) removed outlier: 6.527A pdb=" N VAL A 191 " --> pdb=" O GLY A 148 " (cutoff:3.500A) removed outlier: 5.034A pdb=" N LEU A 150 " --> pdb=" O SER A 189 " (cutoff:3.500A) removed outlier: 6.567A pdb=" N SER A 189 " --> pdb=" O LEU A 150 " (cutoff:3.500A) removed outlier: 4.856A pdb=" N LYS A 152 " --> pdb=" O LEU A 187 " (cutoff:3.500A) removed outlier: 6.727A pdb=" N LEU A 187 " --> pdb=" O LYS A 152 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'A' and resid 160 through 163 removed outlier: 3.668A pdb=" N THR A 214 " --> pdb=" O HIS A 209 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'B' and resid 4 through 7 removed outlier: 5.210A pdb=" N ASP B 17 " --> pdb=" O SER B 77 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 45 through 46 Processing sheet with id=AA8, first strand: chain 'B' and resid 117 through 121 removed outlier: 4.296A pdb=" N SER B 134 " --> pdb=" O PHE B 121 " (cutoff:3.500A) removed outlier: 3.781A pdb=" N SER B 177 " --> pdb=" O THR B 167 " (cutoff:3.500A) removed outlier: 3.790A pdb=" N SER B 179 " --> pdb=" O SER B 165 " (cutoff:3.500A) removed outlier: 4.090A pdb=" N SER B 165 " --> pdb=" O SER B 179 " (cutoff:3.500A) removed outlier: 3.503A pdb=" N GLN B 163 " --> pdb=" O THR B 181 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'B' and resid 156 through 157 Processing sheet with id=AB1, first strand: chain 'C' and resid 162 through 163 149 hydrogen bonds defined for protein. 390 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: 3.34 Time building geometry restraints manager: 3.02 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 0.93 - 1.11: 4263 1.11 - 1.29: 731 1.29 - 1.47: 1478 1.47 - 1.64: 2230 1.64 - 1.82: 13 Bond restraints: 8715 Sorted by residual: bond pdb=" N LEU C 261 " pdb=" H LEU C 261 " ideal model delta sigma weight residual 0.860 1.010 -0.150 2.00e-02 2.50e+03 5.66e+01 bond pdb=" CD1 TYR A 94 " pdb=" HD1 TYR A 94 " ideal model delta sigma weight residual 0.930 1.080 -0.150 2.00e-02 2.50e+03 5.66e+01 bond pdb=" CD2 TYR A 95 " pdb=" HD2 TYR A 95 " ideal model delta sigma weight residual 0.930 1.080 -0.150 2.00e-02 2.50e+03 5.65e+01 bond pdb=" CE1 TYR B 176 " pdb=" HE1 TYR B 176 " ideal model delta sigma weight residual 0.930 1.080 -0.150 2.00e-02 2.50e+03 5.65e+01 bond pdb=" CE1 TYR A 94 " pdb=" HE1 TYR A 94 " ideal model delta sigma weight residual 0.930 1.080 -0.150 2.00e-02 2.50e+03 5.65e+01 ... (remaining 8710 not shown) Histogram of bond angle deviations from ideal: 0.00 - 2.60: 13177 2.60 - 5.21: 1945 5.21 - 7.81: 548 7.81 - 10.41: 50 10.41 - 13.02: 11 Bond angle restraints: 15731 Sorted by residual: angle pdb=" N VAL B 83 " pdb=" CA VAL B 83 " pdb=" C VAL B 83 " ideal model delta sigma weight residual 112.17 119.58 -7.41 9.50e-01 1.11e+00 6.08e+01 angle pdb=" C TRP A 47 " pdb=" N VAL A 48 " pdb=" CA VAL A 48 " ideal model delta sigma weight residual 120.60 129.62 -9.02 1.29e+00 6.01e-01 4.89e+01 angle pdb=" CA PHE B 101 " pdb=" CB PHE B 101 " pdb=" CG PHE B 101 " ideal model delta sigma weight residual 113.80 120.28 -6.48 1.00e+00 1.00e+00 4.20e+01 angle pdb=" N TYR B 94 " pdb=" CA TYR B 94 " pdb=" C TYR B 94 " ideal model delta sigma weight residual 108.60 116.67 -8.07 1.46e+00 4.69e-01 3.06e+01 angle pdb=" N VAL A 48 " pdb=" CA VAL A 48 " pdb=" C VAL A 48 " ideal model delta sigma weight residual 110.36 116.06 -5.70 1.05e+00 9.07e-01 2.95e+01 ... (remaining 15726 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.70: 3766 17.70 - 35.41: 277 35.41 - 53.11: 90 53.11 - 70.81: 17 70.81 - 88.51: 8 Dihedral angle restraints: 4158 sinusoidal: 2253 harmonic: 1905 Sorted by residual: dihedral pdb=" CA ASP A 110 " pdb=" C ASP A 110 " pdb=" N SER A 111 " pdb=" CA SER A 111 " ideal model delta harmonic sigma weight residual -180.00 -150.82 -29.18 0 5.00e+00 4.00e-02 3.41e+01 dihedral pdb=" CA CYS B 88 " pdb=" C CYS B 88 " pdb=" N GLN B 89 " pdb=" CA GLN B 89 " ideal model delta harmonic sigma weight residual 180.00 150.84 29.16 0 5.00e+00 4.00e-02 3.40e+01 dihedral pdb=" CA PHE B 142 " pdb=" C PHE B 142 " pdb=" N TYR B 143 " pdb=" CA TYR B 143 " ideal model delta harmonic sigma weight residual 180.00 150.85 29.15 0 5.00e+00 4.00e-02 3.40e+01 ... (remaining 4155 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.087: 401 0.087 - 0.174: 234 0.174 - 0.262: 36 0.262 - 0.349: 7 0.349 - 0.436: 3 Chirality restraints: 681 Sorted by residual: chirality pdb=" CA HIS B 96 " pdb=" N HIS B 96 " pdb=" C HIS B 96 " pdb=" CB HIS B 96 " both_signs ideal model delta sigma weight residual False 2.51 2.07 0.44 2.00e-01 2.50e+01 4.75e+00 chirality pdb=" CA LYS B 50 " pdb=" N LYS B 50 " pdb=" C LYS B 50 " pdb=" CB LYS B 50 " both_signs ideal model delta sigma weight residual False 2.51 2.11 0.40 2.00e-01 2.50e+01 4.09e+00 chirality pdb=" CA LYS C 258 " pdb=" N LYS C 258 " pdb=" C LYS C 258 " pdb=" CB LYS C 258 " both_signs ideal model delta sigma weight residual False 2.51 2.15 0.36 2.00e-01 2.50e+01 3.30e+00 ... (remaining 678 not shown) Planarity restraints: 1317 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB TYR B 90 " 0.203 2.00e-02 2.50e+03 7.80e-02 1.83e+02 pdb=" CG TYR B 90 " -0.071 2.00e-02 2.50e+03 pdb=" CD1 TYR B 90 " -0.072 2.00e-02 2.50e+03 pdb=" CD2 TYR B 90 " -0.068 2.00e-02 2.50e+03 pdb=" CE1 TYR B 90 " -0.029 2.00e-02 2.50e+03 pdb=" CE2 TYR B 90 " -0.030 2.00e-02 2.50e+03 pdb=" CZ TYR B 90 " 0.013 2.00e-02 2.50e+03 pdb=" OH TYR B 90 " 0.113 2.00e-02 2.50e+03 pdb=" HD1 TYR B 90 " -0.036 2.00e-02 2.50e+03 pdb=" HD2 TYR B 90 " -0.035 2.00e-02 2.50e+03 pdb=" HE1 TYR B 90 " 0.007 2.00e-02 2.50e+03 pdb=" HE2 TYR B 90 " 0.005 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TYR A 59 " 0.195 2.00e-02 2.50e+03 7.42e-02 1.65e+02 pdb=" CG TYR A 59 " -0.062 2.00e-02 2.50e+03 pdb=" CD1 TYR A 59 " -0.066 2.00e-02 2.50e+03 pdb=" CD2 TYR A 59 " -0.065 2.00e-02 2.50e+03 pdb=" CE1 TYR A 59 " -0.021 2.00e-02 2.50e+03 pdb=" CE2 TYR A 59 " -0.030 2.00e-02 2.50e+03 pdb=" CZ TYR A 59 " 0.015 2.00e-02 2.50e+03 pdb=" OH TYR A 59 " 0.104 2.00e-02 2.50e+03 pdb=" HD1 TYR A 59 " -0.050 2.00e-02 2.50e+03 pdb=" HD2 TYR A 59 " -0.031 2.00e-02 2.50e+03 pdb=" HE1 TYR A 59 " 0.011 2.00e-02 2.50e+03 pdb=" HE2 TYR A 59 " -0.000 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TRP B 32 " -0.139 2.00e-02 2.50e+03 6.39e-02 1.64e+02 pdb=" CG TRP B 32 " -0.048 2.00e-02 2.50e+03 pdb=" CD1 TRP B 32 " 0.033 2.00e-02 2.50e+03 pdb=" CD2 TRP B 32 " -0.015 2.00e-02 2.50e+03 pdb=" NE1 TRP B 32 " 0.039 2.00e-02 2.50e+03 pdb=" CE2 TRP B 32 " -0.017 2.00e-02 2.50e+03 pdb=" CE3 TRP B 32 " 0.016 2.00e-02 2.50e+03 pdb=" CZ2 TRP B 32 " -0.055 2.00e-02 2.50e+03 pdb=" CZ3 TRP B 32 " 0.027 2.00e-02 2.50e+03 pdb=" CH2 TRP B 32 " -0.045 2.00e-02 2.50e+03 pdb=" HD1 TRP B 32 " 0.064 2.00e-02 2.50e+03 pdb=" HE1 TRP B 32 " 0.089 2.00e-02 2.50e+03 pdb=" HE3 TRP B 32 " 0.048 2.00e-02 2.50e+03 pdb=" HZ2 TRP B 32 " -0.068 2.00e-02 2.50e+03 pdb=" HZ3 TRP B 32 " 0.117 2.00e-02 2.50e+03 pdb=" HH2 TRP B 32 " -0.046 2.00e-02 2.50e+03 ... (remaining 1314 not shown) Histogram of nonbonded interaction distances: 1.59 - 2.19: 627 2.19 - 2.80: 15346 2.80 - 3.40: 23113 3.40 - 4.00: 29818 4.00 - 4.60: 43383 Nonbonded interactions: 112287 Sorted by model distance: nonbonded pdb=" H ILE B 29 " pdb=" H GLY B 30 " model vdw 1.593 2.100 nonbonded pdb=" HG SER B 185 " pdb=" OD2 ASP B 188 " model vdw 1.615 2.450 nonbonded pdb=" HG SER C 206 " pdb=" OD2 ASP C 209 " model vdw 1.629 2.450 nonbonded pdb=" HG SER B 31 " pdb=" OE1 GLU B 66 " model vdw 1.629 2.450 nonbonded pdb=" HG SER A 186 " pdb=" OE1 GLN B 163 " model vdw 1.654 2.450 ... (remaining 112282 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.00 ========== 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: 0.000 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.000 Extract box with map and model: 0.000 Check model and map are aligned: 0.000 Set scattering table: 0.000 Process input model: 25.760 Find NCS groups from input model: 0.080 Set up NCS constraints: 0.030 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.010 Load rotamer database and sin/cos tables:3.130 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 29.010 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6233 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.013 0.047 4443 Z= 0.859 Angle : 1.974 10.053 6039 Z= 1.288 Chirality : 0.103 0.436 681 Planarity : 0.016 0.116 779 Dihedral : 14.418 88.513 1608 Min Nonbonded Distance : 2.568 Molprobity Statistics. All-atom Clashscore : 0.46 Ramachandran Plot: Outliers : 1.42 % Allowed : 9.56 % Favored : 89.03 % Rotamer: Outliers : 2.87 % Allowed : 3.69 % Favored : 93.44 % Cbeta Deviations : 0.57 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.67 (0.32), residues: 565 helix: -1.12 (0.47), residues: 91 sheet: -0.78 (0.46), residues: 119 loop : -2.44 (0.30), residues: 355 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.097 0.026 TRP A 36 HIS 0.017 0.004 HIS A 209 PHE 0.061 0.014 PHE B 212 TYR 0.182 0.034 TYR B 90 ARG 0.009 0.001 ARG B 61 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 151 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 14 poor density : 137 time to evaluate : 0.741 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 17 SER cc_start: 0.8297 (p) cc_final: 0.8090 (m) REVERT: A 163 TRP cc_start: 0.6230 (m100) cc_final: 0.5555 (m100) REVERT: B 71 PHE cc_start: 0.6120 (t80) cc_final: 0.5349 (t80) REVERT: B 118 VAL cc_start: 0.8393 (t) cc_final: 0.8186 (p) REVERT: B 120 ILE cc_start: 0.3867 (mt) cc_final: 0.3626 (mt) outliers start: 14 outliers final: 8 residues processed: 147 average time/residue: 0.3887 time to fit residues: 70.4902 Evaluate side-chains 110 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 8 poor density : 102 time to evaluate : 0.720 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 56 random chunks: chunk 47 optimal weight: 3.9990 chunk 42 optimal weight: 5.9990 chunk 23 optimal weight: 2.9990 chunk 14 optimal weight: 5.9990 chunk 28 optimal weight: 5.9990 chunk 22 optimal weight: 5.9990 chunk 43 optimal weight: 6.9990 chunk 16 optimal weight: 6.9990 chunk 26 optimal weight: 0.9990 chunk 32 optimal weight: 0.9990 chunk 50 optimal weight: 0.9990 overall best weight: 1.9990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 77 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6190 moved from start: 0.3658 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.037 4443 Z= 0.258 Angle : 0.806 6.875 6039 Z= 0.436 Chirality : 0.047 0.192 681 Planarity : 0.006 0.046 779 Dihedral : 8.788 40.125 611 Min Nonbonded Distance : 2.387 Molprobity Statistics. All-atom Clashscore : 7.08 Ramachandran Plot: Outliers : 0.71 % Allowed : 10.44 % Favored : 88.85 % Rotamer: Outliers : 0.00 % Allowed : 2.25 % Favored : 97.75 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.56 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.52 (0.34), residues: 565 helix: -0.21 (0.53), residues: 98 sheet: -1.35 (0.42), residues: 147 loop : -2.39 (0.32), residues: 320 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.002 TRP A 105 HIS 0.005 0.002 HIS A 173 PHE 0.015 0.002 PHE B 142 TYR 0.031 0.003 TYR B 90 ARG 0.006 0.001 ARG A 44 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 112 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 112 time to evaluate : 0.718 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 163 TRP cc_start: 0.5938 (m100) cc_final: 0.5663 (m100) REVERT: B 71 PHE cc_start: 0.5887 (t80) cc_final: 0.5447 (t80) REVERT: B 200 THR cc_start: 0.8701 (m) cc_final: 0.8415 (p) REVERT: C 220 GLU cc_start: 0.6653 (tp30) cc_final: 0.6198 (pm20) outliers start: 0 outliers final: 0 residues processed: 112 average time/residue: 0.3951 time to fit residues: 54.4794 Evaluate side-chains 95 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 95 time to evaluate : 0.714 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 56 random chunks: chunk 28 optimal weight: 8.9990 chunk 15 optimal weight: 4.9990 chunk 42 optimal weight: 0.8980 chunk 34 optimal weight: 5.9990 chunk 14 optimal weight: 10.0000 chunk 50 optimal weight: 3.9990 chunk 55 optimal weight: 8.9990 chunk 45 optimal weight: 8.9990 chunk 17 optimal weight: 4.9990 chunk 40 optimal weight: 7.9990 chunk 38 optimal weight: 0.9990 overall best weight: 3.1788 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6318 moved from start: 0.4448 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.054 4443 Z= 0.296 Angle : 0.774 6.633 6039 Z= 0.420 Chirality : 0.046 0.159 681 Planarity : 0.007 0.053 779 Dihedral : 8.193 37.220 611 Min Nonbonded Distance : 2.194 Molprobity Statistics. All-atom Clashscore : 8.00 Ramachandran Plot: Outliers : 0.53 % Allowed : 13.98 % Favored : 85.49 % Rotamer: Outliers : 0.00 % Allowed : 2.66 % Favored : 97.34 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.19 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.92 (0.34), residues: 565 helix: -0.14 (0.52), residues: 99 sheet: -1.56 (0.41), residues: 151 loop : -2.84 (0.32), residues: 315 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.015 0.002 TRP C 186 HIS 0.008 0.002 HIS A 173 PHE 0.014 0.002 PHE B 142 TYR 0.034 0.003 TYR B 90 ARG 0.008 0.001 ARG A 44 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 101 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 101 time to evaluate : 0.886 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: B 33 LEU cc_start: 0.7403 (pp) cc_final: 0.7116 (pp) REVERT: B 108 ASP cc_start: 0.7689 (t0) cc_final: 0.7269 (m-30) REVERT: B 148 LYS cc_start: 0.7154 (tppt) cc_final: 0.6953 (tppt) REVERT: B 200 THR cc_start: 0.8734 (m) cc_final: 0.8304 (p) outliers start: 0 outliers final: 0 residues processed: 101 average time/residue: 0.4132 time to fit residues: 50.6741 Evaluate side-chains 93 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 93 time to evaluate : 0.738 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 56 random chunks: chunk 26 optimal weight: 0.3980 chunk 5 optimal weight: 3.9990 chunk 24 optimal weight: 9.9990 chunk 34 optimal weight: 5.9990 chunk 51 optimal weight: 2.9990 chunk 54 optimal weight: 0.6980 chunk 48 optimal weight: 0.8980 chunk 14 optimal weight: 9.9990 chunk 45 optimal weight: 2.9990 chunk 30 optimal weight: 2.9990 chunk 0 optimal weight: 8.9990 overall best weight: 1.5984 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6215 moved from start: 0.4731 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.030 4443 Z= 0.199 Angle : 0.656 5.943 6039 Z= 0.352 Chirality : 0.043 0.138 681 Planarity : 0.006 0.060 779 Dihedral : 7.425 37.809 611 Min Nonbonded Distance : 2.201 Molprobity Statistics. All-atom Clashscore : 5.22 Ramachandran Plot: Outliers : 0.35 % Allowed : 12.21 % Favored : 87.43 % Rotamer: Outliers : 0.00 % Allowed : 3.07 % Favored : 96.93 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.19 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.78 (0.33), residues: 565 helix: 0.19 (0.53), residues: 99 sheet: -1.55 (0.40), residues: 151 loop : -2.79 (0.32), residues: 315 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP B 151 HIS 0.004 0.001 HIS A 173 PHE 0.040 0.002 PHE B 142 TYR 0.023 0.002 TYR B 90 ARG 0.004 0.001 ARG A 44 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 100 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 100 time to evaluate : 0.720 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 3 GLN cc_start: 0.7710 (mt0) cc_final: 0.7207 (mt0) REVERT: A 34 MET cc_start: 0.7974 (mtp) cc_final: 0.7710 (ppp) REVERT: B 108 ASP cc_start: 0.7762 (t0) cc_final: 0.7087 (m-30) REVERT: B 200 THR cc_start: 0.8658 (m) cc_final: 0.8385 (p) outliers start: 0 outliers final: 0 residues processed: 100 average time/residue: 0.4183 time to fit residues: 51.0581 Evaluate side-chains 92 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 92 time to evaluate : 0.758 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 56 random chunks: chunk 40 optimal weight: 9.9990 chunk 22 optimal weight: 5.9990 chunk 46 optimal weight: 5.9990 chunk 37 optimal weight: 10.0000 chunk 0 optimal weight: 9.9990 chunk 27 optimal weight: 8.9990 chunk 48 optimal weight: 2.9990 chunk 13 optimal weight: 7.9990 chunk 18 optimal weight: 7.9990 chunk 10 optimal weight: 7.9990 chunk 31 optimal weight: 1.9990 overall best weight: 4.9990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 77 ASN B 92 ASN ** B 96 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6483 moved from start: 0.5124 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.007 0.058 4443 Z= 0.414 Angle : 0.845 7.852 6039 Z= 0.457 Chirality : 0.049 0.173 681 Planarity : 0.007 0.057 779 Dihedral : 8.383 36.184 611 Min Nonbonded Distance : 2.139 Molprobity Statistics. All-atom Clashscore : 14.27 Ramachandran Plot: Outliers : 0.71 % Allowed : 17.17 % Favored : 82.12 % Rotamer: Outliers : 0.00 % Allowed : 2.66 % Favored : 97.34 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.19 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.49 (0.32), residues: 565 helix: -0.17 (0.51), residues: 99 sheet: -1.86 (0.42), residues: 138 loop : -3.40 (0.29), residues: 328 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.026 0.003 TRP A 36 HIS 0.003 0.002 HIS B 96 PHE 0.022 0.003 PHE B 142 TYR 0.040 0.004 TYR B 90 ARG 0.009 0.001 ARG A 44 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 95 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 95 time to evaluate : 0.707 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 3 GLN cc_start: 0.8056 (mt0) cc_final: 0.7545 (mt0) REVERT: A 163 TRP cc_start: 0.6672 (m100) cc_final: 0.6073 (m100) REVERT: B 33 LEU cc_start: 0.7417 (pp) cc_final: 0.6972 (pp) REVERT: B 106 LYS cc_start: 0.7552 (mttt) cc_final: 0.6957 (mptt) REVERT: B 108 ASP cc_start: 0.7768 (t0) cc_final: 0.7274 (m-30) REVERT: B 200 THR cc_start: 0.8810 (m) cc_final: 0.8439 (p) outliers start: 0 outliers final: 0 residues processed: 95 average time/residue: 0.3828 time to fit residues: 44.9810 Evaluate side-chains 86 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 86 time to evaluate : 0.736 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 56 random chunks: chunk 13 optimal weight: 7.9990 chunk 54 optimal weight: 8.9990 chunk 44 optimal weight: 0.7980 chunk 25 optimal weight: 1.9990 chunk 4 optimal weight: 3.9990 chunk 17 optimal weight: 3.9990 chunk 28 optimal weight: 3.9990 chunk 52 optimal weight: 5.9990 chunk 6 optimal weight: 0.5980 chunk 30 optimal weight: 0.9980 chunk 39 optimal weight: 0.9990 overall best weight: 1.0784 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 96 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6286 moved from start: 0.5481 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.026 4443 Z= 0.183 Angle : 0.639 5.478 6039 Z= 0.342 Chirality : 0.043 0.138 681 Planarity : 0.006 0.059 779 Dihedral : 7.377 37.352 611 Min Nonbonded Distance : 2.279 Molprobity Statistics. All-atom Clashscore : 7.42 Ramachandran Plot: Outliers : 0.35 % Allowed : 10.80 % Favored : 88.85 % Rotamer: Outliers : 0.00 % Allowed : 1.02 % Favored : 98.98 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.19 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.22 (0.33), residues: 565 helix: 0.17 (0.52), residues: 99 sheet: -1.80 (0.39), residues: 151 loop : -3.23 (0.30), residues: 315 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.002 TRP B 32 HIS 0.003 0.001 HIS B 201 PHE 0.026 0.002 PHE B 142 TYR 0.035 0.003 TYR C 177 ARG 0.003 0.000 ARG A 44 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 95 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 95 time to evaluate : 0.700 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 3 GLN cc_start: 0.7833 (mt0) cc_final: 0.7247 (mt0) REVERT: A 163 TRP cc_start: 0.6628 (m100) cc_final: 0.5960 (m100) REVERT: B 33 LEU cc_start: 0.7167 (pp) cc_final: 0.6834 (pp) REVERT: B 106 LYS cc_start: 0.7352 (mttt) cc_final: 0.6775 (mptt) REVERT: B 108 ASP cc_start: 0.7681 (t0) cc_final: 0.7127 (m-30) REVERT: B 148 LYS cc_start: 0.7106 (tppt) cc_final: 0.6804 (tppt) REVERT: B 200 THR cc_start: 0.8759 (m) cc_final: 0.8435 (p) REVERT: C 254 LYS cc_start: 0.6248 (mttt) cc_final: 0.5967 (tttt) outliers start: 0 outliers final: 0 residues processed: 95 average time/residue: 0.4012 time to fit residues: 46.8369 Evaluate side-chains 78 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 78 time to evaluate : 0.713 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 56 random chunks: chunk 30 optimal weight: 0.5980 chunk 45 optimal weight: 5.9990 chunk 53 optimal weight: 5.9990 chunk 33 optimal weight: 2.9990 chunk 32 optimal weight: 0.6980 chunk 24 optimal weight: 8.9990 chunk 21 optimal weight: 4.9990 chunk 16 optimal weight: 9.9990 chunk 10 optimal weight: 0.0170 chunk 34 optimal weight: 8.9990 chunk 36 optimal weight: 0.6980 overall best weight: 1.0020 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 173 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6242 moved from start: 0.5656 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.024 4443 Z= 0.171 Angle : 0.586 5.071 6039 Z= 0.312 Chirality : 0.042 0.135 681 Planarity : 0.006 0.054 779 Dihedral : 6.783 35.117 611 Min Nonbonded Distance : 2.319 Molprobity Statistics. All-atom Clashscore : 6.38 Ramachandran Plot: Outliers : 0.35 % Allowed : 13.45 % Favored : 86.19 % Rotamer: Outliers : 0.00 % Allowed : 0.82 % Favored : 99.18 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.19 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.97 (0.33), residues: 565 helix: 0.49 (0.53), residues: 99 sheet: -1.76 (0.39), residues: 153 loop : -3.04 (0.32), residues: 313 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP C 186 HIS 0.006 0.001 HIS A 173 PHE 0.019 0.001 PHE B 142 TYR 0.017 0.002 TYR B 90 ARG 0.002 0.000 ARG A 44 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 91 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 91 time to evaluate : 0.857 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 163 TRP cc_start: 0.6583 (m100) cc_final: 0.6051 (m100) REVERT: B 33 LEU cc_start: 0.7153 (pp) cc_final: 0.6795 (pp) REVERT: B 108 ASP cc_start: 0.7632 (t0) cc_final: 0.6953 (m-30) REVERT: B 148 LYS cc_start: 0.7074 (tppt) cc_final: 0.6775 (tppt) REVERT: B 200 THR cc_start: 0.8694 (m) cc_final: 0.8352 (p) outliers start: 0 outliers final: 0 residues processed: 91 average time/residue: 0.4119 time to fit residues: 45.5879 Evaluate side-chains 75 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 75 time to evaluate : 0.700 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 56 random chunks: chunk 26 optimal weight: 1.9990 chunk 5 optimal weight: 2.9990 chunk 42 optimal weight: 9.9990 chunk 49 optimal weight: 3.9990 chunk 51 optimal weight: 3.9990 chunk 47 optimal weight: 0.6980 chunk 50 optimal weight: 0.5980 chunk 30 optimal weight: 3.9990 chunk 21 optimal weight: 5.9990 chunk 39 optimal weight: 5.9990 chunk 15 optimal weight: 1.9990 overall best weight: 1.6586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 173 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6273 moved from start: 0.5855 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.027 4443 Z= 0.202 Angle : 0.603 5.428 6039 Z= 0.322 Chirality : 0.042 0.134 681 Planarity : 0.006 0.050 779 Dihedral : 6.750 34.510 611 Min Nonbonded Distance : 2.256 Molprobity Statistics. All-atom Clashscore : 6.15 Ramachandran Plot: Outliers : 0.35 % Allowed : 13.81 % Favored : 85.84 % Rotamer: Outliers : 0.00 % Allowed : 0.41 % Favored : 99.59 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.19 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.00 (0.33), residues: 565 helix: 0.49 (0.53), residues: 99 sheet: -1.75 (0.39), residues: 153 loop : -3.08 (0.31), residues: 313 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP C 186 HIS 0.007 0.001 HIS A 173 PHE 0.010 0.001 PHE B 142 TYR 0.024 0.002 TYR A 94 ARG 0.003 0.000 ARG A 44 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 87 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 87 time to evaluate : 0.742 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 163 TRP cc_start: 0.6661 (m100) cc_final: 0.6098 (m100) REVERT: B 33 LEU cc_start: 0.7197 (pp) cc_final: 0.6872 (pp) REVERT: B 108 ASP cc_start: 0.7670 (t0) cc_final: 0.7094 (m-30) REVERT: B 148 LYS cc_start: 0.7092 (tppt) cc_final: 0.6860 (tppt) REVERT: B 200 THR cc_start: 0.8683 (m) cc_final: 0.8361 (p) outliers start: 0 outliers final: 0 residues processed: 87 average time/residue: 0.4270 time to fit residues: 46.1014 Evaluate side-chains 75 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 75 time to evaluate : 0.748 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 56 random chunks: chunk 45 optimal weight: 10.0000 chunk 47 optimal weight: 0.5980 chunk 50 optimal weight: 1.9990 chunk 33 optimal weight: 0.9980 chunk 53 optimal weight: 0.9980 chunk 32 optimal weight: 1.9990 chunk 25 optimal weight: 4.9990 chunk 36 optimal weight: 2.9990 chunk 55 optimal weight: 0.8980 chunk 51 optimal weight: 0.9990 chunk 44 optimal weight: 0.6980 overall best weight: 0.8380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 173 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6206 moved from start: 0.5985 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.029 4443 Z= 0.167 Angle : 0.571 4.810 6039 Z= 0.305 Chirality : 0.042 0.133 681 Planarity : 0.006 0.052 779 Dihedral : 6.447 34.091 611 Min Nonbonded Distance : 2.303 Molprobity Statistics. All-atom Clashscore : 6.15 Ramachandran Plot: Outliers : 0.35 % Allowed : 11.50 % Favored : 88.14 % Rotamer: Outliers : 0.00 % Allowed : 0.61 % Favored : 99.39 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.19 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.89 (0.33), residues: 565 helix: 0.60 (0.54), residues: 99 sheet: -1.70 (0.38), residues: 158 loop : -3.02 (0.32), residues: 308 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP C 186 HIS 0.007 0.001 HIS A 173 PHE 0.013 0.001 PHE B 142 TYR 0.021 0.002 TYR A 94 ARG 0.002 0.000 ARG A 44 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 88 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 88 time to evaluate : 0.736 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 34 MET cc_start: 0.8254 (mtp) cc_final: 0.7831 (ppp) REVERT: A 163 TRP cc_start: 0.6568 (m100) cc_final: 0.5964 (m100) REVERT: B 33 LEU cc_start: 0.7033 (pp) cc_final: 0.6737 (pp) REVERT: B 108 ASP cc_start: 0.7612 (t0) cc_final: 0.7055 (m-30) REVERT: B 148 LYS cc_start: 0.7090 (tppt) cc_final: 0.6825 (tppt) REVERT: B 200 THR cc_start: 0.8611 (m) cc_final: 0.8299 (p) REVERT: C 139 THR cc_start: 0.6443 (m) cc_final: 0.6211 (p) outliers start: 0 outliers final: 0 residues processed: 88 average time/residue: 0.4117 time to fit residues: 44.4054 Evaluate side-chains 74 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 74 time to evaluate : 0.720 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 56 random chunks: chunk 4 optimal weight: 8.9990 chunk 34 optimal weight: 4.9990 chunk 27 optimal weight: 3.9990 chunk 35 optimal weight: 0.9980 chunk 47 optimal weight: 2.9990 chunk 13 optimal weight: 0.9990 chunk 40 optimal weight: 0.9980 chunk 6 optimal weight: 5.9990 chunk 12 optimal weight: 9.9990 chunk 44 optimal weight: 0.8980 chunk 18 optimal weight: 8.9990 overall best weight: 1.3784 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 173 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6238 moved from start: 0.6061 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.023 4443 Z= 0.185 Angle : 0.579 5.140 6039 Z= 0.309 Chirality : 0.042 0.132 681 Planarity : 0.006 0.052 779 Dihedral : 6.413 34.193 611 Min Nonbonded Distance : 2.276 Molprobity Statistics. All-atom Clashscore : 6.84 Ramachandran Plot: Outliers : 0.35 % Allowed : 12.92 % Favored : 86.73 % Rotamer: Outliers : 0.00 % Allowed : 0.61 % Favored : 99.39 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.19 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.87 (0.33), residues: 565 helix: 0.58 (0.54), residues: 99 sheet: -1.61 (0.38), residues: 153 loop : -3.03 (0.32), residues: 313 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP C 186 HIS 0.007 0.001 HIS A 173 PHE 0.009 0.001 PHE B 142 TYR 0.025 0.002 TYR B 90 ARG 0.003 0.000 ARG A 44 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1130 Ramachandran restraints generated. 565 Oldfield, 0 Emsley, 565 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 89 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 89 time to evaluate : 0.687 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 34 MET cc_start: 0.8284 (mtp) cc_final: 0.7774 (ppp) REVERT: A 163 TRP cc_start: 0.6599 (m100) cc_final: 0.6027 (m100) REVERT: B 33 LEU cc_start: 0.7115 (pp) cc_final: 0.6819 (pp) REVERT: B 108 ASP cc_start: 0.7646 (t0) cc_final: 0.7050 (m-30) REVERT: B 148 LYS cc_start: 0.7078 (tppt) cc_final: 0.6811 (tppt) REVERT: B 200 THR cc_start: 0.8654 (m) cc_final: 0.8338 (p) REVERT: C 139 THR cc_start: 0.6479 (m) cc_final: 0.6217 (p) outliers start: 0 outliers final: 0 residues processed: 89 average time/residue: 0.3544 time to fit residues: 39.7731 Evaluate side-chains 77 residues out of total 488 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 77 time to evaluate : 0.712 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 56 random chunks: chunk 45 optimal weight: 9.9990 chunk 5 optimal weight: 3.9990 chunk 8 optimal weight: 2.9990 chunk 39 optimal weight: 3.9990 chunk 2 optimal weight: 4.9990 chunk 32 optimal weight: 3.9990 chunk 50 optimal weight: 0.7980 chunk 29 optimal weight: 2.9990 chunk 38 optimal weight: 0.9990 chunk 1 optimal weight: 0.7980 chunk 35 optimal weight: 4.9990 overall best weight: 1.7186 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 173 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4527 r_free = 0.4527 target = 0.220298 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 39)----------------| | r_work = 0.4310 r_free = 0.4310 target = 0.194882 restraints weight = 14612.824| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 36)----------------| | r_work = 0.4319 r_free = 0.4319 target = 0.196236 restraints weight = 12278.744| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 28)----------------| | r_work = 0.4346 r_free = 0.4346 target = 0.199706 restraints weight = 11095.554| |-----------------------------------------------------------------------------| r_work (final): 0.4344 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5909 moved from start: 0.6139 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.023 4443 Z= 0.200 Angle : 0.594 5.707 6039 Z= 0.316 Chirality : 0.042 0.133 681 Planarity : 0.007 0.097 779 Dihedral : 6.508 34.721 611 Min Nonbonded Distance : 2.262 Molprobity Statistics. All-atom Clashscore : 7.08 Ramachandran Plot: Outliers : 0.35 % Allowed : 13.45 % Favored : 86.19 % Rotamer: Outliers : 0.00 % Allowed : 0.41 % Favored : 99.59 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 16.67 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.19 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.90 (0.33), residues: 565 helix: 0.75 (0.54), residues: 97 sheet: -1.58 (0.38), residues: 153 loop : -3.13 (0.31), residues: 315 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.016 0.002 TRP C 186 HIS 0.007 0.002 HIS A 173 PHE 0.008 0.001 PHE B 142 TYR 0.022 0.002 TYR A 94 ARG 0.003 0.000 ARG A 44 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 2657.66 seconds wall clock time: 48 minutes 19.56 seconds (2899.56 seconds total)