Starting phenix.real_space_refine (version: 1.21rc1) on Thu Apr 20 00:30:41 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6jnc_9847/04_2023/6jnc_9847_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6jnc_9847/04_2023/6jnc_9847.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.7 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6jnc_9847/04_2023/6jnc_9847.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6jnc_9847/04_2023/6jnc_9847.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6jnc_9847/04_2023/6jnc_9847_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6jnc_9847/04_2023/6jnc_9847_updated.pdb" } resolution = 3.7 write_initial_geo_file = False refinement { macro_cycles = 10 } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= 0.002 sd= 0.014 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 15 5.16 5 C 2103 2.51 5 N 561 2.21 5 O 615 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A GLU 78": "OE1" <-> "OE2" Residue "A TYR 250": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ASP 290": "OD1" <-> "OD2" Residue "A TYR 363": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 381": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 396": "NH1" <-> "NH2" Time to flip residues: 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-4938/modules/chem_data/mon_lib" Total number of atoms: 3294 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 3294 Number of conformers: 1 Conformer: "" Number of residues, atoms: 420, 3294 Unusual residues: {'TRS': 1} Inner-chain residues flagged as termini: ['pdbres="LYS A 419 "'] Classifications: {'peptide': 419, 'undetermined': 1} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 20, 'TRANS': 398, None: 1} Not linked: pdbres="LYS A 419 " pdbres="TRS A 901 " Unresolved non-hydrogen bonds: 3 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 2 Planarities with less than four sites: {'ASP:plan': 1} Unresolved non-hydrogen planarities: 3 Time building chain proxies: 2.30, per 1000 atoms: 0.70 Number of scatterers: 3294 At special positions: 0 Unit cell: (71.34, 75.69, 74.82, 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 15 16.00 O 615 8.00 N 561 7.00 C 2103 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.54 Conformation dependent library (CDL) restraints added in 526.3 milliseconds 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 764 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 18 helices and 4 sheets defined 41.8% alpha, 11.9% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.40 Creating SS restraints... Processing helix chain 'A' and resid 6 through 18 removed outlier: 3.611A pdb=" N GLN A 13 " --> pdb=" O MET A 9 " (cutoff:3.500A) removed outlier: 4.227A pdb=" N LEU A 14 " --> pdb=" O ALA A 10 " (cutoff:3.500A) Processing helix chain 'A' and resid 25 through 28 No H-bonds generated for 'chain 'A' and resid 25 through 28' Processing helix chain 'A' and resid 81 through 90 removed outlier: 3.902A pdb=" N LEU A 86 " --> pdb=" O THR A 82 " (cutoff:3.500A) removed outlier: 3.534A pdb=" N THR A 88 " --> pdb=" O LYS A 84 " (cutoff:3.500A) 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 128 removed outlier: 3.834A pdb=" N ILE A 128 " --> pdb=" O ALA A 124 " (cutoff:3.500A) Processing helix chain 'A' and resid 132 through 134 No H-bonds generated for 'chain 'A' and resid 132 through 134' Processing helix chain 'A' and resid 152 through 164 removed outlier: 3.617A pdb=" N MET A 157 " --> pdb=" O MET A 153 " (cutoff:3.500A) Processing helix chain 'A' and resid 170 through 172 No H-bonds generated for 'chain 'A' and resid 170 through 172' Processing helix chain 'A' and resid 179 through 181 No H-bonds generated for 'chain 'A' and resid 179 through 181' Processing helix chain 'A' and resid 190 through 206 removed outlier: 3.982A pdb=" N ALA A 195 " --> pdb=" O THR A 191 " (cutoff:3.500A) Processing helix chain 'A' and resid 223 through 233 removed outlier: 3.530A pdb=" N THR A 227 " --> pdb=" O ALA A 223 " (cutoff:3.500A) removed outlier: 3.619A pdb=" N ALA A 231 " --> pdb=" O THR A 227 " (cutoff:3.500A) removed outlier: 3.803A pdb=" N LYS A 232 " --> pdb=" O ALA A 228 " (cutoff:3.500A) removed outlier: 3.565A pdb=" N GLU A 233 " --> pdb=" O LYS A 229 " (cutoff:3.500A) Processing helix chain 'A' and resid 257 through 266 removed outlier: 4.105A pdb=" N LEU A 261 " --> pdb=" O PRO A 257 " (cutoff:3.500A) removed outlier: 3.646A pdb=" N LYS A 262 " --> pdb=" O ASP A 258 " (cutoff:3.500A) Processing helix chain 'A' and resid 282 through 284 No H-bonds generated for 'chain 'A' and resid 282 through 284' Processing helix chain 'A' and resid 324 through 332 removed outlier: 3.961A pdb=" N ASP A 328 " --> pdb=" O PRO A 324 " (cutoff:3.500A) removed outlier: 3.784A pdb=" N ILE A 329 " --> pdb=" O GLU A 325 " (cutoff:3.500A) removed outlier: 3.746A pdb=" N GLU A 332 " --> pdb=" O ASP A 328 " (cutoff:3.500A) Processing helix chain 'A' and resid 340 through 343 No H-bonds generated for 'chain 'A' and resid 340 through 343' Processing helix chain 'A' and resid 346 through 360 removed outlier: 3.542A pdb=" N SER A 351 " --> pdb=" O GLY A 347 " (cutoff:3.500A) removed outlier: 4.072A pdb=" N TRP A 355 " --> pdb=" O SER A 351 " (cutoff:3.500A) removed outlier: 3.560A pdb=" N GLN A 357 " --> pdb=" O PHE A 353 " (cutoff:3.500A) removed outlier: 3.904A pdb=" N ASN A 358 " --> pdb=" O GLU A 354 " (cutoff:3.500A) Processing helix chain 'A' and resid 367 through 391 removed outlier: 3.534A pdb=" N LYS A 372 " --> pdb=" O GLU A 368 " (cutoff:3.500A) removed outlier: 4.001A pdb=" N LYS A 375 " --> pdb=" O GLU A 371 " (cutoff:3.500A) removed outlier: 3.618A pdb=" N LYS A 376 " --> pdb=" O LYS A 372 " (cutoff:3.500A) removed outlier: 3.534A pdb=" N LYS A 379 " --> pdb=" O LYS A 375 " (cutoff:3.500A) removed outlier: 3.626A pdb=" N VAL A 384 " --> pdb=" O ALA A 380 " (cutoff:3.500A) removed outlier: 4.018A pdb=" N ASN A 386 " --> pdb=" O TRP A 382 " (cutoff:3.500A) removed outlier: 3.518A pdb=" N THR A 387 " --> pdb=" O GLU A 383 " (cutoff:3.500A) removed outlier: 3.543A pdb=" N ASP A 390 " --> pdb=" O ASN A 386 " (cutoff:3.500A) removed outlier: 3.511A pdb=" N LYS A 391 " --> pdb=" O THR A 387 " (cutoff:3.500A) Processing helix chain 'A' and resid 395 through 413 removed outlier: 3.575A pdb=" N TYR A 400 " --> pdb=" O ARG A 396 " (cutoff:3.500A) removed outlier: 3.604A pdb=" N VAL A 401 " --> pdb=" O ASP A 397 " (cutoff:3.500A) removed outlier: 3.582A pdb=" N VAL A 402 " --> pdb=" O ALA A 398 " (cutoff:3.500A) removed outlier: 3.526A pdb=" N ARG A 406 " --> pdb=" O VAL A 402 " (cutoff:3.500A) removed outlier: 3.848A pdb=" N VAL A 407 " --> pdb=" O ALA A 403 " (cutoff:3.500A) Processing sheet with id= A, first strand: chain 'A' and resid 36 through 45 removed outlier: 3.521A pdb=" N ARG A 36 " --> pdb=" O GLN A 60 " (cutoff:3.500A) removed outlier: 6.599A pdb=" N GLY A 103 " --> pdb=" O HIS A 61 " (cutoff:3.500A) removed outlier: 7.556A pdb=" N GLY A 104 " --> pdb=" O PRO A 67 " (cutoff:3.500A) removed outlier: 5.726A pdb=" N LYS A 69 " --> pdb=" O GLY A 104 " (cutoff:3.500A) removed outlier: 8.161A pdb=" N GLY A 106 " --> pdb=" O LYS A 69 " (cutoff:3.500A) Processing sheet with id= B, first strand: chain 'A' and resid 72 through 75 removed outlier: 6.654A pdb=" N GLY A 106 " --> pdb=" O ARG A 73 " (cutoff:3.500A) removed outlier: 7.771A pdb=" N HIS A 75 " --> pdb=" O GLY A 106 " (cutoff:3.500A) removed outlier: 5.693A pdb=" N ILE A 108 " --> pdb=" O HIS A 75 " (cutoff:3.500A) No H-bonds generated for sheet with id= B Processing sheet with id= C, first strand: chain 'A' and resid 291 through 294 removed outlier: 6.493A pdb=" N THR A 238 " --> pdb=" O ILE A 215 " (cutoff:3.500A) removed outlier: 8.164A pdb=" N VAL A 217 " --> pdb=" O THR A 238 " (cutoff:3.500A) removed outlier: 8.234A pdb=" N VAL A 240 " --> pdb=" O VAL A 217 " (cutoff:3.500A) Processing sheet with id= D, first strand: chain 'A' and resid 313 through 315 removed outlier: 6.716A pdb=" N LEU A 336 " --> pdb=" O VAL A 314 " (cutoff:3.500A) No H-bonds generated for sheet with id= D 92 hydrogen bonds defined for protein. 273 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.83 Time building geometry restraints manager: 1.51 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.34: 1061 1.34 - 1.45: 435 1.45 - 1.57: 1841 1.57 - 1.69: 0 1.69 - 1.80: 29 Bond restraints: 3366 Sorted by residual: bond pdb=" CA PHE A 341 " pdb=" C PHE A 341 " ideal model delta sigma weight residual 1.523 1.467 0.056 1.48e-02 4.57e+03 1.45e+01 bond pdb=" N TRS A 901 " pdb=" C TRS A 901 " ideal model delta sigma weight residual 1.501 1.462 0.039 2.00e-02 2.50e+03 3.79e+00 bond pdb=" CB VAL A 356 " pdb=" CG2 VAL A 356 " ideal model delta sigma weight residual 1.521 1.468 0.053 3.30e-02 9.18e+02 2.59e+00 bond pdb=" CB TRP A 382 " pdb=" CG TRP A 382 " ideal model delta sigma weight residual 1.498 1.449 0.049 3.10e-02 1.04e+03 2.49e+00 bond pdb=" C PHE A 341 " pdb=" N LEU A 342 " ideal model delta sigma weight residual 1.335 1.315 0.020 1.42e-02 4.96e+03 1.96e+00 ... (remaining 3361 not shown) Histogram of bond angle deviations from ideal: 99.04 - 106.08: 88 106.08 - 113.13: 1844 113.13 - 120.17: 1228 120.17 - 127.21: 1352 127.21 - 134.25: 50 Bond angle restraints: 4562 Sorted by residual: angle pdb=" N ILE A 181 " pdb=" CA ILE A 181 " pdb=" C ILE A 181 " ideal model delta sigma weight residual 113.20 107.26 5.94 9.60e-01 1.09e+00 3.83e+01 angle pdb=" C ASP A 340 " pdb=" N PHE A 341 " pdb=" CA PHE A 341 " ideal model delta sigma weight residual 122.65 114.88 7.77 1.60e+00 3.91e-01 2.36e+01 angle pdb=" N VAL A 97 " pdb=" CA VAL A 97 " pdb=" C VAL A 97 " ideal model delta sigma weight residual 112.17 108.44 3.73 9.50e-01 1.11e+00 1.54e+01 angle pdb=" C LYS A 271 " pdb=" N ASP A 272 " pdb=" CA ASP A 272 " ideal model delta sigma weight residual 122.19 127.16 -4.97 1.41e+00 5.03e-01 1.24e+01 angle pdb=" C VAL A 146 " pdb=" N TYR A 147 " pdb=" CA TYR A 147 " ideal model delta sigma weight residual 123.05 128.57 -5.52 1.57e+00 4.06e-01 1.24e+01 ... (remaining 4557 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 16.63: 1792 16.63 - 33.27: 164 33.27 - 49.90: 47 49.90 - 66.53: 11 66.53 - 83.17: 2 Dihedral angle restraints: 2016 sinusoidal: 818 harmonic: 1198 Sorted by residual: dihedral pdb=" CA TRP A 382 " pdb=" C TRP A 382 " pdb=" N GLU A 383 " pdb=" CA GLU A 383 " ideal model delta harmonic sigma weight residual 180.00 162.00 18.00 0 5.00e+00 4.00e-02 1.30e+01 dihedral pdb=" CA ASP A 290 " pdb=" CB ASP A 290 " pdb=" CG ASP A 290 " pdb=" OD1 ASP A 290 " ideal model delta sinusoidal sigma weight residual -30.00 -90.08 60.08 1 2.00e+01 2.50e-03 1.20e+01 dihedral pdb=" CA TYR A 147 " pdb=" C TYR A 147 " pdb=" N THR A 148 " pdb=" CA THR A 148 " ideal model delta harmonic sigma weight residual 180.00 163.58 16.42 0 5.00e+00 4.00e-02 1.08e+01 ... (remaining 2013 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.035: 259 0.035 - 0.070: 161 0.070 - 0.105: 48 0.105 - 0.139: 22 0.139 - 0.174: 7 Chirality restraints: 497 Sorted by residual: chirality pdb=" CB ILE A 23 " pdb=" CA ILE A 23 " pdb=" CG1 ILE A 23 " pdb=" CG2 ILE A 23 " both_signs ideal model delta sigma weight residual False 2.64 2.47 0.17 2.00e-01 2.50e+01 7.58e-01 chirality pdb=" CA PHE A 341 " pdb=" N PHE A 341 " pdb=" C PHE A 341 " pdb=" CB PHE A 341 " both_signs ideal model delta sigma weight residual False 2.51 2.34 0.17 2.00e-01 2.50e+01 7.21e-01 chirality pdb=" CA TYR A 147 " pdb=" N TYR A 147 " pdb=" C TYR A 147 " pdb=" CB TYR A 147 " both_signs ideal model delta sigma weight residual False 2.51 2.35 0.16 2.00e-01 2.50e+01 6.40e-01 ... (remaining 494 not shown) Planarity restraints: 584 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C ASN A 149 " 0.044 5.00e-02 4.00e+02 6.53e-02 6.83e+00 pdb=" N PRO A 150 " -0.113 5.00e-02 4.00e+02 pdb=" CA PRO A 150 " 0.033 5.00e-02 4.00e+02 pdb=" CD PRO A 150 " 0.037 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C ILE A 141 " 0.042 5.00e-02 4.00e+02 6.33e-02 6.42e+00 pdb=" N PRO A 142 " -0.110 5.00e-02 4.00e+02 pdb=" CA PRO A 142 " 0.033 5.00e-02 4.00e+02 pdb=" CD PRO A 142 " 0.035 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C LEU A 99 " 0.040 5.00e-02 4.00e+02 6.03e-02 5.83e+00 pdb=" N PRO A 100 " -0.104 5.00e-02 4.00e+02 pdb=" CA PRO A 100 " 0.031 5.00e-02 4.00e+02 pdb=" CD PRO A 100 " 0.033 5.00e-02 4.00e+02 ... (remaining 581 not shown) Histogram of nonbonded interaction distances: 2.30 - 2.82: 1079 2.82 - 3.34: 2659 3.34 - 3.86: 5082 3.86 - 4.38: 6004 4.38 - 4.90: 10445 Nonbonded interactions: 25269 Sorted by model distance: nonbonded pdb=" N ASP A 340 " pdb=" OD1 ASP A 340 " model vdw 2.304 2.520 nonbonded pdb=" OD2 ASP A 145 " pdb=" O2 TRS A 901 " model vdw 2.368 2.440 nonbonded pdb=" N ASP A 256 " pdb=" OD1 ASP A 256 " model vdw 2.370 2.520 nonbonded pdb=" O PHE A 353 " pdb=" NE2 GLN A 357 " model vdw 2.372 2.520 nonbonded pdb=" N ASN A 149 " pdb=" OD1 ASN A 149 " model vdw 2.377 2.520 ... (remaining 25264 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 group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.480 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: 2.630 Check model and map are aligned: 0.050 Set scattering table: 0.030 Process input model: 13.850 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.000 Load rotamer database and sin/cos tables:3.010 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 20.160 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7626 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.009 0.062 3366 Z= 0.541 Angle : 0.872 7.773 4562 Z= 0.560 Chirality : 0.052 0.174 497 Planarity : 0.008 0.065 584 Dihedral : 14.724 83.169 1252 Min Nonbonded Distance : 2.304 Molprobity Statistics. All-atom Clashscore : 3.18 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.80 % Favored : 95.20 % Rotamer Outliers : 9.09 % 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: -5.69 (0.30), residues: 417 helix: -4.85 (0.13), residues: 188 sheet: -1.75 (0.68), residues: 61 loop : -2.99 (0.40), residues: 168 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 114 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 31 poor density : 83 time to evaluate : 0.414 Fit side-chains revert: symmetry clash outliers start: 31 outliers final: 5 residues processed: 105 average time/residue: 0.1801 time to fit residues: 22.3075 Evaluate side-chains 55 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 50 time to evaluate : 0.397 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.0890 time to fit residues: 1.1197 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 42 random chunks: chunk 35 optimal weight: 0.9990 chunk 31 optimal weight: 7.9990 chunk 17 optimal weight: 2.9990 chunk 10 optimal weight: 0.8980 chunk 21 optimal weight: 10.0000 chunk 32 optimal weight: 20.0000 chunk 12 optimal weight: 0.9990 chunk 20 optimal weight: 5.9990 chunk 24 optimal weight: 9.9990 chunk 38 optimal weight: 1.9990 chunk 11 optimal weight: 6.9990 overall best weight: 1.5788 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 60 GLN A 193 GLN A 386 ASN A 394 HIS Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8225 moved from start: 0.3948 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.030 3366 Z= 0.189 Angle : 0.538 5.084 4562 Z= 0.289 Chirality : 0.046 0.141 497 Planarity : 0.005 0.041 584 Dihedral : 4.373 14.618 458 Min Nonbonded Distance : 2.241 Molprobity Statistics. All-atom Clashscore : 6.51 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.40 % Favored : 97.60 % Rotamer Outliers : 4.11 % 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.89 (0.35), residues: 417 helix: -3.04 (0.28), residues: 184 sheet: -1.50 (0.66), residues: 57 loop : -2.23 (0.42), residues: 176 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 66 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 52 time to evaluate : 0.402 Fit side-chains outliers start: 14 outliers final: 3 residues processed: 64 average time/residue: 0.1427 time to fit residues: 11.5540 Evaluate side-chains 46 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 43 time to evaluate : 0.409 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0359 time to fit residues: 0.7509 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 42 random chunks: chunk 31 optimal weight: 8.9990 chunk 25 optimal weight: 9.9990 chunk 10 optimal weight: 8.9990 chunk 38 optimal weight: 0.9990 chunk 41 optimal weight: 0.0980 chunk 34 optimal weight: 5.9990 chunk 37 optimal weight: 3.9990 chunk 13 optimal weight: 5.9990 chunk 30 optimal weight: 10.0000 chunk 28 optimal weight: 0.9980 chunk 19 optimal weight: 1.9990 overall best weight: 1.6186 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 60 GLN A 193 GLN ** A 357 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 388 HIS Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8448 moved from start: 0.5320 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.026 3366 Z= 0.173 Angle : 0.469 3.890 4562 Z= 0.253 Chirality : 0.044 0.135 497 Planarity : 0.004 0.035 584 Dihedral : 3.992 15.461 458 Min Nonbonded Distance : 2.290 Molprobity Statistics. All-atom Clashscore : 5.00 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.12 % Favored : 96.88 % Rotamer Outliers : 2.64 % 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.90 (0.37), residues: 417 helix: -1.90 (0.34), residues: 188 sheet: -1.04 (0.68), residues: 57 loop : -2.10 (0.43), residues: 172 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 54 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 45 time to evaluate : 0.420 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 9 outliers final: 5 residues processed: 52 average time/residue: 0.1587 time to fit residues: 10.3340 Evaluate side-chains 46 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 41 time to evaluate : 0.364 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.0338 time to fit residues: 0.8338 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 42 random chunks: chunk 4 optimal weight: 10.0000 chunk 18 optimal weight: 9.9990 chunk 25 optimal weight: 9.9990 chunk 38 optimal weight: 5.9990 chunk 40 optimal weight: 7.9990 chunk 20 optimal weight: 0.0070 chunk 36 optimal weight: 5.9990 chunk 10 optimal weight: 5.9990 chunk 33 optimal weight: 8.9990 chunk 23 optimal weight: 10.0000 chunk 0 optimal weight: 10.0000 overall best weight: 5.2006 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 357 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 392 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8641 moved from start: 0.7020 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.007 0.041 3366 Z= 0.419 Angle : 0.610 6.882 4562 Z= 0.321 Chirality : 0.050 0.144 497 Planarity : 0.005 0.031 584 Dihedral : 4.866 29.199 458 Min Nonbonded Distance : 2.237 Molprobity Statistics. All-atom Clashscore : 4.39 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.60 % Favored : 96.40 % Rotamer Outliers : 1.47 % 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.57 (0.39), residues: 417 helix: -1.50 (0.36), residues: 182 sheet: -0.98 (0.72), residues: 57 loop : -2.04 (0.44), residues: 178 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 41 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 36 time to evaluate : 0.427 Fit side-chains revert: symmetry clash outliers start: 5 outliers final: 4 residues processed: 40 average time/residue: 0.1620 time to fit residues: 8.2306 Evaluate side-chains 37 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 33 time to evaluate : 0.416 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.0385 time to fit residues: 0.8253 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 42 random chunks: chunk 30 optimal weight: 9.9990 chunk 16 optimal weight: 2.9990 chunk 34 optimal weight: 1.9990 chunk 28 optimal weight: 0.0030 chunk 0 optimal weight: 10.0000 chunk 20 optimal weight: 10.0000 chunk 36 optimal weight: 3.9990 chunk 10 optimal weight: 6.9990 chunk 13 optimal weight: 2.9990 chunk 8 optimal weight: 3.9990 chunk 23 optimal weight: 10.0000 overall best weight: 2.3998 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 62 ASN ** A 357 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8592 moved from start: 0.7251 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.020 3366 Z= 0.205 Angle : 0.491 7.004 4562 Z= 0.259 Chirality : 0.046 0.151 497 Planarity : 0.004 0.031 584 Dihedral : 4.317 28.040 458 Min Nonbonded Distance : 2.273 Molprobity Statistics. All-atom Clashscore : 4.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.88 % Favored : 97.12 % Rotamer Outliers : 1.17 % 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.17 (0.40), residues: 417 helix: -1.07 (0.37), residues: 187 sheet: -0.91 (0.72), residues: 57 loop : -1.92 (0.44), residues: 173 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 43 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 39 time to evaluate : 0.413 Fit side-chains outliers start: 4 outliers final: 1 residues processed: 41 average time/residue: 0.1649 time to fit residues: 8.6734 Evaluate side-chains 36 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 35 time to evaluate : 0.419 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0619 time to fit residues: 0.6587 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 42 random chunks: chunk 10 optimal weight: 6.9990 chunk 40 optimal weight: 2.9990 chunk 33 optimal weight: 0.9980 chunk 18 optimal weight: 10.0000 chunk 3 optimal weight: 5.9990 chunk 13 optimal weight: 0.6980 chunk 21 optimal weight: 5.9990 chunk 39 optimal weight: 7.9990 chunk 4 optimal weight: 0.8980 chunk 23 optimal weight: 6.9990 chunk 29 optimal weight: 0.0980 overall best weight: 1.1382 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 357 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8547 moved from start: 0.7383 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.020 3366 Z= 0.131 Angle : 0.451 5.328 4562 Z= 0.233 Chirality : 0.045 0.154 497 Planarity : 0.003 0.029 584 Dihedral : 3.974 22.352 458 Min Nonbonded Distance : 2.265 Molprobity Statistics. All-atom Clashscore : 3.94 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.88 % Favored : 97.12 % Rotamer Outliers : 1.47 % 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.56 (0.40), residues: 417 helix: -0.45 (0.40), residues: 181 sheet: -0.57 (0.72), residues: 56 loop : -1.74 (0.43), residues: 180 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 42 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 37 time to evaluate : 0.415 Fit side-chains outliers start: 5 outliers final: 3 residues processed: 39 average time/residue: 0.1654 time to fit residues: 8.1960 Evaluate side-chains 35 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 32 time to evaluate : 0.411 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0363 time to fit residues: 0.7860 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 42 random chunks: chunk 22 optimal weight: 3.9990 chunk 34 optimal weight: 4.9990 chunk 40 optimal weight: 1.9990 chunk 25 optimal weight: 0.9980 chunk 24 optimal weight: 20.0000 chunk 18 optimal weight: 9.9990 chunk 16 optimal weight: 3.9990 chunk 12 optimal weight: 4.9990 chunk 7 optimal weight: 6.9990 chunk 27 optimal weight: 9.9990 chunk 20 optimal weight: 9.9990 overall best weight: 3.1988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 357 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 360 ASN A 392 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8613 moved from start: 0.7645 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.031 3366 Z= 0.267 Angle : 0.504 4.881 4562 Z= 0.262 Chirality : 0.047 0.167 497 Planarity : 0.004 0.031 584 Dihedral : 4.277 24.569 458 Min Nonbonded Distance : 2.278 Molprobity Statistics. All-atom Clashscore : 4.39 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.84 % Favored : 96.16 % Rotamer Outliers : 0.29 % 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.41), residues: 417 helix: -0.62 (0.39), residues: 187 sheet: -0.78 (0.70), residues: 57 loop : -1.78 (0.45), residues: 173 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 37 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 36 time to evaluate : 0.623 Fit side-chains outliers start: 1 outliers final: 0 residues processed: 36 average time/residue: 0.1707 time to fit residues: 7.9009 Evaluate side-chains 33 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 33 time to evaluate : 0.417 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5697 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 42 random chunks: chunk 3 optimal weight: 3.9990 chunk 31 optimal weight: 10.0000 chunk 36 optimal weight: 4.9990 chunk 38 optimal weight: 6.9990 chunk 35 optimal weight: 0.0170 chunk 37 optimal weight: 4.9990 chunk 22 optimal weight: 0.0470 chunk 16 optimal weight: 6.9990 chunk 29 optimal weight: 6.9990 chunk 11 optimal weight: 6.9990 chunk 34 optimal weight: 1.9990 overall best weight: 2.2122 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 357 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8587 moved from start: 0.7725 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.023 3366 Z= 0.194 Angle : 0.474 5.503 4562 Z= 0.244 Chirality : 0.045 0.146 497 Planarity : 0.003 0.031 584 Dihedral : 4.113 23.322 458 Min Nonbonded Distance : 2.263 Molprobity Statistics. All-atom Clashscore : 4.39 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.88 % Favored : 97.12 % Rotamer Outliers : 0.88 % 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.52 (0.41), residues: 417 helix: -0.43 (0.39), residues: 187 sheet: -0.65 (0.71), residues: 57 loop : -1.69 (0.46), residues: 173 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 36 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 33 time to evaluate : 0.434 Fit side-chains outliers start: 3 outliers final: 1 residues processed: 35 average time/residue: 0.1740 time to fit residues: 7.7283 Evaluate side-chains 34 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 33 time to evaluate : 0.387 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0418 time to fit residues: 0.6579 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 42 random chunks: chunk 35 optimal weight: 0.7980 chunk 37 optimal weight: 5.9990 chunk 24 optimal weight: 7.9990 chunk 39 optimal weight: 8.9990 chunk 18 optimal weight: 0.9980 chunk 27 optimal weight: 0.9980 chunk 41 optimal weight: 5.9990 chunk 38 optimal weight: 0.8980 chunk 33 optimal weight: 0.8980 chunk 3 optimal weight: 9.9990 chunk 25 optimal weight: 0.9990 overall best weight: 0.9180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 357 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8532 moved from start: 0.7869 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.016 3366 Z= 0.119 Angle : 0.426 5.156 4562 Z= 0.219 Chirality : 0.044 0.138 497 Planarity : 0.003 0.030 584 Dihedral : 3.641 15.938 458 Min Nonbonded Distance : 2.285 Molprobity Statistics. All-atom Clashscore : 4.09 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.12 % Favored : 96.88 % Rotamer Outliers : 0.59 % 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.76 (0.41), residues: 417 helix: 0.34 (0.41), residues: 176 sheet: -0.15 (0.73), residues: 56 loop : -1.42 (0.44), residues: 185 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 37 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 35 time to evaluate : 0.450 Fit side-chains outliers start: 2 outliers final: 0 residues processed: 36 average time/residue: 0.1700 time to fit residues: 7.8131 Evaluate side-chains 33 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 33 time to evaluate : 0.417 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5705 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 42 random chunks: chunk 20 optimal weight: 9.9990 chunk 26 optimal weight: 10.0000 chunk 35 optimal weight: 0.9990 chunk 10 optimal weight: 6.9990 chunk 30 optimal weight: 6.9990 chunk 4 optimal weight: 10.0000 chunk 9 optimal weight: 4.9990 chunk 33 optimal weight: 5.9990 chunk 13 optimal weight: 5.9990 chunk 34 optimal weight: 0.7980 chunk 6 optimal weight: 1.9990 overall best weight: 2.9588 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 19 GLN ** A 357 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8601 moved from start: 0.7936 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.027 3366 Z= 0.246 Angle : 0.492 5.551 4562 Z= 0.252 Chirality : 0.046 0.164 497 Planarity : 0.003 0.032 584 Dihedral : 4.119 24.168 458 Min Nonbonded Distance : 2.250 Molprobity Statistics. All-atom Clashscore : 4.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.08 % Favored : 95.92 % Rotamer Outliers : 0.59 % 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.95 (0.41), residues: 417 helix: 0.09 (0.40), residues: 182 sheet: -0.08 (0.75), residues: 56 loop : -1.53 (0.45), residues: 179 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 834 Ramachandran restraints generated. 417 Oldfield, 0 Emsley, 417 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 36 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 34 time to evaluate : 0.428 Fit side-chains outliers start: 2 outliers final: 2 residues processed: 35 average time/residue: 0.1618 time to fit residues: 7.3133 Evaluate side-chains 35 residues out of total 342 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 33 time to evaluate : 0.407 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0343 time to fit residues: 0.6226 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 42 random chunks: chunk 29 optimal weight: 3.9990 chunk 1 optimal weight: 0.9980 chunk 24 optimal weight: 10.0000 chunk 38 optimal weight: 6.9990 chunk 22 optimal weight: 5.9990 chunk 28 optimal weight: 6.9990 chunk 26 optimal weight: 5.9990 chunk 25 optimal weight: 0.6980 chunk 16 optimal weight: 6.9990 chunk 15 optimal weight: 3.9990 chunk 41 optimal weight: 8.9990 overall best weight: 3.1386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 357 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3327 r_free = 0.3327 target = 0.116894 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 51)----------------| | r_work = 0.3015 r_free = 0.3015 target = 0.094450 restraints weight = 5177.992| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 27)----------------| | r_work = 0.2990 r_free = 0.2990 target = 0.092916 restraints weight = 6899.183| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 28)----------------| | r_work = 0.3001 r_free = 0.3001 target = 0.093746 restraints weight = 7437.488| |-----------------------------------------------------------------------------| r_work (final): 0.2938 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8569 moved from start: 0.8020 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.032 3366 Z= 0.256 Angle : 0.510 5.784 4562 Z= 0.264 Chirality : 0.047 0.151 497 Planarity : 0.004 0.032 584 Dihedral : 4.331 27.030 458 Min Nonbonded Distance : 2.240 Molprobity Statistics. All-atom Clashscore : 5.61 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.84 % Favored : 96.16 % Rotamer Outliers : 0.29 % 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.17 (0.42), residues: 417 helix: -0.12 (0.39), residues: 188 sheet: -0.14 (0.76), residues: 57 loop : -1.67 (0.46), residues: 172 =============================================================================== Job complete usr+sys time: 986.72 seconds wall clock time: 18 minutes 36.19 seconds (1116.19 seconds total)