Starting phenix.real_space_refine on Tue Feb 13 19:59:21 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gkh_40184/02_2024/8gkh_40184_neut.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gkh_40184/02_2024/8gkh_40184.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.7 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gkh_40184/02_2024/8gkh_40184.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gkh_40184/02_2024/8gkh_40184.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gkh_40184/02_2024/8gkh_40184_neut.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8gkh_40184/02_2024/8gkh_40184_neut.pdb" } resolution = 2.7 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.001 sd= 0.026 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: 7 Type Number sf(0) Gaussians Zn 1 6.06 5 P 125 5.49 5 Mg 2 5.21 5 S 32 5.16 5 C 4337 2.51 5 N 1358 2.21 5 O 1715 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped None Time to flip residues: 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5238/modules/chem_data/mon_lib" Total number of atoms: 7570 Number of models: 1 Model: "" Number of chains: 7 Chain: "N" Number of atoms: 308 Number of conformers: 1 Conformer: "" Number of residues, atoms: 15, 308 Classifications: {'DNA': 15} Link IDs: {'rna3p': 14} Chain: "P" Number of atoms: 4946 Number of conformers: 1 Conformer: "" Number of residues, atoms: 600, 4946 Classifications: {'peptide': 600} Link IDs: {'PTRANS': 25, 'TRANS': 574} Chain breaks: 1 Chain: "T" Number of atoms: 695 Number of conformers: 1 Conformer: "" Number of residues, atoms: 34, 695 Classifications: {'DNA': 34} Link IDs: {'rna3p': 33} Chain breaks: 1 Chain: "W" Number of atoms: 1615 Number of conformers: 1 Conformer: "" Number of residues, atoms: 76, 1615 Classifications: {'RNA': 76} Modifications used: {'rna2p_pur': 6, 'rna2p_pyr': 5, 'rna3p_pur': 30, 'rna3p_pyr': 35} Link IDs: {'rna2p': 11, 'rna3p': 64} Chain breaks: 1 Chain: "P" Number of atoms: 4 Number of conformers: 1 Conformer: "" Number of residues, atoms: 4, 4 Unusual residues: {' MG': 2, ' ZN': 1} Classifications: {'undetermined': 3, 'water': 1} Link IDs: {None: 3} Chain: "T" Number of atoms: 1 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 1 Classifications: {'water': 1} Chain: "W" Number of atoms: 1 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 1 Classifications: {'water': 1} List of CYS excluded from plausible disulfide bonds: (reason: may participate in coordination) ATOM 4558 SG CYS P 548 32.953 45.640 21.668 1.00 51.02 S ATOM 4582 SG CYS P 551 30.062 44.068 19.775 1.00 57.06 S ATOM 4875 SG CYS P 586 33.492 42.516 19.547 1.00 52.44 S ATOM 4911 SG CYS P 591 32.800 45.814 17.814 1.00 44.29 S Time building chain proxies: 4.92, per 1000 atoms: 0.65 Number of scatterers: 7570 At special positions: 0 Unit cell: (75.9, 107.25, 139.425, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 7 Type Number sf(0) Zn 1 29.99 S 32 16.00 P 125 15.00 Mg 2 11.99 O 1715 8.00 N 1358 7.00 C 4337 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: 3.86 Conformation dependent library (CDL) restraints added in 1.1 seconds Dynamic metal coordination Zn2+ tetrahedral coordination pdb=" ZN P 703 " pdb="ZN ZN P 703 " - pdb=" SG CYS P 551 " pdb="ZN ZN P 703 " - pdb=" SG CYS P 586 " pdb="ZN ZN P 703 " - pdb=" SG CYS P 591 " pdb="ZN ZN P 703 " - pdb=" SG CYS P 548 " Number of angles added : 6 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1156 Finding SS restraints... Secondary structure from input PDB file: 30 helices and 4 sheets defined 61.0% alpha, 7.7% beta 52 base pairs and 76 stacking pairs defined. Time for finding SS restraints: 3.29 Creating SS restraints... Processing helix chain 'P' and resid 28 through 32 Processing helix chain 'P' and resid 35 through 68 Processing helix chain 'P' and resid 76 through 90 Processing helix chain 'P' and resid 91 through 93 No H-bonds generated for 'chain 'P' and resid 91 through 93' Processing helix chain 'P' and resid 97 through 119 Proline residue: P 109 - end of helix removed outlier: 4.072A pdb=" N GLN P 112 " --> pdb=" O LEU P 108 " (cutoff:3.500A) removed outlier: 4.205A pdb=" N ARG P 113 " --> pdb=" O PRO P 109 " (cutoff:3.500A) Processing helix chain 'P' and resid 128 through 151 Processing helix chain 'P' and resid 151 through 164 Processing helix chain 'P' and resid 164 through 172 Processing helix chain 'P' and resid 176 through 194 removed outlier: 3.630A pdb=" N PHE P 182 " --> pdb=" O ALA P 178 " (cutoff:3.500A) removed outlier: 3.815A pdb=" N SER P 191 " --> pdb=" O ARG P 187 " (cutoff:3.500A) removed outlier: 3.567A pdb=" N VAL P 192 " --> pdb=" O TYR P 188 " (cutoff:3.500A) Processing helix chain 'P' and resid 200 through 207 Processing helix chain 'P' and resid 208 through 220 Processing helix chain 'P' and resid 230 through 238 Processing helix chain 'P' and resid 239 through 241 No H-bonds generated for 'chain 'P' and resid 239 through 241' Processing helix chain 'P' and resid 242 through 257 Processing helix chain 'P' and resid 279 through 286 Processing helix chain 'P' and resid 290 through 296 removed outlier: 3.847A pdb=" N THR P 296 " --> pdb=" O ARG P 292 " (cutoff:3.500A) Processing helix chain 'P' and resid 298 through 305 removed outlier: 3.601A pdb=" N ARG P 305 " --> pdb=" O ASP P 301 " (cutoff:3.500A) Processing helix chain 'P' and resid 311 through 315 Processing helix chain 'P' and resid 365 through 368 Processing helix chain 'P' and resid 369 through 376 Processing helix chain 'P' and resid 403 through 412 Processing helix chain 'P' and resid 412 through 428 Processing helix chain 'P' and resid 428 through 435 removed outlier: 3.730A pdb=" N ARG P 435 " --> pdb=" O LEU P 431 " (cutoff:3.500A) Processing helix chain 'P' and resid 443 through 457 removed outlier: 3.741A pdb=" N ASP P 457 " --> pdb=" O VAL P 453 " (cutoff:3.500A) Processing helix chain 'P' and resid 457 through 466 Processing helix chain 'P' and resid 468 through 497 Processing helix chain 'P' and resid 523 through 532 removed outlier: 4.036A pdb=" N THR P 527 " --> pdb=" O LYS P 523 " (cutoff:3.500A) Processing helix chain 'P' and resid 568 through 575 Processing helix chain 'P' and resid 588 through 592 Processing helix chain 'P' and resid 605 through 622 Processing sheet with id=AA1, first strand: chain 'P' and resid 21 through 27 removed outlier: 4.929A pdb=" N SER P 325 " --> pdb=" O TYR P 336 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'P' and resid 400 through 402 Processing sheet with id=AA3, first strand: chain 'P' and resid 557 through 559 removed outlier: 3.642A pdb=" N GLY P 585 " --> pdb=" O GLU P 558 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'P' and resid 563 through 565 271 hydrogen bonds defined for protein. 786 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 132 hydrogen bonds 248 hydrogen bond angles 0 basepair planarities 52 basepair parallelities 76 stacking parallelities Total time for adding SS restraints: 3.04 Time building geometry restraints manager: 4.09 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: 1900 1.34 - 1.45: 2153 1.45 - 1.57: 3645 1.57 - 1.69: 245 1.69 - 1.81: 48 Bond restraints: 7991 Sorted by residual: bond pdb=" O3' DG T -17 " pdb=" P DA T -16 " ideal model delta sigma weight residual 1.607 1.625 -0.018 1.50e-02 4.44e+03 1.37e+00 bond pdb=" C1' A W 16 " pdb=" N9 A W 16 " ideal model delta sigma weight residual 1.465 1.480 -0.015 1.50e-02 4.44e+03 1.07e+00 bond pdb=" C1' U W 15 " pdb=" N1 U W 15 " ideal model delta sigma weight residual 1.470 1.484 -0.014 1.50e-02 4.44e+03 8.77e-01 bond pdb=" C3' DC T -6 " pdb=" O3' DC T -6 " ideal model delta sigma weight residual 1.422 1.450 -0.028 3.00e-02 1.11e+03 8.52e-01 bond pdb=" C3' DA T -8 " pdb=" O3' DA T -8 " ideal model delta sigma weight residual 1.422 1.449 -0.027 3.00e-02 1.11e+03 7.82e-01 ... (remaining 7986 not shown) Histogram of bond angle deviations from ideal: 99.16 - 106.13: 791 106.13 - 113.10: 4332 113.10 - 120.06: 2908 120.06 - 127.03: 3008 127.03 - 134.00: 331 Bond angle restraints: 11370 Sorted by residual: angle pdb=" O3' DG T -17 " pdb=" P DA T -16 " pdb=" OP1 DA T -16 " ideal model delta sigma weight residual 108.00 131.22 -23.22 3.00e+00 1.11e-01 5.99e+01 angle pdb=" OP1 DA T -16 " pdb=" P DA T -16 " pdb=" OP2 DA T -16 " ideal model delta sigma weight residual 120.00 103.13 16.87 3.00e+00 1.11e-01 3.16e+01 angle pdb=" O3' C W 19 " pdb=" C3' C W 19 " pdb=" C2' C W 19 " ideal model delta sigma weight residual 109.50 116.88 -7.38 1.50e+00 4.44e-01 2.42e+01 angle pdb=" C4' U W 76 " pdb=" C3' U W 76 " pdb=" O3' U W 76 " ideal model delta sigma weight residual 109.40 116.64 -7.24 1.50e+00 4.44e-01 2.33e+01 angle pdb=" C4' U W 15 " pdb=" C3' U W 15 " pdb=" O3' U W 15 " ideal model delta sigma weight residual 109.40 115.86 -6.46 1.50e+00 4.44e-01 1.85e+01 ... (remaining 11365 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 35.94: 4417 35.94 - 71.88: 334 71.88 - 107.82: 31 107.82 - 143.76: 8 143.76 - 179.70: 11 Dihedral angle restraints: 4801 sinusoidal: 3023 harmonic: 1778 Sorted by residual: dihedral pdb=" C4' U W 76 " pdb=" C3' U W 76 " pdb=" C2' U W 76 " pdb=" C1' U W 76 " ideal model delta sinusoidal sigma weight residual -35.00 36.12 -71.12 1 8.00e+00 1.56e-02 1.01e+02 dihedral pdb=" C4' A W 16 " pdb=" C3' A W 16 " pdb=" C2' A W 16 " pdb=" C1' A W 16 " ideal model delta sinusoidal sigma weight residual -35.00 35.70 -70.70 1 8.00e+00 1.56e-02 1.00e+02 dihedral pdb=" C4' C W 19 " pdb=" C3' C W 19 " pdb=" C2' C W 19 " pdb=" C1' C W 19 " ideal model delta sinusoidal sigma weight residual -35.00 35.27 -70.27 1 8.00e+00 1.56e-02 9.94e+01 ... (remaining 4798 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.149: 1215 0.149 - 0.299: 39 0.299 - 0.448: 52 0.448 - 0.598: 0 0.598 - 0.747: 1 Chirality restraints: 1307 Sorted by residual: chirality pdb=" P DA T -16 " pdb=" OP1 DA T -16 " pdb=" OP2 DA T -16 " pdb=" O5' DA T -16 " both_signs ideal model delta sigma weight residual True 2.35 -3.09 -0.75 2.00e-01 2.50e+01 1.40e+01 chirality pdb=" P DC T 7 " pdb=" OP1 DC T 7 " pdb=" OP2 DC T 7 " pdb=" O5' DC T 7 " both_signs ideal model delta sigma weight residual True 2.35 -2.79 -0.44 2.00e-01 2.50e+01 4.91e+00 chirality pdb=" P DA T -4 " pdb=" OP1 DA T -4 " pdb=" OP2 DA T -4 " pdb=" O5' DA T -4 " both_signs ideal model delta sigma weight residual True 2.35 -2.76 -0.42 2.00e-01 2.50e+01 4.32e+00 ... (remaining 1304 not shown) Planarity restraints: 1000 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C1' G W 20 " 0.014 2.00e-02 2.50e+03 7.39e-03 1.64e+00 pdb=" N9 G W 20 " -0.020 2.00e-02 2.50e+03 pdb=" C8 G W 20 " 0.000 2.00e-02 2.50e+03 pdb=" N7 G W 20 " 0.000 2.00e-02 2.50e+03 pdb=" C5 G W 20 " 0.002 2.00e-02 2.50e+03 pdb=" C6 G W 20 " -0.003 2.00e-02 2.50e+03 pdb=" O6 G W 20 " 0.006 2.00e-02 2.50e+03 pdb=" N1 G W 20 " 0.000 2.00e-02 2.50e+03 pdb=" C2 G W 20 " 0.001 2.00e-02 2.50e+03 pdb=" N2 G W 20 " -0.003 2.00e-02 2.50e+03 pdb=" N3 G W 20 " 0.001 2.00e-02 2.50e+03 pdb=" C4 G W 20 " 0.002 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C1' G W 70 " -0.016 2.00e-02 2.50e+03 7.29e-03 1.59e+00 pdb=" N9 G W 70 " 0.018 2.00e-02 2.50e+03 pdb=" C8 G W 70 " 0.001 2.00e-02 2.50e+03 pdb=" N7 G W 70 " -0.001 2.00e-02 2.50e+03 pdb=" C5 G W 70 " -0.001 2.00e-02 2.50e+03 pdb=" C6 G W 70 " 0.001 2.00e-02 2.50e+03 pdb=" O6 G W 70 " -0.006 2.00e-02 2.50e+03 pdb=" N1 G W 70 " -0.001 2.00e-02 2.50e+03 pdb=" C2 G W 70 " 0.001 2.00e-02 2.50e+03 pdb=" N2 G W 70 " 0.001 2.00e-02 2.50e+03 pdb=" N3 G W 70 " 0.001 2.00e-02 2.50e+03 pdb=" C4 G W 70 " 0.002 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C1' U W 27 " -0.015 2.00e-02 2.50e+03 7.91e-03 1.41e+00 pdb=" N1 U W 27 " 0.017 2.00e-02 2.50e+03 pdb=" C2 U W 27 " 0.001 2.00e-02 2.50e+03 pdb=" O2 U W 27 " 0.002 2.00e-02 2.50e+03 pdb=" N3 U W 27 " -0.001 2.00e-02 2.50e+03 pdb=" C4 U W 27 " -0.002 2.00e-02 2.50e+03 pdb=" O4 U W 27 " -0.004 2.00e-02 2.50e+03 pdb=" C5 U W 27 " 0.000 2.00e-02 2.50e+03 pdb=" C6 U W 27 " 0.003 2.00e-02 2.50e+03 ... (remaining 997 not shown) Histogram of nonbonded interaction distances: 1.88 - 2.48: 63 2.48 - 3.09: 5215 3.09 - 3.69: 12649 3.69 - 4.30: 18239 4.30 - 4.90: 28379 Nonbonded interactions: 64545 Sorted by model distance: nonbonded pdb="MG MG P 702 " pdb=" O HOH P 801 " model vdw 1.881 2.170 nonbonded pdb=" OD1 ASN P 385 " pdb="MG MG P 701 " model vdw 1.886 2.170 nonbonded pdb=" OP1 DT T -19 " pdb="MG MG P 701 " model vdw 1.900 2.170 nonbonded pdb=" OP2 DT T -19 " pdb="MG MG P 702 " model vdw 1.902 2.170 nonbonded pdb=" O2' G W 20 " pdb=" OP1 C W 21 " model vdw 1.902 2.440 ... (remaining 64540 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.600 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.000 Extract box with map and model: 5.770 Check model and map are aligned: 0.120 Set scattering table: 0.070 Process input model: 29.170 Find NCS groups from input model: 0.170 Set up NCS constraints: 0.030 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:11.820 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 48.760 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6744 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.032 7991 Z= 0.130 Angle : 1.275 23.219 11370 Z= 0.504 Chirality : 0.096 0.747 1307 Planarity : 0.002 0.025 1000 Dihedral : 22.525 179.701 3645 Min Nonbonded Distance : 1.881 Molprobity Statistics. All-atom Clashscore : 5.49 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.17 % Favored : 98.83 % Rotamer: Outliers : 0.00 % Allowed : 0.74 % Favored : 99.26 % 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.05 (0.36), residues: 596 helix: 2.95 (0.29), residues: 331 sheet: 0.24 (0.83), residues: 41 loop : 1.06 (0.44), residues: 224 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP P 93 HIS 0.002 0.001 HIS P 55 PHE 0.005 0.001 PHE P 151 TYR 0.016 0.001 TYR P 283 ARG 0.002 0.000 ARG P 461 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 128 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 128 time to evaluate : 0.666 Fit side-chains revert: symmetry clash revert: symmetry clash REVERT: P 238 ASP cc_start: 0.7088 (t0) cc_final: 0.6786 (t0) REVERT: P 487 ASP cc_start: 0.7283 (m-30) cc_final: 0.7050 (m-30) REVERT: P 574 GLU cc_start: 0.7050 (tp30) cc_final: 0.6738 (tp30) outliers start: 0 outliers final: 1 residues processed: 128 average time/residue: 1.3775 time to fit residues: 184.8301 Evaluate side-chains 121 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 120 time to evaluate : 0.594 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain P residue 577 ILE Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 71 random chunks: chunk 59 optimal weight: 1.9990 chunk 53 optimal weight: 0.9980 chunk 29 optimal weight: 1.9990 chunk 18 optimal weight: 1.9990 chunk 36 optimal weight: 0.0670 chunk 28 optimal weight: 1.9990 chunk 55 optimal weight: 2.9990 chunk 21 optimal weight: 1.9990 chunk 33 optimal weight: 0.0770 chunk 41 optimal weight: 2.9990 chunk 64 optimal weight: 2.9990 overall best weight: 1.0280 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: P 35 ASN P 308 ASN P 378 ASN P 393 GLN P 396 ASN P 482 GLN P 517 GLN P 581 HIS Total number of N/Q/H flips: 8 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7047 moved from start: 0.2346 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.040 7991 Z= 0.316 Angle : 0.689 8.022 11370 Z= 0.372 Chirality : 0.043 0.255 1307 Planarity : 0.005 0.039 1000 Dihedral : 25.520 176.573 2411 Min Nonbonded Distance : 1.838 Molprobity Statistics. All-atom Clashscore : 9.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.34 % Favored : 98.66 % Rotamer: Outliers : 2.58 % Allowed : 9.76 % Favored : 87.66 % 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.66 (0.34), residues: 596 helix: 1.64 (0.28), residues: 340 sheet: -0.28 (0.76), residues: 46 loop : 0.74 (0.43), residues: 210 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.002 TRP P 570 HIS 0.008 0.002 HIS P 566 PHE 0.022 0.002 PHE P 62 TYR 0.019 0.002 TYR P 283 ARG 0.010 0.001 ARG P 564 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 137 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 123 time to evaluate : 0.624 Fit side-chains REVERT: P 62 PHE cc_start: 0.5649 (t80) cc_final: 0.3891 (t80) REVERT: P 238 ASP cc_start: 0.7254 (t0) cc_final: 0.7018 (t0) REVERT: P 362 LYS cc_start: 0.7340 (OUTLIER) cc_final: 0.7099 (mtmt) REVERT: P 388 ASP cc_start: 0.7833 (m-30) cc_final: 0.7594 (m-30) REVERT: P 445 MET cc_start: 0.4329 (OUTLIER) cc_final: 0.3869 (tmm) REVERT: P 483 LYS cc_start: 0.8086 (tppt) cc_final: 0.7389 (tppt) REVERT: P 487 ASP cc_start: 0.7347 (m-30) cc_final: 0.7090 (m-30) REVERT: P 602 TYR cc_start: 0.8234 (m-80) cc_final: 0.7921 (m-80) REVERT: P 606 ASP cc_start: 0.7747 (m-30) cc_final: 0.7537 (m-30) outliers start: 14 outliers final: 4 residues processed: 124 average time/residue: 1.2703 time to fit residues: 165.7554 Evaluate side-chains 124 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 118 time to evaluate : 0.636 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain P residue 327 THR Chi-restraints excluded: chain P residue 362 LYS Chi-restraints excluded: chain P residue 445 MET Chi-restraints excluded: chain P residue 545 SER Chi-restraints excluded: chain P residue 577 ILE Chi-restraints excluded: chain P residue 597 THR Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 71 random chunks: chunk 35 optimal weight: 0.9990 chunk 20 optimal weight: 0.8980 chunk 53 optimal weight: 0.7980 chunk 43 optimal weight: 1.9990 chunk 17 optimal weight: 2.9990 chunk 64 optimal weight: 1.9990 chunk 69 optimal weight: 3.9990 chunk 57 optimal weight: 0.9980 chunk 22 optimal weight: 0.9980 chunk 51 optimal weight: 0.2980 chunk 63 optimal weight: 0.9990 overall best weight: 0.7980 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: P 35 ASN P 378 ASN P 396 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7053 moved from start: 0.2792 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.032 7991 Z= 0.261 Angle : 0.623 8.193 11370 Z= 0.338 Chirality : 0.040 0.234 1307 Planarity : 0.005 0.042 1000 Dihedral : 25.165 179.006 2411 Min Nonbonded Distance : 1.853 Molprobity Statistics. All-atom Clashscore : 10.61 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.67 % Favored : 99.33 % Rotamer: Outliers : 3.13 % Allowed : 12.34 % Favored : 84.53 % 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.29 (0.34), residues: 596 helix: 1.33 (0.27), residues: 339 sheet: -0.39 (0.75), residues: 46 loop : 0.57 (0.43), residues: 211 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP P 570 HIS 0.006 0.001 HIS P 566 PHE 0.018 0.002 PHE P 630 TYR 0.018 0.002 TYR P 283 ARG 0.008 0.001 ARG P 374 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 137 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 17 poor density : 120 time to evaluate : 0.685 Fit side-chains revert: symmetry clash REVERT: P 62 PHE cc_start: 0.5799 (t80) cc_final: 0.3982 (t80) REVERT: P 238 ASP cc_start: 0.7357 (t0) cc_final: 0.7138 (t0) REVERT: P 442 MET cc_start: 0.4022 (OUTLIER) cc_final: 0.0334 (mpt) REVERT: P 483 LYS cc_start: 0.8066 (tppt) cc_final: 0.7436 (tppt) REVERT: P 487 ASP cc_start: 0.7342 (m-30) cc_final: 0.7075 (m-30) REVERT: P 602 TYR cc_start: 0.8239 (m-80) cc_final: 0.8004 (m-80) outliers start: 17 outliers final: 5 residues processed: 122 average time/residue: 1.3261 time to fit residues: 170.3974 Evaluate side-chains 116 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 110 time to evaluate : 0.575 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain P residue 327 THR Chi-restraints excluded: chain P residue 395 SER Chi-restraints excluded: chain P residue 442 MET Chi-restraints excluded: chain P residue 545 SER Chi-restraints excluded: chain P residue 577 ILE Chi-restraints excluded: chain P residue 597 THR Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 71 random chunks: chunk 48 optimal weight: 0.6980 chunk 33 optimal weight: 1.9990 chunk 7 optimal weight: 0.7980 chunk 30 optimal weight: 0.2980 chunk 43 optimal weight: 0.9980 chunk 64 optimal weight: 1.9990 chunk 68 optimal weight: 9.9990 chunk 61 optimal weight: 1.9990 chunk 18 optimal weight: 0.7980 chunk 57 optimal weight: 0.4980 chunk 38 optimal weight: 1.9990 overall best weight: 0.6180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: P 35 ASN P 58 HIS P 129 GLN P 321 HIS P 378 ASN P 393 GLN P 396 ASN P 517 GLN ** P 590 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 8 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7043 moved from start: 0.3010 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.030 7991 Z= 0.221 Angle : 0.600 10.032 11370 Z= 0.323 Chirality : 0.038 0.219 1307 Planarity : 0.004 0.043 1000 Dihedral : 25.034 178.226 2411 Min Nonbonded Distance : 1.858 Molprobity Statistics. All-atom Clashscore : 9.60 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.17 % Favored : 98.83 % Rotamer: Outliers : 3.13 % Allowed : 12.34 % Favored : 84.53 % 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.32 (0.34), residues: 596 helix: 1.41 (0.27), residues: 337 sheet: -0.46 (0.74), residues: 46 loop : 0.51 (0.42), residues: 213 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.001 TRP P 570 HIS 0.005 0.001 HIS P 566 PHE 0.019 0.002 PHE P 630 TYR 0.018 0.002 TYR P 283 ARG 0.009 0.001 ARG P 374 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 128 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 17 poor density : 111 time to evaluate : 0.682 Fit side-chains REVERT: P 91 GLU cc_start: 0.6967 (mp0) cc_final: 0.6351 (mp0) REVERT: P 388 ASP cc_start: 0.7875 (m-30) cc_final: 0.7616 (m-30) REVERT: P 442 MET cc_start: 0.3926 (OUTLIER) cc_final: 0.0238 (mpt) REVERT: P 483 LYS cc_start: 0.8063 (tppt) cc_final: 0.7427 (tppt) REVERT: P 487 ASP cc_start: 0.7346 (m-30) cc_final: 0.7063 (m-30) REVERT: P 602 TYR cc_start: 0.8232 (m-80) cc_final: 0.7957 (m-80) outliers start: 17 outliers final: 9 residues processed: 117 average time/residue: 1.3130 time to fit residues: 161.7188 Evaluate side-chains 117 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 107 time to evaluate : 0.616 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain P residue 141 THR Chi-restraints excluded: chain P residue 323 GLU Chi-restraints excluded: chain P residue 327 THR Chi-restraints excluded: chain P residue 395 SER Chi-restraints excluded: chain P residue 442 MET Chi-restraints excluded: chain P residue 450 TYR Chi-restraints excluded: chain P residue 520 SER Chi-restraints excluded: chain P residue 545 SER Chi-restraints excluded: chain P residue 577 ILE Chi-restraints excluded: chain P residue 597 THR Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 71 random chunks: chunk 0 optimal weight: 2.9990 chunk 51 optimal weight: 0.5980 chunk 28 optimal weight: 1.9990 chunk 58 optimal weight: 1.9990 chunk 47 optimal weight: 0.5980 chunk 35 optimal weight: 0.6980 chunk 61 optimal weight: 1.9990 chunk 17 optimal weight: 0.6980 chunk 23 optimal weight: 0.9990 chunk 13 optimal weight: 0.9990 chunk 40 optimal weight: 0.9990 overall best weight: 0.7182 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: P 35 ASN P 124 ASN P 321 HIS P 378 ASN P 396 ASN ** P 590 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7074 moved from start: 0.3275 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.030 7991 Z= 0.247 Angle : 0.625 10.247 11370 Z= 0.332 Chirality : 0.039 0.242 1307 Planarity : 0.005 0.046 1000 Dihedral : 24.804 178.451 2411 Min Nonbonded Distance : 1.853 Molprobity Statistics. All-atom Clashscore : 9.67 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.17 % Favored : 98.83 % Rotamer: Outliers : 3.13 % Allowed : 12.71 % Favored : 84.16 % 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.08 (0.33), residues: 596 helix: 1.23 (0.27), residues: 339 sheet: -0.34 (0.76), residues: 46 loop : 0.32 (0.42), residues: 211 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP P 570 HIS 0.005 0.001 HIS P 566 PHE 0.022 0.002 PHE P 630 TYR 0.018 0.002 TYR P 283 ARG 0.010 0.001 ARG P 374 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 130 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 17 poor density : 113 time to evaluate : 0.678 Fit side-chains revert: symmetry clash REVERT: P 442 MET cc_start: 0.3985 (mpt) cc_final: 0.0290 (mpt) REVERT: P 483 LYS cc_start: 0.8085 (tppt) cc_final: 0.7437 (tppt) REVERT: P 487 ASP cc_start: 0.7358 (m-30) cc_final: 0.7066 (m-30) REVERT: P 574 GLU cc_start: 0.7189 (tp30) cc_final: 0.6982 (tp30) outliers start: 17 outliers final: 8 residues processed: 120 average time/residue: 1.3164 time to fit residues: 166.2134 Evaluate side-chains 115 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 107 time to evaluate : 0.632 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain P residue 327 THR Chi-restraints excluded: chain P residue 395 SER Chi-restraints excluded: chain P residue 450 TYR Chi-restraints excluded: chain P residue 520 SER Chi-restraints excluded: chain P residue 545 SER Chi-restraints excluded: chain P residue 577 ILE Chi-restraints excluded: chain P residue 597 THR Chi-restraints excluded: chain P residue 629 VAL Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 71 random chunks: chunk 16 optimal weight: 6.9990 chunk 68 optimal weight: 9.9990 chunk 57 optimal weight: 0.6980 chunk 31 optimal weight: 0.0870 chunk 5 optimal weight: 1.9990 chunk 22 optimal weight: 0.8980 chunk 36 optimal weight: 0.7980 chunk 66 optimal weight: 2.9990 chunk 7 optimal weight: 0.7980 chunk 39 optimal weight: 0.5980 chunk 50 optimal weight: 0.9980 overall best weight: 0.5958 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: P 35 ASN P 321 HIS P 378 ASN P 393 GLN P 396 ASN P 517 GLN ** P 590 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 6 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7053 moved from start: 0.3393 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.030 7991 Z= 0.218 Angle : 0.588 8.192 11370 Z= 0.317 Chirality : 0.038 0.210 1307 Planarity : 0.004 0.047 1000 Dihedral : 24.625 177.044 2411 Min Nonbonded Distance : 1.866 Molprobity Statistics. All-atom Clashscore : 9.75 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.34 % Favored : 98.66 % Rotamer: Outliers : 2.76 % Allowed : 13.63 % Favored : 83.61 % 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.19 (0.34), residues: 596 helix: 1.33 (0.28), residues: 337 sheet: -0.37 (0.70), residues: 56 loop : 0.42 (0.43), residues: 203 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP P 570 HIS 0.005 0.001 HIS P 566 PHE 0.022 0.002 PHE P 630 TYR 0.018 0.002 TYR P 283 ARG 0.011 0.001 ARG P 374 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 124 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 109 time to evaluate : 0.589 Fit side-chains REVERT: P 388 ASP cc_start: 0.7896 (m-30) cc_final: 0.7623 (m-30) REVERT: P 442 MET cc_start: 0.3946 (mpt) cc_final: 0.0162 (mpt) REVERT: P 483 LYS cc_start: 0.8113 (tppt) cc_final: 0.7458 (tppt) REVERT: P 487 ASP cc_start: 0.7362 (m-30) cc_final: 0.7073 (m-30) outliers start: 15 outliers final: 9 residues processed: 114 average time/residue: 1.2418 time to fit residues: 149.2289 Evaluate side-chains 118 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 109 time to evaluate : 0.579 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain P residue 141 THR Chi-restraints excluded: chain P residue 327 THR Chi-restraints excluded: chain P residue 395 SER Chi-restraints excluded: chain P residue 450 TYR Chi-restraints excluded: chain P residue 486 SER Chi-restraints excluded: chain P residue 520 SER Chi-restraints excluded: chain P residue 545 SER Chi-restraints excluded: chain P residue 559 LYS Chi-restraints excluded: chain P residue 597 THR Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 71 random chunks: chunk 38 optimal weight: 1.9990 chunk 57 optimal weight: 0.6980 chunk 68 optimal weight: 20.0000 chunk 42 optimal weight: 0.9990 chunk 41 optimal weight: 1.9990 chunk 31 optimal weight: 0.9980 chunk 27 optimal weight: 0.7980 chunk 40 optimal weight: 0.6980 chunk 20 optimal weight: 3.9990 chunk 13 optimal weight: 0.9990 chunk 43 optimal weight: 0.9980 overall best weight: 0.8380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: P 35 ASN P 55 HIS P 321 HIS P 378 ASN P 396 ASN ** P 590 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7095 moved from start: 0.3601 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.042 7991 Z= 0.277 Angle : 0.632 7.926 11370 Z= 0.339 Chirality : 0.040 0.242 1307 Planarity : 0.005 0.050 1000 Dihedral : 24.707 178.322 2409 Min Nonbonded Distance : 1.847 Molprobity Statistics. All-atom Clashscore : 11.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.51 % Favored : 98.49 % Rotamer: Outliers : 2.95 % Allowed : 14.36 % Favored : 82.69 % 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.83 (0.34), residues: 596 helix: 1.01 (0.27), residues: 338 sheet: -0.33 (0.71), residues: 56 loop : 0.25 (0.43), residues: 202 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP P 570 HIS 0.005 0.001 HIS P 566 PHE 0.022 0.002 PHE P 630 TYR 0.020 0.002 TYR P 450 ARG 0.010 0.001 ARG P 374 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 129 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 16 poor density : 113 time to evaluate : 0.712 Fit side-chains revert: symmetry clash REVERT: P 388 ASP cc_start: 0.7899 (m-30) cc_final: 0.7624 (m-30) REVERT: P 396 ASN cc_start: 0.7788 (m-40) cc_final: 0.7556 (m110) REVERT: P 442 MET cc_start: 0.4008 (OUTLIER) cc_final: 0.0316 (mpt) REVERT: P 483 LYS cc_start: 0.8110 (tppt) cc_final: 0.7450 (tppt) REVERT: P 487 ASP cc_start: 0.7377 (m-30) cc_final: 0.7092 (m-30) outliers start: 16 outliers final: 11 residues processed: 120 average time/residue: 1.2855 time to fit residues: 162.4776 Evaluate side-chains 119 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 107 time to evaluate : 0.640 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain P residue 141 THR Chi-restraints excluded: chain P residue 327 THR Chi-restraints excluded: chain P residue 395 SER Chi-restraints excluded: chain P residue 419 ARG Chi-restraints excluded: chain P residue 442 MET Chi-restraints excluded: chain P residue 450 TYR Chi-restraints excluded: chain P residue 486 SER Chi-restraints excluded: chain P residue 520 SER Chi-restraints excluded: chain P residue 545 SER Chi-restraints excluded: chain P residue 577 ILE Chi-restraints excluded: chain P residue 597 THR Chi-restraints excluded: chain P residue 629 VAL Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 71 random chunks: chunk 46 optimal weight: 1.9990 chunk 33 optimal weight: 0.8980 chunk 6 optimal weight: 0.0870 chunk 53 optimal weight: 0.7980 chunk 62 optimal weight: 1.9990 chunk 65 optimal weight: 1.9990 chunk 59 optimal weight: 0.7980 chunk 63 optimal weight: 0.9990 chunk 38 optimal weight: 1.9990 chunk 27 optimal weight: 0.8980 chunk 50 optimal weight: 0.6980 overall best weight: 0.6558 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: P 35 ASN P 321 HIS P 378 ASN P 393 GLN P 517 GLN ** P 590 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7071 moved from start: 0.3628 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.040 7991 Z= 0.235 Angle : 0.610 8.414 11370 Z= 0.326 Chirality : 0.038 0.214 1307 Planarity : 0.005 0.053 1000 Dihedral : 24.655 177.293 2409 Min Nonbonded Distance : 1.865 Molprobity Statistics. All-atom Clashscore : 10.97 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.34 % Favored : 98.66 % Rotamer: Outliers : 2.21 % Allowed : 15.84 % Favored : 81.95 % 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.83 (0.33), residues: 596 helix: 1.09 (0.27), residues: 338 sheet: -0.46 (0.76), residues: 46 loop : 0.10 (0.42), residues: 212 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP P 570 HIS 0.005 0.001 HIS P 566 PHE 0.022 0.002 PHE P 630 TYR 0.018 0.002 TYR P 283 ARG 0.011 0.001 ARG P 374 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 127 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 115 time to evaluate : 0.571 Fit side-chains revert: symmetry clash REVERT: P 388 ASP cc_start: 0.7892 (m-30) cc_final: 0.7624 (m-30) REVERT: P 442 MET cc_start: 0.3948 (OUTLIER) cc_final: 0.0211 (mpt) REVERT: P 483 LYS cc_start: 0.8115 (tppt) cc_final: 0.7467 (tppt) REVERT: P 487 ASP cc_start: 0.7378 (m-30) cc_final: 0.7087 (m-30) REVERT: P 602 TYR cc_start: 0.8275 (m-80) cc_final: 0.8023 (m-80) outliers start: 12 outliers final: 9 residues processed: 120 average time/residue: 1.2910 time to fit residues: 163.0298 Evaluate side-chains 121 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 111 time to evaluate : 0.658 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain P residue 141 THR Chi-restraints excluded: chain P residue 327 THR Chi-restraints excluded: chain P residue 419 ARG Chi-restraints excluded: chain P residue 442 MET Chi-restraints excluded: chain P residue 450 TYR Chi-restraints excluded: chain P residue 486 SER Chi-restraints excluded: chain P residue 545 SER Chi-restraints excluded: chain P residue 559 LYS Chi-restraints excluded: chain P residue 577 ILE Chi-restraints excluded: chain P residue 597 THR Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 71 random chunks: chunk 19 optimal weight: 2.9990 chunk 57 optimal weight: 0.5980 chunk 60 optimal weight: 0.0060 chunk 63 optimal weight: 1.9990 chunk 41 optimal weight: 0.7980 chunk 67 optimal weight: 2.9990 chunk 31 optimal weight: 0.6980 chunk 46 optimal weight: 0.6980 chunk 70 optimal weight: 1.9990 chunk 65 optimal weight: 1.9990 chunk 56 optimal weight: 0.8980 overall best weight: 0.5596 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: P 35 ASN P 146 ASN P 321 HIS P 378 ASN P 393 GLN P 517 GLN ** P 590 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 6 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7052 moved from start: 0.3658 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.035 7991 Z= 0.214 Angle : 0.598 8.057 11370 Z= 0.320 Chirality : 0.038 0.202 1307 Planarity : 0.005 0.052 1000 Dihedral : 24.640 176.462 2409 Min Nonbonded Distance : 1.872 Molprobity Statistics. All-atom Clashscore : 11.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.51 % Favored : 98.49 % Rotamer: Outliers : 1.66 % Allowed : 16.94 % Favored : 81.40 % 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.90 (0.34), residues: 596 helix: 1.18 (0.27), residues: 338 sheet: -0.56 (0.77), residues: 46 loop : 0.07 (0.42), residues: 212 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.001 TRP P 570 HIS 0.005 0.001 HIS P 566 PHE 0.017 0.001 PHE P 630 TYR 0.018 0.002 TYR P 134 ARG 0.011 0.000 ARG P 374 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 125 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 116 time to evaluate : 0.622 Fit side-chains revert: symmetry clash REVERT: P 388 ASP cc_start: 0.7890 (m-30) cc_final: 0.7618 (m-30) REVERT: P 442 MET cc_start: 0.3858 (OUTLIER) cc_final: 0.0070 (mpt) REVERT: P 483 LYS cc_start: 0.8108 (tppt) cc_final: 0.7469 (tppt) REVERT: P 487 ASP cc_start: 0.7377 (m-30) cc_final: 0.7090 (m-30) REVERT: P 574 GLU cc_start: 0.7131 (tp30) cc_final: 0.6875 (tp30) REVERT: P 602 TYR cc_start: 0.8266 (m-80) cc_final: 0.8000 (m-80) outliers start: 9 outliers final: 9 residues processed: 121 average time/residue: 1.2379 time to fit residues: 157.8537 Evaluate side-chains 119 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 109 time to evaluate : 0.614 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain P residue 141 THR Chi-restraints excluded: chain P residue 327 THR Chi-restraints excluded: chain P residue 419 ARG Chi-restraints excluded: chain P residue 442 MET Chi-restraints excluded: chain P residue 450 TYR Chi-restraints excluded: chain P residue 486 SER Chi-restraints excluded: chain P residue 545 SER Chi-restraints excluded: chain P residue 559 LYS Chi-restraints excluded: chain P residue 577 ILE Chi-restraints excluded: chain P residue 597 THR Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 71 random chunks: chunk 5 optimal weight: 0.7980 chunk 43 optimal weight: 3.9990 chunk 34 optimal weight: 0.5980 chunk 44 optimal weight: 1.9990 chunk 60 optimal weight: 0.0060 chunk 17 optimal weight: 1.9990 chunk 51 optimal weight: 0.9990 chunk 8 optimal weight: 0.0070 chunk 15 optimal weight: 2.9990 chunk 56 optimal weight: 1.9990 chunk 23 optimal weight: 0.6980 overall best weight: 0.4214 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: P 35 ASN P 321 HIS P 378 ASN P 393 GLN P 517 GLN ** P 590 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7023 moved from start: 0.3650 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.036 7991 Z= 0.193 Angle : 0.596 8.685 11370 Z= 0.318 Chirality : 0.037 0.198 1307 Planarity : 0.005 0.050 1000 Dihedral : 24.620 175.445 2409 Min Nonbonded Distance : 1.894 Molprobity Statistics. All-atom Clashscore : 11.70 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.85 % Favored : 98.15 % Rotamer: Outliers : 1.47 % Allowed : 17.50 % Favored : 81.03 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.02 (0.34), residues: 596 helix: 1.33 (0.28), residues: 337 sheet: -0.55 (0.77), residues: 46 loop : 0.07 (0.42), residues: 213 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP P 570 HIS 0.005 0.001 HIS P 321 PHE 0.020 0.001 PHE P 630 TYR 0.017 0.002 TYR P 283 ARG 0.011 0.000 ARG P 374 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1192 Ramachandran restraints generated. 596 Oldfield, 0 Emsley, 596 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 114 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 106 time to evaluate : 0.520 Fit side-chains revert: symmetry clash REVERT: P 62 PHE cc_start: 0.5525 (t80) cc_final: 0.5325 (t80) REVERT: P 238 ASP cc_start: 0.7267 (t0) cc_final: 0.7044 (t0) REVERT: P 442 MET cc_start: 0.3802 (OUTLIER) cc_final: -0.0036 (mpt) REVERT: P 483 LYS cc_start: 0.8110 (tppt) cc_final: 0.7486 (tppt) REVERT: P 487 ASP cc_start: 0.7370 (m-30) cc_final: 0.7083 (m-30) REVERT: P 602 TYR cc_start: 0.8251 (m-80) cc_final: 0.7992 (m-80) outliers start: 8 outliers final: 6 residues processed: 109 average time/residue: 1.2309 time to fit residues: 141.6305 Evaluate side-chains 117 residues out of total 543 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 110 time to evaluate : 0.676 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain P residue 141 THR Chi-restraints excluded: chain P residue 419 ARG Chi-restraints excluded: chain P residue 442 MET Chi-restraints excluded: chain P residue 450 TYR Chi-restraints excluded: chain P residue 545 SER Chi-restraints excluded: chain P residue 559 LYS Chi-restraints excluded: chain P residue 597 THR Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 71 random chunks: chunk 57 optimal weight: 1.9990 chunk 7 optimal weight: 0.6980 chunk 10 optimal weight: 0.9980 chunk 49 optimal weight: 0.5980 chunk 3 optimal weight: 0.4980 chunk 40 optimal weight: 0.6980 chunk 64 optimal weight: 2.9990 chunk 37 optimal weight: 0.2980 chunk 48 optimal weight: 0.8980 chunk 1 optimal weight: 2.9990 chunk 45 optimal weight: 0.9990 overall best weight: 0.5580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: P 35 ASN P 378 ASN P 396 ASN ** P 590 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4380 r_free = 0.4380 target = 0.187458 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 48)----------------| | r_work = 0.4133 r_free = 0.4133 target = 0.163940 restraints weight = 11168.588| |-----------------------------------------------------------------------------| r_work (start): 0.4123 rms_B_bonded: 1.26 r_work: 0.4033 rms_B_bonded: 2.18 restraints_weight: 0.5000 r_work: 0.3924 rms_B_bonded: 3.81 restraints_weight: 0.2500 r_work (final): 0.3924 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7118 moved from start: 0.3715 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.041 7991 Z= 0.219 Angle : 0.606 8.618 11370 Z= 0.324 Chirality : 0.038 0.199 1307 Planarity : 0.005 0.052 1000 Dihedral : 24.638 175.576 2409 Min Nonbonded Distance : 1.875 Molprobity Statistics. All-atom Clashscore : 11.55 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.51 % Favored : 98.49 % Rotamer: Outliers : 1.66 % Allowed : 16.94 % Favored : 81.40 % 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.34), residues: 596 helix: 1.26 (0.27), residues: 337 sheet: -0.48 (0.78), residues: 46 loop : 0.04 (0.42), residues: 213 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP P 570 HIS 0.005 0.001 HIS P 566 PHE 0.018 0.002 PHE P 630 TYR 0.018 0.002 TYR P 134 ARG 0.012 0.001 ARG P 374 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 3276.85 seconds wall clock time: 58 minutes 40.16 seconds (3520.16 seconds total)