Starting phenix.real_space_refine (version: 1.21rc1) on Mon Jul 3 20:05:52 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wr7_32717/07_2023/7wr7_32717_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wr7_32717/07_2023/7wr7_32717.map Found restraint, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wr7_32717/07_2023/7wr7_32717_ligands.cif Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.1 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { restraint_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wr7_32717/07_2023/7wr7_32717_ligands.cif" default_restraint = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wr7_32717/07_2023/7wr7_32717_ligands.cif" real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wr7_32717/07_2023/7wr7_32717.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wr7_32717/07_2023/7wr7_32717.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wr7_32717/07_2023/7wr7_32717_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wr7_32717/07_2023/7wr7_32717_updated.pdb" } resolution = 3.1 write_initial_geo_file = False refinement { macro_cycles = 10 } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= -0.001 sd= 0.054 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Anisotropic ADP refinement not supported. Converting 34 atoms to isotropic. 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 29 5.16 5 C 2509 2.51 5 N 681 2.21 5 O 776 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped None Time to flip residues: 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5013/modules/chem_data/mon_lib" Total number of atoms: 3995 Number of models: 1 Model: "" Number of chains: 4 Chain: "A" Number of atoms: 2013 Number of conformers: 1 Conformer: "" Number of residues, atoms: 261, 2013 Classifications: {'peptide': 261} Link IDs: {'PTRANS': 15, 'TRANS': 245} Chain: "B" Number of atoms: 1844 Number of conformers: 1 Conformer: "" Number of residues, atoms: 235, 1844 Classifications: {'peptide': 235} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 14, 'TRANS': 220} Chain: "A" Number of atoms: 70 Number of conformers: 1 Conformer: "" Number of residues, atoms: 5, 70 Unusual residues: {'NAG': 5} Classifications: {'undetermined': 5} Link IDs: {None: 4} Unresolved non-hydrogen bonds: 5 Unresolved non-hydrogen angles: 10 Unresolved non-hydrogen chiralities: 5 Chain: "B" Number of atoms: 68 Number of conformers: 1 Conformer: "" Number of residues, atoms: 5, 68 Unusual residues: {'GBS': 1, 'NAG': 4} Classifications: {'undetermined': 5} Link IDs: {None: 4} Unresolved non-hydrogen bonds: 5 Unresolved non-hydrogen angles: 10 Unresolved non-hydrogen chiralities: 4 Planarities with less than four sites: {'GBS:plan-3': 1} Unresolved non-hydrogen planarities: 1 Time building chain proxies: 2.51, per 1000 atoms: 0.63 Number of scatterers: 3995 At special positions: 0 Unit cell: (64.852, 80.542, 95.186, 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 29 16.00 O 776 8.00 N 681 7.00 C 2509 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=8, symmetry=0 Simple disulfide: pdb=" SG CYS A 650 " - pdb=" SG CYS A 668 " distance=2.03 Simple disulfide: pdb=" SG CYS A 662 " - pdb=" SG CYS A 677 " distance=2.03 Simple disulfide: pdb=" SG CYS A 716 " - pdb=" SG CYS A 767 " distance=2.03 Simple disulfide: pdb=" SG CYS A 772 " - pdb=" SG CYS B 896 " distance=2.03 Simple disulfide: pdb=" SG CYS B 810 " - pdb=" SG CYS B 826 " distance=2.03 Simple disulfide: pdb=" SG CYS B 910 " - pdb=" SG CYS B 977 " distance=2.03 Simple disulfide: pdb=" SG CYS B 941 " - pdb=" SG CYS B 956 " distance=2.03 Simple disulfide: pdb=" SG CYS B 967 " - pdb=" SG CYS B 995 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=1, symmetry=0 Links applied NAG-ASN " NAG A 801 " - " ASN A 534 " " NAG A 802 " - " ASN A 630 " " NAG A 803 " - " ASN A 682 " " NAG A 804 " - " ASN A 706 " " NAG A 805 " - " ASN A 725 " " NAG B1101 " - " ASN B 848 " " NAG B1102 " - " ASN B 887 " " NAG B1103 " - " ASN B 909 " " NAG B1104 " - " ASN B 949 " Number of additional bonds: simple=1, symmetry=0 Coordination: Other bonds: Time building additional restraints: 1.63 Conformation dependent library (CDL) restraints added in 771.8 milliseconds 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 904 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 5 helices and 7 sheets defined 6.2% alpha, 18.3% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.35 Creating SS restraints... Processing helix chain 'A' and resid 672 through 674 No H-bonds generated for 'chain 'A' and resid 672 through 674' Processing helix chain 'A' and resid 709 through 718 Processing helix chain 'B' and resid 834 through 836 No H-bonds generated for 'chain 'B' and resid 834 through 836' Processing helix chain 'B' and resid 938 through 944 Processing helix chain 'B' and resid 1011 through 1018 removed outlier: 4.071A pdb=" N PHE B1017 " --> pdb=" O TRP B1013 " (cutoff:3.500A) Processing sheet with id= A, first strand: chain 'A' and resid 536 through 538 removed outlier: 3.744A pdb=" N PHE A 537 " --> pdb=" O ALA A 629 " (cutoff:3.500A) removed outlier: 3.656A pdb=" N ALA A 629 " --> pdb=" O PHE A 537 " (cutoff:3.500A) removed outlier: 3.550A pdb=" N ASP A 572 " --> pdb=" O GLY A 626 " (cutoff:3.500A) removed outlier: 6.975A pdb=" N ASN A 630 " --> pdb=" O PHE A 568 " (cutoff:3.500A) removed outlier: 5.732A pdb=" N PHE A 568 " --> pdb=" O ASN A 630 " (cutoff:3.500A) removed outlier: 7.754A pdb=" N THR A 632 " --> pdb=" O LEU A 566 " (cutoff:3.500A) removed outlier: 5.972A pdb=" N LEU A 566 " --> pdb=" O THR A 632 " (cutoff:3.500A) removed outlier: 7.306A pdb=" N GLY A 634 " --> pdb=" O ILE A 564 " (cutoff:3.500A) removed outlier: 5.992A pdb=" N ILE A 564 " --> pdb=" O GLY A 634 " (cutoff:3.500A) Processing sheet with id= B, first strand: chain 'A' and resid 551 through 557 removed outlier: 3.757A pdb=" N CYS A 552 " --> pdb=" O LEU A 616 " (cutoff:3.500A) removed outlier: 3.545A pdb=" N ALA A 594 " --> pdb=" O ILE A 582 " (cutoff:3.500A) Processing sheet with id= C, first strand: chain 'B' and resid 855 through 858 removed outlier: 3.806A pdb=" N LEU B 802 " --> pdb=" O CYS B 810 " (cutoff:3.500A) removed outlier: 4.142A pdb=" N CYS B 810 " --> pdb=" O LEU B 802 " (cutoff:3.500A) removed outlier: 6.324A pdb=" N TYR B 804 " --> pdb=" O LEU B 808 " (cutoff:3.500A) removed outlier: 5.813A pdb=" N LEU B 808 " --> pdb=" O TYR B 804 " (cutoff:3.500A) Processing sheet with id= D, first strand: chain 'B' and resid 819 through 822 removed outlier: 3.973A pdb=" N GLU B 861 " --> pdb=" O HIS B 881 " (cutoff:3.500A) Processing sheet with id= E, first strand: chain 'B' and resid 909 through 914 removed outlier: 3.622A pdb=" N GLY B 914 " --> pdb=" O GLN B 929 " (cutoff:3.500A) removed outlier: 3.856A pdb=" N GLN B 929 " --> pdb=" O GLY B 914 " (cutoff:3.500A) Processing sheet with id= F, first strand: chain 'B' and resid 954 through 957 removed outlier: 5.851A pdb=" N ALA B1005 " --> pdb=" O THR B 989 " (cutoff:3.500A) removed outlier: 5.978A pdb=" N THR B 989 " --> pdb=" O ALA B1005 " (cutoff:3.500A) Processing sheet with id= G, first strand: chain 'B' and resid 976 through 979 removed outlier: 3.738A pdb=" N CYS B 977 " --> pdb=" O PHE B 984 " (cutoff:3.500A) 59 hydrogen bonds defined for protein. 159 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.06 Time building geometry restraints manager: 1.79 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.21 - 1.33: 1015 1.33 - 1.45: 823 1.45 - 1.57: 2218 1.57 - 1.69: 0 1.69 - 1.81: 36 Bond restraints: 4092 Sorted by residual: bond pdb=" C18 GBS B1105 " pdb=" N3 GBS B1105 " ideal model delta sigma weight residual 1.452 1.333 0.119 2.00e-02 2.50e+03 3.51e+01 bond pdb=" C18 GBS B1105 " pdb=" N2 GBS B1105 " ideal model delta sigma weight residual 1.454 1.341 0.113 2.00e-02 2.50e+03 3.20e+01 bond pdb=" C18 GBS B1105 " pdb=" N4 GBS B1105 " ideal model delta sigma weight residual 1.250 1.336 -0.086 2.00e-02 2.50e+03 1.85e+01 bond pdb=" N GLN B 968 " pdb=" CA GLN B 968 " ideal model delta sigma weight residual 1.456 1.487 -0.032 1.25e-02 6.40e+03 6.52e+00 bond pdb=" O GBS B1105 " pdb=" C6 GBS B1105 " ideal model delta sigma weight residual 1.260 1.210 0.050 2.00e-02 2.50e+03 6.17e+00 ... (remaining 4087 not shown) Histogram of bond angle deviations from ideal: 98.83 - 105.86: 108 105.86 - 112.89: 2193 112.89 - 119.92: 1303 119.92 - 126.94: 1897 126.94 - 133.97: 68 Bond angle restraints: 5569 Sorted by residual: angle pdb=" N ASP B 970 " pdb=" CA ASP B 970 " pdb=" C ASP B 970 " ideal model delta sigma weight residual 114.56 108.41 6.15 1.27e+00 6.20e-01 2.35e+01 angle pdb=" N PRO B 974 " pdb=" CA PRO B 974 " pdb=" CB PRO B 974 " ideal model delta sigma weight residual 103.25 98.83 4.42 1.05e+00 9.07e-01 1.77e+01 angle pdb=" CA GLN B 968 " pdb=" C GLN B 968 " pdb=" O GLN B 968 " ideal model delta sigma weight residual 121.55 118.32 3.23 1.06e+00 8.90e-01 9.31e+00 angle pdb=" N GLY B 972 " pdb=" CA GLY B 972 " pdb=" C GLY B 972 " ideal model delta sigma weight residual 115.32 111.24 4.08 1.38e+00 5.25e-01 8.75e+00 angle pdb=" N PHE A 542 " pdb=" CA PHE A 542 " pdb=" C PHE A 542 " ideal model delta sigma weight residual 109.81 115.70 -5.89 2.21e+00 2.05e-01 7.10e+00 ... (remaining 5564 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 16.22: 2232 16.22 - 32.45: 97 32.45 - 48.67: 20 48.67 - 64.89: 6 64.89 - 81.12: 1 Dihedral angle restraints: 2356 sinusoidal: 944 harmonic: 1412 Sorted by residual: dihedral pdb=" CA PRO A 543 " pdb=" C PRO A 543 " pdb=" N ASN A 544 " pdb=" CA ASN A 544 " ideal model delta harmonic sigma weight residual 180.00 159.94 20.06 0 5.00e+00 4.00e-02 1.61e+01 dihedral pdb=" CB CYS A 772 " pdb=" SG CYS A 772 " pdb=" SG CYS B 896 " pdb=" CB CYS B 896 " ideal model delta sinusoidal sigma weight residual -86.00 -113.70 27.70 1 1.00e+01 1.00e-02 1.10e+01 dihedral pdb=" CA SER A 539 " pdb=" C SER A 539 " pdb=" N THR A 540 " pdb=" CA THR A 540 " ideal model delta harmonic sigma weight residual 180.00 163.44 16.56 0 5.00e+00 4.00e-02 1.10e+01 ... (remaining 2353 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.033: 404 0.033 - 0.067: 140 0.067 - 0.100: 55 0.100 - 0.134: 28 0.134 - 0.167: 1 Chirality restraints: 628 Sorted by residual: chirality pdb=" C2 NAG A 804 " pdb=" C1 NAG A 804 " pdb=" C3 NAG A 804 " pdb=" N2 NAG A 804 " both_signs ideal model delta sigma weight residual False -2.49 -2.32 -0.17 2.00e-01 2.50e+01 6.98e-01 chirality pdb=" CA VAL B 988 " pdb=" N VAL B 988 " pdb=" C VAL B 988 " pdb=" CB VAL B 988 " both_signs ideal model delta sigma weight residual False 2.44 2.57 -0.12 2.00e-01 2.50e+01 3.84e-01 chirality pdb=" CB VAL A 691 " pdb=" CA VAL A 691 " pdb=" CG1 VAL A 691 " pdb=" CG2 VAL A 691 " both_signs ideal model delta sigma weight residual False -2.63 -2.51 -0.12 2.00e-01 2.50e+01 3.71e-01 ... (remaining 625 not shown) Planarity restraints: 733 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C18 GBS B1105 " 0.013 2.00e-02 2.50e+03 3.50e-02 1.53e+01 pdb=" C2 GBS B1105 " -0.041 2.00e-02 2.50e+03 pdb=" N2 GBS B1105 " 0.055 2.00e-02 2.50e+03 pdb=" N3 GBS B1105 " -0.034 2.00e-02 2.50e+03 pdb=" N4 GBS B1105 " 0.007 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C GBS B1105 " -0.026 2.00e-02 2.50e+03 2.49e-02 1.24e+01 pdb=" C1 GBS B1105 " 0.009 2.00e-02 2.50e+03 pdb=" C2 GBS B1105 " 0.028 2.00e-02 2.50e+03 pdb=" C3 GBS B1105 " 0.008 2.00e-02 2.50e+03 pdb=" C4 GBS B1105 " -0.026 2.00e-02 2.50e+03 pdb=" C5 GBS B1105 " -0.031 2.00e-02 2.50e+03 pdb=" C6 GBS B1105 " 0.041 2.00e-02 2.50e+03 pdb=" N2 GBS B1105 " -0.004 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C GLU A 532 " -0.026 5.00e-02 4.00e+02 3.89e-02 2.42e+00 pdb=" N PRO A 533 " 0.067 5.00e-02 4.00e+02 pdb=" CA PRO A 533 " -0.020 5.00e-02 4.00e+02 pdb=" CD PRO A 533 " -0.022 5.00e-02 4.00e+02 ... (remaining 730 not shown) Histogram of nonbonded interaction distances: 2.31 - 2.83: 1246 2.83 - 3.35: 3712 3.35 - 3.86: 6420 3.86 - 4.38: 7358 4.38 - 4.90: 12603 Nonbonded interactions: 31339 Sorted by model distance: nonbonded pdb=" OD2 ASP B 874 " pdb=" ND2 ASN B 875 " model vdw 2.309 2.520 nonbonded pdb=" NH2 ARG B 908 " pdb=" O GLN B 978 " model vdw 2.350 2.520 nonbonded pdb=" O GLY B 811 " pdb=" O GLY B 972 " model vdw 2.357 3.040 nonbonded pdb=" NZ LYS A 740 " pdb=" O SER A 754 " model vdw 2.359 2.520 nonbonded pdb=" O CYS B 826 " pdb=" NH1 ARG B 830 " model vdw 2.367 2.520 ... (remaining 31334 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.850 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.020 Construct map_model_manager: 0.010 Extract box with map and model: 3.770 Check model and map are aligned: 0.070 Set scattering table: 0.050 Process input model: 15.190 Find NCS groups from input model: 0.120 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.590 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 22.690 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6962 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.119 4092 Z= 0.257 Angle : 0.602 6.152 5569 Z= 0.330 Chirality : 0.044 0.167 628 Planarity : 0.003 0.039 724 Dihedral : 9.677 81.116 1428 Min Nonbonded Distance : 2.309 Molprobity Statistics. All-atom Clashscore : 10.13 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.15 % Favored : 90.85 % Rotamer Outliers : 3.29 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.53 (0.30), residues: 492 helix: -2.64 (0.51), residues: 36 sheet: -2.64 (0.42), residues: 105 loop : -3.42 (0.28), residues: 351 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 125 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 111 time to evaluate : 0.461 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 14 outliers final: 5 residues processed: 121 average time/residue: 0.2436 time to fit residues: 34.2687 Evaluate side-chains 76 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 71 time to evaluate : 0.476 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.1341 time to fit residues: 1.4477 ------------------------------------------------------------------------------- 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: 3.9990 chunk 37 optimal weight: 2.9990 chunk 20 optimal weight: 2.9990 chunk 12 optimal weight: 9.9990 chunk 25 optimal weight: 1.9990 chunk 19 optimal weight: 3.9990 chunk 38 optimal weight: 0.9990 chunk 14 optimal weight: 9.9990 chunk 23 optimal weight: 0.7980 chunk 28 optimal weight: 0.3980 chunk 44 optimal weight: 0.9990 overall best weight: 1.0386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 649 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 713 ASN A 717 GLN B 844 HIS ** B 875 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 921 GLN B 938 ASN B 943 GLN B 978 GLN Total number of N/Q/H flips: 7 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7216 moved from start: 0.2654 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.037 4092 Z= 0.319 Angle : 0.631 8.874 5569 Z= 0.317 Chirality : 0.047 0.263 628 Planarity : 0.004 0.041 724 Dihedral : 4.612 47.648 529 Min Nonbonded Distance : 2.172 Molprobity Statistics. All-atom Clashscore : 15.13 Ramachandran Plot: Outliers : 0.00 % Allowed : 11.18 % Favored : 88.82 % Rotamer Outliers : 3.06 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.87 (0.33), residues: 492 helix: -2.22 (0.68), residues: 38 sheet: -1.99 (0.42), residues: 105 loop : -3.03 (0.30), residues: 349 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 80 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 67 time to evaluate : 0.353 Fit side-chains revert: symmetry clash outliers start: 13 outliers final: 8 residues processed: 77 average time/residue: 0.1843 time to fit residues: 16.8722 Evaluate side-chains 76 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 68 time to evaluate : 0.446 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 0 residues processed: 8 average time/residue: 0.0328 time to fit residues: 1.0023 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 24 optimal weight: 4.9990 chunk 13 optimal weight: 7.9990 chunk 37 optimal weight: 1.9990 chunk 30 optimal weight: 0.9990 chunk 12 optimal weight: 2.9990 chunk 44 optimal weight: 2.9990 chunk 48 optimal weight: 0.4980 chunk 39 optimal weight: 0.0170 chunk 15 optimal weight: 4.9990 chunk 35 optimal weight: 0.9980 chunk 33 optimal weight: 0.4980 overall best weight: 0.6020 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 563 ASN ** A 688 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 875 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 938 ASN B 978 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7198 moved from start: 0.2854 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.029 4092 Z= 0.211 Angle : 0.570 9.240 5569 Z= 0.284 Chirality : 0.046 0.248 628 Planarity : 0.004 0.045 724 Dihedral : 4.559 45.907 529 Min Nonbonded Distance : 2.169 Molprobity Statistics. All-atom Clashscore : 13.08 Ramachandran Plot: Outliers : 0.00 % Allowed : 10.37 % Favored : 89.63 % Rotamer Outliers : 2.59 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.56 (0.34), residues: 492 helix: -1.63 (0.74), residues: 38 sheet: -1.84 (0.41), residues: 118 loop : -2.83 (0.32), residues: 336 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 79 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 68 time to evaluate : 0.483 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 11 outliers final: 5 residues processed: 76 average time/residue: 0.2258 time to fit residues: 20.4219 Evaluate side-chains 69 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 64 time to evaluate : 0.473 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.0453 time to fit residues: 1.0105 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 23 optimal weight: 6.9990 chunk 4 optimal weight: 0.0980 chunk 21 optimal weight: 7.9990 chunk 29 optimal weight: 0.8980 chunk 44 optimal weight: 0.4980 chunk 47 optimal weight: 2.9990 chunk 42 optimal weight: 0.2980 chunk 12 optimal weight: 10.0000 chunk 39 optimal weight: 0.9990 chunk 26 optimal weight: 0.3980 chunk 0 optimal weight: 2.9990 overall best weight: 0.4380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 688 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7185 moved from start: 0.3019 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.026 4092 Z= 0.185 Angle : 0.570 9.303 5569 Z= 0.285 Chirality : 0.045 0.245 628 Planarity : 0.004 0.044 724 Dihedral : 4.515 44.413 529 Min Nonbonded Distance : 2.181 Molprobity Statistics. All-atom Clashscore : 13.21 Ramachandran Plot: Outliers : 0.00 % Allowed : 10.77 % Favored : 89.23 % Rotamer Outliers : 1.88 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.24 (0.35), residues: 492 helix: -1.33 (0.78), residues: 38 sheet: -1.34 (0.44), residues: 111 loop : -2.74 (0.32), residues: 343 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 71 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 63 time to evaluate : 0.465 Fit side-chains outliers start: 8 outliers final: 5 residues processed: 67 average time/residue: 0.2089 time to fit residues: 16.9291 Evaluate side-chains 63 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 58 time to evaluate : 0.467 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.0428 time to fit residues: 0.9977 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 35 optimal weight: 0.7980 chunk 19 optimal weight: 4.9990 chunk 40 optimal weight: 0.9980 chunk 32 optimal weight: 5.9990 chunk 0 optimal weight: 2.9990 chunk 24 optimal weight: 0.7980 chunk 42 optimal weight: 1.9990 chunk 11 optimal weight: 0.0040 chunk 15 optimal weight: 0.9980 chunk 9 optimal weight: 0.5980 chunk 27 optimal weight: 0.9980 overall best weight: 0.6392 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 565 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 688 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7217 moved from start: 0.3220 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.031 4092 Z= 0.224 Angle : 0.588 9.264 5569 Z= 0.294 Chirality : 0.045 0.241 628 Planarity : 0.004 0.046 724 Dihedral : 4.564 41.560 529 Min Nonbonded Distance : 2.188 Molprobity Statistics. All-atom Clashscore : 12.70 Ramachandran Plot: Outliers : 0.00 % Allowed : 11.59 % Favored : 88.41 % Rotamer Outliers : 1.18 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.15 (0.35), residues: 492 helix: -1.18 (0.79), residues: 38 sheet: -1.13 (0.44), residues: 119 loop : -2.77 (0.32), residues: 335 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 67 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 62 time to evaluate : 0.449 Fit side-chains revert: symmetry clash outliers start: 5 outliers final: 2 residues processed: 66 average time/residue: 0.2146 time to fit residues: 17.0934 Evaluate side-chains 62 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 60 time to evaluate : 0.485 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0406 time to fit residues: 0.7694 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 11 optimal weight: 0.6980 chunk 47 optimal weight: 2.9990 chunk 39 optimal weight: 3.9990 chunk 21 optimal weight: 0.9980 chunk 3 optimal weight: 0.3980 chunk 15 optimal weight: 0.8980 chunk 24 optimal weight: 3.9990 chunk 45 optimal weight: 0.6980 chunk 5 optimal weight: 2.9990 chunk 27 optimal weight: 0.9980 chunk 34 optimal weight: 1.9990 overall best weight: 0.7380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 565 GLN A 686 ASN ** A 688 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 978 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7237 moved from start: 0.3464 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.035 4092 Z= 0.244 Angle : 0.592 9.253 5569 Z= 0.298 Chirality : 0.046 0.238 628 Planarity : 0.004 0.047 724 Dihedral : 4.673 40.561 529 Min Nonbonded Distance : 2.180 Molprobity Statistics. All-atom Clashscore : 13.21 Ramachandran Plot: Outliers : 0.00 % Allowed : 12.40 % Favored : 87.60 % Rotamer Outliers : 0.71 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.11 (0.35), residues: 492 helix: -1.13 (0.78), residues: 38 sheet: -1.15 (0.44), residues: 118 loop : -2.72 (0.32), residues: 336 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 64 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 61 time to evaluate : 0.477 Fit side-chains revert: symmetry clash outliers start: 3 outliers final: 0 residues processed: 63 average time/residue: 0.2267 time to fit residues: 17.2393 Evaluate side-chains 57 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 57 time to evaluate : 0.477 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.6267 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 26 optimal weight: 0.9990 chunk 39 optimal weight: 3.9990 chunk 47 optimal weight: 4.9990 chunk 29 optimal weight: 0.0970 chunk 28 optimal weight: 0.9980 chunk 21 optimal weight: 4.9990 chunk 18 optimal weight: 1.9990 chunk 14 optimal weight: 8.9990 chunk 9 optimal weight: 0.7980 chunk 30 optimal weight: 0.9980 chunk 32 optimal weight: 3.9990 overall best weight: 0.7780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 565 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7250 moved from start: 0.3651 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.034 4092 Z= 0.257 Angle : 0.615 9.241 5569 Z= 0.309 Chirality : 0.047 0.238 628 Planarity : 0.004 0.047 724 Dihedral : 4.803 41.599 529 Min Nonbonded Distance : 2.196 Molprobity Statistics. All-atom Clashscore : 13.59 Ramachandran Plot: Outliers : 0.00 % Allowed : 12.60 % Favored : 87.40 % Rotamer Outliers : 1.88 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.07 (0.35), residues: 492 helix: -1.03 (0.78), residues: 38 sheet: -1.12 (0.44), residues: 118 loop : -2.70 (0.32), residues: 336 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 71 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 63 time to evaluate : 0.474 Fit side-chains outliers start: 8 outliers final: 3 residues processed: 69 average time/residue: 0.2115 time to fit residues: 17.7327 Evaluate side-chains 61 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 58 time to evaluate : 0.457 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0408 time to fit residues: 0.8485 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 23 optimal weight: 7.9990 chunk 4 optimal weight: 0.4980 chunk 37 optimal weight: 0.6980 chunk 42 optimal weight: 1.9990 chunk 45 optimal weight: 0.4980 chunk 41 optimal weight: 0.7980 chunk 44 optimal weight: 1.9990 chunk 26 optimal weight: 0.0770 chunk 19 optimal weight: 0.9990 chunk 34 optimal weight: 0.8980 chunk 13 optimal weight: 10.0000 overall best weight: 0.5138 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 688 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7221 moved from start: 0.3695 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.028 4092 Z= 0.207 Angle : 0.612 9.270 5569 Z= 0.309 Chirality : 0.046 0.237 628 Planarity : 0.004 0.048 724 Dihedral : 4.800 45.319 529 Min Nonbonded Distance : 2.164 Molprobity Statistics. All-atom Clashscore : 13.46 Ramachandran Plot: Outliers : 0.00 % Allowed : 11.38 % Favored : 88.62 % Rotamer Outliers : 1.18 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.90 (0.36), residues: 492 helix: -1.13 (0.79), residues: 38 sheet: -0.84 (0.45), residues: 118 loop : -2.62 (0.32), residues: 336 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 61 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 56 time to evaluate : 0.476 Fit side-chains revert: symmetry clash outliers start: 5 outliers final: 1 residues processed: 61 average time/residue: 0.2158 time to fit residues: 15.8497 Evaluate side-chains 57 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 56 time to evaluate : 0.494 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0442 time to fit residues: 0.7340 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 39 optimal weight: 0.0270 chunk 41 optimal weight: 2.9990 chunk 43 optimal weight: 0.7980 chunk 28 optimal weight: 4.9990 chunk 46 optimal weight: 2.9990 chunk 22 optimal weight: 10.0000 chunk 32 optimal weight: 5.9990 chunk 48 optimal weight: 0.9990 chunk 44 optimal weight: 0.9990 chunk 38 optimal weight: 0.3980 chunk 4 optimal weight: 0.2980 overall best weight: 0.5040 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 688 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 717 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7218 moved from start: 0.3744 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.028 4092 Z= 0.208 Angle : 0.610 9.277 5569 Z= 0.307 Chirality : 0.046 0.236 628 Planarity : 0.004 0.051 724 Dihedral : 4.839 46.425 529 Min Nonbonded Distance : 2.187 Molprobity Statistics. All-atom Clashscore : 13.21 Ramachandran Plot: Outliers : 0.00 % Allowed : 11.99 % Favored : 88.01 % Rotamer Outliers : 0.94 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.74 (0.36), residues: 492 helix: -1.05 (0.79), residues: 38 sheet: -0.66 (0.46), residues: 118 loop : -2.53 (0.33), residues: 336 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 63 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 59 time to evaluate : 0.470 Fit side-chains revert: symmetry clash outliers start: 4 outliers final: 1 residues processed: 62 average time/residue: 0.2222 time to fit residues: 16.6281 Evaluate side-chains 59 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 58 time to evaluate : 0.404 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0412 time to fit residues: 0.6890 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 30 optimal weight: 0.9990 chunk 23 optimal weight: 5.9990 chunk 41 optimal weight: 1.9990 chunk 11 optimal weight: 1.9990 chunk 35 optimal weight: 0.7980 chunk 5 optimal weight: 0.9990 chunk 10 optimal weight: 0.9980 chunk 38 optimal weight: 0.7980 chunk 16 optimal weight: 6.9990 chunk 39 optimal weight: 0.6980 chunk 4 optimal weight: 0.7980 overall best weight: 0.8180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 717 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7259 moved from start: 0.3944 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.031 4092 Z= 0.271 Angle : 0.628 9.217 5569 Z= 0.316 Chirality : 0.047 0.235 628 Planarity : 0.004 0.055 724 Dihedral : 5.032 45.822 529 Min Nonbonded Distance : 2.167 Molprobity Statistics. All-atom Clashscore : 14.49 Ramachandran Plot: Outliers : 0.00 % Allowed : 13.21 % Favored : 86.79 % Rotamer Outliers : 0.24 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.74 (0.37), residues: 492 helix: -1.01 (0.81), residues: 38 sheet: -0.78 (0.46), residues: 118 loop : -2.48 (0.33), residues: 336 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 984 Ramachandran restraints generated. 492 Oldfield, 0 Emsley, 492 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 59 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 58 time to evaluate : 0.473 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 1 outliers final: 0 residues processed: 58 average time/residue: 0.2368 time to fit residues: 16.5060 Evaluate side-chains 56 residues out of total 425 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 56 time to evaluate : 0.472 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.5640 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 49 random chunks: chunk 7 optimal weight: 0.9990 chunk 34 optimal weight: 0.9980 chunk 2 optimal weight: 0.9990 chunk 28 optimal weight: 3.9990 chunk 44 optimal weight: 0.9980 chunk 26 optimal weight: 0.4980 chunk 33 optimal weight: 1.9990 chunk 1 optimal weight: 0.2980 chunk 31 optimal weight: 0.6980 chunk 29 optimal weight: 0.9980 chunk 19 optimal weight: 2.9990 overall best weight: 0.6980 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 717 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 943 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4047 r_free = 0.4047 target = 0.153935 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 41)----------------| | r_work = 0.3467 r_free = 0.3467 target = 0.103555 restraints weight = 5142.456| |-----------------------------------------------------------------------------| r_work (start): 0.3440 rms_B_bonded: 3.76 r_work: 0.3260 rms_B_bonded: 4.00 restraints_weight: 0.5000 r_work (final): 0.3260 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7009 moved from start: 0.3999 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.036 4092 Z= 0.247 Angle : 0.640 9.233 5569 Z= 0.326 Chirality : 0.047 0.235 628 Planarity : 0.004 0.048 724 Dihedral : 4.999 48.040 529 Min Nonbonded Distance : 2.151 Molprobity Statistics. All-atom Clashscore : 14.36 Ramachandran Plot: Outliers : 0.00 % Allowed : 12.60 % Favored : 87.40 % Rotamer Outliers : 0.24 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.78 (0.36), residues: 492 helix: -1.08 (0.79), residues: 37 sheet: -0.69 (0.46), residues: 118 loop : -2.56 (0.33), residues: 337 =============================================================================== Job complete usr+sys time: 1373.66 seconds wall clock time: 25 minutes 41.48 seconds (1541.48 seconds total)