Starting phenix.real_space_refine on Fri Apr 5 22:49:41 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8u4r_41892/04_2024/8u4r_41892_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8u4r_41892/04_2024/8u4r_41892.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.1 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8u4r_41892/04_2024/8u4r_41892.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8u4r_41892/04_2024/8u4r_41892.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8u4r_41892/04_2024/8u4r_41892_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8u4r_41892/04_2024/8u4r_41892_updated.pdb" } resolution = 3.1 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs 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.049 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Anisotropic ADP refinement not supported. Converting 2008 atoms to isotropic. 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 23 5.16 5 C 2680 2.51 5 N 676 2.21 5 O 728 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "R TYR 121": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5288/modules/chem_data/mon_lib" Total number of atoms: 4107 Number of models: 1 Model: "" Number of chains: 4 Chain: "L" Number of atoms: 864 Number of conformers: 1 Conformer: "" Number of residues, atoms: 111, 864 Classifications: {'peptide': 111} Link IDs: {'PCIS': 2, 'PTRANS': 5, 'TRANS': 103} Chain: "H" Number of atoms: 976 Number of conformers: 1 Conformer: "" Number of residues, atoms: 125, 976 Classifications: {'peptide': 125} Link IDs: {'PTRANS': 2, 'TRANS': 122} Chain: "R" Number of atoms: 2239 Number of conformers: 1 Conformer: "" Number of residues, atoms: 277, 2239 Classifications: {'peptide': 277} Link IDs: {'PTRANS': 10, 'TRANS': 266} Chain breaks: 2 Chain: "R" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 28 Unusual residues: {'CLR': 1} Classifications: {'undetermined': 1} Time building chain proxies: 2.75, per 1000 atoms: 0.67 Number of scatterers: 4107 At special positions: 0 Unit cell: (55.25, 65.45, 110.5, 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 23 16.00 O 728 8.00 N 676 7.00 C 2680 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=4, symmetry=0 Simple disulfide: pdb=" SG CYS L 23 " - pdb=" SG CYS L 93 " distance=2.04 Simple disulfide: pdb=" SG CYS H 22 " - pdb=" SG CYS H 96 " distance=2.03 Simple disulfide: pdb=" SG CYS R 28 " - pdb=" SG CYS R 274 " distance=2.03 Simple disulfide: pdb=" SG CYS R 109 " - pdb=" SG CYS R 186 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 1.73 Conformation dependent library (CDL) restraints added in 788.8 milliseconds 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 954 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 13 helices and 5 sheets defined 41.9% alpha, 18.9% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 12.14 Creating SS restraints... Processing helix chain 'H' and resid 29 through 31 No H-bonds generated for 'chain 'H' and resid 29 through 31' Processing helix chain 'H' and resid 88 through 90 No H-bonds generated for 'chain 'H' and resid 88 through 90' Processing helix chain 'R' and resid 37 through 64 Proline residue: R 42 - end of helix Processing helix chain 'R' and resid 72 through 99 removed outlier: 3.623A pdb=" N LEU R 86 " --> pdb=" O VAL R 82 " (cutoff:3.500A) removed outlier: 3.866A pdb=" N PHE R 87 " --> pdb=" O ALA R 83 " (cutoff:3.500A) Proline residue: R 92 - end of helix Processing helix chain 'R' and resid 107 through 139 removed outlier: 4.133A pdb=" N SER R 123 " --> pdb=" O SER R 119 " (cutoff:3.500A) removed outlier: 4.608A pdb=" N VAL R 124 " --> pdb=" O LEU R 120 " (cutoff:3.500A) removed outlier: 4.182A pdb=" N LEU R 125 " --> pdb=" O TYR R 121 " (cutoff:3.500A) removed outlier: 3.915A pdb=" N ILE R 126 " --> pdb=" O SER R 122 " (cutoff:3.500A) Processing helix chain 'R' and resid 145 through 159 removed outlier: 3.880A pdb=" N VAL R 156 " --> pdb=" O ALA R 152 " (cutoff:3.500A) removed outlier: 4.617A pdb=" N TYR R 157 " --> pdb=" O GLU R 153 " (cutoff:3.500A) removed outlier: 3.965A pdb=" N VAL R 158 " --> pdb=" O LYS R 154 " (cutoff:3.500A) Processing helix chain 'R' and resid 161 through 174 removed outlier: 3.865A pdb=" N THR R 168 " --> pdb=" O ALA R 164 " (cutoff:3.500A) removed outlier: 4.489A pdb=" N ILE R 169 " --> pdb=" O LEU R 165 " (cutoff:3.500A) Proline residue: R 170 - end of helix Processing helix chain 'R' and resid 193 through 208 removed outlier: 4.061A pdb=" N VAL R 197 " --> pdb=" O ASP R 193 " (cutoff:3.500A) Processing helix chain 'R' and resid 210 through 225 Processing helix chain 'R' and resid 235 through 251 removed outlier: 3.819A pdb=" N THR R 240 " --> pdb=" O LYS R 236 " (cutoff:3.500A) Processing helix chain 'R' and resid 253 through 266 Processing helix chain 'R' and resid 275 through 295 removed outlier: 4.333A pdb=" N PHE R 293 " --> pdb=" O ALA R 289 " (cutoff:3.500A) removed outlier: 5.046A pdb=" N HIS R 294 " --> pdb=" O LEU R 290 " (cutoff:3.500A) removed outlier: 5.518A pdb=" N CYS R 295 " --> pdb=" O ALA R 291 " (cutoff:3.500A) Processing helix chain 'R' and resid 297 through 303 Processing sheet with id= A, first strand: chain 'L' and resid 4 through 7 Processing sheet with id= B, first strand: chain 'L' and resid 10 through 12 removed outlier: 6.700A pdb=" N LYS L 108 " --> pdb=" O LEU L 11 " (cutoff:3.500A) removed outlier: 3.624A pdb=" N GLY L 89 " --> pdb=" O VAL L 109 " (cutoff:3.500A) removed outlier: 3.761A pdb=" N ILE L 53 " --> pdb=" O TRP L 40 " (cutoff:3.500A) removed outlier: 5.416A pdb=" N GLN L 42 " --> pdb=" O ARG L 51 " (cutoff:3.500A) removed outlier: 5.113A pdb=" N ARG L 51 " --> pdb=" O GLN L 42 " (cutoff:3.500A) Processing sheet with id= C, first strand: chain 'H' and resid 3 through 6 removed outlier: 3.506A pdb=" N THR H 71 " --> pdb=" O TYR H 80 " (cutoff:3.500A) Processing sheet with id= D, first strand: chain 'H' and resid 10 through 12 removed outlier: 5.997A pdb=" N ARG H 38 " --> pdb=" O TRP H 47 " (cutoff:3.500A) removed outlier: 5.183A pdb=" N TRP H 47 " --> pdb=" O ARG H 38 " (cutoff:3.500A) Processing sheet with id= E, first strand: chain 'R' and resid 175 through 180 206 hydrogen bonds defined for protein. 588 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: 1.03 Time building geometry restraints manager: 1.86 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.34: 1254 1.34 - 1.46: 1017 1.46 - 1.58: 1910 1.58 - 1.70: 1 1.70 - 1.82: 33 Bond restraints: 4215 Sorted by residual: bond pdb=" C6 CLR R 701 " pdb=" C7 CLR R 701 " ideal model delta sigma weight residual 1.492 1.542 -0.050 2.00e-02 2.50e+03 6.17e+00 bond pdb=" C14 CLR R 701 " pdb=" C15 CLR R 701 " ideal model delta sigma weight residual 1.523 1.565 -0.042 2.00e-02 2.50e+03 4.40e+00 bond pdb=" C13 CLR R 701 " pdb=" C14 CLR R 701 " ideal model delta sigma weight residual 1.537 1.578 -0.041 2.00e-02 2.50e+03 4.23e+00 bond pdb=" C20 CLR R 701 " pdb=" C22 CLR R 701 " ideal model delta sigma weight residual 1.535 1.575 -0.040 2.00e-02 2.50e+03 4.01e+00 bond pdb=" C15 CLR R 701 " pdb=" C16 CLR R 701 " ideal model delta sigma weight residual 1.541 1.581 -0.040 2.00e-02 2.50e+03 3.97e+00 ... (remaining 4210 not shown) Histogram of bond angle deviations from ideal: 99.15 - 106.16: 92 106.16 - 113.17: 2269 113.17 - 120.18: 1536 120.18 - 127.18: 1788 127.18 - 134.19: 52 Bond angle restraints: 5737 Sorted by residual: angle pdb=" CA TRP H 116 " pdb=" CB TRP H 116 " pdb=" CG TRP H 116 " ideal model delta sigma weight residual 113.60 119.67 -6.07 1.90e+00 2.77e-01 1.02e+01 angle pdb=" C TYR H 111 " pdb=" N GLY H 112 " pdb=" CA GLY H 112 " ideal model delta sigma weight residual 123.46 116.30 7.16 2.42e+00 1.71e-01 8.74e+00 angle pdb=" C ASP R 181 " pdb=" CA ASP R 181 " pdb=" CB ASP R 181 " ideal model delta sigma weight residual 110.42 116.07 -5.65 1.99e+00 2.53e-01 8.07e+00 angle pdb=" N LEU L 52 " pdb=" CA LEU L 52 " pdb=" C LEU L 52 " ideal model delta sigma weight residual 112.12 108.44 3.68 1.34e+00 5.57e-01 7.53e+00 angle pdb=" C GLY H 112 " pdb=" N MET H 113 " pdb=" CA MET H 113 " ideal model delta sigma weight residual 123.93 128.07 -4.14 1.54e+00 4.22e-01 7.24e+00 ... (remaining 5732 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.92: 2254 17.92 - 35.83: 205 35.83 - 53.75: 35 53.75 - 71.67: 4 71.67 - 89.58: 5 Dihedral angle restraints: 2503 sinusoidal: 1011 harmonic: 1492 Sorted by residual: dihedral pdb=" CB CYS R 28 " pdb=" SG CYS R 28 " pdb=" SG CYS R 274 " pdb=" CB CYS R 274 " ideal model delta sinusoidal sigma weight residual -86.00 -161.73 75.73 1 1.00e+01 1.00e-02 7.23e+01 dihedral pdb=" CB CYS L 23 " pdb=" SG CYS L 23 " pdb=" SG CYS L 93 " pdb=" CB CYS L 93 " ideal model delta sinusoidal sigma weight residual 93.00 135.21 -42.21 1 1.00e+01 1.00e-02 2.49e+01 dihedral pdb=" CB GLU L 86 " pdb=" CG GLU L 86 " pdb=" CD GLU L 86 " pdb=" OE1 GLU L 86 " ideal model delta sinusoidal sigma weight residual 0.00 89.58 -89.58 1 3.00e+01 1.11e-03 1.06e+01 ... (remaining 2500 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.074: 588 0.074 - 0.148: 55 0.148 - 0.222: 3 0.222 - 0.296: 2 0.296 - 0.371: 2 Chirality restraints: 650 Sorted by residual: chirality pdb=" C17 CLR R 701 " pdb=" C13 CLR R 701 " pdb=" C16 CLR R 701 " pdb=" C20 CLR R 701 " both_signs ideal model delta sigma weight residual False 2.55 2.92 -0.37 2.00e-01 2.50e+01 3.43e+00 chirality pdb=" C20 CLR R 701 " pdb=" C17 CLR R 701 " pdb=" C21 CLR R 701 " pdb=" C22 CLR R 701 " both_signs ideal model delta sigma weight residual False 2.59 2.90 -0.31 2.00e-01 2.50e+01 2.42e+00 chirality pdb=" C13 CLR R 701 " pdb=" C12 CLR R 701 " pdb=" C14 CLR R 701 " pdb=" C17 CLR R 701 " both_signs ideal model delta sigma weight residual False -2.93 -3.22 0.29 2.00e-01 2.50e+01 2.11e+00 ... (remaining 647 not shown) Planarity restraints: 704 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA LEU R 216 " -0.014 2.00e-02 2.50e+03 2.79e-02 7.81e+00 pdb=" C LEU R 216 " 0.048 2.00e-02 2.50e+03 pdb=" O LEU R 216 " -0.018 2.00e-02 2.50e+03 pdb=" N SER R 217 " -0.016 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C TYR R 190 " -0.042 5.00e-02 4.00e+02 6.31e-02 6.38e+00 pdb=" N PRO R 191 " 0.109 5.00e-02 4.00e+02 pdb=" CA PRO R 191 " -0.032 5.00e-02 4.00e+02 pdb=" CD PRO R 191 " -0.035 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C LEU R 91 " 0.041 5.00e-02 4.00e+02 6.15e-02 6.05e+00 pdb=" N PRO R 92 " -0.106 5.00e-02 4.00e+02 pdb=" CA PRO R 92 " 0.031 5.00e-02 4.00e+02 pdb=" CD PRO R 92 " 0.034 5.00e-02 4.00e+02 ... (remaining 701 not shown) Histogram of nonbonded interaction distances: 2.17 - 2.71: 108 2.71 - 3.26: 3875 3.26 - 3.81: 6346 3.81 - 4.35: 7822 4.35 - 4.90: 14119 Nonbonded interactions: 32270 Sorted by model distance: nonbonded pdb=" OG1 THR R 117 " pdb=" OD2 ASP R 171 " model vdw 2.165 2.440 nonbonded pdb=" NH1 ARG L 51 " pdb=" OD1 ASP H 114 " model vdw 2.221 2.520 nonbonded pdb=" OE1 GLN R 200 " pdb=" OG SER R 263 " model vdw 2.313 2.440 nonbonded pdb=" O ALA R 83 " pdb=" OG SER R 119 " model vdw 2.334 2.440 nonbonded pdb=" NH2 ARG H 105 " pdb=" OE1 GLU R 288 " model vdw 2.404 2.520 ... (remaining 32265 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 1.910 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.030 Construct map_model_manager: 0.010 Extract box with map and model: 5.810 Check model and map are aligned: 0.050 Set scattering table: 0.040 Process input model: 15.960 Find NCS groups from input model: 0.100 Set up NCS constraints: 0.020 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:2.760 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 26.690 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8155 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.050 4215 Z= 0.197 Angle : 0.672 7.746 5737 Z= 0.345 Chirality : 0.050 0.371 650 Planarity : 0.006 0.063 704 Dihedral : 14.014 89.582 1537 Min Nonbonded Distance : 2.165 Molprobity Statistics. All-atom Clashscore : 4.99 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.58 % Favored : 96.42 % Rotamer: Outliers : 0.00 % Allowed : 0.22 % Favored : 99.78 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.38 (0.40), residues: 503 helix: 1.07 (0.38), residues: 218 sheet: 0.04 (0.46), residues: 133 loop : -0.55 (0.55), residues: 152 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.041 0.002 TRP H 116 HIS 0.004 0.001 HIS R 294 PHE 0.012 0.001 PHE R 293 TYR 0.010 0.001 TYR H 107 ARG 0.004 0.001 ARG L 51 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 67 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 67 time to evaluate : 0.425 Fit side-chains REVERT: R 24 MET cc_start: 0.6356 (tmm) cc_final: 0.6053 (tmm) REVERT: R 301 LEU cc_start: 0.8430 (mp) cc_final: 0.8173 (mp) outliers start: 0 outliers final: 0 residues processed: 67 average time/residue: 0.1995 time to fit residues: 16.2193 Evaluate side-chains 59 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 59 time to evaluate : 0.458 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 42 optimal weight: 0.9980 chunk 37 optimal weight: 1.9990 chunk 21 optimal weight: 0.4980 chunk 12 optimal weight: 0.7980 chunk 25 optimal weight: 0.7980 chunk 20 optimal weight: 2.9990 chunk 39 optimal weight: 0.9990 chunk 15 optimal weight: 0.3980 chunk 23 optimal weight: 0.8980 chunk 29 optimal weight: 0.6980 chunk 45 optimal weight: 0.7980 overall best weight: 0.6380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: R 203 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8222 moved from start: 0.0865 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.035 4215 Z= 0.205 Angle : 0.550 5.187 5737 Z= 0.285 Chirality : 0.042 0.159 650 Planarity : 0.005 0.054 704 Dihedral : 4.381 21.893 627 Min Nonbonded Distance : 2.209 Molprobity Statistics. All-atom Clashscore : 4.26 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.78 % Favored : 96.22 % Rotamer: Outliers : 0.90 % Allowed : 4.93 % Favored : 94.17 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.59 (0.39), residues: 503 helix: 1.43 (0.37), residues: 212 sheet: 0.18 (0.47), residues: 133 loop : -0.74 (0.52), residues: 158 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP H 47 HIS 0.003 0.001 HIS R 294 PHE 0.008 0.001 PHE L 41 TYR 0.009 0.001 TYR H 54 ARG 0.003 0.000 ARG R 235 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 66 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 62 time to evaluate : 0.482 Fit side-chains revert: symmetry clash REVERT: H 18 VAL cc_start: 0.8493 (t) cc_final: 0.8258 (t) REVERT: H 113 MET cc_start: 0.8532 (mmm) cc_final: 0.7986 (mmt) REVERT: R 301 LEU cc_start: 0.8413 (mp) cc_final: 0.8167 (mp) outliers start: 4 outliers final: 3 residues processed: 63 average time/residue: 0.2209 time to fit residues: 16.8594 Evaluate side-chains 63 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 60 time to evaluate : 0.399 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain L residue 39 ASN Chi-restraints excluded: chain L residue 75 ASP Chi-restraints excluded: chain R residue 33 ASN Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 25 optimal weight: 0.8980 chunk 14 optimal weight: 2.9990 chunk 37 optimal weight: 0.9990 chunk 30 optimal weight: 0.7980 chunk 12 optimal weight: 0.9980 chunk 45 optimal weight: 0.9990 chunk 49 optimal weight: 0.9990 chunk 40 optimal weight: 1.9990 chunk 15 optimal weight: 0.8980 chunk 36 optimal weight: 0.5980 chunk 44 optimal weight: 0.8980 overall best weight: 0.8180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: R 203 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8264 moved from start: 0.1269 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.039 4215 Z= 0.239 Angle : 0.551 5.232 5737 Z= 0.286 Chirality : 0.043 0.155 650 Planarity : 0.005 0.055 704 Dihedral : 4.422 20.749 627 Min Nonbonded Distance : 2.238 Molprobity Statistics. All-atom Clashscore : 4.38 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.98 % Favored : 97.02 % Rotamer: Outliers : 0.90 % Allowed : 8.07 % Favored : 91.03 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.70 (0.40), residues: 503 helix: 1.55 (0.37), residues: 212 sheet: 0.22 (0.46), residues: 133 loop : -0.71 (0.53), residues: 158 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP H 47 HIS 0.003 0.001 HIS R 203 PHE 0.005 0.001 PHE L 41 TYR 0.010 0.001 TYR H 54 ARG 0.002 0.000 ARG L 51 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 65 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 61 time to evaluate : 0.493 Fit side-chains revert: symmetry clash REVERT: H 18 VAL cc_start: 0.8536 (t) cc_final: 0.8299 (t) REVERT: R 154 LYS cc_start: 0.7700 (mptt) cc_final: 0.6703 (tttm) outliers start: 4 outliers final: 3 residues processed: 62 average time/residue: 0.2210 time to fit residues: 16.6886 Evaluate side-chains 63 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 60 time to evaluate : 0.558 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain L residue 39 ASN Chi-restraints excluded: chain L residue 75 ASP Chi-restraints excluded: chain R residue 33 ASN Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 34 optimal weight: 7.9990 chunk 23 optimal weight: 3.9990 chunk 5 optimal weight: 0.8980 chunk 21 optimal weight: 2.9990 chunk 30 optimal weight: 1.9990 chunk 45 optimal weight: 0.3980 chunk 48 optimal weight: 2.9990 chunk 43 optimal weight: 0.5980 chunk 13 optimal weight: 0.6980 chunk 40 optimal weight: 0.7980 chunk 27 optimal weight: 5.9990 overall best weight: 0.6780 ------------------------------------------------------------------------------- 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.8255 moved from start: 0.1406 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.034 4215 Z= 0.206 Angle : 0.531 5.533 5737 Z= 0.272 Chirality : 0.042 0.153 650 Planarity : 0.004 0.054 704 Dihedral : 4.326 19.977 627 Min Nonbonded Distance : 2.237 Molprobity Statistics. All-atom Clashscore : 4.26 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.18 % Favored : 96.82 % Rotamer: Outliers : 1.12 % Allowed : 10.76 % Favored : 88.12 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.81 (0.40), residues: 503 helix: 1.70 (0.37), residues: 212 sheet: 0.19 (0.46), residues: 133 loop : -0.67 (0.53), residues: 158 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP H 47 HIS 0.002 0.001 HIS R 79 PHE 0.005 0.001 PHE L 41 TYR 0.009 0.001 TYR H 54 ARG 0.002 0.000 ARG L 51 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 62 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 57 time to evaluate : 0.590 Fit side-chains revert: symmetry clash REVERT: R 154 LYS cc_start: 0.7890 (mptt) cc_final: 0.6855 (tttm) outliers start: 5 outliers final: 4 residues processed: 58 average time/residue: 0.2373 time to fit residues: 16.7184 Evaluate side-chains 61 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 57 time to evaluate : 0.472 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain L residue 39 ASN Chi-restraints excluded: chain L residue 75 ASP Chi-restraints excluded: chain R residue 33 ASN Chi-restraints excluded: chain R residue 304 PHE Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 0 optimal weight: 2.9990 chunk 35 optimal weight: 1.9990 chunk 19 optimal weight: 3.9990 chunk 41 optimal weight: 1.9990 chunk 33 optimal weight: 0.9980 chunk 24 optimal weight: 0.2980 chunk 43 optimal weight: 0.0370 chunk 12 optimal weight: 0.7980 chunk 16 optimal weight: 2.9990 chunk 9 optimal weight: 1.9990 chunk 28 optimal weight: 0.6980 overall best weight: 0.5658 ------------------------------------------------------------------------------- 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.8247 moved from start: 0.1500 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 4215 Z= 0.184 Angle : 0.509 5.325 5737 Z= 0.261 Chirality : 0.041 0.147 650 Planarity : 0.004 0.054 704 Dihedral : 4.214 19.254 627 Min Nonbonded Distance : 2.233 Molprobity Statistics. All-atom Clashscore : 4.26 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.98 % Favored : 97.02 % Rotamer: Outliers : 1.79 % Allowed : 11.66 % Favored : 86.55 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.93 (0.40), residues: 503 helix: 1.89 (0.37), residues: 210 sheet: 0.18 (0.46), residues: 133 loop : -0.63 (0.53), residues: 160 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP H 47 HIS 0.002 0.001 HIS R 294 PHE 0.005 0.001 PHE L 41 TYR 0.009 0.001 TYR H 54 ARG 0.001 0.000 ARG L 51 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 66 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 58 time to evaluate : 0.453 Fit side-chains REVERT: R 28 CYS cc_start: 0.6129 (m) cc_final: 0.5776 (m) REVERT: R 61 LEU cc_start: 0.7828 (tt) cc_final: 0.7600 (tp) REVERT: R 154 LYS cc_start: 0.7927 (mptt) cc_final: 0.6826 (tttm) outliers start: 8 outliers final: 5 residues processed: 60 average time/residue: 0.2104 time to fit residues: 15.2644 Evaluate side-chains 63 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 58 time to evaluate : 0.465 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain L residue 39 ASN Chi-restraints excluded: chain L residue 72 SER Chi-restraints excluded: chain L residue 75 ASP Chi-restraints excluded: chain R residue 33 ASN Chi-restraints excluded: chain R residue 304 PHE Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 11 optimal weight: 0.9990 chunk 48 optimal weight: 2.9990 chunk 40 optimal weight: 0.7980 chunk 22 optimal weight: 3.9990 chunk 4 optimal weight: 0.6980 chunk 16 optimal weight: 2.9990 chunk 25 optimal weight: 0.9990 chunk 46 optimal weight: 0.8980 chunk 5 optimal weight: 0.6980 chunk 27 optimal weight: 4.9990 chunk 35 optimal weight: 0.6980 overall best weight: 0.7580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: R 203 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8272 moved from start: 0.1655 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.040 4215 Z= 0.221 Angle : 0.534 5.387 5737 Z= 0.274 Chirality : 0.042 0.145 650 Planarity : 0.004 0.054 704 Dihedral : 4.281 18.773 627 Min Nonbonded Distance : 2.249 Molprobity Statistics. All-atom Clashscore : 4.50 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.18 % Favored : 96.82 % Rotamer: Outliers : 1.57 % Allowed : 13.23 % Favored : 85.20 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.93 (0.40), residues: 503 helix: 1.89 (0.37), residues: 209 sheet: 0.19 (0.46), residues: 133 loop : -0.65 (0.52), residues: 161 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.001 TRP H 47 HIS 0.002 0.001 HIS R 203 PHE 0.005 0.001 PHE L 41 TYR 0.010 0.001 TYR H 54 ARG 0.002 0.000 ARG L 51 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 67 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 60 time to evaluate : 0.415 Fit side-chains REVERT: R 28 CYS cc_start: 0.6264 (m) cc_final: 0.5914 (m) REVERT: R 61 LEU cc_start: 0.7809 (tt) cc_final: 0.7583 (tp) REVERT: R 154 LYS cc_start: 0.8014 (mptt) cc_final: 0.6881 (tttm) outliers start: 7 outliers final: 7 residues processed: 62 average time/residue: 0.2144 time to fit residues: 16.1883 Evaluate side-chains 65 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 58 time to evaluate : 0.450 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain L residue 39 ASN Chi-restraints excluded: chain L residue 72 SER Chi-restraints excluded: chain L residue 75 ASP Chi-restraints excluded: chain R residue 33 ASN Chi-restraints excluded: chain R residue 122 SER Chi-restraints excluded: chain R residue 185 ILE Chi-restraints excluded: chain R residue 304 PHE Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 27 optimal weight: 4.9990 chunk 40 optimal weight: 0.9990 chunk 48 optimal weight: 2.9990 chunk 30 optimal weight: 1.9990 chunk 29 optimal weight: 0.5980 chunk 22 optimal weight: 1.9990 chunk 19 optimal weight: 3.9990 chunk 28 optimal weight: 0.5980 chunk 14 optimal weight: 1.9990 chunk 9 optimal weight: 1.9990 chunk 32 optimal weight: 0.9990 overall best weight: 1.0386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: R 202 GLN R 203 HIS Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8304 moved from start: 0.1856 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.049 4215 Z= 0.277 Angle : 0.569 5.409 5737 Z= 0.292 Chirality : 0.043 0.148 650 Planarity : 0.004 0.055 704 Dihedral : 4.478 18.240 627 Min Nonbonded Distance : 2.268 Molprobity Statistics. All-atom Clashscore : 4.99 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.78 % Favored : 97.22 % Rotamer: Outliers : 2.24 % Allowed : 14.57 % Favored : 83.18 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.83 (0.40), residues: 503 helix: 1.80 (0.38), residues: 209 sheet: 0.14 (0.46), residues: 133 loop : -0.67 (0.53), residues: 161 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.015 0.001 TRP H 47 HIS 0.003 0.001 HIS H 110 PHE 0.007 0.001 PHE R 199 TYR 0.013 0.001 TYR H 54 ARG 0.003 0.000 ARG L 51 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 70 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 60 time to evaluate : 0.598 Fit side-chains REVERT: R 28 CYS cc_start: 0.6293 (m) cc_final: 0.5940 (m) REVERT: R 61 LEU cc_start: 0.7844 (tt) cc_final: 0.7610 (tp) REVERT: R 154 LYS cc_start: 0.8030 (mptt) cc_final: 0.6895 (tttm) outliers start: 10 outliers final: 7 residues processed: 64 average time/residue: 0.2374 time to fit residues: 18.4026 Evaluate side-chains 67 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 60 time to evaluate : 0.444 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain L residue 39 ASN Chi-restraints excluded: chain L residue 72 SER Chi-restraints excluded: chain L residue 75 ASP Chi-restraints excluded: chain R residue 33 ASN Chi-restraints excluded: chain R residue 122 SER Chi-restraints excluded: chain R residue 132 LEU Chi-restraints excluded: chain R residue 150 LEU Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 23 optimal weight: 3.9990 chunk 4 optimal weight: 0.9980 chunk 37 optimal weight: 0.0020 chunk 43 optimal weight: 7.9990 chunk 46 optimal weight: 1.9990 chunk 42 optimal weight: 0.7980 chunk 44 optimal weight: 1.9990 chunk 27 optimal weight: 4.9990 chunk 19 optimal weight: 2.9990 chunk 35 optimal weight: 0.9990 chunk 13 optimal weight: 0.9980 overall best weight: 0.7590 ------------------------------------------------------------------------------- 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.8285 moved from start: 0.1853 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.041 4215 Z= 0.221 Angle : 0.535 5.370 5737 Z= 0.275 Chirality : 0.042 0.145 650 Planarity : 0.004 0.056 704 Dihedral : 4.365 17.515 627 Min Nonbonded Distance : 2.259 Molprobity Statistics. All-atom Clashscore : 5.35 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.38 % Favored : 96.62 % Rotamer: Outliers : 2.24 % Allowed : 14.80 % Favored : 82.96 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.86 (0.40), residues: 503 helix: 1.85 (0.37), residues: 209 sheet: 0.14 (0.46), residues: 134 loop : -0.70 (0.53), residues: 160 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.001 TRP H 47 HIS 0.002 0.001 HIS H 110 PHE 0.005 0.001 PHE L 41 TYR 0.011 0.001 TYR H 54 ARG 0.002 0.000 ARG L 51 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 68 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 58 time to evaluate : 0.463 Fit side-chains REVERT: R 28 CYS cc_start: 0.6246 (m) cc_final: 0.5890 (m) REVERT: R 61 LEU cc_start: 0.7816 (tt) cc_final: 0.7589 (tp) REVERT: R 154 LYS cc_start: 0.8052 (mptt) cc_final: 0.6900 (tttm) REVERT: R 293 PHE cc_start: 0.8068 (OUTLIER) cc_final: 0.7508 (t80) outliers start: 10 outliers final: 8 residues processed: 61 average time/residue: 0.2217 time to fit residues: 16.4143 Evaluate side-chains 67 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 58 time to evaluate : 0.457 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain L residue 39 ASN Chi-restraints excluded: chain L residue 72 SER Chi-restraints excluded: chain L residue 75 ASP Chi-restraints excluded: chain R residue 33 ASN Chi-restraints excluded: chain R residue 132 LEU Chi-restraints excluded: chain R residue 150 LEU Chi-restraints excluded: chain R residue 185 ILE Chi-restraints excluded: chain R residue 293 PHE Chi-restraints excluded: chain R residue 304 PHE Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 40 optimal weight: 0.9990 chunk 42 optimal weight: 0.6980 chunk 44 optimal weight: 2.9990 chunk 29 optimal weight: 0.9990 chunk 47 optimal weight: 0.9980 chunk 28 optimal weight: 0.9990 chunk 22 optimal weight: 1.9990 chunk 33 optimal weight: 0.6980 chunk 49 optimal weight: 0.7980 chunk 45 optimal weight: 0.5980 chunk 39 optimal weight: 0.5980 overall best weight: 0.6780 ------------------------------------------------------------------------------- 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.8276 moved from start: 0.1861 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.037 4215 Z= 0.203 Angle : 0.519 5.454 5737 Z= 0.267 Chirality : 0.041 0.143 650 Planarity : 0.004 0.056 704 Dihedral : 4.244 16.591 627 Min Nonbonded Distance : 2.253 Molprobity Statistics. All-atom Clashscore : 4.99 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.58 % Favored : 96.42 % Rotamer: Outliers : 2.47 % Allowed : 14.80 % Favored : 82.74 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.90 (0.40), residues: 503 helix: 1.89 (0.38), residues: 209 sheet: 0.14 (0.46), residues: 134 loop : -0.66 (0.53), residues: 160 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP H 47 HIS 0.002 0.001 HIS H 110 PHE 0.005 0.001 PHE L 41 TYR 0.010 0.001 TYR H 54 ARG 0.001 0.000 ARG L 51 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 70 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 59 time to evaluate : 0.420 Fit side-chains REVERT: R 28 CYS cc_start: 0.6233 (m) cc_final: 0.5890 (m) REVERT: R 61 LEU cc_start: 0.7816 (tt) cc_final: 0.7596 (tp) REVERT: R 154 LYS cc_start: 0.8076 (mptt) cc_final: 0.6898 (tttm) outliers start: 11 outliers final: 8 residues processed: 64 average time/residue: 0.1988 time to fit residues: 15.4812 Evaluate side-chains 68 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 60 time to evaluate : 0.449 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain L residue 39 ASN Chi-restraints excluded: chain L residue 72 SER Chi-restraints excluded: chain L residue 75 ASP Chi-restraints excluded: chain R residue 33 ASN Chi-restraints excluded: chain R residue 132 LEU Chi-restraints excluded: chain R residue 150 LEU Chi-restraints excluded: chain R residue 185 ILE Chi-restraints excluded: chain R residue 304 PHE Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 4 optimal weight: 0.0980 chunk 30 optimal weight: 1.9990 chunk 24 optimal weight: 0.9990 chunk 31 optimal weight: 1.9990 chunk 42 optimal weight: 0.5980 chunk 12 optimal weight: 0.6980 chunk 36 optimal weight: 3.9990 chunk 5 optimal weight: 0.3980 chunk 11 optimal weight: 0.9980 chunk 39 optimal weight: 0.7980 chunk 16 optimal weight: 2.9990 overall best weight: 0.5180 ------------------------------------------------------------------------------- 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.8259 moved from start: 0.1849 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.034 4215 Z= 0.176 Angle : 0.511 5.544 5737 Z= 0.261 Chirality : 0.041 0.154 650 Planarity : 0.004 0.055 704 Dihedral : 4.148 16.215 627 Min Nonbonded Distance : 2.241 Molprobity Statistics. All-atom Clashscore : 5.11 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.58 % Favored : 96.42 % Rotamer: Outliers : 1.79 % Allowed : 15.25 % Favored : 82.96 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.00 (0.40), residues: 503 helix: 1.97 (0.37), residues: 209 sheet: 0.19 (0.46), residues: 134 loop : -0.61 (0.53), residues: 160 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP H 47 HIS 0.003 0.001 HIS R 203 PHE 0.005 0.001 PHE L 41 TYR 0.008 0.001 TYR H 54 ARG 0.003 0.000 ARG H 38 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1006 Ramachandran restraints generated. 503 Oldfield, 0 Emsley, 503 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 65 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 57 time to evaluate : 0.524 Fit side-chains REVERT: R 28 CYS cc_start: 0.6217 (m) cc_final: 0.5878 (m) REVERT: R 61 LEU cc_start: 0.7780 (tt) cc_final: 0.7564 (tp) REVERT: R 154 LYS cc_start: 0.8086 (mptt) cc_final: 0.6902 (tttm) outliers start: 8 outliers final: 8 residues processed: 60 average time/residue: 0.2133 time to fit residues: 15.5479 Evaluate side-chains 64 residues out of total 446 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 56 time to evaluate : 0.523 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain L residue 39 ASN Chi-restraints excluded: chain L residue 72 SER Chi-restraints excluded: chain L residue 75 ASP Chi-restraints excluded: chain R residue 33 ASN Chi-restraints excluded: chain R residue 132 LEU Chi-restraints excluded: chain R residue 150 LEU Chi-restraints excluded: chain R residue 185 ILE Chi-restraints excluded: chain R residue 304 PHE Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 40 optimal weight: 1.9990 chunk 5 optimal weight: 0.9980 chunk 7 optimal weight: 2.9990 chunk 34 optimal weight: 3.9990 chunk 2 optimal weight: 0.9990 chunk 28 optimal weight: 0.6980 chunk 45 optimal weight: 0.2980 chunk 26 optimal weight: 0.6980 chunk 1 optimal weight: 0.0970 chunk 31 optimal weight: 2.9990 chunk 30 optimal weight: 0.8980 overall best weight: 0.5378 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3659 r_free = 0.3659 target = 0.139954 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 40)----------------| | r_work = 0.3412 r_free = 0.3412 target = 0.120563 restraints weight = 4672.304| |-----------------------------------------------------------------------------| r_work (start): 0.3418 rms_B_bonded: 1.65 r_work: 0.3298 rms_B_bonded: 2.25 restraints_weight: 0.5000 r_work: 0.3156 rms_B_bonded: 3.91 restraints_weight: 0.2500 r_work (final): 0.3156 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8481 moved from start: 0.1864 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.034 4215 Z= 0.180 Angle : 0.509 5.468 5737 Z= 0.260 Chirality : 0.041 0.148 650 Planarity : 0.004 0.055 704 Dihedral : 4.078 15.903 627 Min Nonbonded Distance : 2.242 Molprobity Statistics. All-atom Clashscore : 4.75 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.38 % Favored : 96.62 % Rotamer: Outliers : 1.79 % Allowed : 15.70 % Favored : 82.51 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 10.53 % 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.05 (0.40), residues: 503 helix: 2.00 (0.37), residues: 210 sheet: 0.19 (0.46), residues: 134 loop : -0.55 (0.54), residues: 159 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP H 47 HIS 0.002 0.001 HIS R 203 PHE 0.006 0.001 PHE L 41 TYR 0.008 0.001 TYR H 54 ARG 0.001 0.000 ARG H 38 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 1306.00 seconds wall clock time: 24 minutes 25.86 seconds (1465.86 seconds total)