Starting phenix.real_space_refine (version: dev) on Sun Dec 11 14:08:37 2022 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8ctj_26992/12_2022/8ctj_26992_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8ctj_26992/12_2022/8ctj_26992.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=4.74 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8ctj_26992/12_2022/8ctj_26992.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8ctj_26992/12_2022/8ctj_26992.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8ctj_26992/12_2022/8ctj_26992_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8ctj_26992/12_2022/8ctj_26992_updated.pdb" } resolution = 4.74 write_initial_geo_file = False refinement { macro_cycles = 10 } 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 ------------------------------------------------------------------------------- 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.005 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 Convert atoms to be neutral Process input model Symmetric amino acids flipped Residue "A GLU 72": "OE1" <-> "OE2" Residue "A TYR 83": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 96": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 102": "OE1" <-> "OE2" Residue "A TYR 104": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 125": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 137": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ASP 140": "OD1" <-> "OD2" Residue "A PHE 141": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 206": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 216": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 217": "OE1" <-> "OE2" Residue "A TYR 218": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ASP 254": "OD1" <-> "OD2" Residue "A PHE 267": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 277": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 319": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 404": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ASP 434": "OD1" <-> "OD2" Residue "A GLU 437": "OE1" <-> "OE2" Residue "A PHE 454": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/chem_data/mon_lib" Total number of atoms: 3312 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 3312 Number of conformers: 1 Conformer: "" Number of residues, atoms: 402, 3312 Unusual residues: {'PEE': 1} Classifications: {'peptide': 401, 'undetermined': 1} Link IDs: {'PTRANS': 10, 'TRANS': 390, None: 1} Not linked: pdbres="LEU A 474 " pdbres="PEE A 601 " Chain breaks: 2 Unresolved non-hydrogen bonds: 25 Unresolved non-hydrogen angles: 25 Unresolved non-hydrogen dihedrals: 23 Time building chain proxies: 2.55, per 1000 atoms: 0.77 Number of scatterers: 3312 At special positions: 0 Unit cell: (55.75, 57.98, 105.925, 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 18 16.00 P 1 15.00 O 562 8.00 N 529 7.00 C 2202 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=2, symmetry=0 Simple disulfide: pdb=" SG CYS A 74 " - pdb=" SG CYS A 128 " distance=2.04 Simple disulfide: pdb=" SG CYS A 89 " - pdb=" SG CYS A 431 " 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.39 Conformation dependent library (CDL) restraints added in 481.7 milliseconds 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 762 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 12 helices and 2 sheets defined 56.6% alpha, 10.5% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.35 Creating SS restraints... Processing helix chain 'A' and resid 91 through 94 No H-bonds generated for 'chain 'A' and resid 91 through 94' Processing helix chain 'A' and resid 98 through 110 removed outlier: 3.651A pdb=" N VAL A 101 " --> pdb=" O ALA A 98 " (cutoff:3.500A) removed outlier: 4.780A pdb=" N GLU A 102 " --> pdb=" O GLU A 99 " (cutoff:3.500A) removed outlier: 4.047A pdb=" N LYS A 109 " --> pdb=" O GLU A 106 " (cutoff:3.500A) removed outlier: 3.789A pdb=" N GLU A 110 " --> pdb=" O LYS A 107 " (cutoff:3.500A) Processing helix chain 'A' and resid 128 through 132 removed outlier: 3.729A pdb=" N PHE A 132 " --> pdb=" O CYS A 128 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 128 through 132' Processing helix chain 'A' and resid 221 through 250 Proline residue: A 225 - end of helix removed outlier: 3.618A pdb=" N ILE A 228 " --> pdb=" O TYR A 224 " (cutoff:3.500A) removed outlier: 4.164A pdb=" N CYS A 250 " --> pdb=" O ALA A 246 " (cutoff:3.500A) Processing helix chain 'A' and resid 257 through 284 Processing helix chain 'A' and resid 291 through 317 Processing helix chain 'A' and resid 328 through 352 removed outlier: 4.554A pdb=" N MET A 346 " --> pdb=" O LEU A 342 " (cutoff:3.500A) Processing helix chain 'A' and resid 357 through 392 removed outlier: 4.074A pdb=" N PHE A 364 " --> pdb=" O ALA A 360 " (cutoff:3.500A) removed outlier: 4.687A pdb=" N ILE A 365 " --> pdb=" O SER A 361 " (cutoff:3.500A) Proline residue: A 366 - end of helix Processing helix chain 'A' and resid 395 through 427 removed outlier: 3.747A pdb=" N ARG A 427 " --> pdb=" O THR A 423 " (cutoff:3.500A) Processing helix chain 'A' and resid 435 through 437 No H-bonds generated for 'chain 'A' and resid 435 through 437' Processing helix chain 'A' and resid 444 through 459 Processing helix chain 'A' and resid 464 through 470 removed outlier: 4.323A pdb=" N GLN A 467 " --> pdb=" O ALA A 464 " (cutoff:3.500A) removed outlier: 4.400A pdb=" N ALA A 470 " --> pdb=" O GLN A 467 " (cutoff:3.500A) Processing sheet with id= A, first strand: chain 'A' and resid 42 through 45 Processing sheet with id= B, first strand: chain 'A' and resid 52 through 59 195 hydrogen bonds defined for protein. 549 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 0 stacking parallelities Total time for adding SS restraints: 0.85 Time building geometry restraints manager: 1.51 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: 977 1.34 - 1.47: 929 1.47 - 1.60: 1461 1.60 - 1.73: 1 1.73 - 1.86: 30 Bond restraints: 3398 Sorted by residual: bond pdb=" O4P PEE A 601 " pdb=" P PEE A 601 " ideal model delta sigma weight residual 1.653 1.776 -0.123 2.00e-02 2.50e+03 3.81e+01 bond pdb=" C1 PEE A 601 " pdb=" C2 PEE A 601 " ideal model delta sigma weight residual 1.519 1.555 -0.036 2.00e-02 2.50e+03 3.19e+00 bond pdb=" CA ASP A 75 " pdb=" CB ASP A 75 " ideal model delta sigma weight residual 1.531 1.577 -0.047 3.12e-02 1.03e+03 2.25e+00 bond pdb=" CB CYS A 74 " pdb=" SG CYS A 74 " ideal model delta sigma weight residual 1.808 1.855 -0.047 3.30e-02 9.18e+02 2.05e+00 bond pdb=" C2 PEE A 601 " pdb=" C3 PEE A 601 " ideal model delta sigma weight residual 1.512 1.540 -0.028 2.00e-02 2.50e+03 2.03e+00 ... (remaining 3393 not shown) Histogram of bond angle deviations from ideal: 90.91 - 99.55: 3 99.55 - 108.20: 143 108.20 - 116.84: 2156 116.84 - 125.48: 2226 125.48 - 134.13: 75 Bond angle restraints: 4603 Sorted by residual: angle pdb=" O1P PEE A 601 " pdb=" P PEE A 601 " pdb=" O2P PEE A 601 " ideal model delta sigma weight residual 119.43 133.35 -13.92 3.00e+00 1.11e-01 2.15e+01 angle pdb=" CA LYS A 213 " pdb=" CB LYS A 213 " pdb=" CG LYS A 213 " ideal model delta sigma weight residual 114.10 121.83 -7.73 2.00e+00 2.50e-01 1.49e+01 angle pdb=" N GLN A 135 " pdb=" CA GLN A 135 " pdb=" C GLN A 135 " ideal model delta sigma weight residual 110.44 106.16 4.28 1.20e+00 6.94e-01 1.27e+01 angle pdb=" CA TRP A 252 " pdb=" CB TRP A 252 " pdb=" CG TRP A 252 " ideal model delta sigma weight residual 113.60 120.05 -6.45 1.90e+00 2.77e-01 1.15e+01 angle pdb=" CA MET A 231 " pdb=" CB MET A 231 " pdb=" CG MET A 231 " ideal model delta sigma weight residual 114.10 120.85 -6.75 2.00e+00 2.50e-01 1.14e+01 ... (remaining 4598 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.14: 1800 17.14 - 34.28: 140 34.28 - 51.42: 25 51.42 - 68.56: 7 68.56 - 85.70: 4 Dihedral angle restraints: 1976 sinusoidal: 799 harmonic: 1177 Sorted by residual: dihedral pdb=" CB CYS A 89 " pdb=" SG CYS A 89 " pdb=" SG CYS A 431 " pdb=" CB CYS A 431 " ideal model delta sinusoidal sigma weight residual 93.00 144.73 -51.73 1 1.00e+01 1.00e-02 3.65e+01 dihedral pdb=" CA ASN A 394 " pdb=" C ASN A 394 " pdb=" N ILE A 395 " pdb=" CA ILE A 395 " ideal model delta harmonic sigma weight residual -180.00 -161.55 -18.45 0 5.00e+00 4.00e-02 1.36e+01 dihedral pdb=" CA PHE A 230 " pdb=" C PHE A 230 " pdb=" N MET A 231 " pdb=" CA MET A 231 " ideal model delta harmonic sigma weight residual 180.00 162.98 17.02 0 5.00e+00 4.00e-02 1.16e+01 ... (remaining 1973 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.077: 459 0.077 - 0.155: 54 0.155 - 0.232: 6 0.232 - 0.309: 0 0.309 - 0.387: 1 Chirality restraints: 520 Sorted by residual: chirality pdb=" C2 PEE A 601 " pdb=" C1 PEE A 601 " pdb=" C3 PEE A 601 " pdb=" O2 PEE A 601 " both_signs ideal model delta sigma weight residual False -2.33 -2.71 0.39 2.00e-01 2.50e+01 3.74e+00 chirality pdb=" CB ILE A 44 " pdb=" CA ILE A 44 " pdb=" CG1 ILE A 44 " pdb=" CG2 ILE A 44 " both_signs ideal model delta sigma weight residual False 2.64 2.41 0.23 2.00e-01 2.50e+01 1.32e+00 chirality pdb=" CB THR A 405 " pdb=" CA THR A 405 " pdb=" OG1 THR A 405 " pdb=" CG2 THR A 405 " both_signs ideal model delta sigma weight residual False 2.55 2.33 0.23 2.00e-01 2.50e+01 1.28e+00 ... (remaining 517 not shown) Planarity restraints: 551 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB ASP A 371 " 0.012 2.00e-02 2.50e+03 2.37e-02 5.62e+00 pdb=" CG ASP A 371 " -0.041 2.00e-02 2.50e+03 pdb=" OD1 ASP A 371 " 0.015 2.00e-02 2.50e+03 pdb=" OD2 ASP A 371 " 0.015 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA GLN A 135 " 0.011 2.00e-02 2.50e+03 2.22e-02 4.93e+00 pdb=" C GLN A 135 " -0.038 2.00e-02 2.50e+03 pdb=" O GLN A 135 " 0.015 2.00e-02 2.50e+03 pdb=" N THR A 136 " 0.013 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE A 404 " -0.010 2.00e-02 2.50e+03 1.54e-02 4.17e+00 pdb=" CG PHE A 404 " 0.035 2.00e-02 2.50e+03 pdb=" CD1 PHE A 404 " -0.014 2.00e-02 2.50e+03 pdb=" CD2 PHE A 404 " -0.013 2.00e-02 2.50e+03 pdb=" CE1 PHE A 404 " 0.001 2.00e-02 2.50e+03 pdb=" CE2 PHE A 404 " 0.000 2.00e-02 2.50e+03 pdb=" CZ PHE A 404 " 0.001 2.00e-02 2.50e+03 ... (remaining 548 not shown) Histogram of nonbonded interaction distances: 2.23 - 2.77: 539 2.77 - 3.30: 3286 3.30 - 3.83: 5411 3.83 - 4.37: 5971 4.37 - 4.90: 10264 Nonbonded interactions: 25471 Sorted by model distance: nonbonded pdb=" O TYR A 340 " pdb=" OG SER A 344 " model vdw 2.233 2.440 nonbonded pdb=" OG SER A 77 " pdb=" O SER A 191 " model vdw 2.233 2.440 nonbonded pdb=" O ILE A 419 " pdb=" OG1 THR A 423 " model vdw 2.276 2.440 nonbonded pdb=" O ASN A 200 " pdb=" OG SER A 203 " model vdw 2.312 2.440 nonbonded pdb=" O SER A 399 " pdb=" ND1 HIS A 403 " model vdw 2.381 2.520 ... (remaining 25466 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 Set scattering table Set to: electron Number of scattering types: 5 Type Number sf(0) Gaussians P 1 5.49 5 S 18 5.16 5 C 2202 2.51 5 N 529 2.21 5 O 562 1.98 5 sf(0) = scattering factor at diffraction angle 0. 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.270 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: 1.460 Check model and map are aligned: 0.050 Convert atoms to be neutral: 0.030 Process input model: 14.110 Find NCS groups from input model: 0.080 Set up NCS constraints: 0.010 Set refine NCS operators: 0.000 Set scattering table: 0.010 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.730 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 18.750 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.4792 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.005 0.123 3398 Z= 0.316 Angle : 0.917 13.923 4603 Z= 0.487 Chirality : 0.051 0.387 520 Planarity : 0.006 0.052 551 Dihedral : 13.466 85.696 1208 Min Nonbonded Distance : 2.233 Molprobity Statistics. All-atom Clashscore : 6.13 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.56 % Favored : 95.44 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.72 (0.41), residues: 395 helix: 0.26 (0.32), residues: 229 sheet: -1.55 (0.73), residues: 47 loop : -1.24 (0.56), residues: 119 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 110 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 110 time to evaluate : 0.386 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 110 average time/residue: 0.1913 time to fit residues: 24.5293 Evaluate side-chains 87 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 87 time to evaluate : 0.387 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5496 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 40 random chunks: chunk 33 optimal weight: 3.9990 chunk 30 optimal weight: 0.8980 chunk 16 optimal weight: 1.9990 chunk 10 optimal weight: 5.9990 chunk 20 optimal weight: 0.8980 chunk 31 optimal weight: 1.9990 chunk 12 optimal weight: 8.9990 chunk 19 optimal weight: 0.8980 chunk 23 optimal weight: 0.9980 chunk 36 optimal weight: 0.3980 chunk 11 optimal weight: 1.9990 overall best weight: 0.8180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 91 ASN A 135 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5186 moved from start: 0.3393 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.033 3398 Z= 0.208 Angle : 0.658 9.636 4603 Z= 0.339 Chirality : 0.043 0.188 520 Planarity : 0.006 0.046 551 Dihedral : 7.239 83.233 447 Min Nonbonded Distance : 2.206 Molprobity Statistics. All-atom Clashscore : 14.19 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.05 % Favored : 95.95 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.07 (0.42), residues: 395 helix: 0.93 (0.33), residues: 230 sheet: -1.25 (0.77), residues: 47 loop : -1.01 (0.56), residues: 118 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 101 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 101 time to evaluate : 0.405 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 101 average time/residue: 0.1774 time to fit residues: 21.2094 Evaluate side-chains 77 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 77 time to evaluate : 0.387 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5161 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 40 random chunks: chunk 30 optimal weight: 0.4980 chunk 24 optimal weight: 2.9990 chunk 10 optimal weight: 2.9990 chunk 36 optimal weight: 0.6980 chunk 39 optimal weight: 0.0270 chunk 32 optimal weight: 1.9990 chunk 12 optimal weight: 10.0000 chunk 29 optimal weight: 1.9990 chunk 35 optimal weight: 0.0970 chunk 27 optimal weight: 0.8980 chunk 18 optimal weight: 2.9990 overall best weight: 0.4436 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 119 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 135 GLN A 187 HIS A 356 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5214 moved from start: 0.4329 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.028 3398 Z= 0.174 Angle : 0.576 9.490 4603 Z= 0.298 Chirality : 0.039 0.144 520 Planarity : 0.004 0.029 551 Dihedral : 6.925 80.326 447 Min Nonbonded Distance : 2.249 Molprobity Statistics. All-atom Clashscore : 12.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.05 % Favored : 95.95 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.43 (0.43), residues: 395 helix: 1.21 (0.34), residues: 224 sheet: -0.78 (0.71), residues: 56 loop : -0.87 (0.60), residues: 115 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 99 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 99 time to evaluate : 0.464 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 99 average time/residue: 0.1544 time to fit residues: 18.5485 Evaluate side-chains 76 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 76 time to evaluate : 0.342 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.4929 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 40 random chunks: chunk 4 optimal weight: 3.9990 chunk 17 optimal weight: 1.9990 chunk 24 optimal weight: 3.9990 chunk 36 optimal weight: 2.9990 chunk 38 optimal weight: 3.9990 chunk 19 optimal weight: 0.5980 chunk 34 optimal weight: 1.9990 chunk 10 optimal weight: 2.9990 chunk 32 optimal weight: 0.9990 chunk 21 optimal weight: 4.9990 chunk 0 optimal weight: 4.9990 overall best weight: 1.7188 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 95 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 119 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 135 GLN A 143 HIS ** A 281 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 356 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5700 moved from start: 0.7470 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.038 3398 Z= 0.317 Angle : 0.836 13.792 4603 Z= 0.441 Chirality : 0.049 0.166 520 Planarity : 0.006 0.042 551 Dihedral : 8.022 85.458 447 Min Nonbonded Distance : 2.129 Molprobity Statistics. All-atom Clashscore : 24.65 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.57 % Favored : 94.43 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.50 (0.41), residues: 395 helix: 0.58 (0.34), residues: 212 sheet: -1.04 (0.69), residues: 55 loop : -1.41 (0.54), residues: 128 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 105 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 105 time to evaluate : 0.394 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 105 average time/residue: 0.1633 time to fit residues: 20.5371 Evaluate side-chains 87 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 87 time to evaluate : 0.388 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5052 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 40 random chunks: chunk 28 optimal weight: 0.2980 chunk 15 optimal weight: 0.8980 chunk 32 optimal weight: 0.6980 chunk 26 optimal weight: 2.9990 chunk 0 optimal weight: 4.9990 chunk 19 optimal weight: 1.9990 chunk 34 optimal weight: 2.9990 chunk 9 optimal weight: 2.9990 chunk 13 optimal weight: 6.9990 chunk 7 optimal weight: 0.4980 chunk 22 optimal weight: 0.9980 overall best weight: 0.6780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 95 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 135 GLN A 187 HIS A 356 GLN A 394 ASN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5637 moved from start: 0.8015 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.036 3398 Z= 0.195 Angle : 0.670 11.422 4603 Z= 0.343 Chirality : 0.044 0.243 520 Planarity : 0.005 0.046 551 Dihedral : 7.196 79.866 447 Min Nonbonded Distance : 2.208 Molprobity Statistics. All-atom Clashscore : 16.58 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.81 % Favored : 95.19 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.02 (0.42), residues: 395 helix: 1.10 (0.35), residues: 212 sheet: -1.06 (0.66), residues: 58 loop : -1.23 (0.57), residues: 125 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 102 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 102 time to evaluate : 0.712 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 102 average time/residue: 0.1680 time to fit residues: 21.0054 Evaluate side-chains 87 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 87 time to evaluate : 0.398 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5489 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 40 random chunks: chunk 9 optimal weight: 0.0870 chunk 38 optimal weight: 0.9990 chunk 32 optimal weight: 1.9990 chunk 17 optimal weight: 0.9990 chunk 3 optimal weight: 0.8980 chunk 12 optimal weight: 10.0000 chunk 20 optimal weight: 2.9990 chunk 37 optimal weight: 0.4980 chunk 4 optimal weight: 0.9990 chunk 22 optimal weight: 0.7980 chunk 28 optimal weight: 10.0000 overall best weight: 0.6560 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 135 GLN A 187 HIS A 356 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5644 moved from start: 0.8675 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.032 3398 Z= 0.184 Angle : 0.633 9.125 4603 Z= 0.332 Chirality : 0.042 0.158 520 Planarity : 0.005 0.045 551 Dihedral : 7.001 79.064 447 Min Nonbonded Distance : 2.169 Molprobity Statistics. All-atom Clashscore : 14.94 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.30 % Favored : 95.70 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.18 (0.43), residues: 395 helix: 1.27 (0.35), residues: 212 sheet: -1.24 (0.64), residues: 62 loop : -1.09 (0.58), residues: 121 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 101 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 101 time to evaluate : 0.383 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 101 average time/residue: 0.1645 time to fit residues: 20.1735 Evaluate side-chains 83 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 83 time to evaluate : 0.394 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5553 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 40 random chunks: chunk 21 optimal weight: 3.9990 chunk 32 optimal weight: 0.9990 chunk 38 optimal weight: 0.8980 chunk 24 optimal weight: 0.7980 chunk 23 optimal weight: 0.7980 chunk 17 optimal weight: 0.7980 chunk 15 optimal weight: 2.9990 chunk 11 optimal weight: 0.6980 chunk 7 optimal weight: 4.9990 chunk 26 optimal weight: 0.6980 chunk 19 optimal weight: 0.6980 overall best weight: 0.7380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 135 GLN A 187 HIS A 356 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5681 moved from start: 0.9180 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.023 3398 Z= 0.180 Angle : 0.648 9.230 4603 Z= 0.338 Chirality : 0.042 0.144 520 Planarity : 0.005 0.041 551 Dihedral : 6.912 79.108 447 Min Nonbonded Distance : 2.189 Molprobity Statistics. All-atom Clashscore : 15.09 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.56 % Favored : 95.44 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.25 (0.43), residues: 395 helix: 1.31 (0.35), residues: 212 sheet: -1.22 (0.64), residues: 58 loop : -1.04 (0.58), residues: 125 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 94 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 94 time to evaluate : 0.396 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 94 average time/residue: 0.1469 time to fit residues: 16.9807 Evaluate side-chains 85 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 85 time to evaluate : 0.346 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5240 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 40 random chunks: chunk 3 optimal weight: 0.2980 chunk 30 optimal weight: 0.8980 chunk 35 optimal weight: 0.3980 chunk 36 optimal weight: 0.0020 chunk 33 optimal weight: 3.9990 chunk 21 optimal weight: 0.0980 chunk 15 optimal weight: 2.9990 chunk 28 optimal weight: 1.9990 chunk 11 optimal weight: 0.9980 chunk 32 optimal weight: 0.6980 chunk 23 optimal weight: 3.9990 overall best weight: 0.2988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 51 ASN A 95 ASN A 135 GLN A 356 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5614 moved from start: 0.9469 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.037 3398 Z= 0.163 Angle : 0.608 7.636 4603 Z= 0.319 Chirality : 0.041 0.138 520 Planarity : 0.005 0.048 551 Dihedral : 6.701 76.619 447 Min Nonbonded Distance : 2.209 Molprobity Statistics. All-atom Clashscore : 12.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.80 % Favored : 96.20 % Rotamer Outliers : 0.28 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.50 (0.44), residues: 395 helix: 1.46 (0.36), residues: 213 sheet: -0.87 (0.67), residues: 58 loop : -0.96 (0.59), residues: 124 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 94 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 93 time to evaluate : 0.430 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 1 outliers final: 0 residues processed: 94 average time/residue: 0.1453 time to fit residues: 16.9792 Evaluate side-chains 86 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 86 time to evaluate : 0.527 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.8904 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 40 random chunks: chunk 37 optimal weight: 0.8980 chunk 23 optimal weight: 1.9990 chunk 18 optimal weight: 2.9990 chunk 26 optimal weight: 0.9990 chunk 39 optimal weight: 2.9990 chunk 36 optimal weight: 0.9980 chunk 31 optimal weight: 0.7980 chunk 3 optimal weight: 0.0270 chunk 24 optimal weight: 0.8980 chunk 19 optimal weight: 2.9990 chunk 25 optimal weight: 0.0770 overall best weight: 0.5396 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 135 GLN A 356 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5678 moved from start: 1.0039 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.037 3398 Z= 0.177 Angle : 0.622 7.905 4603 Z= 0.329 Chirality : 0.042 0.239 520 Planarity : 0.005 0.050 551 Dihedral : 6.709 78.089 447 Min Nonbonded Distance : 2.228 Molprobity Statistics. All-atom Clashscore : 13.30 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.05 % Favored : 95.95 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.58 (0.44), residues: 395 helix: 1.52 (0.36), residues: 213 sheet: -0.97 (0.68), residues: 62 loop : -0.83 (0.60), residues: 120 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 94 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 94 time to evaluate : 0.422 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 94 average time/residue: 0.1458 time to fit residues: 16.9513 Evaluate side-chains 87 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 87 time to evaluate : 0.409 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5709 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 40 random chunks: chunk 33 optimal weight: 2.9990 chunk 9 optimal weight: 0.6980 chunk 29 optimal weight: 3.9990 chunk 4 optimal weight: 2.9990 chunk 8 optimal weight: 4.9990 chunk 31 optimal weight: 0.6980 chunk 13 optimal weight: 3.9990 chunk 32 optimal weight: 0.8980 chunk 5 optimal weight: 0.5980 chunk 27 optimal weight: 0.9980 chunk 1 optimal weight: 1.9990 overall best weight: 0.7780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 135 GLN A 356 GLN ** A 465 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5739 moved from start: 1.0663 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.033 3398 Z= 0.196 Angle : 0.683 8.461 4603 Z= 0.361 Chirality : 0.043 0.217 520 Planarity : 0.005 0.037 551 Dihedral : 6.963 79.974 447 Min Nonbonded Distance : 2.147 Molprobity Statistics. All-atom Clashscore : 16.29 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.56 % Favored : 95.44 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.27 (0.44), residues: 395 helix: 1.20 (0.36), residues: 219 sheet: -1.04 (0.67), residues: 61 loop : -0.98 (0.61), residues: 115 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 790 Ramachandran restraints generated. 395 Oldfield, 0 Emsley, 395 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 103 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 103 time to evaluate : 0.395 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 103 average time/residue: 0.1492 time to fit residues: 18.8177 Evaluate side-chains 94 residues out of total 356 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 94 time to evaluate : 0.368 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5264 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 40 random chunks: chunk 22 optimal weight: 0.0770 chunk 36 optimal weight: 0.5980 chunk 21 optimal weight: 1.9990 chunk 27 optimal weight: 2.9990 chunk 1 optimal weight: 0.0470 chunk 25 optimal weight: 1.9990 chunk 24 optimal weight: 0.9990 chunk 23 optimal weight: 0.7980 chunk 15 optimal weight: 1.9990 chunk 14 optimal weight: 0.5980 chunk 39 optimal weight: 0.5980 overall best weight: 0.3836 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 135 GLN A 187 HIS A 356 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4362 r_free = 0.4362 target = 0.202058 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 39)----------------| | r_work = 0.4106 r_free = 0.4106 target = 0.177535 restraints weight = 5765.635| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 41)----------------| | r_work = 0.4179 r_free = 0.4179 target = 0.184285 restraints weight = 3470.879| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 35)----------------| | r_work = 0.4227 r_free = 0.4227 target = 0.188980 restraints weight = 2217.793| |-----------------------------------------------------------------------------| r_work (final): 0.4203 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5721 moved from start: 1.0861 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.033 3398 Z= 0.174 Angle : 0.655 7.758 4603 Z= 0.348 Chirality : 0.043 0.230 520 Planarity : 0.005 0.039 551 Dihedral : 6.702 77.818 447 Min Nonbonded Distance : 2.179 Molprobity Statistics. All-atom Clashscore : 12.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.81 % Favored : 95.19 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.49 (0.44), residues: 395 helix: 1.46 (0.36), residues: 213 sheet: -0.98 (0.69), residues: 60 loop : -0.91 (0.59), residues: 122 =============================================================================== Job complete usr+sys time: 1144.02 seconds wall clock time: 21 minutes 21.65 seconds (1281.65 seconds total)