Starting phenix.real_space_refine on Thu Dec 7 23:07:12 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8qkk_18464/12_2023/8qkk_18464.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8qkk_18464/12_2023/8qkk_18464.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.23 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8qkk_18464/12_2023/8qkk_18464.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8qkk_18464/12_2023/8qkk_18464.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8qkk_18464/12_2023/8qkk_18464.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8qkk_18464/12_2023/8qkk_18464.pdb" } resolution = 3.23 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.072 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 4 Type Number sf(0) Gaussians S 26 5.16 5 C 3622 2.51 5 N 922 2.21 5 O 1011 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A GLU 271": "OE1" <-> "OE2" Time to flip residues: 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/chem_data/mon_lib" Total number of atoms: 5581 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 5581 Number of conformers: 1 Conformer: "" Number of residues, atoms: 730, 5581 Classifications: {'peptide': 730} Link IDs: {'PTRANS': 36, 'TRANS': 693} Chain breaks: 1 Time building chain proxies: 3.32, per 1000 atoms: 0.59 Number of scatterers: 5581 At special positions: 0 Unit cell: (69.6, 71.775, 105.125, 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 26 16.00 O 1011 8.00 N 922 7.00 C 3622 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.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: 2.05 Conformation dependent library (CDL) restraints added in 998.7 milliseconds 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1338 Finding SS restraints... Secondary structure from input PDB file: 38 helices and 3 sheets defined 76.3% alpha, 4.2% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.55 Creating SS restraints... Processing helix chain 'A' and resid 1 through 12 removed outlier: 4.000A pdb=" N TRP A 5 " --> pdb=" O GLY A 1 " (cutoff:3.500A) Processing helix chain 'A' and resid 12 through 32 Processing helix chain 'A' and resid 33 through 36 Processing helix chain 'A' and resid 48 through 62 Processing helix chain 'A' and resid 64 through 67 removed outlier: 4.189A pdb=" N SER A 67 " --> pdb=" O ASP A 64 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 64 through 67' Processing helix chain 'A' and resid 75 through 79 removed outlier: 4.077A pdb=" N ASP A 78 " --> pdb=" O PRO A 75 " (cutoff:3.500A) Processing helix chain 'A' and resid 83 through 101 Processing helix chain 'A' and resid 112 through 114 No H-bonds generated for 'chain 'A' and resid 112 through 114' Processing helix chain 'A' and resid 119 through 124 Processing helix chain 'A' and resid 143 through 159 Proline residue: A 156 - end of helix Processing helix chain 'A' and resid 160 through 164 Processing helix chain 'A' and resid 171 through 205 removed outlier: 3.845A pdb=" N ILE A 182 " --> pdb=" O LEU A 178 " (cutoff:3.500A) removed outlier: 3.626A pdb=" N GLY A 183 " --> pdb=" O THR A 179 " (cutoff:3.500A) removed outlier: 4.705A pdb=" N ILE A 194 " --> pdb=" O GLU A 190 " (cutoff:3.500A) Proline residue: A 195 - end of helix Processing helix chain 'A' and resid 208 through 211 Processing helix chain 'A' and resid 212 through 234 Processing helix chain 'A' and resid 241 through 272 removed outlier: 4.091A pdb=" N LEU A 251 " --> pdb=" O THR A 247 " (cutoff:3.500A) Processing helix chain 'A' and resid 274 through 285 Processing helix chain 'A' and resid 285 through 302 Processing helix chain 'A' and resid 303 through 307 Processing helix chain 'A' and resid 309 through 332 Processing helix chain 'A' and resid 332 through 342 Processing helix chain 'A' and resid 343 through 346 removed outlier: 3.748A pdb=" N ASP A 346 " --> pdb=" O PRO A 343 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 343 through 346' Processing helix chain 'A' and resid 350 through 356 removed outlier: 4.159A pdb=" N ILE A 356 " --> pdb=" O THR A 352 " (cutoff:3.500A) Processing helix chain 'A' and resid 363 through 374 removed outlier: 3.691A pdb=" N ARG A 367 " --> pdb=" O GLN A 363 " (cutoff:3.500A) Processing helix chain 'A' and resid 385 through 400 Processing helix chain 'A' and resid 400 through 416 Proline residue: A 407 - end of helix Processing helix chain 'A' and resid 417 through 422 Processing helix chain 'A' and resid 428 through 432 Processing helix chain 'A' and resid 436 through 449 Processing helix chain 'A' and resid 470 through 484 removed outlier: 3.725A pdb=" N VAL A 484 " --> pdb=" O LYS A 480 " (cutoff:3.500A) Processing helix chain 'A' and resid 493 through 497 removed outlier: 3.857A pdb=" N LYS A 496 " --> pdb=" O ASP A 493 " (cutoff:3.500A) removed outlier: 3.766A pdb=" N MET A 497 " --> pdb=" O PRO A 494 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 493 through 497' Processing helix chain 'A' and resid 525 through 536 removed outlier: 3.572A pdb=" N ASP A 531 " --> pdb=" O ALA A 527 " (cutoff:3.500A) Processing helix chain 'A' and resid 548 through 584 Proline residue: A 565 - end of helix Processing helix chain 'A' and resid 587 through 612 Processing helix chain 'A' and resid 629 through 660 removed outlier: 3.829A pdb=" N GLY A 633 " --> pdb=" O ALA A 629 " (cutoff:3.500A) removed outlier: 3.595A pdb=" N ALA A 637 " --> pdb=" O GLY A 633 " (cutoff:3.500A) removed outlier: 4.320A pdb=" N ARG A 653 " --> pdb=" O PHE A 649 " (cutoff:3.500A) removed outlier: 4.065A pdb=" N MET A 654 " --> pdb=" O LEU A 650 " (cutoff:3.500A) Processing helix chain 'A' and resid 663 through 695 removed outlier: 4.910A pdb=" N ARG A 677 " --> pdb=" O ALA A 673 " (cutoff:3.500A) removed outlier: 4.980A pdb=" N LEU A 678 " --> pdb=" O THR A 674 " (cutoff:3.500A) removed outlier: 3.683A pdb=" N PHE A 693 " --> pdb=" O VAL A 689 " (cutoff:3.500A) removed outlier: 3.727A pdb=" N PHE A 695 " --> pdb=" O GLY A 691 " (cutoff:3.500A) Processing helix chain 'A' and resid 698 through 722 removed outlier: 4.639A pdb=" N ILE A 719 " --> pdb=" O ASP A 715 " (cutoff:3.500A) removed outlier: 4.887A pdb=" N ARG A 720 " --> pdb=" O ALA A 716 " (cutoff:3.500A) removed outlier: 3.752A pdb=" N PHE A 722 " --> pdb=" O ILE A 718 " (cutoff:3.500A) Processing helix chain 'A' and resid 722 through 732 removed outlier: 3.605A pdb=" N ALA A 726 " --> pdb=" O PHE A 722 " (cutoff:3.500A) Processing helix chain 'A' and resid 733 through 737 Processing sheet with id=AA1, first strand: chain 'A' and resid 109 through 110 Processing sheet with id=AA2, first strand: chain 'A' and resid 424 through 425 Processing sheet with id=AA3, first strand: chain 'A' and resid 499 through 500 404 hydrogen bonds defined for protein. 1158 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.65 Time building geometry restraints manager: 2.23 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.23 - 1.34: 1799 1.34 - 1.46: 1211 1.46 - 1.58: 2636 1.58 - 1.69: 0 1.69 - 1.81: 50 Bond restraints: 5696 Sorted by residual: bond pdb=" N GLY A 1 " pdb=" CA GLY A 1 " ideal model delta sigma weight residual 1.451 1.492 -0.041 1.60e-02 3.91e+03 6.44e+00 bond pdb=" N LEU A 248 " pdb=" CA LEU A 248 " ideal model delta sigma weight residual 1.468 1.458 0.010 1.24e-02 6.50e+03 6.49e-01 bond pdb=" CB LYS A 187 " pdb=" CG LYS A 187 " ideal model delta sigma weight residual 1.520 1.539 -0.019 3.00e-02 1.11e+03 4.16e-01 bond pdb=" CA ILE A 194 " pdb=" CB ILE A 194 " ideal model delta sigma weight residual 1.539 1.536 0.003 5.40e-03 3.43e+04 3.66e-01 bond pdb=" CA TRP A 362 " pdb=" CB TRP A 362 " ideal model delta sigma weight residual 1.530 1.540 -0.010 1.69e-02 3.50e+03 3.62e-01 ... (remaining 5691 not shown) Histogram of bond angle deviations from ideal: 98.44 - 105.56: 148 105.56 - 112.67: 3122 112.67 - 119.79: 1806 119.79 - 126.90: 2596 126.90 - 134.02: 82 Bond angle restraints: 7754 Sorted by residual: angle pdb=" N VAL A 586 " pdb=" CA VAL A 586 " pdb=" C VAL A 586 " ideal model delta sigma weight residual 112.96 110.87 2.09 1.00e+00 1.00e+00 4.37e+00 angle pdb=" C THR A 133 " pdb=" N PHE A 134 " pdb=" CA PHE A 134 " ideal model delta sigma weight residual 121.75 118.99 2.76 1.73e+00 3.34e-01 2.55e+00 angle pdb=" C ASN A 361 " pdb=" N TRP A 362 " pdb=" CA TRP A 362 " ideal model delta sigma weight residual 121.54 124.52 -2.98 1.91e+00 2.74e-01 2.43e+00 angle pdb=" CB LYS A 187 " pdb=" CG LYS A 187 " pdb=" CD LYS A 187 " ideal model delta sigma weight residual 111.30 114.68 -3.38 2.30e+00 1.89e-01 2.16e+00 angle pdb=" C VAL A 302 " pdb=" N PRO A 303 " pdb=" CD PRO A 303 " ideal model delta sigma weight residual 125.00 130.92 -5.92 4.10e+00 5.95e-02 2.08e+00 ... (remaining 7749 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.18: 3063 17.18 - 34.36: 260 34.36 - 51.55: 41 51.55 - 68.73: 13 68.73 - 85.91: 11 Dihedral angle restraints: 3388 sinusoidal: 1309 harmonic: 2079 Sorted by residual: dihedral pdb=" CA THR A 717 " pdb=" C THR A 717 " pdb=" N ILE A 718 " pdb=" CA ILE A 718 " ideal model delta harmonic sigma weight residual 180.00 -163.57 -16.43 0 5.00e+00 4.00e-02 1.08e+01 dihedral pdb=" CB GLU A 190 " pdb=" CG GLU A 190 " pdb=" CD GLU A 190 " pdb=" OE1 GLU A 190 " ideal model delta sinusoidal sigma weight residual 0.00 -85.91 85.91 1 3.00e+01 1.11e-03 9.91e+00 dihedral pdb=" CB MET A 701 " pdb=" CG MET A 701 " pdb=" SD MET A 701 " pdb=" CE MET A 701 " ideal model delta sinusoidal sigma weight residual -180.00 -122.31 -57.69 3 1.50e+01 4.44e-03 9.45e+00 ... (remaining 3385 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.025: 546 0.025 - 0.049: 223 0.049 - 0.074: 94 0.074 - 0.099: 39 0.099 - 0.123: 24 Chirality restraints: 926 Sorted by residual: chirality pdb=" CA VAL A 460 " pdb=" N VAL A 460 " pdb=" C VAL A 460 " pdb=" CB VAL A 460 " both_signs ideal model delta sigma weight residual False 2.44 2.56 -0.12 2.00e-01 2.50e+01 3.81e-01 chirality pdb=" CA ILE A 72 " pdb=" N ILE A 72 " pdb=" C ILE A 72 " pdb=" CB ILE A 72 " both_signs ideal model delta sigma weight residual False 2.43 2.56 -0.12 2.00e-01 2.50e+01 3.74e-01 chirality pdb=" CA VAL A 484 " pdb=" N VAL A 484 " pdb=" C VAL A 484 " pdb=" CB VAL A 484 " both_signs ideal model delta sigma weight residual False 2.44 2.56 -0.12 2.00e-01 2.50e+01 3.47e-01 ... (remaining 923 not shown) Planarity restraints: 978 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA ALA A 221 " 0.007 2.00e-02 2.50e+03 1.40e-02 1.96e+00 pdb=" C ALA A 221 " -0.024 2.00e-02 2.50e+03 pdb=" O ALA A 221 " 0.009 2.00e-02 2.50e+03 pdb=" N ILE A 222 " 0.008 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA VAL A 246 " -0.007 2.00e-02 2.50e+03 1.37e-02 1.87e+00 pdb=" C VAL A 246 " 0.024 2.00e-02 2.50e+03 pdb=" O VAL A 246 " -0.009 2.00e-02 2.50e+03 pdb=" N THR A 247 " -0.008 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA ALA A 223 " -0.007 2.00e-02 2.50e+03 1.34e-02 1.79e+00 pdb=" C ALA A 223 " 0.023 2.00e-02 2.50e+03 pdb=" O ALA A 223 " -0.009 2.00e-02 2.50e+03 pdb=" N GLY A 224 " -0.008 2.00e-02 2.50e+03 ... (remaining 975 not shown) Histogram of nonbonded interaction distances: 2.13 - 2.69: 93 2.69 - 3.24: 6018 3.24 - 3.79: 9052 3.79 - 4.35: 11642 4.35 - 4.90: 20149 Nonbonded interactions: 46954 Sorted by model distance: nonbonded pdb=" OG1 THR A 127 " pdb=" OD1 ASP A 129 " model vdw 2.132 2.440 nonbonded pdb=" NH1 ARG A 167 " pdb=" OE1 GLN A 444 " model vdw 2.290 2.520 nonbonded pdb=" OG SER A 293 " pdb=" OD1 ASP A 645 " model vdw 2.296 2.440 nonbonded pdb=" O GLN A 473 " pdb=" OD1 ASP A 476 " model vdw 2.309 3.040 nonbonded pdb=" NH2 ARG A 288 " pdb=" O ALA A 582 " model vdw 2.328 2.520 ... (remaining 46949 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.770 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.010 Extract box with map and model: 7.290 Check model and map are aligned: 0.080 Set scattering table: 0.050 Process input model: 18.180 Find NCS groups from input model: 0.140 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.590 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 30.130 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7565 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.002 0.041 5696 Z= 0.147 Angle : 0.493 5.917 7754 Z= 0.257 Chirality : 0.038 0.123 926 Planarity : 0.003 0.027 978 Dihedral : 14.137 85.910 2050 Min Nonbonded Distance : 2.132 Molprobity Statistics. All-atom Clashscore : 9.00 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.89 % Favored : 97.11 % Rotamer: Outliers : 0.00 % Allowed : 0.17 % Favored : 99.83 % 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: 2.28 (0.34), residues: 726 helix: 1.88 (0.25), residues: 490 sheet: 1.27 (0.85), residues: 41 loop : 0.79 (0.51), residues: 195 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP A 736 HIS 0.002 0.001 HIS A 68 PHE 0.014 0.001 PHE A 134 TYR 0.009 0.001 TYR A 151 ARG 0.002 0.000 ARG A 367 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 92 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 92 time to evaluate : 0.614 Fit side-chains revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 92 average time/residue: 0.1785 time to fit residues: 21.7578 Evaluate side-chains 79 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 79 time to evaluate : 0.753 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.7852 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 73 random chunks: chunk 61 optimal weight: 3.9990 chunk 55 optimal weight: 0.6980 chunk 30 optimal weight: 1.9990 chunk 18 optimal weight: 2.9990 chunk 37 optimal weight: 6.9990 chunk 29 optimal weight: 0.6980 chunk 57 optimal weight: 0.9990 chunk 22 optimal weight: 1.9990 chunk 34 optimal weight: 3.9990 chunk 42 optimal weight: 1.9990 chunk 66 optimal weight: 2.9990 overall best weight: 1.2786 ------------------------------------------------------------------------------- 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.7706 moved from start: 0.1143 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.065 5696 Z= 0.365 Angle : 0.640 7.336 7754 Z= 0.323 Chirality : 0.045 0.183 926 Planarity : 0.004 0.032 978 Dihedral : 4.117 16.011 774 Min Nonbonded Distance : 2.244 Molprobity Statistics. All-atom Clashscore : 8.83 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.62 % Favored : 97.38 % Rotamer: Outliers : 1.51 % Allowed : 8.24 % Favored : 90.25 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.72 (0.32), residues: 726 helix: 1.54 (0.24), residues: 500 sheet: 1.07 (0.85), residues: 41 loop : 0.26 (0.49), residues: 185 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A 736 HIS 0.005 0.001 HIS A 37 PHE 0.021 0.002 PHE A 311 TYR 0.010 0.001 TYR A 151 ARG 0.006 0.001 ARG A 281 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 90 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 81 time to evaluate : 0.624 Fit side-chains revert: symmetry clash outliers start: 9 outliers final: 7 residues processed: 86 average time/residue: 0.1725 time to fit residues: 19.8901 Evaluate side-chains 85 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 78 time to evaluate : 0.633 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.0582 time to fit residues: 1.5925 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 73 random chunks: chunk 36 optimal weight: 0.5980 chunk 20 optimal weight: 0.6980 chunk 55 optimal weight: 0.9980 chunk 45 optimal weight: 0.0060 chunk 18 optimal weight: 0.9990 chunk 66 optimal weight: 2.9990 chunk 71 optimal weight: 1.9990 chunk 59 optimal weight: 0.6980 chunk 65 optimal weight: 3.9990 chunk 22 optimal weight: 0.6980 chunk 53 optimal weight: 0.9990 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 128 GLN A 625 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7652 moved from start: 0.1315 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.041 5696 Z= 0.213 Angle : 0.554 7.189 7754 Z= 0.276 Chirality : 0.041 0.158 926 Planarity : 0.004 0.032 978 Dihedral : 3.989 16.790 774 Min Nonbonded Distance : 2.242 Molprobity Statistics. All-atom Clashscore : 8.12 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.34 % Favored : 97.66 % Rotamer: Outliers : 1.01 % Allowed : 9.92 % Favored : 89.08 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.89 (0.32), residues: 726 helix: 1.73 (0.24), residues: 494 sheet: 0.93 (0.88), residues: 41 loop : 0.21 (0.48), residues: 191 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP A 736 HIS 0.006 0.001 HIS A 37 PHE 0.017 0.001 PHE A 311 TYR 0.010 0.001 TYR A 151 ARG 0.006 0.000 ARG A 281 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 91 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 85 time to evaluate : 0.617 Fit side-chains revert: symmetry clash outliers start: 6 outliers final: 2 residues processed: 91 average time/residue: 0.1732 time to fit residues: 21.0423 Evaluate side-chains 77 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 75 time to evaluate : 0.622 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.2126 time to fit residues: 1.3839 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 73 random chunks: chunk 65 optimal weight: 10.0000 chunk 49 optimal weight: 5.9990 chunk 34 optimal weight: 0.4980 chunk 7 optimal weight: 0.1980 chunk 31 optimal weight: 0.5980 chunk 44 optimal weight: 0.0980 chunk 66 optimal weight: 0.6980 chunk 70 optimal weight: 0.9980 chunk 63 optimal weight: 0.9990 chunk 19 optimal weight: 0.9980 chunk 58 optimal weight: 0.9990 overall best weight: 0.4180 ------------------------------------------------------------------------------- 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.7635 moved from start: 0.1472 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.039 5696 Z= 0.194 Angle : 0.549 6.968 7754 Z= 0.271 Chirality : 0.041 0.297 926 Planarity : 0.004 0.032 978 Dihedral : 3.967 16.806 774 Min Nonbonded Distance : 2.229 Molprobity Statistics. All-atom Clashscore : 8.74 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.75 % Favored : 97.25 % Rotamer: Outliers : 1.34 % Allowed : 10.92 % Favored : 87.73 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.94 (0.32), residues: 726 helix: 1.81 (0.24), residues: 494 sheet: 0.79 (0.89), residues: 41 loop : 0.16 (0.47), residues: 191 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A 736 HIS 0.005 0.001 HIS A 37 PHE 0.017 0.001 PHE A 311 TYR 0.010 0.001 TYR A 151 ARG 0.004 0.000 ARG A 281 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 90 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 82 time to evaluate : 0.603 Fit side-chains outliers start: 8 outliers final: 2 residues processed: 88 average time/residue: 0.1676 time to fit residues: 19.8957 Evaluate side-chains 80 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 78 time to evaluate : 0.611 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0565 time to fit residues: 1.0975 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 73 random chunks: chunk 40 optimal weight: 1.9990 chunk 1 optimal weight: 0.5980 chunk 52 optimal weight: 0.5980 chunk 29 optimal weight: 0.9990 chunk 60 optimal weight: 3.9990 chunk 48 optimal weight: 5.9990 chunk 0 optimal weight: 0.5980 chunk 36 optimal weight: 0.9990 chunk 63 optimal weight: 1.9990 chunk 17 optimal weight: 2.9990 chunk 23 optimal weight: 0.5980 overall best weight: 0.6782 ------------------------------------------------------------------------------- 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.7664 moved from start: 0.1639 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.038 5696 Z= 0.227 Angle : 0.568 8.278 7754 Z= 0.281 Chirality : 0.041 0.216 926 Planarity : 0.004 0.032 978 Dihedral : 3.981 16.594 774 Min Nonbonded Distance : 2.248 Molprobity Statistics. All-atom Clashscore : 8.92 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.48 % Favored : 97.52 % Rotamer: Outliers : 1.01 % Allowed : 12.27 % Favored : 86.72 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.81 (0.32), residues: 726 helix: 1.69 (0.24), residues: 493 sheet: 0.86 (0.89), residues: 41 loop : 0.12 (0.47), residues: 192 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP A 736 HIS 0.005 0.001 HIS A 37 PHE 0.017 0.001 PHE A 311 TYR 0.010 0.001 TYR A 646 ARG 0.004 0.000 ARG A 281 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 86 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 80 time to evaluate : 0.633 Fit side-chains outliers start: 6 outliers final: 3 residues processed: 84 average time/residue: 0.1801 time to fit residues: 20.2696 Evaluate side-chains 80 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 77 time to evaluate : 0.567 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0666 time to fit residues: 1.2039 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 73 random chunks: chunk 63 optimal weight: 0.9980 chunk 13 optimal weight: 0.8980 chunk 41 optimal weight: 0.9990 chunk 17 optimal weight: 1.9990 chunk 70 optimal weight: 0.5980 chunk 58 optimal weight: 0.9990 chunk 32 optimal weight: 0.7980 chunk 5 optimal weight: 2.9990 chunk 23 optimal weight: 0.7980 chunk 37 optimal weight: 5.9990 chunk 68 optimal weight: 0.9980 overall best weight: 0.8180 ------------------------------------------------------------------------------- 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.7687 moved from start: 0.1762 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.039 5696 Z= 0.249 Angle : 0.590 9.112 7754 Z= 0.289 Chirality : 0.043 0.404 926 Planarity : 0.004 0.031 978 Dihedral : 4.022 16.609 774 Min Nonbonded Distance : 2.263 Molprobity Statistics. All-atom Clashscore : 9.09 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.20 % Favored : 97.80 % Rotamer: Outliers : 0.67 % Allowed : 12.94 % Favored : 86.39 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.74 (0.32), residues: 726 helix: 1.65 (0.24), residues: 493 sheet: 0.77 (0.89), residues: 41 loop : 0.09 (0.46), residues: 192 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP A 736 HIS 0.004 0.001 HIS A 37 PHE 0.018 0.001 PHE A 311 TYR 0.012 0.001 TYR A 646 ARG 0.003 0.000 ARG A 266 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 90 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 86 time to evaluate : 0.597 Fit side-chains revert: symmetry clash outliers start: 4 outliers final: 1 residues processed: 90 average time/residue: 0.1775 time to fit residues: 21.2294 Evaluate side-chains 85 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 84 time to evaluate : 0.615 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0542 time to fit residues: 0.9481 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 73 random chunks: chunk 7 optimal weight: 0.6980 chunk 40 optimal weight: 0.9990 chunk 51 optimal weight: 0.6980 chunk 39 optimal weight: 1.9990 chunk 59 optimal weight: 0.8980 chunk 70 optimal weight: 0.9980 chunk 44 optimal weight: 1.9990 chunk 42 optimal weight: 0.4980 chunk 32 optimal weight: 0.7980 chunk 43 optimal weight: 0.9990 chunk 28 optimal weight: 0.1980 overall best weight: 0.5780 ------------------------------------------------------------------------------- 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.7666 moved from start: 0.1845 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.032 5696 Z= 0.208 Angle : 0.585 12.118 7754 Z= 0.285 Chirality : 0.042 0.330 926 Planarity : 0.004 0.032 978 Dihedral : 4.011 17.364 774 Min Nonbonded Distance : 2.220 Molprobity Statistics. All-atom Clashscore : 9.80 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.34 % Favored : 97.66 % Rotamer: Outliers : 1.51 % Allowed : 13.61 % Favored : 84.87 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.82 (0.32), residues: 726 helix: 1.70 (0.24), residues: 493 sheet: 0.76 (0.90), residues: 41 loop : 0.15 (0.46), residues: 192 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.009 0.001 TRP A 736 HIS 0.005 0.001 HIS A 37 PHE 0.017 0.001 PHE A 311 TYR 0.010 0.001 TYR A 646 ARG 0.003 0.000 ARG A 281 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 93 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 84 time to evaluate : 0.595 Fit side-chains outliers start: 9 outliers final: 5 residues processed: 88 average time/residue: 0.1762 time to fit residues: 20.7385 Evaluate side-chains 89 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 84 time to evaluate : 0.623 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.0689 time to fit residues: 1.4611 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 73 random chunks: chunk 42 optimal weight: 0.0270 chunk 21 optimal weight: 0.9980 chunk 13 optimal weight: 0.5980 chunk 44 optimal weight: 1.9990 chunk 47 optimal weight: 0.9980 chunk 34 optimal weight: 3.9990 chunk 6 optimal weight: 0.0370 chunk 55 optimal weight: 0.6980 chunk 64 optimal weight: 0.9980 chunk 67 optimal weight: 0.8980 chunk 61 optimal weight: 0.7980 overall best weight: 0.4316 ------------------------------------------------------------------------------- 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.7644 moved from start: 0.1906 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.030 5696 Z= 0.193 Angle : 0.608 14.999 7754 Z= 0.290 Chirality : 0.042 0.274 926 Planarity : 0.004 0.031 978 Dihedral : 4.002 17.613 774 Min Nonbonded Distance : 2.229 Molprobity Statistics. All-atom Clashscore : 9.71 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.34 % Favored : 97.66 % Rotamer: Outliers : 0.84 % Allowed : 14.62 % Favored : 84.54 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.84 (0.32), residues: 726 helix: 1.72 (0.24), residues: 493 sheet: 0.75 (0.90), residues: 41 loop : 0.16 (0.46), residues: 192 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP A 736 HIS 0.005 0.001 HIS A 37 PHE 0.017 0.001 PHE A 311 TYR 0.012 0.001 TYR A 646 ARG 0.003 0.000 ARG A 281 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 88 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 83 time to evaluate : 0.562 Fit side-chains outliers start: 5 outliers final: 1 residues processed: 86 average time/residue: 0.1705 time to fit residues: 19.7137 Evaluate side-chains 79 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 78 time to evaluate : 0.599 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0668 time to fit residues: 0.9730 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 73 random chunks: chunk 65 optimal weight: 6.9990 chunk 67 optimal weight: 2.9990 chunk 39 optimal weight: 0.8980 chunk 28 optimal weight: 0.5980 chunk 51 optimal weight: 1.9990 chunk 20 optimal weight: 0.9980 chunk 59 optimal weight: 2.9990 chunk 62 optimal weight: 0.8980 chunk 43 optimal weight: 0.7980 chunk 69 optimal weight: 0.9990 chunk 42 optimal weight: 0.0670 overall best weight: 0.6518 ------------------------------------------------------------------------------- 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.7667 moved from start: 0.1960 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 5696 Z= 0.226 Angle : 0.616 14.851 7754 Z= 0.294 Chirality : 0.042 0.268 926 Planarity : 0.004 0.032 978 Dihedral : 4.010 17.399 774 Min Nonbonded Distance : 2.240 Molprobity Statistics. All-atom Clashscore : 10.06 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.20 % Favored : 97.80 % Rotamer: Outliers : 1.01 % Allowed : 14.79 % Favored : 84.20 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.81 (0.32), residues: 726 helix: 1.70 (0.24), residues: 493 sheet: 0.68 (0.89), residues: 41 loop : 0.18 (0.47), residues: 192 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.009 0.001 TRP A 736 HIS 0.005 0.001 HIS A 37 PHE 0.017 0.001 PHE A 311 TYR 0.015 0.001 TYR A 646 ARG 0.002 0.000 ARG A 266 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 92 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 86 time to evaluate : 0.570 Fit side-chains outliers start: 6 outliers final: 3 residues processed: 89 average time/residue: 0.1818 time to fit residues: 21.4535 Evaluate side-chains 88 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 85 time to evaluate : 0.638 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0670 time to fit residues: 1.1884 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 73 random chunks: chunk 32 optimal weight: 0.5980 chunk 48 optimal weight: 9.9990 chunk 72 optimal weight: 0.0980 chunk 66 optimal weight: 4.9990 chunk 57 optimal weight: 0.9980 chunk 6 optimal weight: 0.5980 chunk 44 optimal weight: 0.0070 chunk 35 optimal weight: 4.9990 chunk 46 optimal weight: 0.4980 chunk 61 optimal weight: 0.9980 chunk 17 optimal weight: 0.6980 overall best weight: 0.3598 ------------------------------------------------------------------------------- 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.7631 moved from start: 0.2035 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.028 5696 Z= 0.186 Angle : 0.610 14.731 7754 Z= 0.287 Chirality : 0.041 0.237 926 Planarity : 0.004 0.032 978 Dihedral : 3.969 17.877 774 Min Nonbonded Distance : 2.225 Molprobity Statistics. All-atom Clashscore : 9.71 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.48 % Favored : 97.52 % Rotamer: Outliers : 0.50 % Allowed : 15.46 % Favored : 84.03 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.92 (0.32), residues: 726 helix: 1.77 (0.24), residues: 493 sheet: 0.79 (0.91), residues: 41 loop : 0.20 (0.46), residues: 192 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP A 736 HIS 0.005 0.001 HIS A 37 PHE 0.016 0.001 PHE A 311 TYR 0.012 0.001 TYR A 646 ARG 0.003 0.000 ARG A 281 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1452 Ramachandran restraints generated. 726 Oldfield, 0 Emsley, 726 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 87 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 84 time to evaluate : 0.564 Fit side-chains outliers start: 3 outliers final: 1 residues processed: 86 average time/residue: 0.1720 time to fit residues: 19.8517 Evaluate side-chains 82 residues out of total 595 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 81 time to evaluate : 0.644 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0571 time to fit residues: 1.0082 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 73 random chunks: chunk 53 optimal weight: 1.9990 chunk 8 optimal weight: 0.9980 chunk 16 optimal weight: 1.9990 chunk 58 optimal weight: 0.8980 chunk 24 optimal weight: 0.8980 chunk 59 optimal weight: 0.9980 chunk 7 optimal weight: 0.8980 chunk 10 optimal weight: 3.9990 chunk 50 optimal weight: 2.9990 chunk 3 optimal weight: 0.8980 chunk 41 optimal weight: 0.0000 overall best weight: 0.7184 ------------------------------------------------------------------------------- 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.4295 r_free = 0.4295 target = 0.212525 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 70)----------------| | r_work = 0.4024 r_free = 0.4024 target = 0.181304 restraints weight = 5659.628| |-----------------------------------------------------------------------------| r_work (start): 0.3975 rms_B_bonded: 1.61 r_work: 0.3798 rms_B_bonded: 3.10 restraints_weight: 0.5000 r_work: 0.3658 rms_B_bonded: 4.98 restraints_weight: 0.2500 r_work (final): 0.3658 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7773 moved from start: 0.2071 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.035 5696 Z= 0.237 Angle : 0.634 15.107 7754 Z= 0.300 Chirality : 0.043 0.240 926 Planarity : 0.004 0.031 978 Dihedral : 4.028 17.301 774 Min Nonbonded Distance : 2.254 Molprobity Statistics. All-atom Clashscore : 10.33 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.34 % Favored : 97.66 % Rotamer: Outliers : 0.17 % Allowed : 16.64 % Favored : 83.19 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.79 (0.32), residues: 726 helix: 1.69 (0.24), residues: 493 sheet: 0.70 (0.90), residues: 41 loop : 0.14 (0.46), residues: 192 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.001 TRP A 736 HIS 0.006 0.001 HIS A 37 PHE 0.017 0.001 PHE A 311 TYR 0.014 0.001 TYR A 646 ARG 0.003 0.000 ARG A 367 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 1524.88 seconds wall clock time: 28 minutes 6.62 seconds (1686.62 seconds total)