Starting phenix.real_space_refine (version: dev) on Thu May 12 16:15:47 2022 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/3j9d_6239/05_2022/3j9d_6239.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/3j9d_6239/05_2022/3j9d_6239.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.3 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/3j9d_6239/05_2022/3j9d_6239.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/3j9d_6239/05_2022/3j9d_6239.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/3j9d_6239/05_2022/3j9d_6239.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/3j9d_6239/05_2022/3j9d_6239.pdb" } resolution = 3.3 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.013 sd= 0.220 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 Convert atoms to be neutral Process input model Symmetric amino acids flipped Residue "A ASP 2": "OD1" <-> "OD2" Residue "A GLU 3": "OE1" <-> "OE2" Residue "A GLU 22": "OE1" <-> "OE2" Residue "A GLU 54": "OE1" <-> "OE2" Residue "A TYR 83": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ASP 141": "OD1" <-> "OD2" Residue "A ARG 175": "NH1" <-> "NH2" Residue "A GLU 177": "OE1" <-> "OE2" Residue "A GLU 437": "OE1" <-> "OE2" Residue "A GLU 480": "OE1" <-> "OE2" Residue "A ASP 498": "OD1" <-> "OD2" Residue "A PHE 499": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 500": "OE1" <-> "OE2" Residue "A GLU 509": "OE1" <-> "OE2" Residue "A GLU 535": "OE1" <-> "OE2" Residue "A ASP 557": "OD1" <-> "OD2" Residue "A GLU 560": "OE1" <-> "OE2" Residue "A TYR 570": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 628": "NH1" <-> "NH2" Residue "A GLU 647": "OE1" <-> "OE2" Residue "A GLU 648": "OE1" <-> "OE2" Residue "A GLU 652": "OE1" <-> "OE2" Residue "A ASP 669": "OD1" <-> "OD2" Residue "A ARG 684": "NH1" <-> "NH2" Residue "A GLU 693": "OE1" <-> "OE2" Residue "A ARG 705": "NH1" <-> "NH2" Residue "A GLU 708": "OE1" <-> "OE2" Residue "A GLU 715": "OE1" <-> "OE2" Residue "A GLU 776": "OE1" <-> "OE2" Residue "A ASP 782": "OD1" <-> "OD2" Residue "A TYR 789": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 796": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 802": "NH1" <-> "NH2" Residue "A PHE 803": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 808": "OE1" <-> "OE2" Residue "A GLU 879": "OE1" <-> "OE2" Residue "A ARG 893": "NH1" <-> "NH2" Residue "A ASP 910": "OD1" <-> "OD2" Residue "A GLU 914": "OE1" <-> "OE2" Residue "A TYR 918": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4597/modules/chem_data/mon_lib" Total number of atoms: 6086 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 6086 Number of conformers: 1 Conformer: "" Number of residues, atoms: 747, 6086 Unusual residues: {' ZN': 1} Classifications: {'undetermined': 1, 'peptide': 746} Link IDs: {'PTRANS': 28, 'CIS': 14, 'TRANS': 701, None: 1, 'PCIS': 2} Not linked: pdbres="VAL A 961 " pdbres=" ZN A1001 " Chain breaks: 1 Unresolved chain link angles: 2 Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 2 List of CYS excluded from plausible disulfide bonds: (reason: may participate in coordination) ATOM 1310 SG CYS A 162 20.064 54.313 30.963 1.00 92.09 S ATOM 3290 SG CYS A 617 23.647 56.157 31.648 1.00 70.99 S ATOM 5217 SG CYS A 851 20.362 53.869 34.614 1.00 39.61 S Time building chain proxies: 4.04, per 1000 atoms: 0.66 Number of scatterers: 6086 At special positions: 0 Unit cell: (84.84, 107.06, 81.81, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) Zn 1 29.99 S 32 16.00 O 1105 8.00 N 1057 7.00 C 3891 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 2.25 Conformation dependent library (CDL) restraints added in 973.0 milliseconds Dynamic metal coordination Zn2+ tetrahedral coordination Zn2+ tetrahedral coordination pdb=" ZN A1001 " pdb="ZN ZN A1001 " - pdb=" NE2 HIS A 164 " pdb="ZN ZN A1001 " - pdb=" SG CYS A 162 " pdb="ZN ZN A1001 " - pdb=" SG CYS A 617 " pdb="ZN ZN A1001 " - pdb=" SG CYS A 851 " Number of angles added : 3 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1414 Finding SS restraints... Secondary structure from input PDB file: 33 helices and 8 sheets defined 35.9% alpha, 12.5% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.40 Creating SS restraints... Processing helix chain 'A' and resid 40 through 43 Processing helix chain 'A' and resid 55 through 59 Processing helix chain 'A' and resid 74 through 85 Processing helix chain 'A' and resid 103 through 110 Processing helix chain 'A' and resid 144 through 148 removed outlier: 4.027A pdb=" N THR A 147 " --> pdb=" O LYS A 144 " (cutoff:3.500A) Processing helix chain 'A' and resid 157 through 161 Processing helix chain 'A' and resid 163 through 165 No H-bonds generated for 'chain 'A' and resid 163 through 165' Processing helix chain 'A' and resid 166 through 176 Processing helix chain 'A' and resid 412 through 416 removed outlier: 3.544A pdb=" N ASP A 415 " --> pdb=" O ASN A 412 " (cutoff:3.500A) removed outlier: 3.953A pdb=" N PHE A 416 " --> pdb=" O LEU A 413 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 412 through 416' Processing helix chain 'A' and resid 461 through 476 Processing helix chain 'A' and resid 484 through 488 removed outlier: 3.506A pdb=" N ILE A 487 " --> pdb=" O MET A 484 " (cutoff:3.500A) Processing helix chain 'A' and resid 523 through 527 Processing helix chain 'A' and resid 557 through 562 Processing helix chain 'A' and resid 589 through 596 Processing helix chain 'A' and resid 599 through 607 removed outlier: 3.917A pdb=" N ARG A 606 " --> pdb=" O GLU A 602 " (cutoff:3.500A) removed outlier: 3.656A pdb=" N ARG A 607 " --> pdb=" O ILE A 603 " (cutoff:3.500A) Processing helix chain 'A' and resid 617 through 638 removed outlier: 3.598A pdb=" N THR A 621 " --> pdb=" O CYS A 617 " (cutoff:3.500A) Processing helix chain 'A' and resid 639 through 641 No H-bonds generated for 'chain 'A' and resid 639 through 641' Processing helix chain 'A' and resid 669 through 682 WARNING: missing atoms! Processing helix chain 'A' and resid 690 through 704 removed outlier: 3.817A pdb=" N MET A 704 " --> pdb=" O ARG A 700 " (cutoff:3.500A) Processing helix chain 'A' and resid 705 through 713 removed outlier: 3.851A pdb=" N LEU A 710 " --> pdb=" O GLY A 706 " (cutoff:3.500A) Processing helix chain 'A' and resid 717 through 725 Processing helix chain 'A' and resid 732 through 737 removed outlier: 3.998A pdb=" N ILE A 736 " --> pdb=" O VAL A 732 " (cutoff:3.500A) Processing helix chain 'A' and resid 738 through 741 Processing helix chain 'A' and resid 742 through 747 removed outlier: 4.454A pdb=" N PHE A 746 " --> pdb=" O LEU A 742 " (cutoff:3.500A) Processing helix chain 'A' and resid 786 through 791 removed outlier: 4.407A pdb=" N GLU A 791 " --> pdb=" O ILE A 787 " (cutoff:3.500A) Processing helix chain 'A' and resid 792 through 795 Processing helix chain 'A' and resid 798 through 804 Processing helix chain 'A' and resid 808 through 822 removed outlier: 4.476A pdb=" N LYS A 812 " --> pdb=" O GLU A 808 " (cutoff:3.500A) Processing helix chain 'A' and resid 832 through 837 Processing helix chain 'A' and resid 839 through 851 Processing helix chain 'A' and resid 882 through 892 Processing helix chain 'A' and resid 897 through 899 No H-bonds generated for 'chain 'A' and resid 897 through 899' Processing helix chain 'A' and resid 953 through 960 Processing sheet with id= A, first strand: chain 'A' and resid 130 through 133 removed outlier: 10.664A pdb=" N PHE A 23 " --> pdb=" O GLU A 3 " (cutoff:3.500A) removed outlier: 7.318A pdb=" N GLY A 5 " --> pdb=" O PHE A 23 " (cutoff:3.500A) removed outlier: 8.360A pdb=" N ILE A 25 " --> pdb=" O GLY A 5 " (cutoff:3.500A) removed outlier: 3.663A pdb=" N LEU A 4 " --> pdb=" O LEU A 942 " (cutoff:3.500A) removed outlier: 3.671A pdb=" N ASP A 938 " --> pdb=" O VAL A 8 " (cutoff:3.500A) Processing sheet with id= B, first strand: chain 'A' and resid 33 through 34 removed outlier: 3.864A pdb=" N ARG A 37 " --> pdb=" O VAL A 34 " (cutoff:3.500A) Processing sheet with id= C, first strand: chain 'A' and resid 70 through 71 Processing sheet with id= D, first strand: chain 'A' and resid 87 through 90 Processing sheet with id= E, first strand: chain 'A' and resid 149 through 150 removed outlier: 6.508A pdb=" N ILE A 871 " --> pdb=" O LEU A 901 " (cutoff:3.500A) removed outlier: 3.831A pdb=" N GLY A 903 " --> pdb=" O ILE A 871 " (cutoff:3.500A) Processing sheet with id= F, first strand: chain 'A' and resid 180 through 184 Processing sheet with id= G, first strand: chain 'A' and resid 409 through 410 Processing sheet with id= H, first strand: chain 'A' and resid 760 through 765 197 hydrogen bonds defined for protein. 528 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.80 Time building geometry restraints manager: 2.46 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.19 - 1.54: 6038 1.54 - 1.90: 179 1.90 - 2.25: 3 2.25 - 2.61: 0 2.61 - 2.96: 1 Bond restraints: 6221 Sorted by residual: bond pdb=" SD MET A 771 " pdb=" CE MET A 771 " ideal model delta sigma weight residual 1.791 2.960 -1.169 2.50e-02 1.60e+03 2.19e+03 bond pdb=" CD2 TYR A 465 " pdb=" CE2 TYR A 465 " ideal model delta sigma weight residual 1.382 2.067 -0.685 3.00e-02 1.11e+03 5.21e+02 bond pdb=" CD1 TYR A 465 " pdb=" CE1 TYR A 465 " ideal model delta sigma weight residual 1.382 2.059 -0.677 3.00e-02 1.11e+03 5.10e+02 bond pdb=" CE1 TYR A 465 " pdb=" CZ TYR A 465 " ideal model delta sigma weight residual 1.378 1.812 -0.434 2.40e-02 1.74e+03 3.27e+02 bond pdb=" CE2 TYR A 465 " pdb=" CZ TYR A 465 " ideal model delta sigma weight residual 1.378 1.801 -0.423 2.40e-02 1.74e+03 3.11e+02 ... (remaining 6216 not shown) Histogram of bond angle deviations from ideal: 95.49 - 105.27: 151 105.27 - 115.06: 3692 115.06 - 124.84: 4446 124.84 - 134.63: 128 134.63 - 144.42: 3 Bond angle restraints: 8420 Sorted by residual: angle pdb=" N HIS A 16 " pdb=" CA HIS A 16 " pdb=" C HIS A 16 " ideal model delta sigma weight residual 111.11 95.49 15.62 1.20e+00 6.94e-01 1.70e+02 angle pdb=" CG MET A 771 " pdb=" SD MET A 771 " pdb=" CE MET A 771 " ideal model delta sigma weight residual 100.90 129.48 -28.58 2.20e+00 2.07e-01 1.69e+02 angle pdb=" C HIS A 653 " pdb=" N PRO A 654 " pdb=" CA PRO A 654 " ideal model delta sigma weight residual 119.84 104.48 15.36 1.25e+00 6.40e-01 1.51e+02 angle pdb=" N VAL A 27 " pdb=" CA VAL A 27 " pdb=" C VAL A 27 " ideal model delta sigma weight residual 112.29 122.07 -9.78 9.40e-01 1.13e+00 1.08e+02 angle pdb=" N GLU A 715 " pdb=" CA GLU A 715 " pdb=" C GLU A 715 " ideal model delta sigma weight residual 113.18 99.89 13.29 1.33e+00 5.65e-01 9.99e+01 ... (remaining 8415 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.68: 3047 17.68 - 35.36: 463 35.36 - 53.04: 191 53.04 - 70.72: 47 70.72 - 88.40: 8 Dihedral angle restraints: 3756 sinusoidal: 1560 harmonic: 2196 Sorted by residual: dihedral pdb=" CA VAL A 27 " pdb=" C VAL A 27 " pdb=" N GLY A 28 " pdb=" CA GLY A 28 " ideal model delta harmonic sigma weight residual 0.00 38.21 -38.21 0 5.00e+00 4.00e-02 5.84e+01 dihedral pdb=" CA HIS A 38 " pdb=" C HIS A 38 " pdb=" N ASP A 39 " pdb=" CA ASP A 39 " ideal model delta harmonic sigma weight residual -180.00 -145.72 -34.28 0 5.00e+00 4.00e-02 4.70e+01 dihedral pdb=" CA SER A 55 " pdb=" C SER A 55 " pdb=" N ILE A 56 " pdb=" CA ILE A 56 " ideal model delta harmonic sigma weight residual -180.00 -148.27 -31.73 0 5.00e+00 4.00e-02 4.03e+01 ... (remaining 3753 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.101: 749 0.101 - 0.203: 160 0.203 - 0.304: 17 0.304 - 0.405: 2 0.405 - 0.506: 1 Chirality restraints: 929 Sorted by residual: chirality pdb=" CA ASN A 63 " pdb=" N ASN A 63 " pdb=" C ASN A 63 " pdb=" CB ASN A 63 " both_signs ideal model delta sigma weight residual False 2.51 2.00 0.51 2.00e-01 2.50e+01 6.41e+00 chirality pdb=" CA ASN A 128 " pdb=" N ASN A 128 " pdb=" C ASN A 128 " pdb=" CB ASN A 128 " both_signs ideal model delta sigma weight residual False 2.51 2.19 0.32 2.00e-01 2.50e+01 2.60e+00 chirality pdb=" CB ILE A 670 " pdb=" CA ILE A 670 " pdb=" CG1 ILE A 670 " pdb=" CG2 ILE A 670 " both_signs ideal model delta sigma weight residual False 2.64 2.34 0.31 2.00e-01 2.50e+01 2.36e+00 ... (remaining 926 not shown) Planarity restraints: 1078 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA PRO A 618 " -0.018 2.00e-02 2.50e+03 3.60e-02 1.29e+01 pdb=" C PRO A 618 " 0.062 2.00e-02 2.50e+03 pdb=" O PRO A 618 " -0.023 2.00e-02 2.50e+03 pdb=" N THR A 619 " -0.022 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C HIS A 653 " 0.057 5.00e-02 4.00e+02 8.58e-02 1.18e+01 pdb=" N PRO A 654 " -0.148 5.00e-02 4.00e+02 pdb=" CA PRO A 654 " 0.036 5.00e-02 4.00e+02 pdb=" CD PRO A 654 " 0.055 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CA HIS A 16 " 0.016 2.00e-02 2.50e+03 3.34e-02 1.11e+01 pdb=" C HIS A 16 " -0.058 2.00e-02 2.50e+03 pdb=" O HIS A 16 " 0.022 2.00e-02 2.50e+03 pdb=" N LEU A 17 " 0.019 2.00e-02 2.50e+03 ... (remaining 1075 not shown) Histogram of nonbonded interaction distances: 2.03 - 2.61: 158 2.61 - 3.18: 5710 3.18 - 3.75: 9704 3.75 - 4.33: 13753 4.33 - 4.90: 21601 Nonbonded interactions: 50926 Sorted by model distance: nonbonded pdb=" O ASN A 47 " pdb=" ND2 ASN A 47 " model vdw 2.034 2.520 nonbonded pdb=" O PHE A 803 " pdb=" N LYS A 805 " model vdw 2.131 2.520 nonbonded pdb=" NE ARG A 74 " pdb=" OD2 ASP A 114 " model vdw 2.162 2.520 nonbonded pdb=" N GLU A 184 " pdb=" OE1 GLU A 184 " model vdw 2.234 2.520 nonbonded pdb=" OG1 THR A 506 " pdb=" N THR A 507 " model vdw 2.247 2.520 ... (remaining 50921 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 Set scattering table Set to: electron Number of scattering types: 5 Type Number sf(0) Gaussians Zn 1 6.06 5 S 32 5.16 5 C 3891 2.51 5 N 1057 2.21 5 O 1105 1.98 5 sf(0) = scattering factor at diffraction angle 0. 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.130 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: 1.170 Check model and map are aligned: 0.100 Convert atoms to be neutral: 0.050 Process input model: 20.970 Find NCS groups from input model: 0.150 Set up NCS constraints: 0.030 Set refine NCS operators: 0.000 Set scattering table: 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.690 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 25.300 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8929 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.027 1.169 6221 Z= 1.945 Angle : 2.041 28.578 8420 Z= 1.287 Chirality : 0.083 0.506 929 Planarity : 0.009 0.086 1078 Dihedral : 20.567 88.397 2342 Min Nonbonded Distance : 2.034 Molprobity Statistics. All-atom Clashscore : 23.82 Ramachandran Plot: Outliers : 2.16 % Allowed : 8.22 % Favored : 89.62 % Rotamer Outliers : 24.92 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.82 % Twisted Proline : 0.00 % Twisted General : 0.84 % 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.64 (0.29), residues: 742 helix: -3.11 (0.25), residues: 222 sheet: -2.97 (0.37), residues: 153 loop : -1.54 (0.34), residues: 367 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 310 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 166 poor density : 144 time to evaluate : 0.805 Fit side-chains TARDY: cannot create tardy model for: "SER A 671 " (corrupted residue). Skipping it. revert: symmetry clash outliers start: 166 outliers final: 72 residues processed: 273 average time/residue: 0.2195 time to fit residues: 74.4331 Evaluate side-chains 183 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 72 poor density : 111 time to evaluate : 0.833 Switching outliers to nearest non-outliers outliers start: 72 outliers final: 0 residues processed: 72 average time/residue: 0.0831 time to fit residues: 10.5651 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 74 random chunks: chunk 62 optimal weight: 0.7980 chunk 56 optimal weight: 0.9990 chunk 31 optimal weight: 0.7980 chunk 19 optimal weight: 1.9990 chunk 37 optimal weight: 3.9990 chunk 29 optimal weight: 9.9990 chunk 58 optimal weight: 2.9990 chunk 22 optimal weight: 2.9990 chunk 35 optimal weight: 0.7980 chunk 43 optimal weight: 0.9990 chunk 67 optimal weight: 0.8980 overall best weight: 0.8582 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 131 HIS ** A 180 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 541 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 587 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 758 GLN A 947 HIS Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8761 moved from start: 0.2926 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.022 1.168 6221 Z= 1.646 Angle : 0.954 28.702 8420 Z= 0.476 Chirality : 0.049 0.163 929 Planarity : 0.005 0.046 1078 Dihedral : 7.603 32.826 833 Min Nonbonded Distance : 2.111 Molprobity Statistics. All-atom Clashscore : 18.89 Ramachandran Plot: Outliers : 0.67 % Allowed : 8.36 % Favored : 90.97 % Rotamer Outliers : 4.35 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.96 % Twisted Proline : 0.00 % Twisted General : 0.28 % 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.28), residues: 742 helix: -1.57 (0.31), residues: 221 sheet: -2.27 (0.37), residues: 150 loop : -1.68 (0.31), residues: 371 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 171 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 29 poor density : 142 time to evaluate : 0.724 Fit side-chains TARDY: cannot create tardy model for: "SER A 671 " (corrupted residue). Skipping it. outliers start: 29 outliers final: 11 residues processed: 164 average time/residue: 0.1833 time to fit residues: 39.2018 Evaluate side-chains 117 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 106 time to evaluate : 0.721 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 0 residues processed: 11 average time/residue: 0.0581 time to fit residues: 2.1838 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 74 random chunks: chunk 37 optimal weight: 7.9990 chunk 20 optimal weight: 5.9990 chunk 55 optimal weight: 6.9990 chunk 45 optimal weight: 3.9990 chunk 18 optimal weight: 10.0000 chunk 67 optimal weight: 4.9990 chunk 72 optimal weight: 4.9990 chunk 59 optimal weight: 6.9990 chunk 66 optimal weight: 7.9990 chunk 22 optimal weight: 1.9990 chunk 54 optimal weight: 0.8980 overall best weight: 3.3788 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 180 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 444 ASN ** A 541 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 587 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 947 HIS Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8820 moved from start: 0.3533 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.023 1.168 6221 Z= 1.681 Angle : 0.961 28.561 8420 Z= 0.476 Chirality : 0.051 0.205 929 Planarity : 0.006 0.077 1078 Dihedral : 7.119 34.903 833 Min Nonbonded Distance : 2.149 Molprobity Statistics. All-atom Clashscore : 17.58 Ramachandran Plot: Outliers : 0.67 % Allowed : 8.49 % Favored : 90.84 % Rotamer Outliers : 6.01 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.96 % Twisted Proline : 0.00 % Twisted General : 0.14 % 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.59 (0.28), residues: 742 helix: -1.27 (0.30), residues: 229 sheet: -2.07 (0.36), residues: 144 loop : -1.79 (0.31), residues: 369 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 150 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 40 poor density : 110 time to evaluate : 0.715 Fit side-chains TARDY: cannot create tardy model for: "SER A 671 " (corrupted residue). Skipping it. outliers start: 40 outliers final: 22 residues processed: 139 average time/residue: 0.1824 time to fit residues: 33.4290 Evaluate side-chains 123 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 22 poor density : 101 time to evaluate : 0.753 Switching outliers to nearest non-outliers outliers start: 22 outliers final: 0 residues processed: 22 average time/residue: 0.0722 time to fit residues: 3.6633 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 74 random chunks: chunk 66 optimal weight: 6.9990 chunk 50 optimal weight: 0.9990 chunk 34 optimal weight: 0.9990 chunk 7 optimal weight: 0.9990 chunk 32 optimal weight: 0.9980 chunk 45 optimal weight: 0.7980 chunk 67 optimal weight: 0.7980 chunk 71 optimal weight: 5.9990 chunk 35 optimal weight: 4.9990 chunk 64 optimal weight: 3.9990 chunk 19 optimal weight: 4.9990 overall best weight: 0.9184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 180 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 444 ASN ** A 541 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 587 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8721 moved from start: 0.4146 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.022 1.164 6221 Z= 1.634 Angle : 0.816 28.622 8420 Z= 0.401 Chirality : 0.046 0.149 929 Planarity : 0.004 0.047 1078 Dihedral : 6.363 34.802 833 Min Nonbonded Distance : 2.178 Molprobity Statistics. All-atom Clashscore : 16.35 Ramachandran Plot: Outliers : 0.67 % Allowed : 7.68 % Favored : 91.64 % Rotamer Outliers : 3.45 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.96 % Twisted Proline : 0.00 % Twisted General : 0.14 % 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.43 (0.28), residues: 742 helix: -0.67 (0.33), residues: 228 sheet: -2.19 (0.38), residues: 144 loop : -1.94 (0.30), residues: 370 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 149 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 23 poor density : 126 time to evaluate : 0.742 Fit side-chains TARDY: cannot create tardy model for: "SER A 671 " (corrupted residue). Skipping it. outliers start: 23 outliers final: 7 residues processed: 142 average time/residue: 0.1784 time to fit residues: 33.5703 Evaluate side-chains 114 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 107 time to evaluate : 0.744 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.0603 time to fit residues: 1.8093 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 74 random chunks: chunk 59 optimal weight: 3.9990 chunk 40 optimal weight: 4.9990 chunk 1 optimal weight: 0.0870 chunk 53 optimal weight: 5.9990 chunk 29 optimal weight: 0.8980 chunk 61 optimal weight: 4.9990 chunk 49 optimal weight: 4.9990 chunk 0 optimal weight: 5.9990 chunk 36 optimal weight: 7.9990 chunk 64 optimal weight: 2.9990 chunk 18 optimal weight: 10.0000 overall best weight: 2.5964 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 16 HIS ** A 180 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 444 ASN A 474 ASN ** A 541 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8789 moved from start: 0.4200 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.022 1.165 6221 Z= 1.660 Angle : 0.869 28.507 8420 Z= 0.429 Chirality : 0.048 0.172 929 Planarity : 0.004 0.040 1078 Dihedral : 6.446 34.635 833 Min Nonbonded Distance : 2.212 Molprobity Statistics. All-atom Clashscore : 15.69 Ramachandran Plot: Outliers : 0.67 % Allowed : 8.63 % Favored : 90.70 % Rotamer Outliers : 3.45 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.96 % Twisted Proline : 0.00 % Twisted General : 0.14 % 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.38 (0.29), residues: 742 helix: -0.78 (0.33), residues: 234 sheet: -2.07 (0.38), residues: 144 loop : -1.85 (0.31), residues: 364 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 130 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 23 poor density : 107 time to evaluate : 0.782 Fit side-chains TARDY: cannot create tardy model for: "SER A 671 " (corrupted residue). Skipping it. outliers start: 23 outliers final: 10 residues processed: 121 average time/residue: 0.1984 time to fit residues: 31.2600 Evaluate side-chains 105 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 95 time to evaluate : 0.727 Switching outliers to nearest non-outliers outliers start: 10 outliers final: 0 residues processed: 10 average time/residue: 0.0613 time to fit residues: 2.1053 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 74 random chunks: chunk 24 optimal weight: 4.9990 chunk 64 optimal weight: 0.9990 chunk 14 optimal weight: 2.9990 chunk 41 optimal weight: 2.9990 chunk 17 optimal weight: 3.9990 chunk 71 optimal weight: 6.9990 chunk 59 optimal weight: 0.6980 chunk 33 optimal weight: 6.9990 chunk 5 optimal weight: 6.9990 chunk 23 optimal weight: 4.9990 chunk 37 optimal weight: 5.9990 overall best weight: 2.3388 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 180 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 474 ASN ** A 541 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8782 moved from start: 0.4307 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.022 1.162 6221 Z= 1.657 Angle : 0.840 29.176 8420 Z= 0.414 Chirality : 0.048 0.165 929 Planarity : 0.005 0.047 1078 Dihedral : 6.381 35.819 833 Min Nonbonded Distance : 2.199 Molprobity Statistics. All-atom Clashscore : 15.44 Ramachandran Plot: Outliers : 0.67 % Allowed : 8.49 % Favored : 90.84 % Rotamer Outliers : 2.25 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.96 % Twisted Proline : 0.00 % Twisted General : 0.14 % 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.19 (0.29), residues: 742 helix: -0.74 (0.33), residues: 234 sheet: -1.55 (0.40), residues: 134 loop : -1.86 (0.31), residues: 374 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 119 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 104 time to evaluate : 0.734 Fit side-chains TARDY: cannot create tardy model for: "SER A 671 " (corrupted residue). Skipping it. outliers start: 15 outliers final: 8 residues processed: 115 average time/residue: 0.1853 time to fit residues: 28.5159 Evaluate side-chains 106 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 98 time to evaluate : 0.772 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 0 residues processed: 8 average time/residue: 0.0642 time to fit residues: 1.9374 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 74 random chunks: chunk 69 optimal weight: 0.6980 chunk 8 optimal weight: 4.9990 chunk 40 optimal weight: 1.9990 chunk 52 optimal weight: 3.9990 chunk 60 optimal weight: 1.9990 chunk 39 optimal weight: 2.9990 chunk 71 optimal weight: 6.9990 chunk 44 optimal weight: 0.1980 chunk 43 optimal weight: 6.9990 chunk 32 optimal weight: 9.9990 chunk 28 optimal weight: 0.0020 overall best weight: 0.9792 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 180 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 474 ASN ** A 541 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8714 moved from start: 0.4685 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.022 1.160 6221 Z= 1.634 Angle : 0.786 28.675 8420 Z= 0.384 Chirality : 0.045 0.164 929 Planarity : 0.004 0.041 1078 Dihedral : 5.997 33.986 833 Min Nonbonded Distance : 2.195 Molprobity Statistics. All-atom Clashscore : 16.59 Ramachandran Plot: Outliers : 0.67 % Allowed : 7.82 % Favored : 91.51 % Rotamer Outliers : 2.25 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.96 % Twisted Proline : 0.00 % Twisted General : 0.14 % 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.93 (0.29), residues: 742 helix: -0.35 (0.34), residues: 230 sheet: -1.36 (0.45), residues: 112 loop : -1.84 (0.29), residues: 400 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 128 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 113 time to evaluate : 0.767 Fit side-chains TARDY: cannot create tardy model for: "SER A 671 " (corrupted residue). Skipping it. outliers start: 15 outliers final: 8 residues processed: 121 average time/residue: 0.1961 time to fit residues: 31.3129 Evaluate side-chains 107 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 99 time to evaluate : 0.750 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 0 residues processed: 8 average time/residue: 0.0634 time to fit residues: 2.0113 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 74 random chunks: chunk 42 optimal weight: 5.9990 chunk 21 optimal weight: 8.9990 chunk 14 optimal weight: 0.0170 chunk 13 optimal weight: 3.9990 chunk 45 optimal weight: 5.9990 chunk 48 optimal weight: 1.9990 chunk 35 optimal weight: 3.9990 chunk 6 optimal weight: 0.7980 chunk 56 optimal weight: 0.6980 chunk 64 optimal weight: 0.6980 chunk 68 optimal weight: 4.9990 overall best weight: 0.8420 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 180 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 541 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8693 moved from start: 0.4859 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.022 1.159 6221 Z= 1.634 Angle : 0.787 28.578 8420 Z= 0.386 Chirality : 0.045 0.164 929 Planarity : 0.004 0.041 1078 Dihedral : 5.882 33.836 833 Min Nonbonded Distance : 2.196 Molprobity Statistics. All-atom Clashscore : 16.76 Ramachandran Plot: Outliers : 0.67 % Allowed : 7.68 % Favored : 91.64 % Rotamer Outliers : 0.75 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.96 % Twisted Proline : 0.00 % Twisted General : 0.14 % 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.98 (0.30), residues: 742 helix: -0.24 (0.35), residues: 231 sheet: -1.79 (0.44), residues: 124 loop : -1.82 (0.30), residues: 387 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 111 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 106 time to evaluate : 0.733 Fit side-chains TARDY: cannot create tardy model for: "SER A 671 " (corrupted residue). Skipping it. outliers start: 5 outliers final: 3 residues processed: 109 average time/residue: 0.1869 time to fit residues: 26.9954 Evaluate side-chains 105 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 102 time to evaluate : 0.702 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0587 time to fit residues: 1.3347 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 74 random chunks: chunk 62 optimal weight: 2.9990 chunk 66 optimal weight: 1.9990 chunk 68 optimal weight: 3.9990 chunk 40 optimal weight: 6.9990 chunk 28 optimal weight: 0.8980 chunk 52 optimal weight: 2.9990 chunk 20 optimal weight: 7.9990 chunk 60 optimal weight: 0.7980 chunk 43 optimal weight: 0.8980 chunk 70 optimal weight: 0.8980 chunk 42 optimal weight: 2.9990 overall best weight: 1.0982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 180 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 541 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8696 moved from start: 0.4935 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.022 1.161 6221 Z= 1.637 Angle : 0.781 28.514 8420 Z= 0.384 Chirality : 0.045 0.174 929 Planarity : 0.004 0.042 1078 Dihedral : 5.842 34.172 833 Min Nonbonded Distance : 2.195 Molprobity Statistics. All-atom Clashscore : 16.27 Ramachandran Plot: Outliers : 0.67 % Allowed : 7.28 % Favored : 92.05 % Rotamer Outliers : 0.75 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.96 % Twisted Proline : 0.00 % Twisted General : 0.14 % 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.88 (0.30), residues: 742 helix: -0.14 (0.35), residues: 232 sheet: -1.62 (0.45), residues: 123 loop : -1.82 (0.30), residues: 387 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 110 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 105 time to evaluate : 0.784 Fit side-chains TARDY: cannot create tardy model for: "SER A 671 " (corrupted residue). Skipping it. outliers start: 5 outliers final: 2 residues processed: 109 average time/residue: 0.1794 time to fit residues: 26.2928 Evaluate side-chains 100 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 98 time to evaluate : 0.741 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0596 time to fit residues: 1.2515 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 74 random chunks: chunk 33 optimal weight: 0.4980 chunk 48 optimal weight: 4.9990 chunk 73 optimal weight: 4.9990 chunk 67 optimal weight: 4.9990 chunk 58 optimal weight: 0.0570 chunk 6 optimal weight: 4.9990 chunk 45 optimal weight: 3.9990 chunk 35 optimal weight: 7.9990 chunk 46 optimal weight: 1.9990 chunk 62 optimal weight: 0.9980 chunk 17 optimal weight: 0.3980 overall best weight: 0.7900 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 135 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 180 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 541 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8669 moved from start: 0.5060 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.022 1.160 6221 Z= 1.634 Angle : 0.769 28.529 8420 Z= 0.378 Chirality : 0.044 0.163 929 Planarity : 0.004 0.045 1078 Dihedral : 5.721 33.674 833 Min Nonbonded Distance : 2.196 Molprobity Statistics. All-atom Clashscore : 16.59 Ramachandran Plot: Outliers : 0.67 % Allowed : 7.41 % Favored : 91.91 % Rotamer Outliers : 0.45 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.96 % Twisted Proline : 0.00 % Twisted General : 0.14 % 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.81 (0.30), residues: 742 helix: -0.02 (0.35), residues: 232 sheet: -1.57 (0.46), residues: 122 loop : -1.83 (0.29), residues: 388 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1484 Ramachandran restraints generated. 742 Oldfield, 0 Emsley, 742 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 110 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 107 time to evaluate : 0.668 Fit side-chains revert: symmetry clash TARDY: cannot create tardy model for: "SER A 671 " (corrupted residue). Skipping it. outliers start: 3 outliers final: 3 residues processed: 109 average time/residue: 0.1734 time to fit residues: 25.3613 Evaluate side-chains 100 residues out of total 666 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 97 time to evaluate : 0.715 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0704 time to fit residues: 1.3634 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 74 random chunks: chunk 54 optimal weight: 0.9990 chunk 8 optimal weight: 4.9990 chunk 16 optimal weight: 0.9990 chunk 58 optimal weight: 0.3980 chunk 24 optimal weight: 5.9990 chunk 60 optimal weight: 2.9990 chunk 7 optimal weight: 0.2980 chunk 10 optimal weight: 1.9990 chunk 51 optimal weight: 1.9990 chunk 3 optimal weight: 0.3980 chunk 42 optimal weight: 0.6980 overall best weight: 0.5582 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 135 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 180 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 541 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 562 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3222 r_free = 0.3222 target = 0.113776 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 39)----------------| | r_work = 0.3015 r_free = 0.3015 target = 0.098056 restraints weight = 8604.674| |-----------------------------------------------------------------------------| r_work (start): 0.3018 rms_B_bonded: 1.59 r_work: 0.3006 rms_B_bonded: 1.42 restraints_weight: 0.5000 r_work: 0.2995 rms_B_bonded: 1.40 restraints_weight: 0.2500 r_work: 0.2984 rms_B_bonded: 1.45 restraints_weight: 0.1250 r_work: 0.2973 rms_B_bonded: 1.53 restraints_weight: 0.0625 r_work: 0.2960 rms_B_bonded: 1.66 restraints_weight: 0.0312 r_work: 0.2946 rms_B_bonded: 1.82 restraints_weight: 0.0156 r_work: 0.2931 rms_B_bonded: 2.03 restraints_weight: 0.0078 r_work: 0.2913 rms_B_bonded: 2.28 restraints_weight: 0.0039 r_work: 0.2894 rms_B_bonded: 2.57 restraints_weight: 0.0020 r_work: 0.2872 rms_B_bonded: 2.92 restraints_weight: 0.0010 r_work: 0.2848 rms_B_bonded: 3.34 restraints_weight: 0.0005 r_work: 0.2821 rms_B_bonded: 3.83 restraints_weight: 0.0002 r_work: 0.2789 rms_B_bonded: 4.41 restraints_weight: 0.0001 r_work (final): 0.2789 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8450 moved from start: 0.5216 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.022 1.160 6221 Z= 1.633 Angle : 0.765 28.512 8420 Z= 0.376 Chirality : 0.044 0.178 929 Planarity : 0.004 0.045 1078 Dihedral : 5.573 32.467 833 Min Nonbonded Distance : 2.179 Molprobity Statistics. All-atom Clashscore : 16.27 Ramachandran Plot: Outliers : 0.81 % Allowed : 7.01 % Favored : 92.18 % Rotamer Outliers : 0.15 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 6.67 % Cis-general : 1.96 % Twisted Proline : 0.00 % Twisted General : 0.14 % 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.73 (0.30), residues: 742 helix: 0.02 (0.35), residues: 232 sheet: -1.36 (0.47), residues: 121 loop : -1.82 (0.30), residues: 389 =============================================================================== Job complete usr+sys time: 1514.95 seconds wall clock time: 28 minutes 6.87 seconds (1686.87 seconds total)