Starting phenix.real_space_refine on Wed Mar 5 22:09:06 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/6ouc_20202/03_2025/6ouc_20202.cif Found real_map, /net/cci-nas-00/data/ceres_data/6ouc_20202/03_2025/6ouc_20202.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.8 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/6ouc_20202/03_2025/6ouc_20202.map" default_real_map = "/net/cci-nas-00/data/ceres_data/6ouc_20202/03_2025/6ouc_20202.map" model { file = "/net/cci-nas-00/data/ceres_data/6ouc_20202/03_2025/6ouc_20202.cif" } default_model = "/net/cci-nas-00/data/ceres_data/6ouc_20202/03_2025/6ouc_20202.cif" } resolution = 2.8 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.041 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: 5 Type Number sf(0) Gaussians Zn 1 6.06 5 S 15 5.16 5 C 2421 2.51 5 N 642 2.21 5 O 820 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 4 residue(s): 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5624/modules/chem_data/mon_lib" Total number of atoms: 3899 Number of models: 1 Model: "" Number of chains: 2 Chain: "A" Number of atoms: 3774 Number of conformers: 1 Conformer: "" Number of residues, atoms: 483, 3774 Classifications: {'peptide': 483} Link IDs: {'PTRANS': 46, 'TRANS': 436} Chain breaks: 1 Chain: "A" Number of atoms: 125 Number of conformers: 1 Conformer: "" Number of residues, atoms: 125, 125 Unusual residues: {' ZN': 1} Classifications: {'undetermined': 1, 'water': 124} Link IDs: {None: 124} Time building chain proxies: 2.80, per 1000 atoms: 0.72 Number of scatterers: 3899 At special positions: 0 Unit cell: (95.23, 95.23, 78.11, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) Zn 1 29.99 S 15 16.00 O 820 8.00 N 642 7.00 C 2421 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.00 Amino acid : True - 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: 0.87 Conformation dependent library (CDL) restraints added in 471.9 milliseconds Dynamic metal coordination Zn2+ tetrahedral coordination pdb=" ZN A 601 " pdb="ZN ZN A 601 " - pdb=" ND1 HIS A 299 " pdb="ZN ZN A 601 " - pdb=" NE2 HIS A 293 " 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 900 Finding SS restraints... Secondary structure from input PDB file: 4 helices and 9 sheets defined 6.2% alpha, 28.6% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.46 Creating SS restraints... Processing helix chain 'A' and resid 82 through 85 removed outlier: 3.522A pdb=" N LEU A 85 " --> pdb=" O GLY A 82 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 82 through 85' Processing helix chain 'A' and resid 86 through 93 Processing helix chain 'A' and resid 136 through 141 removed outlier: 3.628A pdb=" N THR A 139 " --> pdb=" O GLN A 136 " (cutoff:3.500A) removed outlier: 3.508A pdb=" N PHE A 141 " --> pdb=" O ILE A 138 " (cutoff:3.500A) Processing helix chain 'A' and resid 456 through 467 Processing sheet with id=AA1, first strand: chain 'A' and resid 57 through 58 Processing sheet with id=AA2, first strand: chain 'A' and resid 57 through 58 removed outlier: 4.053A pdb=" N SER A 199 " --> pdb=" O ALA A 110 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 75 through 80 removed outlier: 6.342A pdb=" N ALA A 180 " --> pdb=" O LEU A 76 " (cutoff:3.500A) removed outlier: 3.931A pdb=" N ASN A 78 " --> pdb=" O ILE A 178 " (cutoff:3.500A) removed outlier: 6.392A pdb=" N ILE A 178 " --> pdb=" O ASN A 78 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'A' and resid 167 through 168 Processing sheet with id=AA5, first strand: chain 'A' and resid 115 through 116 Processing sheet with id=AA6, first strand: chain 'A' and resid 450 through 453 removed outlier: 4.853A pdb=" N MET A 249 " --> pdb=" O VAL A 509 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'A' and resid 450 through 453 removed outlier: 4.743A pdb=" N MET A 499 " --> pdb=" O LEU A 434 " (cutoff:3.500A) removed outlier: 4.237A pdb=" N PHE A 489 " --> pdb=" O TYR A 479 " (cutoff:3.500A) removed outlier: 7.196A pdb=" N VAL A 474 " --> pdb=" O TRP A 522 " (cutoff:3.500A) removed outlier: 4.875A pdb=" N TRP A 522 " --> pdb=" O VAL A 474 " (cutoff:3.500A) removed outlier: 6.940A pdb=" N LEU A 476 " --> pdb=" O ASP A 520 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'A' and resid 265 through 266 removed outlier: 4.565A pdb=" N ARG A 265 " --> pdb=" O GLN A 273 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'A' and resid 352 through 360 removed outlier: 7.700A pdb=" N ARG A 331 " --> pdb=" O LEU A 396 " (cutoff:3.500A) removed outlier: 4.587A pdb=" N LEU A 396 " --> pdb=" O ARG A 331 " (cutoff:3.500A) removed outlier: 6.881A pdb=" N PHE A 333 " --> pdb=" O VAL A 394 " (cutoff:3.500A) removed outlier: 3.953A pdb=" N ARG A 339 " --> pdb=" O PRO A 388 " (cutoff:3.500A) removed outlier: 3.939A pdb=" N PHE A 286 " --> pdb=" O PHE A 391 " (cutoff:3.500A) 100 hydrogen bonds defined for protein. 258 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 0 stacking parallelities Total time for adding SS restraints: 1.07 Time building geometry restraints manager: 1.06 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.21 - 1.33: 969 1.33 - 1.45: 735 1.45 - 1.56: 2154 1.56 - 1.68: 0 1.68 - 1.80: 25 Bond restraints: 3883 Sorted by residual: bond pdb=" CA TYR A 214 " pdb=" C TYR A 214 " ideal model delta sigma weight residual 1.520 1.489 0.032 1.22e-02 6.72e+03 6.77e+00 bond pdb=" CA GLN A 170 " pdb=" C GLN A 170 " ideal model delta sigma weight residual 1.528 1.561 -0.033 1.33e-02 5.65e+03 6.30e+00 bond pdb=" N TRP A 378 " pdb=" CA TRP A 378 " ideal model delta sigma weight residual 1.459 1.430 0.029 1.24e-02 6.50e+03 5.52e+00 bond pdb=" CG LEU A 455 " pdb=" CD1 LEU A 455 " ideal model delta sigma weight residual 1.521 1.447 0.074 3.30e-02 9.18e+02 5.06e+00 bond pdb=" CA GLN A 170 " pdb=" CB GLN A 170 " ideal model delta sigma weight residual 1.527 1.498 0.029 1.44e-02 4.82e+03 3.96e+00 ... (remaining 3878 not shown) Histogram of bond angle deviations from ideal: 0.00 - 3.39: 5207 3.39 - 6.79: 93 6.79 - 10.18: 8 10.18 - 13.57: 6 13.57 - 16.97: 1 Bond angle restraints: 5315 Sorted by residual: angle pdb=" N GLY A 62 " pdb=" CA GLY A 62 " pdb=" C GLY A 62 " ideal model delta sigma weight residual 112.73 129.70 -16.97 1.20e+00 6.94e-01 2.00e+02 angle pdb=" N THR A 213 " pdb=" CA THR A 213 " pdb=" C THR A 213 " ideal model delta sigma weight residual 112.68 125.51 -12.83 1.33e+00 5.65e-01 9.31e+01 angle pdb=" N ASN A 297 " pdb=" CA ASN A 297 " pdb=" C ASN A 297 " ideal model delta sigma weight residual 111.28 101.34 9.94 1.09e+00 8.42e-01 8.32e+01 angle pdb=" N TYR A 301 " pdb=" CA TYR A 301 " pdb=" C TYR A 301 " ideal model delta sigma weight residual 111.02 121.74 -10.72 1.25e+00 6.40e-01 7.36e+01 angle pdb=" N LEU A 300 " pdb=" CA LEU A 300 " pdb=" C LEU A 300 " ideal model delta sigma weight residual 114.04 123.08 -9.04 1.24e+00 6.50e-01 5.31e+01 ... (remaining 5310 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 13.63: 2184 13.63 - 27.25: 101 27.25 - 40.87: 25 40.87 - 54.49: 3 54.49 - 68.11: 2 Dihedral angle restraints: 2315 sinusoidal: 917 harmonic: 1398 Sorted by residual: dihedral pdb=" C LEU A 300 " pdb=" N LEU A 300 " pdb=" CA LEU A 300 " pdb=" CB LEU A 300 " ideal model delta harmonic sigma weight residual -122.60 -138.63 16.03 0 2.50e+00 1.60e-01 4.11e+01 dihedral pdb=" CA ALA A 344 " pdb=" C ALA A 344 " pdb=" N ALA A 345 " pdb=" CA ALA A 345 " ideal model delta harmonic sigma weight residual -180.00 -150.63 -29.37 0 5.00e+00 4.00e-02 3.45e+01 dihedral pdb=" N LEU A 300 " pdb=" C LEU A 300 " pdb=" CA LEU A 300 " pdb=" CB LEU A 300 " ideal model delta harmonic sigma weight residual 122.80 137.25 -14.45 0 2.50e+00 1.60e-01 3.34e+01 ... (remaining 2312 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.161: 574 0.161 - 0.321: 11 0.321 - 0.482: 1 0.482 - 0.642: 1 0.642 - 0.802: 1 Chirality restraints: 588 Sorted by residual: chirality pdb=" CA LEU A 300 " pdb=" N LEU A 300 " pdb=" C LEU A 300 " pdb=" CB LEU A 300 " both_signs ideal model delta sigma weight residual False 2.51 1.71 0.80 2.00e-01 2.50e+01 1.61e+01 chirality pdb=" CA TYR A 301 " pdb=" N TYR A 301 " pdb=" C TYR A 301 " pdb=" CB TYR A 301 " both_signs ideal model delta sigma weight residual False 2.51 2.00 0.51 2.00e-01 2.50e+01 6.45e+00 chirality pdb=" CA THR A 213 " pdb=" N THR A 213 " pdb=" C THR A 213 " pdb=" CB THR A 213 " both_signs ideal model delta sigma weight residual False 2.53 2.09 0.43 2.00e-01 2.50e+01 4.66e+00 ... (remaining 585 not shown) Planarity restraints: 708 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB TRP A 51 " 0.020 2.00e-02 2.50e+03 1.84e-02 8.44e+00 pdb=" CG TRP A 51 " -0.048 2.00e-02 2.50e+03 pdb=" CD1 TRP A 51 " 0.025 2.00e-02 2.50e+03 pdb=" CD2 TRP A 51 " -0.004 2.00e-02 2.50e+03 pdb=" NE1 TRP A 51 " -0.000 2.00e-02 2.50e+03 pdb=" CE2 TRP A 51 " -0.006 2.00e-02 2.50e+03 pdb=" CE3 TRP A 51 " 0.004 2.00e-02 2.50e+03 pdb=" CZ2 TRP A 51 " 0.002 2.00e-02 2.50e+03 pdb=" CZ3 TRP A 51 " 0.004 2.00e-02 2.50e+03 pdb=" CH2 TRP A 51 " 0.004 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TRP A 405 " -0.010 2.00e-02 2.50e+03 1.45e-02 5.23e+00 pdb=" CG TRP A 405 " 0.037 2.00e-02 2.50e+03 pdb=" CD1 TRP A 405 " -0.023 2.00e-02 2.50e+03 pdb=" CD2 TRP A 405 " -0.002 2.00e-02 2.50e+03 pdb=" NE1 TRP A 405 " 0.000 2.00e-02 2.50e+03 pdb=" CE2 TRP A 405 " 0.007 2.00e-02 2.50e+03 pdb=" CE3 TRP A 405 " -0.006 2.00e-02 2.50e+03 pdb=" CZ2 TRP A 405 " 0.001 2.00e-02 2.50e+03 pdb=" CZ3 TRP A 405 " -0.003 2.00e-02 2.50e+03 pdb=" CH2 TRP A 405 " -0.001 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C ASP A 49 " 0.038 5.00e-02 4.00e+02 5.68e-02 5.17e+00 pdb=" N PRO A 50 " -0.098 5.00e-02 4.00e+02 pdb=" CA PRO A 50 " 0.029 5.00e-02 4.00e+02 pdb=" CD PRO A 50 " 0.031 5.00e-02 4.00e+02 ... (remaining 705 not shown) Histogram of nonbonded interaction distances: 2.13 - 2.68: 158 2.68 - 3.24: 3304 3.24 - 3.79: 5839 3.79 - 4.35: 8181 4.35 - 4.90: 13624 Nonbonded interactions: 31106 Sorted by model distance: nonbonded pdb=" O PRO A 243 " pdb=" O HOH A 701 " model vdw 2.131 3.040 nonbonded pdb=" OE1 GLU A 84 " pdb=" O HOH A 702 " model vdw 2.140 3.040 nonbonded pdb=" O ALA A 496 " pdb=" O HOH A 703 " model vdw 2.163 3.040 nonbonded pdb=" O MET A 532 " pdb=" O HOH A 704 " model vdw 2.194 3.040 nonbonded pdb=" OD1 ASP A 211 " pdb=" N PHE A 212 " model vdw 2.197 3.120 ... (remaining 31101 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: 0.170 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: 0.160 Check model and map are aligned: 0.020 Set scattering table: 0.040 Process input model: 14.090 Find NCS groups from input model: 0.040 Set up NCS constraints: 0.020 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:11.980 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 26.530 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8603 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.008 0.074 3883 Z= 0.532 Angle : 1.153 16.966 5315 Z= 0.699 Chirality : 0.081 0.802 588 Planarity : 0.010 0.057 708 Dihedral : 9.030 68.110 1415 Min Nonbonded Distance : 2.131 Molprobity Statistics. All-atom Clashscore : 6.96 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.55 % Favored : 96.45 % Rotamer: Outliers : 0.00 % Allowed : 0.96 % Favored : 99.04 % Cbeta Deviations : 0.22 % 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.32 (0.31), residues: 479 helix: -3.43 (0.57), residues: 27 sheet: -0.80 (0.41), residues: 130 loop : -1.78 (0.27), residues: 322 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.048 0.006 TRP A 51 HIS 0.013 0.002 HIS A 494 PHE 0.035 0.004 PHE A 498 TYR 0.025 0.004 TYR A 410 ARG 0.015 0.002 ARG A 495 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 112 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 112 time to evaluate : 0.435 Fit side-chains REVERT: A 84 GLU cc_start: 0.7264 (mp0) cc_final: 0.7053 (mp0) REVERT: A 155 LEU cc_start: 0.8987 (mt) cc_final: 0.8739 (mt) REVERT: A 356 ASP cc_start: 0.8347 (m-30) cc_final: 0.8122 (m-30) REVERT: A 379 GLN cc_start: 0.7891 (mm-40) cc_final: 0.7605 (mm110) REVERT: A 443 LYS cc_start: 0.9112 (tttm) cc_final: 0.8850 (tttm) outliers start: 0 outliers final: 0 residues processed: 112 average time/residue: 1.0901 time to fit residues: 126.2016 Evaluate side-chains 100 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 100 time to evaluate : 0.395 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 48 random chunks: chunk 40 optimal weight: 0.9980 chunk 36 optimal weight: 0.7980 chunk 20 optimal weight: 1.9990 chunk 12 optimal weight: 4.9990 chunk 24 optimal weight: 3.9990 chunk 19 optimal weight: 3.9990 chunk 37 optimal weight: 0.7980 chunk 14 optimal weight: 0.0670 chunk 22 optimal weight: 3.9990 chunk 28 optimal weight: 4.9990 chunk 43 optimal weight: 0.0010 overall best weight: 0.5324 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 58 GLN A 167 HIS A 189 ASN A 347 GLN A 386 ASN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3167 r_free = 0.3167 target = 0.112015 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 49)----------------| | r_work = 0.2954 r_free = 0.2954 target = 0.096746 restraints weight = 4609.366| |-----------------------------------------------------------------------------| r_work (start): 0.2897 rms_B_bonded: 1.58 r_work: 0.2820 rms_B_bonded: 1.83 restraints_weight: 0.5000 r_work: 0.2719 rms_B_bonded: 3.14 restraints_weight: 0.2500 r_work (final): 0.2719 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8744 moved from start: 0.1743 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.024 3883 Z= 0.136 Angle : 0.539 6.936 5315 Z= 0.283 Chirality : 0.045 0.129 588 Planarity : 0.006 0.037 708 Dihedral : 5.077 20.881 519 Min Nonbonded Distance : 2.506 Molprobity Statistics. All-atom Clashscore : 2.68 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.09 % Favored : 97.91 % Rotamer: Outliers : 1.44 % Allowed : 7.21 % Favored : 91.35 % 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.25 (0.34), residues: 479 helix: -1.60 (0.85), residues: 27 sheet: -0.35 (0.40), residues: 155 loop : -1.01 (0.31), residues: 297 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.015 0.002 TRP A 51 HIS 0.004 0.001 HIS A 494 PHE 0.014 0.001 PHE A 64 TYR 0.009 0.001 TYR A 168 ARG 0.003 0.001 ARG A 495 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 107 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 6 poor density : 101 time to evaluate : 0.473 Fit side-chains REVERT: A 187 ARG cc_start: 0.8020 (mtm180) cc_final: 0.7486 (ttp80) REVERT: A 238 SER cc_start: 0.9158 (m) cc_final: 0.8856 (p) REVERT: A 398 ASP cc_start: 0.8663 (OUTLIER) cc_final: 0.8384 (p0) REVERT: A 443 LYS cc_start: 0.9237 (tttm) cc_final: 0.8920 (tttm) outliers start: 6 outliers final: 0 residues processed: 101 average time/residue: 1.1155 time to fit residues: 116.5458 Evaluate side-chains 96 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 1 poor density : 95 time to evaluate : 0.476 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 398 ASP Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 48 random chunks: chunk 38 optimal weight: 8.9990 chunk 44 optimal weight: 2.9990 chunk 45 optimal weight: 3.9990 chunk 7 optimal weight: 1.9990 chunk 6 optimal weight: 1.9990 chunk 26 optimal weight: 10.0000 chunk 19 optimal weight: 0.9980 chunk 0 optimal weight: 7.9990 chunk 13 optimal weight: 0.5980 chunk 37 optimal weight: 2.9990 chunk 36 optimal weight: 7.9990 overall best weight: 1.7186 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 293 HIS A 338 GLN A 386 ASN A 461 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3106 r_free = 0.3106 target = 0.107723 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 40)----------------| | r_work = 0.2885 r_free = 0.2885 target = 0.092132 restraints weight = 4930.415| |-----------------------------------------------------------------------------| r_work (start): 0.2922 rms_B_bonded: 1.59 r_work: 0.2848 rms_B_bonded: 1.76 restraints_weight: 0.5000 r_work: 0.2753 rms_B_bonded: 3.00 restraints_weight: 0.2500 r_work (final): 0.2753 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8808 moved from start: 0.1991 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.030 3883 Z= 0.243 Angle : 0.579 6.568 5315 Z= 0.296 Chirality : 0.046 0.128 588 Planarity : 0.006 0.037 708 Dihedral : 4.967 20.266 519 Min Nonbonded Distance : 2.529 Molprobity Statistics. All-atom Clashscore : 2.54 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.71 % Favored : 97.29 % Rotamer: Outliers : 2.64 % Allowed : 9.86 % Favored : 87.50 % 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.01 (0.34), residues: 479 helix: -1.55 (0.81), residues: 33 sheet: -0.16 (0.42), residues: 144 loop : -0.81 (0.31), residues: 302 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.002 TRP A 51 HIS 0.005 0.001 HIS A 494 PHE 0.014 0.002 PHE A 64 TYR 0.013 0.001 TYR A 168 ARG 0.004 0.001 ARG A 437 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 103 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 11 poor density : 92 time to evaluate : 0.498 Fit side-chains REVERT: A 187 ARG cc_start: 0.8063 (mtm180) cc_final: 0.7679 (mpp-170) REVERT: A 443 LYS cc_start: 0.9205 (tttm) cc_final: 0.8944 (tttm) REVERT: A 525 GLN cc_start: 0.8357 (OUTLIER) cc_final: 0.7731 (pm20) outliers start: 11 outliers final: 7 residues processed: 96 average time/residue: 1.1409 time to fit residues: 113.2036 Evaluate side-chains 95 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 8 poor density : 87 time to evaluate : 0.372 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 118 LYS Chi-restraints excluded: chain A residue 272 LEU Chi-restraints excluded: chain A residue 337 THR Chi-restraints excluded: chain A residue 389 VAL Chi-restraints excluded: chain A residue 455 LEU Chi-restraints excluded: chain A residue 500 THR Chi-restraints excluded: chain A residue 525 GLN Chi-restraints excluded: chain A residue 528 SER Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 48 random chunks: chunk 3 optimal weight: 5.9990 chunk 23 optimal weight: 0.3980 chunk 26 optimal weight: 6.9990 chunk 44 optimal weight: 2.9990 chunk 34 optimal weight: 0.9980 chunk 30 optimal weight: 1.9990 chunk 21 optimal weight: 4.9990 chunk 27 optimal weight: 0.8980 chunk 11 optimal weight: 4.9990 chunk 37 optimal weight: 3.9990 chunk 38 optimal weight: 1.9990 overall best weight: 1.2584 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 386 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3160 r_free = 0.3160 target = 0.111234 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 43)----------------| | r_work = 0.2940 r_free = 0.2940 target = 0.095554 restraints weight = 4849.319| |-----------------------------------------------------------------------------| r_work (start): 0.2891 rms_B_bonded: 1.64 r_work: 0.2814 rms_B_bonded: 1.80 restraints_weight: 0.5000 r_work: 0.2719 rms_B_bonded: 3.08 restraints_weight: 0.2500 r_work (final): 0.2719 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8804 moved from start: 0.2272 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.025 3883 Z= 0.191 Angle : 0.534 6.202 5315 Z= 0.272 Chirality : 0.045 0.129 588 Planarity : 0.006 0.037 708 Dihedral : 4.660 20.328 519 Min Nonbonded Distance : 2.548 Molprobity Statistics. All-atom Clashscore : 2.54 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.92 % Favored : 97.08 % Rotamer: Outliers : 2.64 % Allowed : 12.74 % Favored : 84.62 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.79 (0.35), residues: 479 helix: -1.18 (0.86), residues: 33 sheet: -0.17 (0.41), residues: 154 loop : -0.59 (0.33), residues: 292 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP A 459 HIS 0.003 0.001 HIS A 494 PHE 0.012 0.001 PHE A 64 TYR 0.010 0.001 TYR A 168 ARG 0.002 0.000 ARG A 437 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 102 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 11 poor density : 91 time to evaluate : 0.405 Fit side-chains REVERT: A 187 ARG cc_start: 0.8125 (mtm180) cc_final: 0.7657 (mtp85) REVERT: A 204 THR cc_start: 0.9308 (p) cc_final: 0.9077 (m) REVERT: A 238 SER cc_start: 0.9190 (m) cc_final: 0.8942 (p) REVERT: A 443 LYS cc_start: 0.9189 (tttm) cc_final: 0.8945 (tttm) REVERT: A 525 GLN cc_start: 0.8406 (OUTLIER) cc_final: 0.7801 (pm20) outliers start: 11 outliers final: 5 residues processed: 97 average time/residue: 1.1378 time to fit residues: 114.0949 Evaluate side-chains 96 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 6 poor density : 90 time to evaluate : 0.427 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 118 LYS Chi-restraints excluded: chain A residue 272 LEU Chi-restraints excluded: chain A residue 455 LEU Chi-restraints excluded: chain A residue 500 THR Chi-restraints excluded: chain A residue 525 GLN Chi-restraints excluded: chain A residue 528 SER Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 48 random chunks: chunk 9 optimal weight: 4.9990 chunk 23 optimal weight: 0.9980 chunk 47 optimal weight: 6.9990 chunk 29 optimal weight: 9.9990 chunk 37 optimal weight: 6.9990 chunk 12 optimal weight: 0.0050 chunk 27 optimal weight: 5.9990 chunk 25 optimal weight: 2.9990 chunk 30 optimal weight: 6.9990 chunk 19 optimal weight: 2.9990 chunk 42 optimal weight: 4.9990 overall best weight: 2.4000 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 374 GLN A 386 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3055 r_free = 0.3055 target = 0.104182 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 39)----------------| | r_work = 0.2827 r_free = 0.2827 target = 0.088413 restraints weight = 5125.593| |-----------------------------------------------------------------------------| r_work (start): 0.2898 rms_B_bonded: 1.61 r_work: 0.2824 rms_B_bonded: 1.78 restraints_weight: 0.5000 r_work: 0.2726 rms_B_bonded: 3.04 restraints_weight: 0.2500 r_work (final): 0.2726 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8835 moved from start: 0.2339 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.035 3883 Z= 0.317 Angle : 0.595 6.411 5315 Z= 0.303 Chirality : 0.048 0.139 588 Planarity : 0.006 0.040 708 Dihedral : 4.934 20.610 519 Min Nonbonded Distance : 2.543 Molprobity Statistics. All-atom Clashscore : 2.54 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.34 % Favored : 96.66 % Rotamer: Outliers : 2.88 % Allowed : 12.74 % Favored : 84.38 % 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.78 (0.35), residues: 479 helix: -1.42 (0.82), residues: 33 sheet: -0.11 (0.41), residues: 154 loop : -0.59 (0.32), residues: 292 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP A 51 HIS 0.005 0.001 HIS A 299 PHE 0.013 0.002 PHE A 64 TYR 0.016 0.002 TYR A 214 ARG 0.004 0.001 ARG A 437 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 106 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 12 poor density : 94 time to evaluate : 0.425 Fit side-chains REVERT: A 204 THR cc_start: 0.9332 (p) cc_final: 0.9100 (m) REVERT: A 398 ASP cc_start: 0.8542 (p0) cc_final: 0.8184 (p0) REVERT: A 443 LYS cc_start: 0.9229 (tttm) cc_final: 0.9007 (tttm) REVERT: A 525 GLN cc_start: 0.8465 (OUTLIER) cc_final: 0.7870 (pm20) outliers start: 12 outliers final: 6 residues processed: 101 average time/residue: 1.1050 time to fit residues: 115.4939 Evaluate side-chains 99 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 7 poor density : 92 time to evaluate : 0.409 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 118 LYS Chi-restraints excluded: chain A residue 272 LEU Chi-restraints excluded: chain A residue 383 LEU Chi-restraints excluded: chain A residue 455 LEU Chi-restraints excluded: chain A residue 500 THR Chi-restraints excluded: chain A residue 525 GLN Chi-restraints excluded: chain A residue 528 SER Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 48 random chunks: chunk 0 optimal weight: 6.9990 chunk 35 optimal weight: 2.9990 chunk 43 optimal weight: 3.9990 chunk 5 optimal weight: 3.9990 chunk 45 optimal weight: 0.0020 chunk 46 optimal weight: 0.8980 chunk 26 optimal weight: 3.9990 chunk 36 optimal weight: 6.9990 chunk 9 optimal weight: 4.9990 chunk 29 optimal weight: 0.3980 chunk 30 optimal weight: 4.9990 overall best weight: 1.6592 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 386 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3088 r_free = 0.3088 target = 0.106340 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 42)----------------| | r_work = 0.2865 r_free = 0.2865 target = 0.090753 restraints weight = 5003.863| |-----------------------------------------------------------------------------| r_work (start): 0.2916 rms_B_bonded: 1.59 r_work: 0.2843 rms_B_bonded: 1.76 restraints_weight: 0.5000 r_work: 0.2746 rms_B_bonded: 3.00 restraints_weight: 0.2500 r_work (final): 0.2746 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8820 moved from start: 0.2488 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.029 3883 Z= 0.232 Angle : 0.545 6.463 5315 Z= 0.278 Chirality : 0.046 0.136 588 Planarity : 0.006 0.041 708 Dihedral : 4.754 20.785 519 Min Nonbonded Distance : 2.547 Molprobity Statistics. All-atom Clashscore : 2.14 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.13 % Favored : 96.87 % Rotamer: Outliers : 2.40 % Allowed : 13.46 % Favored : 84.13 % 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.64 (0.35), residues: 479 helix: -1.14 (0.85), residues: 33 sheet: 0.00 (0.41), residues: 154 loop : -0.52 (0.32), residues: 292 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP A 522 HIS 0.003 0.001 HIS A 494 PHE 0.012 0.001 PHE A 64 TYR 0.011 0.001 TYR A 214 ARG 0.003 0.001 ARG A 437 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 106 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 10 poor density : 96 time to evaluate : 0.416 Fit side-chains REVERT: A 187 ARG cc_start: 0.8084 (mtp85) cc_final: 0.7666 (mtp85) REVERT: A 204 THR cc_start: 0.9313 (p) cc_final: 0.9085 (m) REVERT: A 443 LYS cc_start: 0.9258 (tttm) cc_final: 0.9042 (tttm) REVERT: A 525 GLN cc_start: 0.8364 (OUTLIER) cc_final: 0.7772 (pm20) outliers start: 10 outliers final: 6 residues processed: 101 average time/residue: 1.0918 time to fit residues: 114.2021 Evaluate side-chains 101 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 7 poor density : 94 time to evaluate : 0.414 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 118 LYS Chi-restraints excluded: chain A residue 272 LEU Chi-restraints excluded: chain A residue 383 LEU Chi-restraints excluded: chain A residue 455 LEU Chi-restraints excluded: chain A residue 500 THR Chi-restraints excluded: chain A residue 525 GLN Chi-restraints excluded: chain A residue 528 SER Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 48 random chunks: chunk 44 optimal weight: 3.9990 chunk 18 optimal weight: 3.9990 chunk 13 optimal weight: 0.7980 chunk 35 optimal weight: 3.9990 chunk 6 optimal weight: 1.9990 chunk 4 optimal weight: 3.9990 chunk 3 optimal weight: 2.9990 chunk 47 optimal weight: 4.9990 chunk 16 optimal weight: 0.9980 chunk 33 optimal weight: 0.7980 chunk 27 optimal weight: 0.0060 overall best weight: 0.9198 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 386 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3106 r_free = 0.3106 target = 0.107820 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 36)----------------| | r_work = 0.2886 r_free = 0.2886 target = 0.092264 restraints weight = 4917.590| |-----------------------------------------------------------------------------| r_work (start): 0.2924 rms_B_bonded: 1.59 r_work: 0.2852 rms_B_bonded: 1.75 restraints_weight: 0.5000 r_work: 0.2757 rms_B_bonded: 2.98 restraints_weight: 0.2500 r_work (final): 0.2757 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8789 moved from start: 0.2642 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.022 3883 Z= 0.155 Angle : 0.500 6.338 5315 Z= 0.253 Chirality : 0.044 0.131 588 Planarity : 0.005 0.042 708 Dihedral : 4.394 19.054 519 Min Nonbonded Distance : 2.556 Molprobity Statistics. All-atom Clashscore : 2.41 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.92 % Favored : 97.08 % Rotamer: Outliers : 1.68 % Allowed : 14.66 % Favored : 83.65 % 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.41 (0.36), residues: 479 helix: -0.85 (0.89), residues: 33 sheet: 0.12 (0.42), residues: 154 loop : -0.36 (0.33), residues: 292 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A 522 HIS 0.003 0.001 HIS A 299 PHE 0.012 0.001 PHE A 64 TYR 0.010 0.001 TYR A 362 ARG 0.002 0.000 ARG A 94 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 105 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 7 poor density : 98 time to evaluate : 0.455 Fit side-chains REVERT: A 114 PHE cc_start: 0.8195 (m-80) cc_final: 0.7942 (m-80) REVERT: A 187 ARG cc_start: 0.8061 (mtp85) cc_final: 0.7615 (ttp80) REVERT: A 238 SER cc_start: 0.9159 (m) cc_final: 0.8940 (p) REVERT: A 379 GLN cc_start: 0.8077 (mm-40) cc_final: 0.7780 (mm-40) REVERT: A 398 ASP cc_start: 0.8738 (p0) cc_final: 0.8466 (p0) REVERT: A 443 LYS cc_start: 0.9185 (tttm) cc_final: 0.8983 (tttm) outliers start: 7 outliers final: 5 residues processed: 102 average time/residue: 1.1366 time to fit residues: 119.9141 Evaluate side-chains 99 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 94 time to evaluate : 0.457 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 118 LYS Chi-restraints excluded: chain A residue 383 LEU Chi-restraints excluded: chain A residue 455 LEU Chi-restraints excluded: chain A residue 500 THR Chi-restraints excluded: chain A residue 528 SER Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 48 random chunks: chunk 22 optimal weight: 6.9990 chunk 23 optimal weight: 6.9990 chunk 24 optimal weight: 8.9990 chunk 33 optimal weight: 1.9990 chunk 6 optimal weight: 0.9980 chunk 14 optimal weight: 0.9990 chunk 8 optimal weight: 2.9990 chunk 47 optimal weight: 6.9990 chunk 37 optimal weight: 7.9990 chunk 10 optimal weight: 0.9990 chunk 5 optimal weight: 4.9990 overall best weight: 1.5988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3093 r_free = 0.3093 target = 0.106772 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 40)----------------| | r_work = 0.2871 r_free = 0.2871 target = 0.091191 restraints weight = 4876.495| |-----------------------------------------------------------------------------| r_work (start): 0.2881 rms_B_bonded: 1.57 r_work: 0.2806 rms_B_bonded: 1.71 restraints_weight: 0.5000 r_work: 0.2712 rms_B_bonded: 2.92 restraints_weight: 0.2500 r_work (final): 0.2712 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8777 moved from start: 0.2614 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.028 3883 Z= 0.223 Angle : 0.538 6.220 5315 Z= 0.272 Chirality : 0.046 0.135 588 Planarity : 0.006 0.040 708 Dihedral : 4.540 19.131 519 Min Nonbonded Distance : 2.551 Molprobity Statistics. All-atom Clashscore : 2.01 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.76 % Favored : 96.24 % Rotamer: Outliers : 1.92 % Allowed : 15.38 % 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.42 (0.36), residues: 479 helix: -0.96 (0.87), residues: 33 sheet: 0.15 (0.41), residues: 154 loop : -0.38 (0.33), residues: 292 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.003 0.001 TRP A 522 HIS 0.003 0.001 HIS A 494 PHE 0.012 0.001 PHE A 64 TYR 0.012 0.001 TYR A 214 ARG 0.003 0.001 ARG A 409 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 101 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 8 poor density : 93 time to evaluate : 0.460 Fit side-chains REVERT: A 114 PHE cc_start: 0.8229 (m-80) cc_final: 0.7964 (m-80) REVERT: A 187 ARG cc_start: 0.8006 (mtp85) cc_final: 0.7562 (ttp80) REVERT: A 379 GLN cc_start: 0.8121 (mm-40) cc_final: 0.7783 (mm-40) REVERT: A 443 LYS cc_start: 0.9238 (tttm) cc_final: 0.9035 (tttm) outliers start: 8 outliers final: 6 residues processed: 99 average time/residue: 1.1228 time to fit residues: 115.0017 Evaluate side-chains 99 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 6 poor density : 93 time to evaluate : 0.453 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 118 LYS Chi-restraints excluded: chain A residue 272 LEU Chi-restraints excluded: chain A residue 383 LEU Chi-restraints excluded: chain A residue 455 LEU Chi-restraints excluded: chain A residue 500 THR Chi-restraints excluded: chain A residue 528 SER Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 48 random chunks: chunk 14 optimal weight: 6.9990 chunk 9 optimal weight: 0.5980 chunk 35 optimal weight: 4.9990 chunk 5 optimal weight: 4.9990 chunk 41 optimal weight: 0.0010 chunk 7 optimal weight: 4.9990 chunk 45 optimal weight: 0.9980 chunk 36 optimal weight: 6.9990 chunk 16 optimal weight: 0.9990 chunk 32 optimal weight: 6.9990 chunk 25 optimal weight: 1.9990 overall best weight: 0.9190 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3168 r_free = 0.3168 target = 0.111882 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 32)----------------| | r_work = 0.2952 r_free = 0.2952 target = 0.096270 restraints weight = 4813.504| |-----------------------------------------------------------------------------| r_work (start): 0.2955 rms_B_bonded: 1.63 r_work: 0.2881 rms_B_bonded: 1.82 restraints_weight: 0.5000 r_work: 0.2784 rms_B_bonded: 3.12 restraints_weight: 0.2500 r_work (final): 0.2784 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8779 moved from start: 0.2694 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.025 3883 Z= 0.157 Angle : 0.510 6.234 5315 Z= 0.256 Chirality : 0.044 0.129 588 Planarity : 0.005 0.041 708 Dihedral : 4.337 18.912 519 Min Nonbonded Distance : 2.556 Molprobity Statistics. All-atom Clashscore : 2.68 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.92 % Favored : 97.08 % Rotamer: Outliers : 1.44 % Allowed : 15.38 % Favored : 83.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: -0.30 (0.36), residues: 479 helix: -0.75 (0.91), residues: 33 sheet: 0.16 (0.42), residues: 154 loop : -0.27 (0.33), residues: 292 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A 522 HIS 0.002 0.001 HIS A 299 PHE 0.011 0.001 PHE A 64 TYR 0.009 0.001 TYR A 168 ARG 0.005 0.000 ARG A 409 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 105 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 6 poor density : 99 time to evaluate : 0.468 Fit side-chains REVERT: A 187 ARG cc_start: 0.8134 (mtp85) cc_final: 0.7691 (ttp80) REVERT: A 238 SER cc_start: 0.9185 (m) cc_final: 0.8942 (p) REVERT: A 379 GLN cc_start: 0.8068 (mm-40) cc_final: 0.7733 (mm-40) REVERT: A 384 LYS cc_start: 0.8385 (mttt) cc_final: 0.8071 (mttt) REVERT: A 398 ASP cc_start: 0.8813 (p0) cc_final: 0.8517 (p0) REVERT: A 443 LYS cc_start: 0.9221 (tttm) cc_final: 0.9019 (tttm) REVERT: A 452 ASP cc_start: 0.8837 (m-30) cc_final: 0.8427 (m-30) outliers start: 6 outliers final: 5 residues processed: 103 average time/residue: 1.0614 time to fit residues: 113.2634 Evaluate side-chains 97 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 92 time to evaluate : 0.462 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 118 LYS Chi-restraints excluded: chain A residue 383 LEU Chi-restraints excluded: chain A residue 455 LEU Chi-restraints excluded: chain A residue 500 THR Chi-restraints excluded: chain A residue 528 SER Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 48 random chunks: chunk 46 optimal weight: 0.1980 chunk 2 optimal weight: 3.9990 chunk 44 optimal weight: 2.9990 chunk 23 optimal weight: 0.8980 chunk 36 optimal weight: 5.9990 chunk 20 optimal weight: 6.9990 chunk 24 optimal weight: 3.9990 chunk 18 optimal weight: 2.9990 chunk 26 optimal weight: 0.8980 chunk 12 optimal weight: 0.0970 chunk 1 optimal weight: 0.0980 overall best weight: 0.4378 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3159 r_free = 0.3159 target = 0.111535 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 41)----------------| | r_work = 0.2944 r_free = 0.2944 target = 0.095972 restraints weight = 4929.961| |-----------------------------------------------------------------------------| r_work (start): 0.2977 rms_B_bonded: 1.60 r_work: 0.2904 rms_B_bonded: 1.77 restraints_weight: 0.5000 r_work: 0.2811 rms_B_bonded: 3.03 restraints_weight: 0.2500 r_work (final): 0.2811 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8752 moved from start: 0.2883 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.024 3883 Z= 0.115 Angle : 0.488 6.207 5315 Z= 0.242 Chirality : 0.044 0.124 588 Planarity : 0.005 0.040 708 Dihedral : 4.052 18.604 519 Min Nonbonded Distance : 2.570 Molprobity Statistics. All-atom Clashscore : 2.95 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.51 % Favored : 97.49 % Rotamer: Outliers : 1.20 % Allowed : 15.62 % Favored : 83.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: -0.11 (0.36), residues: 479 helix: -0.58 (0.94), residues: 33 sheet: 0.30 (0.43), residues: 145 loop : -0.12 (0.33), residues: 301 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.007 0.001 TRP A 522 HIS 0.002 0.001 HIS A 299 PHE 0.010 0.001 PHE A 64 TYR 0.008 0.001 TYR A 168 ARG 0.003 0.000 ARG A 69 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 958 Ramachandran restraints generated. 479 Oldfield, 0 Emsley, 479 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 98 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 93 time to evaluate : 0.532 Fit side-chains REVERT: A 187 ARG cc_start: 0.8216 (mtp85) cc_final: 0.7802 (ttp80) REVERT: A 384 LYS cc_start: 0.8320 (mttt) cc_final: 0.8115 (mttt) REVERT: A 409 ARG cc_start: 0.8421 (mmm-85) cc_final: 0.8167 (mtp180) REVERT: A 443 LYS cc_start: 0.9125 (tttm) cc_final: 0.8919 (tttm) outliers start: 5 outliers final: 4 residues processed: 97 average time/residue: 1.1901 time to fit residues: 119.2511 Evaluate side-chains 94 residues out of total 416 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 4 poor density : 90 time to evaluate : 0.413 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 118 LYS Chi-restraints excluded: chain A residue 383 LEU Chi-restraints excluded: chain A residue 455 LEU Chi-restraints excluded: chain A residue 528 SER Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 48 random chunks: chunk 21 optimal weight: 4.9990 chunk 44 optimal weight: 3.9990 chunk 12 optimal weight: 4.9990 chunk 31 optimal weight: 0.6980 chunk 43 optimal weight: 2.9990 chunk 14 optimal weight: 5.9990 chunk 30 optimal weight: 0.0470 chunk 15 optimal weight: 5.9990 chunk 3 optimal weight: 0.5980 chunk 40 optimal weight: 5.9990 chunk 11 optimal weight: 1.9990 overall best weight: 1.2682 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3127 r_free = 0.3127 target = 0.109098 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 38)----------------| | r_work = 0.2909 r_free = 0.2909 target = 0.093591 restraints weight = 4858.304| |-----------------------------------------------------------------------------| r_work (start): 0.2914 rms_B_bonded: 1.58 r_work: 0.2841 rms_B_bonded: 1.74 restraints_weight: 0.5000 r_work: 0.2747 rms_B_bonded: 2.98 restraints_weight: 0.2500 r_work (final): 0.2747 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8791 moved from start: 0.2797 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.024 3883 Z= 0.191 Angle : 0.530 6.092 5315 Z= 0.264 Chirality : 0.045 0.123 588 Planarity : 0.005 0.039 708 Dihedral : 4.305 18.231 519 Min Nonbonded Distance : 2.567 Molprobity Statistics. All-atom Clashscore : 2.81 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.55 % Favored : 96.45 % Rotamer: Outliers : 0.96 % Allowed : 16.35 % 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.13 (0.36), residues: 479 helix: -0.68 (0.91), residues: 33 sheet: 0.26 (0.42), residues: 153 loop : -0.12 (0.33), residues: 293 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.004 0.001 TRP A 522 HIS 0.003 0.001 HIS A 494 PHE 0.011 0.001 PHE A 64 TYR 0.010 0.001 TYR A 214 ARG 0.003 0.000 ARG A 69 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 3673.17 seconds wall clock time: 63 minutes 24.10 seconds (3804.10 seconds total)