Starting phenix.real_space_refine on Mon Jan 13 20:50:37 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/8s5u_19750/01_2025/8s5u_19750.cif Found real_map, /net/cci-nas-00/data/ceres_data/8s5u_19750/01_2025/8s5u_19750.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.28 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/8s5u_19750/01_2025/8s5u_19750.map" default_real_map = "/net/cci-nas-00/data/ceres_data/8s5u_19750/01_2025/8s5u_19750.map" model { file = "/net/cci-nas-00/data/ceres_data/8s5u_19750/01_2025/8s5u_19750.cif" } default_model = "/net/cci-nas-00/data/ceres_data/8s5u_19750/01_2025/8s5u_19750.cif" } resolution = 3.28 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.050 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 20 5.16 5 C 2652 2.51 5 N 620 2.21 5 O 706 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/chem_data/mon_lib" Total number of atoms: 3998 Number of models: 1 Model: "" Number of chains: 3 Chain: "B" Number of atoms: 982 Number of conformers: 1 Conformer: "" Number of residues, atoms: 128, 982 Classifications: {'peptide': 128} Link IDs: {'PTRANS': 3, 'TRANS': 124} Chain: "A" Number of atoms: 2984 Number of conformers: 1 Conformer: "" Number of residues, atoms: 374, 2984 Classifications: {'peptide': 374} Link IDs: {'PTRANS': 13, 'TRANS': 360} Chain breaks: 1 Chain: "A" Number of atoms: 32 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 32 Unusual residues: {'NAG': 1, 'VIB': 1} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen dihedrals: 3 Unresolved non-hydrogen chiralities: 1 Time building chain proxies: 3.97, per 1000 atoms: 0.99 Number of scatterers: 3998 At special positions: 0 Unit cell: (107.1, 72.25, 64.6, 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 20 16.00 O 706 8.00 N 620 7.00 C 2652 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=1, symmetry=0 Simple disulfide: pdb=" SG CYS B 22 " - pdb=" SG CYS B 96 " 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 Links applied NAG-ASN " NAG A 601 " - " ASN A 45 " Time building additional restraints: 0.90 Conformation dependent library (CDL) restraints added in 456.3 milliseconds 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 940 Finding SS restraints... Secondary structure from input PDB file: 22 helices and 3 sheets defined 66.9% alpha, 12.2% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.43 Creating SS restraints... Processing helix chain 'B' and resid 87 through 91 removed outlier: 3.728A pdb=" N THR B 91 " --> pdb=" O PRO B 88 " (cutoff:3.500A) Processing helix chain 'B' and resid 110 through 113 removed outlier: 3.503A pdb=" N THR B 113 " --> pdb=" O ASN B 110 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 110 through 113' Processing helix chain 'A' and resid 12 through 28 removed outlier: 3.721A pdb=" N THR A 16 " --> pdb=" O TRP A 12 " (cutoff:3.500A) Processing helix chain 'A' and resid 31 through 33 No H-bonds generated for 'chain 'A' and resid 31 through 33' Processing helix chain 'A' and resid 34 through 41 Processing helix chain 'A' and resid 47 through 54 Processing helix chain 'A' and resid 55 through 74 removed outlier: 3.884A pdb=" N TRP A 59 " --> pdb=" O ILE A 55 " (cutoff:3.500A) removed outlier: 3.962A pdb=" N THR A 60 " --> pdb=" O PHE A 56 " (cutoff:3.500A) removed outlier: 3.525A pdb=" N TYR A 61 " --> pdb=" O PRO A 57 " (cutoff:3.500A) removed outlier: 3.948A pdb=" N LEU A 68 " --> pdb=" O LEU A 64 " (cutoff:3.500A) Proline residue: A 69 - end of helix Processing helix chain 'A' and resid 78 through 99 removed outlier: 4.364A pdb=" N VAL A 82 " --> pdb=" O ARG A 78 " (cutoff:3.500A) Processing helix chain 'A' and resid 103 through 119 Processing helix chain 'A' and resid 120 through 127 removed outlier: 3.504A pdb=" N TYR A 124 " --> pdb=" O GLU A 120 " (cutoff:3.500A) removed outlier: 4.482A pdb=" N ILE A 127 " --> pdb=" O TYR A 123 " (cutoff:3.500A) Processing helix chain 'A' and resid 128 through 130 No H-bonds generated for 'chain 'A' and resid 128 through 130' Processing helix chain 'A' and resid 132 through 134 No H-bonds generated for 'chain 'A' and resid 132 through 134' Processing helix chain 'A' and resid 135 through 166 removed outlier: 3.540A pdb=" N ARG A 144 " --> pdb=" O SER A 140 " (cutoff:3.500A) removed outlier: 3.742A pdb=" N THR A 147 " --> pdb=" O CYS A 143 " (cutoff:3.500A) Processing helix chain 'A' and resid 168 through 188 removed outlier: 3.805A pdb=" N LEU A 188 " --> pdb=" O PHE A 184 " (cutoff:3.500A) Processing helix chain 'A' and resid 271 through 274 removed outlier: 4.120A pdb=" N SER A 274 " --> pdb=" O GLU A 271 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 271 through 274' Processing helix chain 'A' and resid 275 through 296 removed outlier: 3.505A pdb=" N GLN A 294 " --> pdb=" O ALA A 290 " (cutoff:3.500A) Processing helix chain 'A' and resid 298 through 307 Processing helix chain 'A' and resid 311 through 335 removed outlier: 3.571A pdb=" N TYR A 315 " --> pdb=" O ASP A 311 " (cutoff:3.500A) removed outlier: 3.913A pdb=" N GLY A 335 " --> pdb=" O ALA A 331 " (cutoff:3.500A) Processing helix chain 'A' and resid 340 through 366 removed outlier: 3.757A pdb=" N LEU A 344 " --> pdb=" O ASN A 340 " (cutoff:3.500A) removed outlier: 5.775A pdb=" N GLU A 346 " --> pdb=" O ASP A 342 " (cutoff:3.500A) removed outlier: 4.691A pdb=" N LEU A 347 " --> pdb=" O LEU A 343 " (cutoff:3.500A) removed outlier: 3.557A pdb=" N VAL A 355 " --> pdb=" O VAL A 351 " (cutoff:3.500A) removed outlier: 3.599A pdb=" N LEU A 360 " --> pdb=" O ASN A 356 " (cutoff:3.500A) Processing helix chain 'A' and resid 368 through 397 removed outlier: 3.529A pdb=" N CYS A 372 " --> pdb=" O ASN A 368 " (cutoff:3.500A) removed outlier: 3.587A pdb=" N THR A 388 " --> pdb=" O MET A 384 " (cutoff:3.500A) removed outlier: 3.518A pdb=" N ILE A 389 " --> pdb=" O LEU A 385 " (cutoff:3.500A) Processing helix chain 'A' and resid 402 through 426 Processing helix chain 'A' and resid 434 through 455 removed outlier: 4.267A pdb=" N PHE A 439 " --> pdb=" O VAL A 435 " (cutoff:3.500A) removed outlier: 5.528A pdb=" N LEU A 440 " --> pdb=" O SER A 436 " (cutoff:3.500A) removed outlier: 3.776A pdb=" N VAL A 441 " --> pdb=" O ILE A 437 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'B' and resid 3 through 7 removed outlier: 3.896A pdb=" N VAL B 5 " --> pdb=" O ALA B 23 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'B' and resid 11 through 13 removed outlier: 6.371A pdb=" N VAL B 12 " --> pdb=" O SER B 127 " (cutoff:3.500A) removed outlier: 3.524A pdb=" N MET B 34 " --> pdb=" O ILE B 51 " (cutoff:3.500A) removed outlier: 5.781A pdb=" N ARG B 38 " --> pdb=" O PHE B 47 " (cutoff:3.500A) removed outlier: 5.644A pdb=" N PHE B 47 " --> pdb=" O ARG B 38 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'B' and resid 11 through 13 removed outlier: 6.371A pdb=" N VAL B 12 " --> pdb=" O SER B 127 " (cutoff:3.500A) 259 hydrogen bonds defined for protein. 753 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.05 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.20 - 1.32: 630 1.32 - 1.44: 1220 1.44 - 1.56: 2225 1.56 - 1.68: 1 1.68 - 1.81: 32 Bond restraints: 4108 Sorted by residual: bond pdb=" C5 VIB A 602 " pdb=" S1 VIB A 602 " ideal model delta sigma weight residual 1.721 1.508 0.213 2.00e-02 2.50e+03 1.13e+02 bond pdb=" C4 VIB A 602 " pdb=" N3 VIB A 602 " ideal model delta sigma weight residual 1.385 1.572 -0.187 2.00e-02 2.50e+03 8.79e+01 bond pdb=" C2A VIB A 602 " pdb=" N3A VIB A 602 " ideal model delta sigma weight residual 1.323 1.478 -0.155 2.00e-02 2.50e+03 6.02e+01 bond pdb=" C2A VIB A 602 " pdb=" N1A VIB A 602 " ideal model delta sigma weight residual 1.331 1.481 -0.150 2.00e-02 2.50e+03 5.63e+01 bond pdb=" C2 VIB A 602 " pdb=" S1 VIB A 602 " ideal model delta sigma weight residual 1.670 1.539 0.131 2.00e-02 2.50e+03 4.30e+01 ... (remaining 4103 not shown) Histogram of bond angle deviations from ideal: 0.00 - 3.33: 5562 3.33 - 6.67: 40 6.67 - 10.00: 1 10.00 - 13.34: 0 13.34 - 16.67: 1 Bond angle restraints: 5604 Sorted by residual: angle pdb=" N TYR A 79 " pdb=" CA TYR A 79 " pdb=" C TYR A 79 " ideal model delta sigma weight residual 112.34 120.30 -7.96 1.30e+00 5.92e-01 3.75e+01 angle pdb=" C2 VIB A 602 " pdb=" S1 VIB A 602 " pdb=" C5 VIB A 602 " ideal model delta sigma weight residual 91.25 107.92 -16.67 3.00e+00 1.11e-01 3.09e+01 angle pdb=" N LEU A 95 " pdb=" CA LEU A 95 " pdb=" C LEU A 95 " ideal model delta sigma weight residual 113.23 107.31 5.92 1.24e+00 6.50e-01 2.28e+01 angle pdb=" N SER A 353 " pdb=" CA SER A 353 " pdb=" C SER A 353 " ideal model delta sigma weight residual 112.90 107.37 5.53 1.31e+00 5.83e-01 1.78e+01 angle pdb=" N LEU A 97 " pdb=" CA LEU A 97 " pdb=" C LEU A 97 " ideal model delta sigma weight residual 112.23 106.92 5.31 1.26e+00 6.30e-01 1.77e+01 ... (remaining 5599 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 15.89: 2015 15.89 - 31.78: 247 31.78 - 47.68: 73 47.68 - 63.57: 14 63.57 - 79.46: 6 Dihedral angle restraints: 2355 sinusoidal: 880 harmonic: 1475 Sorted by residual: dihedral pdb=" CA TYR A 126 " pdb=" C TYR A 126 " pdb=" N ILE A 127 " pdb=" CA ILE A 127 " ideal model delta harmonic sigma weight residual -180.00 -163.32 -16.68 0 5.00e+00 4.00e-02 1.11e+01 dihedral pdb=" CG ARG A 402 " pdb=" CD ARG A 402 " pdb=" NE ARG A 402 " pdb=" CZ ARG A 402 " ideal model delta sinusoidal sigma weight residual 90.00 133.07 -43.07 2 1.50e+01 4.44e-03 9.95e+00 dihedral pdb=" CA LYS A 306 " pdb=" CB LYS A 306 " pdb=" CG LYS A 306 " pdb=" CD LYS A 306 " ideal model delta sinusoidal sigma weight residual -180.00 -121.08 -58.92 3 1.50e+01 4.44e-03 9.47e+00 ... (remaining 2352 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.044: 478 0.044 - 0.088: 133 0.088 - 0.133: 23 0.133 - 0.177: 5 0.177 - 0.221: 3 Chirality restraints: 642 Sorted by residual: chirality pdb=" CA TYR A 273 " pdb=" N TYR A 273 " pdb=" C TYR A 273 " pdb=" CB TYR A 273 " both_signs ideal model delta sigma weight residual False 2.51 2.29 0.22 2.00e-01 2.50e+01 1.22e+00 chirality pdb=" CA TYR A 79 " pdb=" N TYR A 79 " pdb=" C TYR A 79 " pdb=" CB TYR A 79 " both_signs ideal model delta sigma weight residual False 2.51 2.31 0.20 2.00e-01 2.50e+01 1.03e+00 chirality pdb=" C1 NAG A 601 " pdb=" ND2 ASN A 45 " pdb=" C2 NAG A 601 " pdb=" O5 NAG A 601 " both_signs ideal model delta sigma weight residual False -2.40 -2.58 0.18 2.00e-01 2.50e+01 8.04e-01 ... (remaining 639 not shown) Planarity restraints: 682 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C2 NAG A 601 " -0.359 2.00e-02 2.50e+03 3.06e-01 1.17e+03 pdb=" C7 NAG A 601 " 0.091 2.00e-02 2.50e+03 pdb=" C8 NAG A 601 " -0.163 2.00e-02 2.50e+03 pdb=" N2 NAG A 601 " 0.541 2.00e-02 2.50e+03 pdb=" O7 NAG A 601 " -0.110 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C LYS A 80 " -0.029 5.00e-02 4.00e+02 4.35e-02 3.02e+00 pdb=" N PRO A 81 " 0.075 5.00e-02 4.00e+02 pdb=" CA PRO A 81 " -0.022 5.00e-02 4.00e+02 pdb=" CD PRO A 81 " -0.024 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CA PHE A 352 " 0.008 2.00e-02 2.50e+03 1.54e-02 2.36e+00 pdb=" C PHE A 352 " -0.027 2.00e-02 2.50e+03 pdb=" O PHE A 352 " 0.010 2.00e-02 2.50e+03 pdb=" N SER A 353 " 0.009 2.00e-02 2.50e+03 ... (remaining 679 not shown) Histogram of nonbonded interaction distances: 2.34 - 2.85: 1493 2.85 - 3.36: 3830 3.36 - 3.88: 6459 3.88 - 4.39: 7731 4.39 - 4.90: 13382 Nonbonded interactions: 32895 Sorted by model distance: nonbonded pdb=" OG1 THR B 106 " pdb=" OD2 ASP A 304 " model vdw 2.338 3.040 nonbonded pdb=" O LEU A 377 " pdb=" OG SER A 381 " model vdw 2.394 3.040 nonbonded pdb=" OG SER A 309 " pdb=" OD1 ASP A 311 " model vdw 2.402 3.040 nonbonded pdb=" O ARG B 31 " pdb=" O PRO B 102 " model vdw 2.404 3.040 nonbonded pdb=" O SER B 117 " pdb=" OG SER B 117 " model vdw 2.408 3.040 ... (remaining 32890 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 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: 3.010 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.180 Check model and map are aligned: 0.030 Set scattering table: 0.040 Process input model: 14.880 Find NCS groups from input model: 0.030 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:2.840 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 21.040 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8036 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.009 0.213 4108 Z= 0.547 Angle : 0.708 16.670 5604 Z= 0.414 Chirality : 0.043 0.221 642 Planarity : 0.012 0.306 681 Dihedral : 15.772 79.460 1412 Min Nonbonded Distance : 2.338 Molprobity Statistics. All-atom Clashscore : 5.38 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.63 % Favored : 96.37 % Rotamer: Outliers : 0.47 % Allowed : 17.73 % Favored : 81.80 % 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.65 (0.38), residues: 496 helix: 0.13 (0.30), residues: 290 sheet: -1.03 (0.78), residues: 46 loop : -1.07 (0.49), residues: 160 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP A 94 HIS 0.002 0.001 HIS A 364 PHE 0.026 0.001 PHE A 71 TYR 0.008 0.001 TYR A 171 ARG 0.002 0.000 ARG A 456 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 71 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 2 poor density : 69 time to evaluate : 0.426 Fit side-chains REVERT: A 75 ASP cc_start: 0.8172 (m-30) cc_final: 0.7907 (m-30) REVERT: A 421 MET cc_start: 0.7598 (ttp) cc_final: 0.7123 (ttp) outliers start: 2 outliers final: 0 residues processed: 69 average time/residue: 0.1670 time to fit residues: 14.5863 Evaluate side-chains 60 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 60 time to evaluate : 0.447 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 41 optimal weight: 1.9990 chunk 37 optimal weight: 0.8980 chunk 20 optimal weight: 0.0970 chunk 12 optimal weight: 0.9980 chunk 25 optimal weight: 0.5980 chunk 19 optimal weight: 0.0970 chunk 38 optimal weight: 0.8980 chunk 14 optimal weight: 0.7980 chunk 23 optimal weight: 0.9980 chunk 28 optimal weight: 0.5980 chunk 44 optimal weight: 0.0970 overall best weight: 0.2974 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 110 ASN A 340 ASN A 356 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3357 r_free = 0.3357 target = 0.128683 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 48)----------------| | r_work = 0.3201 r_free = 0.3201 target = 0.115968 restraints weight = 5173.820| |-----------------------------------------------------------------------------| r_work (start): 0.3197 rms_B_bonded: 1.27 r_work: 0.3103 rms_B_bonded: 1.85 restraints_weight: 0.5000 r_work: 0.2994 rms_B_bonded: 3.26 restraints_weight: 0.2500 r_work (final): 0.2994 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8391 moved from start: 0.1501 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.046 4108 Z= 0.198 Angle : 0.536 5.560 5604 Z= 0.291 Chirality : 0.041 0.140 642 Planarity : 0.005 0.063 681 Dihedral : 5.981 59.072 575 Min Nonbonded Distance : 2.585 Molprobity Statistics. All-atom Clashscore : 2.63 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.42 % Favored : 97.58 % Rotamer: Outliers : 1.89 % Allowed : 18.68 % Favored : 79.43 % 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.61 (0.38), residues: 496 helix: 1.18 (0.30), residues: 300 sheet: -0.70 (0.75), residues: 46 loop : -0.66 (0.52), residues: 150 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.019 0.001 TRP A 94 HIS 0.001 0.000 HIS A 364 PHE 0.018 0.001 PHE A 90 TYR 0.012 0.001 TYR A 445 ARG 0.002 0.000 ARG B 72 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 73 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 8 poor density : 65 time to evaluate : 0.485 Fit side-chains outliers start: 8 outliers final: 3 residues processed: 70 average time/residue: 0.1214 time to fit residues: 11.5524 Evaluate side-chains 55 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 52 time to evaluate : 0.454 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 81 LEU Chi-restraints excluded: chain A residue 356 ASN Chi-restraints excluded: chain A residue 455 MET Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 30 optimal weight: 0.0770 chunk 47 optimal weight: 0.9990 chunk 11 optimal weight: 2.9990 chunk 0 optimal weight: 1.9990 chunk 41 optimal weight: 1.9990 chunk 29 optimal weight: 0.9980 chunk 6 optimal weight: 0.9990 chunk 12 optimal weight: 0.7980 chunk 46 optimal weight: 0.8980 chunk 28 optimal weight: 0.9980 chunk 27 optimal weight: 0.5980 overall best weight: 0.6738 ------------------------------------------------------------------------------- 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.3316 r_free = 0.3316 target = 0.125725 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 42)----------------| | r_work = 0.3159 r_free = 0.3159 target = 0.113159 restraints weight = 5166.421| |-----------------------------------------------------------------------------| r_work (start): 0.3154 rms_B_bonded: 1.26 r_work: 0.3064 rms_B_bonded: 1.80 restraints_weight: 0.5000 r_work: 0.2956 rms_B_bonded: 3.15 restraints_weight: 0.2500 r_work (final): 0.2956 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8441 moved from start: 0.1709 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.043 4108 Z= 0.313 Angle : 0.551 5.989 5604 Z= 0.299 Chirality : 0.043 0.148 642 Planarity : 0.005 0.063 681 Dihedral : 5.617 58.392 575 Min Nonbonded Distance : 2.494 Molprobity Statistics. All-atom Clashscore : 3.13 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.42 % Favored : 97.58 % Rotamer: Outliers : 2.60 % Allowed : 20.57 % Favored : 76.83 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.93 (0.38), residues: 496 helix: 1.43 (0.30), residues: 300 sheet: -0.55 (0.73), residues: 46 loop : -0.56 (0.52), residues: 150 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.020 0.002 TRP A 94 HIS 0.001 0.001 HIS A 135 PHE 0.025 0.002 PHE A 71 TYR 0.014 0.002 TYR A 14 ARG 0.002 0.000 ARG A 402 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 69 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 11 poor density : 58 time to evaluate : 0.425 Fit side-chains REVERT: A 75 ASP cc_start: 0.8208 (m-30) cc_final: 0.7751 (m-30) REVERT: A 455 MET cc_start: 0.7022 (OUTLIER) cc_final: 0.6230 (ttp) outliers start: 11 outliers final: 5 residues processed: 65 average time/residue: 0.1245 time to fit residues: 10.9273 Evaluate side-chains 62 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 6 poor density : 56 time to evaluate : 0.427 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 69 THR Chi-restraints excluded: chain B residue 81 LEU Chi-restraints excluded: chain A residue 145 SER Chi-restraints excluded: chain A residue 283 LEU Chi-restraints excluded: chain A residue 294 GLN Chi-restraints excluded: chain A residue 455 MET Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 8 optimal weight: 0.9990 chunk 20 optimal weight: 0.8980 chunk 45 optimal weight: 0.6980 chunk 0 optimal weight: 1.9990 chunk 7 optimal weight: 0.9990 chunk 32 optimal weight: 0.9990 chunk 35 optimal weight: 1.9990 chunk 17 optimal weight: 0.0770 chunk 3 optimal weight: 0.4980 chunk 27 optimal weight: 1.9990 chunk 48 optimal weight: 0.9990 overall best weight: 0.6340 ------------------------------------------------------------------------------- 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.3317 r_free = 0.3317 target = 0.125836 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 42)----------------| | r_work = 0.3161 r_free = 0.3161 target = 0.113309 restraints weight = 5128.866| |-----------------------------------------------------------------------------| r_work (start): 0.3158 rms_B_bonded: 1.25 r_work: 0.3067 rms_B_bonded: 1.80 restraints_weight: 0.5000 r_work: 0.2963 rms_B_bonded: 3.11 restraints_weight: 0.2500 r_work (final): 0.2963 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8438 moved from start: 0.1898 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.039 4108 Z= 0.295 Angle : 0.537 5.990 5604 Z= 0.292 Chirality : 0.043 0.146 642 Planarity : 0.005 0.061 681 Dihedral : 5.417 53.648 575 Min Nonbonded Distance : 2.537 Molprobity Statistics. All-atom Clashscore : 4.50 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.82 % Favored : 97.18 % Rotamer: Outliers : 2.13 % Allowed : 21.04 % Favored : 76.83 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.12 (0.38), residues: 496 helix: 1.54 (0.30), residues: 300 sheet: -0.44 (0.73), residues: 46 loop : -0.39 (0.53), residues: 150 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.020 0.002 TRP A 94 HIS 0.001 0.001 HIS A 135 PHE 0.024 0.002 PHE A 71 TYR 0.023 0.002 TYR A 14 ARG 0.002 0.000 ARG A 429 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 69 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 9 poor density : 60 time to evaluate : 0.460 Fit side-chains REVERT: A 455 MET cc_start: 0.6995 (OUTLIER) cc_final: 0.6195 (ttp) outliers start: 9 outliers final: 6 residues processed: 67 average time/residue: 0.1201 time to fit residues: 11.0351 Evaluate side-chains 65 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 7 poor density : 58 time to evaluate : 0.420 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 12 VAL Chi-restraints excluded: chain B residue 69 THR Chi-restraints excluded: chain B residue 71 SER Chi-restraints excluded: chain B residue 81 LEU Chi-restraints excluded: chain A residue 145 SER Chi-restraints excluded: chain A residue 283 LEU Chi-restraints excluded: chain A residue 455 MET Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 45 optimal weight: 0.5980 chunk 32 optimal weight: 0.9980 chunk 22 optimal weight: 0.9980 chunk 3 optimal weight: 0.7980 chunk 1 optimal weight: 2.9990 chunk 9 optimal weight: 0.9990 chunk 10 optimal weight: 0.6980 chunk 23 optimal weight: 0.9980 chunk 38 optimal weight: 0.9990 chunk 5 optimal weight: 0.2980 chunk 34 optimal weight: 0.6980 overall best weight: 0.6180 ------------------------------------------------------------------------------- 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.3322 r_free = 0.3322 target = 0.126128 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 46)----------------| | r_work = 0.3163 r_free = 0.3163 target = 0.113383 restraints weight = 5231.398| |-----------------------------------------------------------------------------| r_work (start): 0.3157 rms_B_bonded: 1.29 r_work: 0.3065 rms_B_bonded: 1.84 restraints_weight: 0.5000 r_work: 0.2957 rms_B_bonded: 3.23 restraints_weight: 0.2500 r_work (final): 0.2957 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8444 moved from start: 0.2009 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.036 4108 Z= 0.291 Angle : 0.538 6.012 5604 Z= 0.292 Chirality : 0.043 0.147 642 Planarity : 0.004 0.057 681 Dihedral : 5.318 49.279 575 Min Nonbonded Distance : 2.559 Molprobity Statistics. All-atom Clashscore : 4.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.62 % Favored : 97.38 % Rotamer: Outliers : 3.55 % Allowed : 20.80 % Favored : 75.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: 1.17 (0.38), residues: 496 helix: 1.57 (0.30), residues: 300 sheet: -0.43 (0.73), residues: 46 loop : -0.36 (0.54), residues: 150 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.019 0.002 TRP A 94 HIS 0.001 0.001 HIS A 135 PHE 0.023 0.002 PHE A 71 TYR 0.013 0.001 TYR A 445 ARG 0.001 0.000 ARG A 138 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 71 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 15 poor density : 56 time to evaluate : 0.471 Fit side-chains REVERT: A 363 MET cc_start: 0.9006 (mtm) cc_final: 0.8787 (ttp) REVERT: A 455 MET cc_start: 0.6975 (OUTLIER) cc_final: 0.6163 (ttp) outliers start: 15 outliers final: 9 residues processed: 68 average time/residue: 0.1326 time to fit residues: 12.0031 Evaluate side-chains 67 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 10 poor density : 57 time to evaluate : 0.416 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 12 VAL Chi-restraints excluded: chain B residue 51 ILE Chi-restraints excluded: chain B residue 69 THR Chi-restraints excluded: chain B residue 81 LEU Chi-restraints excluded: chain A residue 145 SER Chi-restraints excluded: chain A residue 283 LEU Chi-restraints excluded: chain A residue 294 GLN Chi-restraints excluded: chain A residue 366 THR Chi-restraints excluded: chain A residue 454 LEU Chi-restraints excluded: chain A residue 455 MET Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 31 optimal weight: 0.9990 chunk 48 optimal weight: 0.9990 chunk 39 optimal weight: 0.9980 chunk 5 optimal weight: 0.4980 chunk 45 optimal weight: 0.9990 chunk 14 optimal weight: 0.7980 chunk 19 optimal weight: 0.0970 chunk 47 optimal weight: 0.6980 chunk 28 optimal weight: 0.9980 chunk 2 optimal weight: 0.9980 chunk 42 optimal weight: 0.6980 overall best weight: 0.5578 ------------------------------------------------------------------------------- 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.3328 r_free = 0.3328 target = 0.126522 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 36)----------------| | r_work = 0.3167 r_free = 0.3167 target = 0.113691 restraints weight = 5198.078| |-----------------------------------------------------------------------------| r_work (start): 0.3166 rms_B_bonded: 1.41 r_work: 0.3064 rms_B_bonded: 1.92 restraints_weight: 0.5000 r_work: 0.2954 rms_B_bonded: 3.35 restraints_weight: 0.2500 r_work (final): 0.2954 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8450 moved from start: 0.2169 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.031 4108 Z= 0.268 Angle : 0.538 6.002 5604 Z= 0.289 Chirality : 0.042 0.144 642 Planarity : 0.004 0.056 681 Dihedral : 5.241 47.271 575 Min Nonbonded Distance : 2.583 Molprobity Statistics. All-atom Clashscore : 4.13 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.62 % Favored : 97.38 % Rotamer: Outliers : 3.55 % Allowed : 20.33 % Favored : 76.12 % 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.39), residues: 496 helix: 1.64 (0.30), residues: 300 sheet: -0.36 (0.74), residues: 46 loop : -0.33 (0.54), residues: 150 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.019 0.002 TRP A 94 HIS 0.001 0.000 HIS A 135 PHE 0.021 0.002 PHE A 71 TYR 0.014 0.001 TYR A 280 ARG 0.002 0.000 ARG B 38 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 70 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 15 poor density : 55 time to evaluate : 0.449 Fit side-chains REVERT: A 363 MET cc_start: 0.8997 (mtm) cc_final: 0.8792 (ttp) REVERT: A 455 MET cc_start: 0.7035 (OUTLIER) cc_final: 0.6319 (ttp) outliers start: 15 outliers final: 12 residues processed: 65 average time/residue: 0.1149 time to fit residues: 10.1716 Evaluate side-chains 66 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 13 poor density : 53 time to evaluate : 0.466 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 12 VAL Chi-restraints excluded: chain B residue 51 ILE Chi-restraints excluded: chain B residue 69 THR Chi-restraints excluded: chain B residue 81 LEU Chi-restraints excluded: chain A residue 18 ILE Chi-restraints excluded: chain A residue 145 SER Chi-restraints excluded: chain A residue 283 LEU Chi-restraints excluded: chain A residue 294 GLN Chi-restraints excluded: chain A residue 366 THR Chi-restraints excluded: chain A residue 415 LEU Chi-restraints excluded: chain A residue 454 LEU Chi-restraints excluded: chain A residue 455 MET Chi-restraints excluded: chain A residue 458 MET Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 7 optimal weight: 0.4980 chunk 33 optimal weight: 0.0980 chunk 17 optimal weight: 0.0570 chunk 31 optimal weight: 0.9990 chunk 39 optimal weight: 0.9980 chunk 13 optimal weight: 0.2980 chunk 16 optimal weight: 0.7980 chunk 38 optimal weight: 0.6980 chunk 11 optimal weight: 2.9990 chunk 12 optimal weight: 0.0980 chunk 6 optimal weight: 0.0060 overall best weight: 0.1114 ------------------------------------------------------------------------------- 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.3397 r_free = 0.3397 target = 0.132204 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 53)----------------| | r_work = 0.3251 r_free = 0.3251 target = 0.120160 restraints weight = 5081.666| |-----------------------------------------------------------------------------| r_work (start): 0.3246 rms_B_bonded: 1.22 r_work: 0.3157 rms_B_bonded: 1.78 restraints_weight: 0.5000 r_work: 0.3051 rms_B_bonded: 3.12 restraints_weight: 0.2500 r_work (final): 0.3051 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8351 moved from start: 0.2470 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.027 4108 Z= 0.137 Angle : 0.488 7.022 5604 Z= 0.263 Chirality : 0.040 0.132 642 Planarity : 0.004 0.053 681 Dihedral : 4.952 47.524 575 Min Nonbonded Distance : 2.613 Molprobity Statistics. All-atom Clashscore : 4.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.22 % Favored : 97.78 % Rotamer: Outliers : 1.89 % Allowed : 21.99 % Favored : 76.12 % 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.67 (0.39), residues: 496 helix: 1.96 (0.30), residues: 301 sheet: -0.11 (0.75), residues: 46 loop : -0.19 (0.54), residues: 149 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.020 0.001 TRP A 94 HIS 0.002 0.001 HIS A 364 PHE 0.012 0.001 PHE A 71 TYR 0.012 0.001 TYR A 445 ARG 0.002 0.000 ARG B 19 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 70 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 8 poor density : 62 time to evaluate : 0.479 Fit side-chains REVERT: A 74 THR cc_start: 0.8374 (m) cc_final: 0.8031 (m) REVERT: A 455 MET cc_start: 0.6841 (OUTLIER) cc_final: 0.6134 (ttp) outliers start: 8 outliers final: 4 residues processed: 66 average time/residue: 0.1207 time to fit residues: 10.8916 Evaluate side-chains 63 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 58 time to evaluate : 0.407 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 12 VAL Chi-restraints excluded: chain B residue 81 LEU Chi-restraints excluded: chain A residue 294 GLN Chi-restraints excluded: chain A residue 455 MET Chi-restraints excluded: chain A residue 458 MET Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 25 optimal weight: 0.9980 chunk 46 optimal weight: 0.2980 chunk 47 optimal weight: 0.9990 chunk 15 optimal weight: 0.6980 chunk 33 optimal weight: 0.7980 chunk 18 optimal weight: 0.9980 chunk 2 optimal weight: 0.9980 chunk 3 optimal weight: 0.9990 chunk 6 optimal weight: 0.7980 chunk 32 optimal weight: 0.5980 chunk 38 optimal weight: 0.6980 overall best weight: 0.6180 ------------------------------------------------------------------------------- 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.3333 r_free = 0.3333 target = 0.126975 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 36)----------------| | r_work = 0.3171 r_free = 0.3171 target = 0.114195 restraints weight = 5124.849| |-----------------------------------------------------------------------------| r_work (start): 0.3168 rms_B_bonded: 1.41 r_work: 0.3070 rms_B_bonded: 1.88 restraints_weight: 0.5000 r_work: 0.2959 rms_B_bonded: 3.29 restraints_weight: 0.2500 r_work (final): 0.2959 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8443 moved from start: 0.2330 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.043 4108 Z= 0.290 Angle : 0.537 6.046 5604 Z= 0.290 Chirality : 0.043 0.141 642 Planarity : 0.004 0.053 681 Dihedral : 5.110 47.263 575 Min Nonbonded Distance : 2.599 Molprobity Statistics. All-atom Clashscore : 4.50 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.62 % Favored : 97.38 % Rotamer: Outliers : 1.65 % Allowed : 22.46 % Favored : 75.89 % 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.50 (0.39), residues: 496 helix: 1.86 (0.30), residues: 300 sheet: -0.25 (0.73), residues: 46 loop : -0.29 (0.54), residues: 150 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.017 0.002 TRP A 94 HIS 0.001 0.000 HIS A 135 PHE 0.022 0.002 PHE A 71 TYR 0.013 0.001 TYR A 445 ARG 0.002 0.000 ARG A 402 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 62 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 7 poor density : 55 time to evaluate : 0.490 Fit side-chains REVERT: A 455 MET cc_start: 0.6853 (OUTLIER) cc_final: 0.6084 (ttp) outliers start: 7 outliers final: 4 residues processed: 59 average time/residue: 0.1305 time to fit residues: 10.4905 Evaluate side-chains 60 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 55 time to evaluate : 0.454 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 12 VAL Chi-restraints excluded: chain B residue 51 ILE Chi-restraints excluded: chain B residue 81 LEU Chi-restraints excluded: chain A residue 294 GLN Chi-restraints excluded: chain A residue 455 MET Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 40 optimal weight: 0.9980 chunk 28 optimal weight: 0.2980 chunk 23 optimal weight: 0.9990 chunk 4 optimal weight: 0.9990 chunk 38 optimal weight: 0.9980 chunk 35 optimal weight: 1.9990 chunk 31 optimal weight: 0.5980 chunk 34 optimal weight: 1.9990 chunk 7 optimal weight: 0.9980 chunk 39 optimal weight: 0.6980 chunk 29 optimal weight: 0.1980 overall best weight: 0.5580 ------------------------------------------------------------------------------- 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.3333 r_free = 0.3333 target = 0.127023 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 36)----------------| | r_work = 0.3170 r_free = 0.3170 target = 0.114086 restraints weight = 5145.447| |-----------------------------------------------------------------------------| r_work (start): 0.3169 rms_B_bonded: 1.42 r_work: 0.3068 rms_B_bonded: 1.91 restraints_weight: 0.5000 r_work: 0.2956 rms_B_bonded: 3.35 restraints_weight: 0.2500 r_work (final): 0.2956 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8445 moved from start: 0.2377 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.039 4108 Z= 0.271 Angle : 0.544 6.759 5604 Z= 0.293 Chirality : 0.043 0.142 642 Planarity : 0.004 0.054 681 Dihedral : 5.162 47.341 575 Min Nonbonded Distance : 2.593 Molprobity Statistics. All-atom Clashscore : 4.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.62 % Favored : 97.38 % Rotamer: Outliers : 1.65 % Allowed : 23.64 % Favored : 74.70 % 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.48 (0.39), residues: 496 helix: 1.84 (0.30), residues: 300 sheet: -0.34 (0.73), residues: 46 loop : -0.28 (0.54), residues: 150 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.018 0.002 TRP A 94 HIS 0.001 0.000 HIS A 135 PHE 0.021 0.002 PHE A 71 TYR 0.012 0.001 TYR A 445 ARG 0.002 0.000 ARG A 456 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 61 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 7 poor density : 54 time to evaluate : 0.349 Fit side-chains REVERT: A 455 MET cc_start: 0.6791 (OUTLIER) cc_final: 0.5987 (ttp) outliers start: 7 outliers final: 4 residues processed: 58 average time/residue: 0.1318 time to fit residues: 10.5803 Evaluate side-chains 59 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 54 time to evaluate : 0.453 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 12 VAL Chi-restraints excluded: chain B residue 51 ILE Chi-restraints excluded: chain B residue 81 LEU Chi-restraints excluded: chain A residue 294 GLN Chi-restraints excluded: chain A residue 455 MET Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 8 optimal weight: 0.9990 chunk 29 optimal weight: 0.2980 chunk 7 optimal weight: 0.3980 chunk 33 optimal weight: 0.3980 chunk 1 optimal weight: 1.9990 chunk 28 optimal weight: 0.5980 chunk 3 optimal weight: 0.0870 chunk 10 optimal weight: 0.6980 chunk 40 optimal weight: 0.1980 chunk 27 optimal weight: 0.9990 chunk 9 optimal weight: 0.7980 overall best weight: 0.2758 ------------------------------------------------------------------------------- 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.3375 r_free = 0.3375 target = 0.130614 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 34)----------------| | r_work = 0.3218 r_free = 0.3218 target = 0.117762 restraints weight = 5069.195| |-----------------------------------------------------------------------------| r_work (start): 0.3217 rms_B_bonded: 1.39 r_work: 0.3121 rms_B_bonded: 1.87 restraints_weight: 0.5000 r_work: 0.3011 rms_B_bonded: 3.30 restraints_weight: 0.2500 r_work (final): 0.3011 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8386 moved from start: 0.2584 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.043 4108 Z= 0.179 Angle : 0.492 7.429 5604 Z= 0.265 Chirality : 0.040 0.135 642 Planarity : 0.004 0.053 681 Dihedral : 4.957 47.219 575 Min Nonbonded Distance : 2.620 Molprobity Statistics. All-atom Clashscore : 4.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.22 % Favored : 97.78 % Rotamer: Outliers : 1.42 % Allowed : 23.17 % Favored : 75.41 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.74 (0.39), residues: 496 helix: 2.04 (0.31), residues: 300 sheet: -0.22 (0.74), residues: 46 loop : -0.15 (0.55), residues: 150 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.020 0.002 TRP A 94 HIS 0.001 0.001 HIS A 364 PHE 0.016 0.001 PHE A 71 TYR 0.012 0.001 TYR A 445 ARG 0.002 0.000 ARG A 456 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 992 Ramachandran restraints generated. 496 Oldfield, 0 Emsley, 496 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 61 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 6 poor density : 55 time to evaluate : 0.489 Fit side-chains REVERT: A 74 THR cc_start: 0.8411 (m) cc_final: 0.8063 (m) REVERT: A 455 MET cc_start: 0.6667 (OUTLIER) cc_final: 0.5848 (ttp) outliers start: 6 outliers final: 4 residues processed: 58 average time/residue: 0.1301 time to fit residues: 10.2134 Evaluate side-chains 60 residues out of total 423 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 55 time to evaluate : 0.422 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 12 VAL Chi-restraints excluded: chain B residue 81 LEU Chi-restraints excluded: chain A residue 294 GLN Chi-restraints excluded: chain A residue 454 LEU Chi-restraints excluded: chain A residue 455 MET Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 23 optimal weight: 0.9980 chunk 6 optimal weight: 0.7980 chunk 10 optimal weight: 0.6980 chunk 14 optimal weight: 0.9990 chunk 27 optimal weight: 0.6980 chunk 25 optimal weight: 0.9990 chunk 1 optimal weight: 2.9990 chunk 16 optimal weight: 0.4980 chunk 2 optimal weight: 0.9990 chunk 24 optimal weight: 0.5980 chunk 48 optimal weight: 0.9980 overall best weight: 0.6580 ------------------------------------------------------------------------------- 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.3325 r_free = 0.3325 target = 0.126439 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 34)----------------| | r_work = 0.3162 r_free = 0.3162 target = 0.113486 restraints weight = 5104.149| |-----------------------------------------------------------------------------| r_work (start): 0.3158 rms_B_bonded: 1.43 r_work: 0.3057 rms_B_bonded: 1.92 restraints_weight: 0.5000 r_work: 0.2946 rms_B_bonded: 3.36 restraints_weight: 0.2500 r_work (final): 0.2946 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8457 moved from start: 0.2408 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.042 4108 Z= 0.305 Angle : 0.559 7.167 5604 Z= 0.300 Chirality : 0.043 0.145 642 Planarity : 0.004 0.053 681 Dihedral : 5.193 47.900 575 Min Nonbonded Distance : 2.594 Molprobity Statistics. All-atom Clashscore : 4.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.62 % Favored : 97.38 % Rotamer: Outliers : 1.89 % Allowed : 23.40 % Favored : 74.70 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.52 (0.39), residues: 496 helix: 1.84 (0.31), residues: 300 sheet: -0.15 (0.75), residues: 44 loop : -0.22 (0.54), residues: 152 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.017 0.002 TRP A 94 HIS 0.001 0.001 HIS A 135 PHE 0.022 0.002 PHE A 71 TYR 0.013 0.001 TYR A 445 ARG 0.002 0.000 ARG B 38 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 2202.85 seconds wall clock time: 39 minutes 49.77 seconds (2389.77 seconds total)