Starting phenix.real_space_refine on Sun Apr 27 06:51:23 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/7yr5_34046/04_2025/7yr5_34046.cif Found real_map, /net/cci-nas-00/data/ceres_data/7yr5_34046/04_2025/7yr5_34046.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.63 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/7yr5_34046/04_2025/7yr5_34046.map" default_real_map = "/net/cci-nas-00/data/ceres_data/7yr5_34046/04_2025/7yr5_34046.map" model { file = "/net/cci-nas-00/data/ceres_data/7yr5_34046/04_2025/7yr5_34046.cif" } default_model = "/net/cci-nas-00/data/ceres_data/7yr5_34046/04_2025/7yr5_34046.cif" } resolution = 3.63 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= -0.000 sd= 0.015 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 4 Type Number sf(0) Gaussians S 35 5.16 5 C 3164 2.51 5 N 775 2.21 5 O 854 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 7 residue(s): 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/chem_data/mon_lib" Total number of atoms: 4828 Number of models: 1 Model: "" Number of chains: 2 Chain: "B" Number of atoms: 1899 Number of conformers: 1 Conformer: "" Number of residues, atoms: 237, 1899 Classifications: {'peptide': 237} Link IDs: {'PTRANS': 8, 'TRANS': 228} Chain: "A" Number of atoms: 2929 Number of conformers: 1 Conformer: "" Number of residues, atoms: 385, 2929 Classifications: {'peptide': 385} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 17, 'TRANS': 367} Chain breaks: 1 Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 1 Unresolved non-hydrogen dihedrals: 1 Time building chain proxies: 4.12, per 1000 atoms: 0.85 Number of scatterers: 4828 At special positions: 0 Unit cell: (73.08, 103.32, 123.48, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 4 Type Number sf(0) S 35 16.00 O 854 8.00 N 775 7.00 C 3164 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=2, symmetry=0 Simple disulfide: pdb=" SG CYS B 88 " - pdb=" SG CYS B 142 " distance=2.04 Simple disulfide: pdb=" SG CYS B 180 " - pdb=" SG CYS B 237 " distance=2.03 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: 1.24 Conformation dependent library (CDL) restraints added in 629.1 milliseconds 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1126 Finding SS restraints... Secondary structure from input PDB file: 22 helices and 5 sheets defined 61.7% alpha, 7.4% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.55 Creating SS restraints... Processing helix chain 'B' and resid 257 through 290 removed outlier: 4.103A pdb=" N LYS B 262 " --> pdb=" O LEU B 258 " (cutoff:3.500A) Proline residue: B 263 - end of helix Processing helix chain 'A' and resid 19 through 40 removed outlier: 4.237A pdb=" N PHE A 36 " --> pdb=" O PHE A 32 " (cutoff:3.500A) Proline residue: A 37 - end of helix Processing helix chain 'A' and resid 40 through 51 removed outlier: 4.150A pdb=" N LYS A 45 " --> pdb=" O THR A 41 " (cutoff:3.500A) removed outlier: 4.608A pdb=" N GLU A 46 " --> pdb=" O VAL A 42 " (cutoff:3.500A) removed outlier: 4.322A pdb=" N GLY A 49 " --> pdb=" O LYS A 45 " (cutoff:3.500A) removed outlier: 3.517A pdb=" N ILE A 50 " --> pdb=" O GLU A 46 " (cutoff:3.500A) Processing helix chain 'A' and resid 54 through 72 removed outlier: 5.085A pdb=" N ILE A 61 " --> pdb=" O GLU A 57 " (cutoff:3.500A) removed outlier: 4.715A pdb=" N SER A 62 " --> pdb=" O VAL A 58 " (cutoff:3.500A) Processing helix chain 'A' and resid 72 through 83 Processing helix chain 'A' and resid 84 through 104 Processing helix chain 'A' and resid 111 through 139 removed outlier: 3.690A pdb=" N CYS A 115 " --> pdb=" O GLN A 111 " (cutoff:3.500A) Proline residue: A 130 - end of helix Processing helix chain 'A' and resid 142 through 152 Processing helix chain 'A' and resid 152 through 162 removed outlier: 4.313A pdb=" N VAL A 156 " --> pdb=" O ALA A 152 " (cutoff:3.500A) Processing helix chain 'A' and resid 163 through 172 removed outlier: 5.924A pdb=" N PHE A 169 " --> pdb=" O ASN A 165 " (cutoff:3.500A) removed outlier: 3.720A pdb=" N GLY A 170 " --> pdb=" O GLN A 166 " (cutoff:3.500A) removed outlier: 3.602A pdb=" N ILE A 171 " --> pdb=" O VAL A 167 " (cutoff:3.500A) Processing helix chain 'A' and resid 175 through 188 removed outlier: 3.519A pdb=" N ILE A 180 " --> pdb=" O GLY A 176 " (cutoff:3.500A) removed outlier: 3.760A pdb=" N LEU A 181 " --> pdb=" O SER A 177 " (cutoff:3.500A) Processing helix chain 'A' and resid 189 through 194 Processing helix chain 'A' and resid 254 through 259 removed outlier: 3.551A pdb=" N LEU A 257 " --> pdb=" O ASP A 254 " (cutoff:3.500A) Processing helix chain 'A' and resid 260 through 291 removed outlier: 4.491A pdb=" N PHE A 275 " --> pdb=" O VAL A 271 " (cutoff:3.500A) removed outlier: 4.633A pdb=" N GLY A 276 " --> pdb=" O ILE A 272 " (cutoff:3.500A) Proline residue: A 280 - end of helix removed outlier: 5.279A pdb=" N SER A 285 " --> pdb=" O LEU A 281 " (cutoff:3.500A) removed outlier: 5.017A pdb=" N SER A 286 " --> pdb=" O VAL A 282 " (cutoff:3.500A) Processing helix chain 'A' and resid 294 through 322 removed outlier: 3.590A pdb=" N LEU A 302 " --> pdb=" O SER A 298 " (cutoff:3.500A) Proline residue: A 314 - end of helix Processing helix chain 'A' and resid 323 through 326 Processing helix chain 'A' and resid 328 through 347 removed outlier: 3.607A pdb=" N VAL A 338 " --> pdb=" O ALA A 334 " (cutoff:3.500A) removed outlier: 3.698A pdb=" N ALA A 339 " --> pdb=" O ALA A 335 " (cutoff:3.500A) removed outlier: 3.661A pdb=" N ASN A 340 " --> pdb=" O SER A 336 " (cutoff:3.500A) Processing helix chain 'A' and resid 352 through 374 removed outlier: 3.859A pdb=" N TRP A 369 " --> pdb=" O PHE A 365 " (cutoff:3.500A) removed outlier: 4.302A pdb=" N LEU A 370 " --> pdb=" O ALA A 366 " (cutoff:3.500A) Processing helix chain 'A' and resid 375 through 382 Processing helix chain 'A' and resid 383 through 412 removed outlier: 3.983A pdb=" N SER A 388 " --> pdb=" O PRO A 384 " (cutoff:3.500A) removed outlier: 5.138A pdb=" N SER A 389 " --> pdb=" O GLN A 385 " (cutoff:3.500A) removed outlier: 4.108A pdb=" N GLY A 392 " --> pdb=" O SER A 388 " (cutoff:3.500A) removed outlier: 3.634A pdb=" N LEU A 393 " --> pdb=" O SER A 389 " (cutoff:3.500A) removed outlier: 3.742A pdb=" N ILE A 396 " --> pdb=" O GLY A 392 " (cutoff:3.500A) removed outlier: 3.606A pdb=" N VAL A 397 " --> pdb=" O LEU A 393 " (cutoff:3.500A) removed outlier: 4.562A pdb=" N CYS A 400 " --> pdb=" O ILE A 396 " (cutoff:3.500A) Proline residue: A 401 - end of helix Proline residue: A 406 - end of helix Processing helix chain 'A' and resid 418 through 420 No H-bonds generated for 'chain 'A' and resid 418 through 420' Processing helix chain 'A' and resid 421 through 448 removed outlier: 3.934A pdb=" N ALA A 425 " --> pdb=" O TYR A 421 " (cutoff:3.500A) removed outlier: 3.837A pdb=" N CYS A 426 " --> pdb=" O THR A 422 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'B' and resid 61 through 64 removed outlier: 5.621A pdb=" N ILE B 62 " --> pdb=" O GLN B 89 " (cutoff:3.500A) removed outlier: 7.364A pdb=" N THR B 91 " --> pdb=" O ILE B 62 " (cutoff:3.500A) removed outlier: 6.641A pdb=" N LEU B 64 " --> pdb=" O THR B 91 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'B' and resid 72 through 74 removed outlier: 3.715A pdb=" N GLN B 150 " --> pdb=" O PHE B 144 " (cutoff:3.500A) removed outlier: 4.314A pdb=" N TRP B 103 " --> pdb=" O LEU B 110 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'B' and resid 169 through 171 Processing sheet with id=AA4, first strand: chain 'B' and resid 191 through 194 Processing sheet with id=AA5, first strand: chain 'B' and resid 210 through 213 278 hydrogen bonds defined for protein. 798 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.41 Time building geometry restraints manager: 1.44 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: 1472 1.34 - 1.46: 1090 1.46 - 1.58: 2336 1.58 - 1.70: 0 1.70 - 1.81: 51 Bond restraints: 4949 Sorted by residual: bond pdb=" N VAL B 84 " pdb=" CA VAL B 84 " ideal model delta sigma weight residual 1.456 1.493 -0.038 1.11e-02 8.12e+03 1.16e+01 bond pdb=" N ILE B 76 " pdb=" CA ILE B 76 " ideal model delta sigma weight residual 1.458 1.495 -0.037 1.19e-02 7.06e+03 9.69e+00 bond pdb=" N VAL B 116 " pdb=" CA VAL B 116 " ideal model delta sigma weight residual 1.457 1.493 -0.036 1.22e-02 6.72e+03 8.68e+00 bond pdb=" N VAL B 72 " pdb=" CA VAL B 72 " ideal model delta sigma weight residual 1.459 1.494 -0.036 1.24e-02 6.50e+03 8.36e+00 bond pdb=" N LYS B 104 " pdb=" CA LYS B 104 " ideal model delta sigma weight residual 1.454 1.489 -0.035 1.27e-02 6.20e+03 7.61e+00 ... (remaining 4944 not shown) Histogram of bond angle deviations from ideal: 0.00 - 1.38: 6220 1.38 - 2.75: 367 2.75 - 4.13: 89 4.13 - 5.51: 33 5.51 - 6.88: 6 Bond angle restraints: 6715 Sorted by residual: angle pdb=" C ARG B 128 " pdb=" N PHE B 129 " pdb=" CA PHE B 129 " ideal model delta sigma weight residual 122.77 117.48 5.29 1.33e+00 5.65e-01 1.58e+01 angle pdb=" CA PHE B 129 " pdb=" CB PHE B 129 " pdb=" CG PHE B 129 " ideal model delta sigma weight residual 113.80 117.77 -3.97 1.00e+00 1.00e+00 1.58e+01 angle pdb=" N PHE B 154 " pdb=" CA PHE B 154 " pdb=" C PHE B 154 " ideal model delta sigma weight residual 109.46 103.18 6.28 1.66e+00 3.63e-01 1.43e+01 angle pdb=" C CYS A 106 " pdb=" N ASN A 107 " pdb=" CA ASN A 107 " ideal model delta sigma weight residual 122.61 128.02 -5.41 1.56e+00 4.11e-01 1.20e+01 angle pdb=" C CYS B 88 " pdb=" N GLN B 89 " pdb=" CA GLN B 89 " ideal model delta sigma weight residual 123.05 118.23 4.82 1.40e+00 5.10e-01 1.18e+01 ... (remaining 6710 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 14.54: 2481 14.54 - 29.07: 292 29.07 - 43.61: 82 43.61 - 58.14: 13 58.14 - 72.68: 7 Dihedral angle restraints: 2875 sinusoidal: 1101 harmonic: 1774 Sorted by residual: dihedral pdb=" CA GLN A 166 " pdb=" C GLN A 166 " pdb=" N VAL A 167 " pdb=" CA VAL A 167 " ideal model delta harmonic sigma weight residual 180.00 153.92 26.08 0 5.00e+00 4.00e-02 2.72e+01 dihedral pdb=" CA LEU A 164 " pdb=" C LEU A 164 " pdb=" N ASN A 165 " pdb=" CA ASN A 165 " ideal model delta harmonic sigma weight residual 180.00 158.35 21.65 0 5.00e+00 4.00e-02 1.88e+01 dihedral pdb=" CA GLN A 110 " pdb=" C GLN A 110 " pdb=" N GLN A 111 " pdb=" CA GLN A 111 " ideal model delta harmonic sigma weight residual 180.00 159.98 20.02 0 5.00e+00 4.00e-02 1.60e+01 ... (remaining 2872 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.052: 571 0.052 - 0.104: 141 0.104 - 0.156: 34 0.156 - 0.208: 12 0.208 - 0.261: 4 Chirality restraints: 762 Sorted by residual: chirality pdb=" CA PHE B 154 " pdb=" N PHE B 154 " pdb=" C PHE B 154 " pdb=" CB PHE B 154 " both_signs ideal model delta sigma weight residual False 2.51 2.77 -0.26 2.00e-01 2.50e+01 1.70e+00 chirality pdb=" CA CYS B 88 " pdb=" N CYS B 88 " pdb=" C CYS B 88 " pdb=" CB CYS B 88 " both_signs ideal model delta sigma weight residual False 2.51 2.75 -0.24 2.00e-01 2.50e+01 1.39e+00 chirality pdb=" CA ILE B 62 " pdb=" N ILE B 62 " pdb=" C ILE B 62 " pdb=" CB ILE B 62 " both_signs ideal model delta sigma weight residual False 2.43 2.66 -0.22 2.00e-01 2.50e+01 1.26e+00 ... (remaining 759 not shown) Planarity restraints: 834 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA VAL B 253 " 0.021 2.00e-02 2.50e+03 4.20e-02 1.76e+01 pdb=" C VAL B 253 " -0.073 2.00e-02 2.50e+03 pdb=" O VAL B 253 " 0.027 2.00e-02 2.50e+03 pdb=" N VAL B 254 " 0.024 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA VAL B 99 " -0.012 2.00e-02 2.50e+03 2.48e-02 6.15e+00 pdb=" C VAL B 99 " 0.043 2.00e-02 2.50e+03 pdb=" O VAL B 99 " -0.016 2.00e-02 2.50e+03 pdb=" N ASN B 100 " -0.014 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA VAL A 114 " -0.009 2.00e-02 2.50e+03 1.85e-02 3.42e+00 pdb=" C VAL A 114 " 0.032 2.00e-02 2.50e+03 pdb=" O VAL A 114 " -0.012 2.00e-02 2.50e+03 pdb=" N CYS A 115 " -0.011 2.00e-02 2.50e+03 ... (remaining 831 not shown) Histogram of nonbonded interaction distances: 2.25 - 2.78: 1062 2.78 - 3.31: 4795 3.31 - 3.84: 7906 3.84 - 4.37: 9100 4.37 - 4.90: 15390 Nonbonded interactions: 38253 Sorted by model distance: nonbonded pdb=" OG SER A 63 " pdb=" O LEU A 403 " model vdw 2.250 3.040 nonbonded pdb=" OG SER A 29 " pdb=" O ASN A 127 " model vdw 2.276 3.040 nonbonded pdb=" OH TYR A 138 " pdb=" O ARG A 196 " model vdw 2.290 3.040 nonbonded pdb=" OH TYR B 210 " pdb=" OH TYR B 235 " model vdw 2.292 3.040 nonbonded pdb=" OG SER B 134 " pdb=" O GLU B 160 " model vdw 2.318 3.040 ... (remaining 38248 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.00 ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 2.890 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.190 Check model and map are aligned: 0.030 Set scattering table: 0.060 Process input model: 16.790 Find NCS groups from input model: 0.050 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:3.090 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 23.130 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7635 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.005 0.042 4951 Z= 0.340 Angle : 0.807 6.881 6719 Z= 0.501 Chirality : 0.054 0.261 762 Planarity : 0.005 0.042 834 Dihedral : 13.919 72.678 1743 Min Nonbonded Distance : 2.250 Molprobity Statistics. All-atom Clashscore : 8.97 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.52 % Favored : 94.48 % Rotamer: Outliers : 0.00 % Allowed : 1.52 % Favored : 98.48 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.05 (0.34), residues: 616 helix: 0.65 (0.29), residues: 362 sheet: 0.32 (0.53), residues: 102 loop : -1.81 (0.42), residues: 152 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP B 236 HIS 0.002 0.000 HIS B 249 PHE 0.017 0.001 PHE B 129 TYR 0.013 0.001 TYR A 445 ARG 0.003 0.000 ARG B 128 Details of bonding type rmsd hydrogen bonds : bond 0.16864 ( 278) hydrogen bonds : angle 6.21610 ( 798) SS BOND : bond 0.00626 ( 2) SS BOND : angle 0.62855 ( 4) covalent geometry : bond 0.00491 ( 4949) covalent geometry : angle 0.80756 ( 6715) *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 112 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 112 time to evaluate : 0.527 Fit side-chains REVERT: A 151 MET cc_start: 0.8056 (ttm) cc_final: 0.7434 (ttp) REVERT: A 416 TYR cc_start: 0.7273 (t80) cc_final: 0.6929 (t80) REVERT: A 418 ASP cc_start: 0.7536 (t70) cc_final: 0.7060 (p0) REVERT: A 446 ARG cc_start: 0.8201 (mmp80) cc_final: 0.7579 (mtp180) outliers start: 0 outliers final: 0 residues processed: 112 average time/residue: 0.1589 time to fit residues: 23.6512 Evaluate side-chains 84 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 84 time to evaluate : 0.520 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 61 random chunks: chunk 51 optimal weight: 2.9990 chunk 46 optimal weight: 1.9990 chunk 25 optimal weight: 1.9990 chunk 15 optimal weight: 10.0000 chunk 31 optimal weight: 3.9990 chunk 24 optimal weight: 2.9990 chunk 47 optimal weight: 2.9990 chunk 18 optimal weight: 1.9990 chunk 29 optimal weight: 2.9990 chunk 35 optimal weight: 3.9990 chunk 55 optimal weight: 0.9990 overall best weight: 1.9990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 89 GLN A 165 ASN A 260 HIS Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4129 r_free = 0.4129 target = 0.162265 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 42)----------------| | r_work = 0.3559 r_free = 0.3559 target = 0.118176 restraints weight = 7829.141| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 47)----------------| | r_work = 0.3607 r_free = 0.3607 target = 0.121505 restraints weight = 3758.848| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 46)----------------| | r_work = 0.3637 r_free = 0.3637 target = 0.123556 restraints weight = 2307.192| |-----------------------------------------------------------------------------| r_work (final): 0.3621 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7852 moved from start: 0.1859 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.043 4951 Z= 0.269 Angle : 0.748 10.678 6719 Z= 0.396 Chirality : 0.047 0.168 762 Planarity : 0.005 0.049 834 Dihedral : 5.310 26.618 665 Min Nonbonded Distance : 2.538 Molprobity Statistics. All-atom Clashscore : 12.17 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.68 % Favored : 94.32 % Rotamer: Outliers : 2.85 % Allowed : 14.42 % Favored : 82.73 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.12 (0.34), residues: 616 helix: 0.54 (0.28), residues: 365 sheet: 0.12 (0.53), residues: 100 loop : -1.54 (0.44), residues: 151 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.016 0.002 TRP A 369 HIS 0.006 0.001 HIS A 260 PHE 0.025 0.002 PHE B 129 TYR 0.014 0.002 TYR A 140 ARG 0.002 0.000 ARG A 86 Details of bonding type rmsd hydrogen bonds : bond 0.05642 ( 278) hydrogen bonds : angle 5.52299 ( 798) SS BOND : bond 0.00212 ( 2) SS BOND : angle 0.82750 ( 4) covalent geometry : bond 0.00623 ( 4949) covalent geometry : angle 0.74764 ( 6715) *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 102 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 15 poor density : 87 time to evaluate : 0.543 Fit side-chains REVERT: A 48 GLU cc_start: 0.7984 (mm-30) cc_final: 0.7438 (mm-30) REVERT: A 83 TYR cc_start: 0.6992 (m-10) cc_final: 0.6773 (m-10) REVERT: A 416 TYR cc_start: 0.7602 (t80) cc_final: 0.7237 (t80) REVERT: A 418 ASP cc_start: 0.7718 (t70) cc_final: 0.7355 (p0) outliers start: 15 outliers final: 8 residues processed: 95 average time/residue: 0.1267 time to fit residues: 16.6704 Evaluate side-chains 90 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 8 poor density : 82 time to evaluate : 0.519 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 179 THR Chi-restraints excluded: chain B residue 285 THR Chi-restraints excluded: chain A residue 133 THR Chi-restraints excluded: chain A residue 329 ILE Chi-restraints excluded: chain A residue 346 LEU Chi-restraints excluded: chain A residue 351 THR Chi-restraints excluded: chain A residue 371 SER Chi-restraints excluded: chain A residue 372 SER Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 61 random chunks: chunk 43 optimal weight: 0.9980 chunk 2 optimal weight: 0.8980 chunk 29 optimal weight: 2.9990 chunk 25 optimal weight: 0.9990 chunk 46 optimal weight: 1.9990 chunk 59 optimal weight: 0.6980 chunk 1 optimal weight: 0.7980 chunk 56 optimal weight: 1.9990 chunk 44 optimal weight: 0.6980 chunk 23 optimal weight: 2.9990 chunk 14 optimal weight: 4.9990 overall best weight: 0.8180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 89 GLN A 129 ASN A 187 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4189 r_free = 0.4189 target = 0.168569 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 45)----------------| | r_work = 0.3617 r_free = 0.3617 target = 0.123802 restraints weight = 7894.035| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 36)----------------| | r_work = 0.3668 r_free = 0.3668 target = 0.127330 restraints weight = 3752.975| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 42)----------------| | r_work = 0.3698 r_free = 0.3698 target = 0.129462 restraints weight = 2254.380| |-----------------------------------------------------------------------------| r_work (final): 0.3687 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7766 moved from start: 0.1951 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.038 4951 Z= 0.154 Angle : 0.634 8.783 6719 Z= 0.332 Chirality : 0.043 0.153 762 Planarity : 0.005 0.045 834 Dihedral : 5.025 22.824 665 Min Nonbonded Distance : 2.348 Molprobity Statistics. All-atom Clashscore : 10.00 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.22 % Favored : 95.78 % Rotamer: Outliers : 2.47 % Allowed : 19.92 % Favored : 77.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: 0.12 (0.35), residues: 616 helix: 0.79 (0.29), residues: 365 sheet: -0.01 (0.53), residues: 104 loop : -1.45 (0.46), residues: 147 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP B 236 HIS 0.003 0.001 HIS B 249 PHE 0.019 0.001 PHE B 129 TYR 0.009 0.001 TYR B 210 ARG 0.002 0.000 ARG B 238 Details of bonding type rmsd hydrogen bonds : bond 0.04603 ( 278) hydrogen bonds : angle 5.06059 ( 798) SS BOND : bond 0.00189 ( 2) SS BOND : angle 0.78727 ( 4) covalent geometry : bond 0.00348 ( 4949) covalent geometry : angle 0.63386 ( 6715) *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 104 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 13 poor density : 91 time to evaluate : 0.557 Fit side-chains REVERT: A 82 LYS cc_start: 0.8314 (ttpp) cc_final: 0.7361 (tmmt) REVERT: A 416 TYR cc_start: 0.7481 (t80) cc_final: 0.7152 (t80) REVERT: A 418 ASP cc_start: 0.7545 (t70) cc_final: 0.7340 (p0) REVERT: A 420 LYS cc_start: 0.9034 (tmmt) cc_final: 0.8690 (tmmt) REVERT: A 438 PHE cc_start: 0.8059 (OUTLIER) cc_final: 0.7651 (t80) outliers start: 13 outliers final: 8 residues processed: 98 average time/residue: 0.1428 time to fit residues: 18.9509 Evaluate side-chains 94 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 9 poor density : 85 time to evaluate : 0.504 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 143 PHE Chi-restraints excluded: chain B residue 201 VAL Chi-restraints excluded: chain A residue 256 THR Chi-restraints excluded: chain A residue 264 LEU Chi-restraints excluded: chain A residue 281 LEU Chi-restraints excluded: chain A residue 329 ILE Chi-restraints excluded: chain A residue 371 SER Chi-restraints excluded: chain A residue 372 SER Chi-restraints excluded: chain A residue 438 PHE Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 61 random chunks: chunk 33 optimal weight: 1.9990 chunk 17 optimal weight: 3.9990 chunk 1 optimal weight: 1.9990 chunk 5 optimal weight: 9.9990 chunk 16 optimal weight: 1.9990 chunk 53 optimal weight: 2.9990 chunk 36 optimal weight: 0.9990 chunk 41 optimal weight: 0.0470 chunk 58 optimal weight: 2.9990 chunk 38 optimal weight: 0.6980 chunk 24 optimal weight: 1.9990 overall best weight: 1.1484 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 89 GLN A 129 ASN A 291 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4171 r_free = 0.4171 target = 0.166316 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 38)----------------| | r_work = 0.3594 r_free = 0.3594 target = 0.121521 restraints weight = 7923.430| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 37)----------------| | r_work = 0.3644 r_free = 0.3644 target = 0.125002 restraints weight = 3851.141| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 46)----------------| | r_work = 0.3675 r_free = 0.3675 target = 0.127135 restraints weight = 2344.177| |-----------------------------------------------------------------------------| r_work (final): 0.3652 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7778 moved from start: 0.2216 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.038 4951 Z= 0.179 Angle : 0.658 8.078 6719 Z= 0.342 Chirality : 0.043 0.153 762 Planarity : 0.005 0.042 834 Dihedral : 5.059 24.691 665 Min Nonbonded Distance : 2.560 Molprobity Statistics. All-atom Clashscore : 10.93 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.84 % Favored : 94.16 % Rotamer: Outliers : 3.98 % Allowed : 21.44 % Favored : 74.57 % 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.35), residues: 616 helix: 0.76 (0.28), residues: 368 sheet: -0.03 (0.53), residues: 104 loop : -1.42 (0.48), residues: 144 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP A 369 HIS 0.003 0.001 HIS B 249 PHE 0.019 0.002 PHE B 129 TYR 0.010 0.001 TYR B 210 ARG 0.001 0.000 ARG B 151 Details of bonding type rmsd hydrogen bonds : bond 0.04712 ( 278) hydrogen bonds : angle 5.06080 ( 798) SS BOND : bond 0.00136 ( 2) SS BOND : angle 0.60914 ( 4) covalent geometry : bond 0.00412 ( 4949) covalent geometry : angle 0.65847 ( 6715) *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 108 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 21 poor density : 87 time to evaluate : 1.059 Fit side-chains REVERT: A 48 GLU cc_start: 0.7910 (mm-30) cc_final: 0.7340 (mm-30) REVERT: A 416 TYR cc_start: 0.7568 (t80) cc_final: 0.7156 (t80) REVERT: A 420 LYS cc_start: 0.9039 (tmmt) cc_final: 0.8311 (ttpp) REVERT: A 438 PHE cc_start: 0.8130 (OUTLIER) cc_final: 0.7681 (t80) outliers start: 21 outliers final: 12 residues processed: 100 average time/residue: 0.2171 time to fit residues: 30.4294 Evaluate side-chains 95 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 13 poor density : 82 time to evaluate : 0.763 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 116 VAL Chi-restraints excluded: chain B residue 143 PHE Chi-restraints excluded: chain B residue 201 VAL Chi-restraints excluded: chain B residue 212 ILE Chi-restraints excluded: chain A residue 189 CYS Chi-restraints excluded: chain A residue 256 THR Chi-restraints excluded: chain A residue 281 LEU Chi-restraints excluded: chain A residue 329 ILE Chi-restraints excluded: chain A residue 349 LEU Chi-restraints excluded: chain A residue 351 THR Chi-restraints excluded: chain A residue 371 SER Chi-restraints excluded: chain A residue 372 SER Chi-restraints excluded: chain A residue 438 PHE Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 61 random chunks: chunk 14 optimal weight: 2.9990 chunk 56 optimal weight: 0.5980 chunk 49 optimal weight: 2.9990 chunk 15 optimal weight: 4.9990 chunk 8 optimal weight: 10.0000 chunk 55 optimal weight: 0.5980 chunk 21 optimal weight: 1.9990 chunk 20 optimal weight: 0.9990 chunk 44 optimal weight: 0.9980 chunk 29 optimal weight: 2.9990 chunk 47 optimal weight: 2.9990 overall best weight: 1.0384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 89 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 113 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 129 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4180 r_free = 0.4180 target = 0.167456 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 42)----------------| | r_work = 0.3603 r_free = 0.3603 target = 0.122638 restraints weight = 7962.989| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 37)----------------| | r_work = 0.3654 r_free = 0.3654 target = 0.126154 restraints weight = 3848.497| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 43)----------------| | r_work = 0.3686 r_free = 0.3686 target = 0.128377 restraints weight = 2333.973| |-----------------------------------------------------------------------------| r_work (final): 0.3674 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7791 moved from start: 0.2383 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.041 4951 Z= 0.171 Angle : 0.653 8.191 6719 Z= 0.342 Chirality : 0.043 0.147 762 Planarity : 0.005 0.043 834 Dihedral : 5.067 23.515 665 Min Nonbonded Distance : 2.559 Molprobity Statistics. All-atom Clashscore : 10.41 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.68 % Favored : 94.32 % Rotamer: Outliers : 3.61 % Allowed : 23.72 % Favored : 72.68 % 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.15 (0.35), residues: 616 helix: 0.84 (0.29), residues: 365 sheet: -0.08 (0.54), residues: 104 loop : -1.45 (0.46), residues: 147 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP B 236 HIS 0.003 0.001 HIS B 249 PHE 0.019 0.002 PHE B 129 TYR 0.018 0.001 TYR A 287 ARG 0.001 0.000 ARG B 238 Details of bonding type rmsd hydrogen bonds : bond 0.04566 ( 278) hydrogen bonds : angle 4.97937 ( 798) SS BOND : bond 0.00321 ( 2) SS BOND : angle 1.59513 ( 4) covalent geometry : bond 0.00394 ( 4949) covalent geometry : angle 0.65192 ( 6715) *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 108 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 19 poor density : 89 time to evaluate : 0.613 Fit side-chains REVERT: A 404 LEU cc_start: 0.8299 (pp) cc_final: 0.8054 (mt) REVERT: A 420 LYS cc_start: 0.9066 (tmmt) cc_final: 0.8350 (ttpp) REVERT: A 438 PHE cc_start: 0.8149 (OUTLIER) cc_final: 0.7688 (t80) outliers start: 19 outliers final: 15 residues processed: 99 average time/residue: 0.1823 time to fit residues: 24.6549 Evaluate side-chains 96 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 16 poor density : 80 time to evaluate : 0.549 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 116 VAL Chi-restraints excluded: chain B residue 158 VAL Chi-restraints excluded: chain B residue 201 VAL Chi-restraints excluded: chain B residue 212 ILE Chi-restraints excluded: chain B residue 274 LEU Chi-restraints excluded: chain B residue 285 THR Chi-restraints excluded: chain A residue 81 ASN Chi-restraints excluded: chain A residue 156 VAL Chi-restraints excluded: chain A residue 189 CYS Chi-restraints excluded: chain A residue 256 THR Chi-restraints excluded: chain A residue 281 LEU Chi-restraints excluded: chain A residue 329 ILE Chi-restraints excluded: chain A residue 351 THR Chi-restraints excluded: chain A residue 371 SER Chi-restraints excluded: chain A residue 372 SER Chi-restraints excluded: chain A residue 438 PHE Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 61 random chunks: chunk 16 optimal weight: 0.9980 chunk 19 optimal weight: 1.9990 chunk 46 optimal weight: 2.9990 chunk 1 optimal weight: 3.9990 chunk 36 optimal weight: 1.9990 chunk 23 optimal weight: 3.9990 chunk 57 optimal weight: 0.8980 chunk 26 optimal weight: 0.5980 chunk 52 optimal weight: 0.9980 chunk 17 optimal weight: 0.8980 chunk 30 optimal weight: 2.9990 overall best weight: 0.8780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 89 GLN ** B 113 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 127 ASN A 129 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4190 r_free = 0.4190 target = 0.168058 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 48)----------------| | r_work = 0.3616 r_free = 0.3616 target = 0.123284 restraints weight = 7912.210| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 41)----------------| | r_work = 0.3668 r_free = 0.3668 target = 0.126855 restraints weight = 3760.113| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 44)----------------| | r_work = 0.3701 r_free = 0.3701 target = 0.129098 restraints weight = 2279.540| |-----------------------------------------------------------------------------| r_work (final): 0.3690 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7763 moved from start: 0.2456 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.037 4951 Z= 0.155 Angle : 0.635 7.394 6719 Z= 0.332 Chirality : 0.042 0.145 762 Planarity : 0.005 0.042 834 Dihedral : 4.991 23.815 665 Min Nonbonded Distance : 2.551 Molprobity Statistics. All-atom Clashscore : 10.11 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.36 % Favored : 94.64 % Rotamer: Outliers : 3.98 % Allowed : 24.48 % Favored : 71.54 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.17 (0.35), residues: 616 helix: 0.90 (0.29), residues: 365 sheet: -0.23 (0.53), residues: 105 loop : -1.45 (0.46), residues: 146 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP B 236 HIS 0.002 0.001 HIS B 249 PHE 0.019 0.002 PHE A 307 TYR 0.018 0.001 TYR A 287 ARG 0.002 0.000 ARG B 238 Details of bonding type rmsd hydrogen bonds : bond 0.04371 ( 278) hydrogen bonds : angle 4.90649 ( 798) SS BOND : bond 0.00238 ( 2) SS BOND : angle 1.14931 ( 4) covalent geometry : bond 0.00356 ( 4949) covalent geometry : angle 0.63451 ( 6715) *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 107 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 21 poor density : 86 time to evaluate : 0.604 Fit side-chains REVERT: A 48 GLU cc_start: 0.7723 (mm-30) cc_final: 0.7196 (mm-30) REVERT: A 65 MET cc_start: 0.8463 (tpp) cc_final: 0.8221 (ttm) REVERT: A 420 LYS cc_start: 0.9053 (tmmt) cc_final: 0.8341 (ttpp) REVERT: A 438 PHE cc_start: 0.8155 (OUTLIER) cc_final: 0.7680 (t80) outliers start: 21 outliers final: 17 residues processed: 99 average time/residue: 0.1485 time to fit residues: 20.0816 Evaluate side-chains 101 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 18 poor density : 83 time to evaluate : 0.535 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 89 GLN Chi-restraints excluded: chain B residue 116 VAL Chi-restraints excluded: chain B residue 158 VAL Chi-restraints excluded: chain B residue 201 VAL Chi-restraints excluded: chain B residue 212 ILE Chi-restraints excluded: chain B residue 274 LEU Chi-restraints excluded: chain B residue 285 THR Chi-restraints excluded: chain A residue 156 VAL Chi-restraints excluded: chain A residue 189 CYS Chi-restraints excluded: chain A residue 281 LEU Chi-restraints excluded: chain A residue 329 ILE Chi-restraints excluded: chain A residue 346 LEU Chi-restraints excluded: chain A residue 349 LEU Chi-restraints excluded: chain A residue 351 THR Chi-restraints excluded: chain A residue 371 SER Chi-restraints excluded: chain A residue 372 SER Chi-restraints excluded: chain A residue 374 LEU Chi-restraints excluded: chain A residue 438 PHE Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 61 random chunks: chunk 57 optimal weight: 0.0570 chunk 25 optimal weight: 0.6980 chunk 21 optimal weight: 0.9990 chunk 19 optimal weight: 0.9980 chunk 5 optimal weight: 9.9990 chunk 59 optimal weight: 2.9990 chunk 31 optimal weight: 1.9990 chunk 56 optimal weight: 0.7980 chunk 3 optimal weight: 4.9990 chunk 16 optimal weight: 4.9990 chunk 11 optimal weight: 0.9990 overall best weight: 0.7100 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 89 GLN ** B 113 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 129 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4204 r_free = 0.4204 target = 0.169433 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 51)----------------| | r_work = 0.3633 r_free = 0.3633 target = 0.124585 restraints weight = 7751.154| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 39)----------------| | r_work = 0.3684 r_free = 0.3684 target = 0.128171 restraints weight = 3732.890| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 52)----------------| | r_work = 0.3716 r_free = 0.3716 target = 0.130390 restraints weight = 2267.218| |-----------------------------------------------------------------------------| r_work (final): 0.3698 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7743 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.003 0.037 4951 Z= 0.147 Angle : 0.631 8.996 6719 Z= 0.328 Chirality : 0.042 0.143 762 Planarity : 0.005 0.042 834 Dihedral : 4.946 22.643 665 Min Nonbonded Distance : 2.553 Molprobity Statistics. All-atom Clashscore : 9.18 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.87 % Favored : 95.13 % Rotamer: Outliers : 3.98 % Allowed : 25.62 % Favored : 70.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.20 (0.35), residues: 616 helix: 0.94 (0.29), residues: 365 sheet: -0.30 (0.53), residues: 105 loop : -1.40 (0.47), residues: 146 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP B 236 HIS 0.002 0.000 HIS B 249 PHE 0.017 0.001 PHE B 129 TYR 0.018 0.001 TYR A 287 ARG 0.002 0.000 ARG B 238 Details of bonding type rmsd hydrogen bonds : bond 0.04232 ( 278) hydrogen bonds : angle 4.84563 ( 798) SS BOND : bond 0.00392 ( 2) SS BOND : angle 1.45876 ( 4) covalent geometry : bond 0.00337 ( 4949) covalent geometry : angle 0.62990 ( 6715) *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1232 Ramachandran restraints generated. 616 Oldfield, 0 Emsley, 616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 108 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 21 poor density : 87 time to evaluate : 0.535 Fit side-chains REVERT: A 418 ASP cc_start: 0.7406 (p0) cc_final: 0.6721 (m-30) REVERT: A 420 LYS cc_start: 0.9049 (tmmt) cc_final: 0.8412 (ttpp) REVERT: A 438 PHE cc_start: 0.8154 (OUTLIER) cc_final: 0.7685 (t80) outliers start: 21 outliers final: 16 residues processed: 100 average time/residue: 0.1375 time to fit residues: 18.5728 Evaluate side-chains 99 residues out of total 528 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 17 poor density : 82 time to evaluate : 0.534 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 89 GLN Chi-restraints excluded: chain B residue 116 VAL Chi-restraints excluded: chain B residue 158 VAL Chi-restraints excluded: chain B residue 201 VAL Chi-restraints excluded: chain B residue 212 ILE Chi-restraints excluded: chain B residue 254 VAL Chi-restraints excluded: chain B residue 274 LEU Chi-restraints excluded: chain B residue 285 THR Chi-restraints excluded: chain A residue 81 ASN Chi-restraints excluded: chain A residue 156 VAL Chi-restraints excluded: chain A residue 189 CYS Chi-restraints excluded: chain A residue 264 LEU Chi-restraints excluded: chain A residue 281 LEU Chi-restraints excluded: chain A residue 329 ILE Chi-restraints excluded: chain A residue 346 LEU Chi-restraints excluded: chain A residue 372 SER Chi-restraints excluded: chain A residue 438 PHE Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 8, in run_program(real_space_refine.Program) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/iotbx/cli_parser.py", line 994, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/phenix/phenix/programs/real_space_refine.py", line 210, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 306, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 702, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 767, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 51, in __init__ gradients_method = gradients_method) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/mmtbx/refinement/real_space/weight.py", line 87, in __init__ rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 368, in refine rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 112, in refine states_collector = self.states_accumulator) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 174, in __init__ exception_handling_params=lbfgs_exception_handling_params) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 131, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 245, in compute_functional_and_gradients compute_gradients=True) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1525, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1427, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5669/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1308, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 53.7234 > 50: distance: 24 - 30: 32.556 distance: 30 - 31: 53.175 distance: 32 - 33: 3.670 distance: 32 - 34: 18.442 distance: 33 - 59: 12.652 distance: 34 - 35: 4.425 distance: 35 - 36: 4.214 distance: 36 - 37: 7.321 distance: 36 - 40: 4.706 distance: 37 - 66: 13.241 distance: 38 - 39: 6.628 distance: 40 - 41: 4.810 distance: 41 - 44: 5.880 distance: 42 - 43: 3.938 distance: 42 - 51: 3.086 distance: 43 - 74: 16.635 distance: 44 - 45: 6.202 distance: 47 - 48: 3.834 distance: 51 - 52: 5.718 distance: 52 - 53: 4.930 distance: 52 - 55: 4.849 distance: 53 - 59: 4.553 distance: 54 - 82: 9.873 distance: 55 - 57: 10.642 distance: 56 - 58: 7.701 distance: 59 - 60: 4.837 distance: 60 - 61: 3.713 distance: 60 - 63: 5.136 distance: 61 - 62: 3.446 distance: 62 - 89: 4.854 distance: 63 - 65: 5.989 distance: 66 - 67: 5.445 distance: 67 - 68: 3.893 distance: 68 - 69: 3.415 distance: 69 - 93: 4.591 distance: 70 - 71: 4.263 distance: 71 - 72: 5.886 distance: 72 - 73: 6.963 distance: 74 - 75: 3.791 distance: 75 - 76: 5.184 distance: 75 - 78: 5.558 distance: 76 - 82: 7.744 distance: 77 - 97: 7.186 distance: 78 - 79: 5.633 distance: 78 - 80: 6.729 distance: 79 - 81: 4.268 distance: 82 - 83: 4.715 distance: 83 - 84: 7.446 distance: 83 - 86: 9.573 distance: 85 - 103: 15.942 distance: 86 - 87: 8.290 distance: 86 - 88: 3.284 distance: 89 - 90: 4.619 distance: 90 - 91: 6.508 distance: 91 - 92: 5.146 distance: 91 - 93: 6.581 distance: 92 - 111: 13.347 distance: 93 - 94: 4.968 distance: 94 - 95: 4.218 distance: 95 - 97: 5.695 distance: 96 - 116: 5.022 distance: 164 - 175: 3.027