Starting phenix.real_space_refine on Mon Nov 13 16:55:34 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gof_34177/11_2023/8gof_34177_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gof_34177/11_2023/8gof_34177.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.0 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gof_34177/11_2023/8gof_34177.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gof_34177/11_2023/8gof_34177.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gof_34177/11_2023/8gof_34177_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gof_34177/11_2023/8gof_34177_updated.pdb" } resolution = 3.0 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.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 Set scattering table Set to: electron Number of scattering types: 4 Type Number sf(0) Gaussians S 10 5.16 5 C 2192 2.51 5 N 545 2.21 5 O 552 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A GLU 67": "OE1" <-> "OE2" Residue "A ARG 256": "NH1" <-> "NH2" Residue "A GLU 257": "OE1" <-> "OE2" Residue "A GLU 294": "OE1" <-> "OE2" Time to flip residues: 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5144/modules/chem_data/mon_lib" Total number of atoms: 3299 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 3299 Number of conformers: 1 Conformer: "" Number of residues, atoms: 408, 3299 Unusual residues: {'LYA': 1} Classifications: {'peptide': 407, 'undetermined': 1} Link IDs: {'PTRANS': 14, 'TRANS': 392, None: 1} Not linked: pdbres="HIS A 457 " pdbres="LYA A 601 " Chain breaks: 1 Time building chain proxies: 2.34, per 1000 atoms: 0.71 Number of scatterers: 3299 At special positions: 0 Unit cell: (80, 72, 66, 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 10 16.00 O 552 8.00 N 545 7.00 C 2192 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: 1.66 Conformation dependent library (CDL) restraints added in 816.9 milliseconds 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 770 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 18 helices and 0 sheets defined 74.0% alpha, 0.0% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.34 Creating SS restraints... Processing helix chain 'A' and resid 26 through 40 removed outlier: 3.520A pdb=" N CYS A 30 " --> pdb=" O ARG A 26 " (cutoff:3.500A) removed outlier: 3.931A pdb=" N TYR A 31 " --> pdb=" O HIS A 27 " (cutoff:3.500A) removed outlier: 4.274A pdb=" N GLN A 40 " --> pdb=" O GLY A 36 " (cutoff:3.500A) Processing helix chain 'A' and resid 49 through 52 No H-bonds generated for 'chain 'A' and resid 49 through 52' Processing helix chain 'A' and resid 61 through 67 removed outlier: 3.517A pdb=" N ASN A 66 " --> pdb=" O GLU A 62 " (cutoff:3.500A) Processing helix chain 'A' and resid 69 through 89 removed outlier: 3.587A pdb=" N TYR A 74 " --> pdb=" O VAL A 71 " (cutoff:3.500A) removed outlier: 4.577A pdb=" N SER A 75 " --> pdb=" O LEU A 72 " (cutoff:3.500A) removed outlier: 3.597A pdb=" N TYR A 76 " --> pdb=" O SER A 73 " (cutoff:3.500A) removed outlier: 4.135A pdb=" N VAL A 79 " --> pdb=" O TYR A 76 " (cutoff:3.500A) removed outlier: 4.040A pdb=" N LEU A 80 " --> pdb=" O LEU A 77 " (cutoff:3.500A) removed outlier: 4.777A pdb=" N VAL A 81 " --> pdb=" O ALA A 78 " (cutoff:3.500A) Proline residue: A 82 - end of helix removed outlier: 4.001A pdb=" N THR A 87 " --> pdb=" O PHE A 84 " (cutoff:3.500A) removed outlier: 3.739A pdb=" N TYR A 89 " --> pdb=" O LEU A 86 " (cutoff:3.500A) Processing helix chain 'A' and resid 94 through 112 removed outlier: 3.524A pdb=" N PHE A 103 " --> pdb=" O GLN A 99 " (cutoff:3.500A) Processing helix chain 'A' and resid 116 through 131 removed outlier: 3.517A pdb=" N THR A 129 " --> pdb=" O PHE A 125 " (cutoff:3.500A) Processing helix chain 'A' and resid 134 through 143 removed outlier: 3.868A pdb=" N TYR A 139 " --> pdb=" O ALA A 135 " (cutoff:3.500A) removed outlier: 3.684A pdb=" N PHE A 141 " --> pdb=" O SER A 137 " (cutoff:3.500A) removed outlier: 3.586A pdb=" N SER A 142 " --> pdb=" O SER A 138 " (cutoff:3.500A) Processing helix chain 'A' and resid 149 through 175 removed outlier: 3.824A pdb=" N SER A 168 " --> pdb=" O VAL A 164 " (cutoff:3.500A) removed outlier: 3.533A pdb=" N VAL A 175 " --> pdb=" O GLY A 171 " (cutoff:3.500A) Processing helix chain 'A' and resid 182 through 201 Processing helix chain 'A' and resid 250 through 262 removed outlier: 3.886A pdb=" N ASP A 260 " --> pdb=" O ARG A 256 " (cutoff:3.500A) Processing helix chain 'A' and resid 267 through 276 Processing helix chain 'A' and resid 279 through 295 removed outlier: 3.558A pdb=" N TYR A 287 " --> pdb=" O LEU A 283 " (cutoff:3.500A) removed outlier: 4.970A pdb=" N HIS A 289 " --> pdb=" O VAL A 285 " (cutoff:3.500A) removed outlier: 5.210A pdb=" N ILE A 290 " --> pdb=" O TYR A 286 " (cutoff:3.500A) Processing helix chain 'A' and resid 301 through 303 No H-bonds generated for 'chain 'A' and resid 301 through 303' Processing helix chain 'A' and resid 306 through 326 removed outlier: 3.698A pdb=" N THR A 314 " --> pdb=" O ASP A 310 " (cutoff:3.500A) removed outlier: 3.846A pdb=" N GLY A 325 " --> pdb=" O SER A 321 " (cutoff:3.500A) removed outlier: 4.411A pdb=" N PHE A 326 " --> pdb=" O PHE A 322 " (cutoff:3.500A) Processing helix chain 'A' and resid 336 through 354 removed outlier: 3.549A pdb=" N THR A 345 " --> pdb=" O GLY A 341 " (cutoff:3.500A) removed outlier: 3.711A pdb=" N VAL A 350 " --> pdb=" O GLN A 346 " (cutoff:3.500A) Processing helix chain 'A' and resid 359 through 387 removed outlier: 3.670A pdb=" N ILE A 362 " --> pdb=" O PRO A 359 " (cutoff:3.500A) removed outlier: 4.357A pdb=" N TRP A 363 " --> pdb=" O SER A 360 " (cutoff:3.500A) removed outlier: 3.623A pdb=" N TYR A 366 " --> pdb=" O TRP A 363 " (cutoff:3.500A) removed outlier: 3.730A pdb=" N ALA A 368 " --> pdb=" O CYS A 365 " (cutoff:3.500A) removed outlier: 3.734A pdb=" N VAL A 370 " --> pdb=" O ALA A 367 " (cutoff:3.500A) removed outlier: 4.932A pdb=" N LEU A 371 " --> pdb=" O ALA A 368 " (cutoff:3.500A) removed outlier: 3.629A pdb=" N GLY A 374 " --> pdb=" O LEU A 371 " (cutoff:3.500A) removed outlier: 3.759A pdb=" N GLN A 377 " --> pdb=" O GLY A 374 " (cutoff:3.500A) removed outlier: 3.549A pdb=" N PHE A 378 " --> pdb=" O SER A 375 " (cutoff:3.500A) removed outlier: 4.292A pdb=" N VAL A 380 " --> pdb=" O GLN A 377 " (cutoff:3.500A) Proline residue: A 381 - end of helix removed outlier: 3.848A pdb=" N THR A 384 " --> pdb=" O PRO A 381 " (cutoff:3.500A) removed outlier: 3.556A pdb=" N PHE A 385 " --> pdb=" O ILE A 382 " (cutoff:3.500A) removed outlier: 3.575A pdb=" N GLN A 386 " --> pdb=" O ALA A 383 " (cutoff:3.500A) removed outlier: 4.267A pdb=" N ILE A 387 " --> pdb=" O THR A 384 " (cutoff:3.500A) Processing helix chain 'A' and resid 393 through 417 removed outlier: 3.520A pdb=" N ALA A 397 " --> pdb=" O LYS A 393 " (cutoff:3.500A) Processing helix chain 'A' and resid 428 through 452 removed outlier: 3.504A pdb=" N GLN A 433 " --> pdb=" O ARG A 429 " (cutoff:3.500A) 182 hydrogen bonds defined for protein. 498 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.40 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: 987 1.34 - 1.46: 765 1.46 - 1.57: 1623 1.57 - 1.69: 0 1.69 - 1.81: 16 Bond restraints: 3391 Sorted by residual: bond pdb=" N MET A 130 " pdb=" CA MET A 130 " ideal model delta sigma weight residual 1.457 1.486 -0.029 1.29e-02 6.01e+03 5.04e+00 bond pdb=" C12 LYA A 601 " pdb=" N11 LYA A 601 " ideal model delta sigma weight residual 1.340 1.370 -0.030 2.00e-02 2.50e+03 2.20e+00 bond pdb=" C10 LYA A 601 " pdb=" N11 LYA A 601 " ideal model delta sigma weight residual 1.384 1.355 0.029 2.00e-02 2.50e+03 2.16e+00 bond pdb=" C VAL A 81 " pdb=" O VAL A 81 " ideal model delta sigma weight residual 1.240 1.223 0.017 1.26e-02 6.30e+03 1.82e+00 bond pdb=" C17 LYA A 601 " pdb=" N16 LYA A 601 " ideal model delta sigma weight residual 1.305 1.331 -0.026 2.00e-02 2.50e+03 1.72e+00 ... (remaining 3386 not shown) Histogram of bond angle deviations from ideal: 98.30 - 105.57: 55 105.57 - 112.84: 1783 112.84 - 120.12: 1390 120.12 - 127.39: 1355 127.39 - 134.66: 39 Bond angle restraints: 4622 Sorted by residual: angle pdb=" N PHE A 103 " pdb=" CA PHE A 103 " pdb=" C PHE A 103 " ideal model delta sigma weight residual 111.40 105.48 5.92 1.22e+00 6.72e-01 2.35e+01 angle pdb=" C TYR A 92 " pdb=" CA TYR A 92 " pdb=" CB TYR A 92 " ideal model delta sigma weight residual 110.42 102.32 8.10 1.99e+00 2.53e-01 1.66e+01 angle pdb=" N LEU A 162 " pdb=" CA LEU A 162 " pdb=" C LEU A 162 " ideal model delta sigma weight residual 112.23 107.69 4.54 1.26e+00 6.30e-01 1.30e+01 angle pdb=" N ILE A 329 " pdb=" CA ILE A 329 " pdb=" C ILE A 329 " ideal model delta sigma weight residual 112.96 109.62 3.34 1.00e+00 1.00e+00 1.12e+01 angle pdb=" N VAL A 428 " pdb=" CA VAL A 428 " pdb=" C VAL A 428 " ideal model delta sigma weight residual 111.17 115.69 -4.52 1.38e+00 5.25e-01 1.07e+01 ... (remaining 4617 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 15.09: 1753 15.09 - 30.17: 144 30.17 - 45.25: 27 45.25 - 60.34: 13 60.34 - 75.42: 3 Dihedral angle restraints: 1940 sinusoidal: 742 harmonic: 1198 Sorted by residual: dihedral pdb=" CA ARG A 91 " pdb=" C ARG A 91 " pdb=" N TYR A 92 " pdb=" CA TYR A 92 " ideal model delta harmonic sigma weight residual 180.00 154.19 25.81 0 5.00e+00 4.00e-02 2.66e+01 dihedral pdb=" CB ARG A 373 " pdb=" CG ARG A 373 " pdb=" CD ARG A 373 " pdb=" NE ARG A 373 " ideal model delta sinusoidal sigma weight residual 60.00 118.74 -58.74 3 1.50e+01 4.44e-03 9.47e+00 dihedral pdb=" CA VAL A 418 " pdb=" C VAL A 418 " pdb=" N SER A 419 " pdb=" CA SER A 419 " ideal model delta harmonic sigma weight residual 180.00 -164.62 -15.38 0 5.00e+00 4.00e-02 9.46e+00 ... (remaining 1937 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.041: 343 0.041 - 0.081: 157 0.081 - 0.122: 26 0.122 - 0.163: 7 0.163 - 0.203: 1 Chirality restraints: 534 Sorted by residual: chirality pdb=" CA ARG A 42 " pdb=" N ARG A 42 " pdb=" C ARG A 42 " pdb=" CB ARG A 42 " both_signs ideal model delta sigma weight residual False 2.51 2.71 -0.20 2.00e-01 2.50e+01 1.03e+00 chirality pdb=" CB VAL A 288 " pdb=" CA VAL A 288 " pdb=" CG1 VAL A 288 " pdb=" CG2 VAL A 288 " both_signs ideal model delta sigma weight residual False -2.63 -2.48 -0.15 2.00e-01 2.50e+01 5.37e-01 chirality pdb=" CA LYS A 57 " pdb=" N LYS A 57 " pdb=" C LYS A 57 " pdb=" CB LYS A 57 " both_signs ideal model delta sigma weight residual False 2.51 2.37 0.14 2.00e-01 2.50e+01 4.88e-01 ... (remaining 531 not shown) Planarity restraints: 559 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C20 LYA A 601 " -0.041 2.00e-02 2.50e+03 1.04e-01 1.36e+02 pdb=" C23 LYA A 601 " 0.127 2.00e-02 2.50e+03 pdb=" C4 LYA A 601 " 0.097 2.00e-02 2.50e+03 pdb=" N22 LYA A 601 " -0.163 2.00e-02 2.50e+03 pdb=" O21 LYA A 601 " -0.021 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C THR A 69 " 0.044 5.00e-02 4.00e+02 6.50e-02 6.75e+00 pdb=" N PRO A 70 " -0.112 5.00e-02 4.00e+02 pdb=" CA PRO A 70 " 0.032 5.00e-02 4.00e+02 pdb=" CD PRO A 70 " 0.037 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CA ILE A 41 " 0.011 2.00e-02 2.50e+03 2.23e-02 4.96e+00 pdb=" C ILE A 41 " -0.039 2.00e-02 2.50e+03 pdb=" O ILE A 41 " 0.015 2.00e-02 2.50e+03 pdb=" N ARG A 42 " 0.013 2.00e-02 2.50e+03 ... (remaining 556 not shown) Histogram of nonbonded interaction distances: 2.27 - 2.79: 991 2.79 - 3.32: 2962 3.32 - 3.85: 5677 3.85 - 4.37: 7168 4.37 - 4.90: 11662 Nonbonded interactions: 28460 Sorted by model distance: nonbonded pdb=" O VAL A 83 " pdb=" OG1 THR A 87 " model vdw 2.269 2.440 nonbonded pdb=" OD1 ASP A 88 " pdb=" OH TYR A 92 " model vdw 2.271 2.440 nonbonded pdb=" O VAL A 304 " pdb=" OH TYR A 366 " model vdw 2.278 2.440 nonbonded pdb=" O LEU A 32 " pdb=" OH TYR A 136 " model vdw 2.298 2.440 nonbonded pdb=" O VAL A 380 " pdb=" OG1 THR A 384 " model vdw 2.324 2.440 ... (remaining 28455 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.380 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.040 Check model and map are aligned: 0.050 Set scattering table: 0.030 Process input model: 14.260 Find NCS groups from input model: 0.090 Set up NCS constraints: 0.010 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.010 Load rotamer database and sin/cos tables:12.140 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 35.020 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7749 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.006 0.056 3391 Z= 0.387 Angle : 0.799 8.103 4622 Z= 0.453 Chirality : 0.046 0.203 534 Planarity : 0.007 0.104 559 Dihedral : 12.708 75.420 1170 Min Nonbonded Distance : 2.269 Molprobity Statistics. All-atom Clashscore : 7.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.48 % Favored : 97.52 % Rotamer: Outliers : 0.00 % Allowed : 8.26 % Favored : 91.74 % 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: -6.52 (0.22), residues: 403 helix: -4.45 (0.11), residues: 322 sheet: None (None), residues: 0 loop : -2.16 (0.54), residues: 81 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 41 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 41 time to evaluate : 0.398 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 41 average time/residue: 1.3941 time to fit residues: 58.8515 Evaluate side-chains 28 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 28 time to evaluate : 0.362 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5132 ------------------------------------------------------------------------------- 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: 0.7980 chunk 30 optimal weight: 0.9980 chunk 16 optimal weight: 0.9990 chunk 10 optimal weight: 0.5980 chunk 20 optimal weight: 0.9980 chunk 31 optimal weight: 0.6980 chunk 12 optimal weight: 0.6980 chunk 19 optimal weight: 0.9980 chunk 23 optimal weight: 2.9990 chunk 36 optimal weight: 0.6980 chunk 11 optimal weight: 0.3980 overall best weight: 0.6180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 114 HIS ** A 172 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 293 ASN ** A 346 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7701 moved from start: 0.1787 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.035 3391 Z= 0.198 Angle : 0.513 5.679 4622 Z= 0.273 Chirality : 0.038 0.116 534 Planarity : 0.005 0.037 559 Dihedral : 6.054 55.624 461 Min Nonbonded Distance : 2.242 Molprobity Statistics. All-atom Clashscore : 8.55 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.49 % Favored : 98.51 % Rotamer: Outliers : 1.42 % Allowed : 11.40 % Favored : 87.18 % 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: -3.33 (0.32), residues: 403 helix: -2.11 (0.22), residues: 321 sheet: None (None), residues: 0 loop : -1.54 (0.56), residues: 82 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 44 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 39 time to evaluate : 0.351 Fit side-chains outliers start: 5 outliers final: 0 residues processed: 42 average time/residue: 1.2423 time to fit residues: 53.8532 Evaluate side-chains 36 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 36 time to evaluate : 0.373 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5104 ------------------------------------------------------------------------------- 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.6980 chunk 24 optimal weight: 0.6980 chunk 10 optimal weight: 0.5980 chunk 36 optimal weight: 0.5980 chunk 39 optimal weight: 1.9990 chunk 32 optimal weight: 0.6980 chunk 12 optimal weight: 0.5980 chunk 29 optimal weight: 0.5980 chunk 35 optimal weight: 0.1980 chunk 27 optimal weight: 3.9990 chunk 18 optimal weight: 0.8980 overall best weight: 0.5180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 114 HIS A 172 GLN ** A 346 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 358 HIS Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7683 moved from start: 0.2070 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.032 3391 Z= 0.174 Angle : 0.459 5.392 4622 Z= 0.241 Chirality : 0.037 0.114 534 Planarity : 0.003 0.029 559 Dihedral : 5.894 59.822 461 Min Nonbonded Distance : 2.231 Molprobity Statistics. All-atom Clashscore : 8.40 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.49 % Favored : 98.51 % Rotamer: Outliers : 1.71 % Allowed : 13.68 % Favored : 84.62 % 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.14 (0.38), residues: 403 helix: -0.48 (0.27), residues: 321 sheet: None (None), residues: 0 loop : -1.25 (0.58), residues: 82 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 41 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 35 time to evaluate : 0.392 Fit side-chains outliers start: 6 outliers final: 1 residues processed: 39 average time/residue: 1.3033 time to fit residues: 52.5054 Evaluate side-chains 34 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 33 time to evaluate : 0.366 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 1 residues processed: 0 time to fit residues: 0.4989 ------------------------------------------------------------------------------- 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: 0.9980 chunk 17 optimal weight: 0.6980 chunk 24 optimal weight: 1.9990 chunk 36 optimal weight: 0.7980 chunk 38 optimal weight: 0.8980 chunk 19 optimal weight: 3.9990 chunk 34 optimal weight: 0.5980 chunk 10 optimal weight: 0.5980 chunk 32 optimal weight: 1.9990 chunk 21 optimal weight: 1.9990 chunk 0 optimal weight: 3.9990 overall best weight: 0.7180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 114 HIS A 346 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7696 moved from start: 0.2214 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.034 3391 Z= 0.211 Angle : 0.466 5.338 4622 Z= 0.243 Chirality : 0.038 0.118 534 Planarity : 0.003 0.028 559 Dihedral : 5.796 56.648 461 Min Nonbonded Distance : 2.165 Molprobity Statistics. All-atom Clashscore : 8.55 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.99 % Favored : 98.01 % Rotamer: Outliers : 1.71 % Allowed : 15.38 % Favored : 82.91 % 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.12 (0.41), residues: 403 helix: 0.26 (0.29), residues: 321 sheet: None (None), residues: 0 loop : -1.02 (0.59), residues: 82 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 38 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 32 time to evaluate : 0.377 Fit side-chains outliers start: 6 outliers final: 2 residues processed: 37 average time/residue: 1.1915 time to fit residues: 45.6752 Evaluate side-chains 34 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 32 time to evaluate : 0.352 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 2 residues processed: 0 time to fit residues: 0.5044 ------------------------------------------------------------------------------- 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.9990 chunk 15 optimal weight: 1.9990 chunk 32 optimal weight: 0.9980 chunk 26 optimal weight: 3.9990 chunk 0 optimal weight: 3.9990 chunk 19 optimal weight: 3.9990 chunk 34 optimal weight: 0.6980 chunk 9 optimal weight: 0.9980 chunk 13 optimal weight: 0.7980 chunk 7 optimal weight: 2.9990 chunk 22 optimal weight: 1.9990 overall best weight: 0.8982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 114 HIS A 213 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7709 moved from start: 0.2232 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.037 3391 Z= 0.249 Angle : 0.493 6.570 4622 Z= 0.252 Chirality : 0.039 0.120 534 Planarity : 0.003 0.029 559 Dihedral : 5.791 53.583 461 Min Nonbonded Distance : 2.147 Molprobity Statistics. All-atom Clashscore : 8.10 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.99 % Favored : 98.01 % Rotamer: Outliers : 1.99 % Allowed : 14.53 % Favored : 83.48 % 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.21 (0.41), residues: 403 helix: 0.49 (0.29), residues: 321 sheet: None (None), residues: 0 loop : -0.88 (0.60), residues: 82 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 41 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 34 time to evaluate : 0.399 Fit side-chains outliers start: 7 outliers final: 2 residues processed: 39 average time/residue: 1.3137 time to fit residues: 52.9006 Evaluate side-chains 35 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 33 time to evaluate : 0.369 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 2 residues processed: 0 time to fit residues: 0.5042 ------------------------------------------------------------------------------- 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.9990 chunk 38 optimal weight: 0.5980 chunk 32 optimal weight: 1.9990 chunk 17 optimal weight: 0.6980 chunk 3 optimal weight: 0.5980 chunk 12 optimal weight: 0.0670 chunk 20 optimal weight: 1.9990 chunk 37 optimal weight: 0.6980 chunk 4 optimal weight: 0.8980 chunk 22 optimal weight: 1.9990 chunk 28 optimal weight: 1.9990 overall best weight: 0.5318 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 114 HIS A 277 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7678 moved from start: 0.2388 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.031 3391 Z= 0.173 Angle : 0.454 5.453 4622 Z= 0.235 Chirality : 0.037 0.117 534 Planarity : 0.003 0.027 559 Dihedral : 5.525 45.995 461 Min Nonbonded Distance : 2.214 Molprobity Statistics. All-atom Clashscore : 9.30 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.49 % Favored : 98.51 % Rotamer: Outliers : 1.71 % Allowed : 15.67 % Favored : 82.62 % 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.42), residues: 403 helix: 0.76 (0.30), residues: 321 sheet: None (None), residues: 0 loop : -0.80 (0.60), residues: 82 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 41 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 35 time to evaluate : 0.371 Fit side-chains outliers start: 6 outliers final: 2 residues processed: 41 average time/residue: 1.1523 time to fit residues: 48.9205 Evaluate side-chains 34 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 32 time to evaluate : 0.376 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 2 residues processed: 0 time to fit residues: 0.5063 ------------------------------------------------------------------------------- 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: 2.9990 chunk 32 optimal weight: 0.6980 chunk 38 optimal weight: 1.9990 chunk 24 optimal weight: 0.8980 chunk 23 optimal weight: 0.7980 chunk 17 optimal weight: 1.9990 chunk 15 optimal weight: 0.8980 chunk 11 optimal weight: 0.6980 chunk 7 optimal weight: 0.4980 chunk 26 optimal weight: 0.1980 chunk 19 optimal weight: 0.9980 overall best weight: 0.5780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 114 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7684 moved from start: 0.2431 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.031 3391 Z= 0.180 Angle : 0.448 5.222 4622 Z= 0.233 Chirality : 0.037 0.116 534 Planarity : 0.003 0.027 559 Dihedral : 5.653 49.269 461 Min Nonbonded Distance : 2.235 Molprobity Statistics. All-atom Clashscore : 8.10 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.74 % Favored : 98.26 % Rotamer: Outliers : 1.42 % Allowed : 16.24 % Favored : 82.34 % 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.85 (0.42), residues: 403 helix: 0.95 (0.30), residues: 321 sheet: None (None), residues: 0 loop : -0.71 (0.61), residues: 82 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 38 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 33 time to evaluate : 0.416 Fit side-chains outliers start: 5 outliers final: 2 residues processed: 38 average time/residue: 1.2550 time to fit residues: 49.3343 Evaluate side-chains 35 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 33 time to evaluate : 0.296 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 2 residues processed: 0 time to fit residues: 0.5135 ------------------------------------------------------------------------------- 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: 1.9990 chunk 30 optimal weight: 2.9990 chunk 35 optimal weight: 0.7980 chunk 36 optimal weight: 0.7980 chunk 33 optimal weight: 0.8980 chunk 21 optimal weight: 0.7980 chunk 15 optimal weight: 0.6980 chunk 28 optimal weight: 0.6980 chunk 11 optimal weight: 0.6980 chunk 32 optimal weight: 0.9980 chunk 23 optimal weight: 0.5980 overall best weight: 0.6980 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 114 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7692 moved from start: 0.2445 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 3391 Z= 0.204 Angle : 0.459 5.258 4622 Z= 0.238 Chirality : 0.038 0.117 534 Planarity : 0.003 0.028 559 Dihedral : 5.592 45.157 461 Min Nonbonded Distance : 2.236 Molprobity Statistics. All-atom Clashscore : 8.40 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.99 % Favored : 98.01 % Rotamer: Outliers : 1.71 % Allowed : 15.95 % Favored : 82.34 % 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.93 (0.42), residues: 403 helix: 1.02 (0.30), residues: 321 sheet: None (None), residues: 0 loop : -0.72 (0.61), residues: 82 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 41 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 35 time to evaluate : 0.353 Fit side-chains outliers start: 6 outliers final: 2 residues processed: 39 average time/residue: 1.2284 time to fit residues: 49.5771 Evaluate side-chains 37 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 35 time to evaluate : 0.373 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 2 residues processed: 0 time to fit residues: 0.5072 ------------------------------------------------------------------------------- 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.6980 chunk 23 optimal weight: 2.9990 chunk 18 optimal weight: 1.9990 chunk 26 optimal weight: 0.7980 chunk 39 optimal weight: 0.9990 chunk 36 optimal weight: 0.5980 chunk 31 optimal weight: 0.7980 chunk 3 optimal weight: 1.9990 chunk 24 optimal weight: 1.9990 chunk 19 optimal weight: 0.2980 chunk 25 optimal weight: 0.8980 overall best weight: 0.6380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 114 HIS A 277 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7686 moved from start: 0.2463 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.031 3391 Z= 0.193 Angle : 0.455 5.234 4622 Z= 0.237 Chirality : 0.038 0.116 534 Planarity : 0.003 0.027 559 Dihedral : 5.528 42.951 461 Min Nonbonded Distance : 2.240 Molprobity Statistics. All-atom Clashscore : 8.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.99 % Favored : 98.01 % Rotamer: Outliers : 1.42 % Allowed : 16.52 % Favored : 82.05 % 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.07 (0.42), residues: 403 helix: 1.12 (0.30), residues: 321 sheet: None (None), residues: 0 loop : -0.70 (0.61), residues: 82 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 38 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 33 time to evaluate : 0.373 Fit side-chains outliers start: 5 outliers final: 2 residues processed: 38 average time/residue: 1.1935 time to fit residues: 46.9495 Evaluate side-chains 35 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 33 time to evaluate : 0.372 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 2 residues processed: 0 time to fit residues: 0.5069 ------------------------------------------------------------------------------- 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: 0.1980 chunk 9 optimal weight: 0.0980 chunk 29 optimal weight: 1.9990 chunk 4 optimal weight: 0.4980 chunk 8 optimal weight: 0.5980 chunk 31 optimal weight: 1.9990 chunk 13 optimal weight: 1.9990 chunk 32 optimal weight: 0.5980 chunk 5 optimal weight: 0.7980 chunk 27 optimal weight: 0.7980 chunk 1 optimal weight: 0.6980 overall best weight: 0.3980 ------------------------------------------------------------------------------- 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.7659 moved from start: 0.2583 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.028 3391 Z= 0.149 Angle : 0.433 5.177 4622 Z= 0.226 Chirality : 0.036 0.112 534 Planarity : 0.003 0.026 559 Dihedral : 5.328 38.560 461 Min Nonbonded Distance : 2.255 Molprobity Statistics. All-atom Clashscore : 8.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.99 % Favored : 98.01 % Rotamer: Outliers : 0.85 % Allowed : 17.38 % Favored : 81.77 % 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.31 (0.42), residues: 403 helix: 1.30 (0.30), residues: 321 sheet: None (None), residues: 0 loop : -0.65 (0.61), residues: 82 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 806 Ramachandran restraints generated. 403 Oldfield, 0 Emsley, 403 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 32 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 29 time to evaluate : 0.366 Fit side-chains outliers start: 3 outliers final: 2 residues processed: 31 average time/residue: 1.3437 time to fit residues: 43.1417 Evaluate side-chains 31 residues out of total 351 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 29 time to evaluate : 0.372 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 2 residues processed: 0 time to fit residues: 0.5068 ------------------------------------------------------------------------------- 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.7980 chunk 36 optimal weight: 2.9990 chunk 21 optimal weight: 1.9990 chunk 27 optimal weight: 4.9990 chunk 1 optimal weight: 0.9980 chunk 25 optimal weight: 2.9990 chunk 24 optimal weight: 1.9990 chunk 23 optimal weight: 0.9980 chunk 15 optimal weight: 0.9990 chunk 14 optimal weight: 0.7980 chunk 39 optimal weight: 0.8980 overall best weight: 0.8980 ------------------------------------------------------------------------------- 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.3141 r_free = 0.3141 target = 0.102362 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 36)----------------| | r_work = 0.2920 r_free = 0.2920 target = 0.086760 restraints weight = 4133.194| |-----------------------------------------------------------------------------| r_work (start): 0.2930 rms_B_bonded: 1.65 r_work: 0.2806 rms_B_bonded: 2.36 restraints_weight: 0.5000 r_work: 0.2671 rms_B_bonded: 3.94 restraints_weight: 0.2500 r_work (final): 0.2671 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7864 moved from start: 0.2465 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.034 3391 Z= 0.244 Angle : 0.477 5.303 4622 Z= 0.247 Chirality : 0.039 0.120 534 Planarity : 0.003 0.028 559 Dihedral : 5.634 44.458 461 Min Nonbonded Distance : 2.219 Molprobity Statistics. All-atom Clashscore : 8.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.48 % Favored : 97.52 % Rotamer: Outliers : 0.85 % Allowed : 17.38 % Favored : 81.77 % 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.12 (0.42), residues: 403 helix: 1.16 (0.30), residues: 321 sheet: None (None), residues: 0 loop : -0.69 (0.61), residues: 82 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 1541.73 seconds wall clock time: 28 minutes 15.94 seconds (1695.94 seconds total)