Starting phenix.real_space_refine (version: dev) on Thu Jun 2 21:11:02 2022 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xd1_33132/06_2022/7xd1_33132.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xd1_33132/06_2022/7xd1_33132.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.2 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xd1_33132/06_2022/7xd1_33132.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xd1_33132/06_2022/7xd1_33132.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xd1_33132/06_2022/7xd1_33132.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7xd1_33132/06_2022/7xd1_33132.pdb" } resolution = 3.2 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.000 sd= 0.003 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 "C PHE 25": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C TYR 50": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "E PHE 84": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "G PHE 25": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "G TYR 50": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4617/modules/chem_data/mon_lib" Total number of atoms: 12048 Number of models: 1 Model: "" Number of chains: 10 Chain: "A" Number of atoms: 810 Number of conformers: 1 Conformer: "" Number of residues, atoms: 98, 810 Classifications: {'peptide': 98} Link IDs: {'PTRANS': 4, 'TRANS': 93} Chain: "B" Number of atoms: 641 Number of conformers: 1 Conformer: "" Number of residues, atoms: 80, 641 Classifications: {'peptide': 80} Link IDs: {'PTRANS': 1, 'TRANS': 78} Chain: "C" Number of atoms: 836 Number of conformers: 1 Conformer: "" Number of residues, atoms: 109, 836 Classifications: {'peptide': 109} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 5, 'TRANS': 103} Unresolved non-hydrogen bonds: 4 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 4 Chain: "E" Number of atoms: 801 Number of conformers: 1 Conformer: "" Number of residues, atoms: 97, 801 Classifications: {'peptide': 97} Link IDs: {'PTRANS': 3, 'TRANS': 93} Chain: "F" Number of atoms: 627 Number of conformers: 1 Conformer: "" Number of residues, atoms: 79, 627 Classifications: {'peptide': 79} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 1, 'TRANS': 77} Unresolved non-hydrogen bonds: 6 Unresolved non-hydrogen angles: 7 Unresolved non-hydrogen dihedrals: 5 Planarities with less than four sites: {'ARG:plan': 1} Unresolved non-hydrogen planarities: 5 Chain: "G" Number of atoms: 836 Number of conformers: 1 Conformer: "" Number of residues, atoms: 109, 836 Classifications: {'peptide': 109} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 5, 'TRANS': 103} Unresolved non-hydrogen bonds: 4 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 4 Chain: "I" Number of atoms: 3031 Number of conformers: 1 Conformer: "" Number of residues, atoms: 147, 3031 Classifications: {'DNA': 147} Link IDs: {'rna3p': 146} Chain: "J" Number of atoms: 2996 Number of conformers: 1 Conformer: "" Number of residues, atoms: 147, 2996 Classifications: {'DNA': 147} Link IDs: {'rna3p': 146} Chain: "D" Number of atoms: 735 Number of conformers: 1 Conformer: "" Number of residues, atoms: 94, 735 Classifications: {'peptide': 94} Link IDs: {'PTRANS': 2, 'TRANS': 91} Chain: "H" Number of atoms: 735 Number of conformers: 1 Conformer: "" Number of residues, atoms: 94, 735 Classifications: {'peptide': 94} Link IDs: {'PTRANS': 2, 'TRANS': 91} Time building chain proxies: 7.54, per 1000 atoms: 0.63 Number of scatterers: 12048 At special positions: 0 Unit cell: (119.214, 127.806, 91.29, 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) S 14 16.00 P 294 15.00 O 2824 8.00 N 2268 7.00 C 6648 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: 7.38 Conformation dependent library (CDL) restraints added in 1.6 seconds 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1420 Finding SS restraints... Secondary structure from input PDB file: 35 helices and 10 sheets defined 68.8% alpha, 3.2% beta 129 base pairs and 276 stacking pairs defined. Time for finding SS restraints: 6.08 Creating SS restraints... Processing helix chain 'A' and resid 44 through 55 removed outlier: 3.927A pdb=" N TYR A 54 " --> pdb=" O GLU A 50 " (cutoff:3.500A) Processing helix chain 'A' and resid 63 through 76 Processing helix chain 'A' and resid 85 through 114 removed outlier: 3.525A pdb=" N ALA A 98 " --> pdb=" O GLU A 94 " (cutoff:3.500A) removed outlier: 3.531A pdb=" N LEU A 109 " --> pdb=" O GLU A 105 " (cutoff:3.500A) Processing helix chain 'A' and resid 120 through 131 removed outlier: 3.847A pdb=" N ILE A 124 " --> pdb=" O MET A 120 " (cutoff:3.500A) removed outlier: 3.755A pdb=" N GLN A 125 " --> pdb=" O PRO A 121 " (cutoff:3.500A) removed outlier: 3.538A pdb=" N ARG A 131 " --> pdb=" O ALA A 127 " (cutoff:3.500A) Processing helix chain 'B' and resid 25 through 29 Processing helix chain 'B' and resid 30 through 41 removed outlier: 3.581A pdb=" N GLY B 41 " --> pdb=" O LEU B 37 " (cutoff:3.500A) Processing helix chain 'B' and resid 49 through 76 removed outlier: 4.132A pdb=" N GLU B 53 " --> pdb=" O LEU B 49 " (cutoff:3.500A) removed outlier: 4.137A pdb=" N ARG B 67 " --> pdb=" O GLU B 63 " (cutoff:3.500A) removed outlier: 3.563A pdb=" N ASP B 68 " --> pdb=" O ASN B 64 " (cutoff:3.500A) removed outlier: 3.646A pdb=" N THR B 71 " --> pdb=" O ARG B 67 " (cutoff:3.500A) removed outlier: 3.677A pdb=" N ALA B 76 " --> pdb=" O TYR B 72 " (cutoff:3.500A) Processing helix chain 'B' and resid 82 through 92 removed outlier: 3.637A pdb=" N VAL B 86 " --> pdb=" O THR B 82 " (cutoff:3.500A) Processing helix chain 'C' and resid 16 through 21 removed outlier: 3.641A pdb=" N ALA C 21 " --> pdb=" O ARG C 17 " (cutoff:3.500A) Processing helix chain 'C' and resid 26 through 36 removed outlier: 3.737A pdb=" N VAL C 30 " --> pdb=" O PRO C 26 " (cutoff:3.500A) Processing helix chain 'C' and resid 46 through 72 removed outlier: 3.535A pdb=" N ALA C 53 " --> pdb=" O VAL C 49 " (cutoff:3.500A) Processing helix chain 'C' and resid 79 through 90 removed outlier: 3.646A pdb=" N LEU C 83 " --> pdb=" O ILE C 79 " (cutoff:3.500A) Processing helix chain 'C' and resid 90 through 97 removed outlier: 3.639A pdb=" N LYS C 95 " --> pdb=" O GLU C 91 " (cutoff:3.500A) Processing helix chain 'E' and resid 45 through 55 Processing helix chain 'E' and resid 63 through 75 removed outlier: 3.588A pdb=" N PHE E 67 " --> pdb=" O ARG E 63 " (cutoff:3.500A) removed outlier: 3.610A pdb=" N LEU E 70 " --> pdb=" O PRO E 66 " (cutoff:3.500A) removed outlier: 3.526A pdb=" N GLU E 73 " --> pdb=" O ARG E 69 " (cutoff:3.500A) Processing helix chain 'E' and resid 76 through 78 No H-bonds generated for 'chain 'E' and resid 76 through 78' Processing helix chain 'E' and resid 85 through 114 removed outlier: 3.565A pdb=" N ALA E 111 " --> pdb=" O THR E 107 " (cutoff:3.500A) Processing helix chain 'E' and resid 120 through 131 removed outlier: 3.845A pdb=" N ILE E 124 " --> pdb=" O MET E 120 " (cutoff:3.500A) Processing helix chain 'F' and resid 25 through 29 Processing helix chain 'F' and resid 30 through 42 Processing helix chain 'F' and resid 50 through 76 removed outlier: 3.999A pdb=" N ARG F 55 " --> pdb=" O TYR F 51 " (cutoff:3.500A) removed outlier: 3.960A pdb=" N GLY F 56 " --> pdb=" O GLU F 52 " (cutoff:3.500A) removed outlier: 3.531A pdb=" N LYS F 59 " --> pdb=" O ARG F 55 " (cutoff:3.500A) removed outlier: 3.848A pdb=" N ARG F 67 " --> pdb=" O GLU F 63 " (cutoff:3.500A) removed outlier: 3.756A pdb=" N ASP F 68 " --> pdb=" O ASN F 64 " (cutoff:3.500A) removed outlier: 3.789A pdb=" N THR F 71 " --> pdb=" O ARG F 67 " (cutoff:3.500A) removed outlier: 3.551A pdb=" N HIS F 75 " --> pdb=" O THR F 71 " (cutoff:3.500A) Processing helix chain 'F' and resid 82 through 92 removed outlier: 3.749A pdb=" N VAL F 86 " --> pdb=" O THR F 82 " (cutoff:3.500A) removed outlier: 3.557A pdb=" N ARG F 92 " --> pdb=" O TYR F 88 " (cutoff:3.500A) Processing helix chain 'G' and resid 16 through 21 removed outlier: 3.641A pdb=" N ALA G 21 " --> pdb=" O ARG G 17 " (cutoff:3.500A) Processing helix chain 'G' and resid 26 through 36 removed outlier: 3.737A pdb=" N VAL G 30 " --> pdb=" O PRO G 26 " (cutoff:3.500A) Processing helix chain 'G' and resid 46 through 72 removed outlier: 3.535A pdb=" N ALA G 53 " --> pdb=" O VAL G 49 " (cutoff:3.500A) Processing helix chain 'G' and resid 79 through 90 removed outlier: 3.646A pdb=" N LEU G 83 " --> pdb=" O ILE G 79 " (cutoff:3.500A) Processing helix chain 'G' and resid 90 through 97 removed outlier: 3.639A pdb=" N LYS G 95 " --> pdb=" O GLU G 91 " (cutoff:3.500A) Processing helix chain 'D' and resid 37 through 49 removed outlier: 3.572A pdb=" N LYS D 46 " --> pdb=" O TYR D 42 " (cutoff:3.500A) Processing helix chain 'D' and resid 58 through 84 removed outlier: 3.640A pdb=" N ARG D 79 " --> pdb=" O GLY D 75 " (cutoff:3.500A) Processing helix chain 'D' and resid 91 through 102 Processing helix chain 'D' and resid 104 through 124 removed outlier: 3.580A pdb=" N LYS D 108 " --> pdb=" O GLY D 104 " (cutoff:3.500A) removed outlier: 3.599A pdb=" N ALA D 110 " --> pdb=" O LEU D 106 " (cutoff:3.500A) removed outlier: 3.541A pdb=" N SER D 112 " --> pdb=" O LYS D 108 " (cutoff:3.500A) removed outlier: 3.663A pdb=" N ALA D 117 " --> pdb=" O GLU D 113 " (cutoff:3.500A) Processing helix chain 'H' and resid 37 through 49 removed outlier: 3.572A pdb=" N LYS H 46 " --> pdb=" O TYR H 42 " (cutoff:3.500A) Processing helix chain 'H' and resid 58 through 84 removed outlier: 3.640A pdb=" N ARG H 79 " --> pdb=" O GLY H 75 " (cutoff:3.500A) Processing helix chain 'H' and resid 91 through 102 Processing helix chain 'H' and resid 104 through 124 removed outlier: 3.580A pdb=" N LYS H 108 " --> pdb=" O GLY H 104 " (cutoff:3.500A) removed outlier: 3.599A pdb=" N ALA H 110 " --> pdb=" O LEU H 106 " (cutoff:3.500A) removed outlier: 3.541A pdb=" N SER H 112 " --> pdb=" O LYS H 108 " (cutoff:3.500A) removed outlier: 3.663A pdb=" N ALA H 117 " --> pdb=" O GLU H 113 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 83 through 84 removed outlier: 6.907A pdb=" N ARG A 83 " --> pdb=" O VAL B 81 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA1 Processing sheet with id=AA2, first strand: chain 'A' and resid 118 through 119 Processing sheet with id=AA3, first strand: chain 'B' and resid 96 through 98 removed outlier: 6.597A pdb=" N THR B 96 " --> pdb=" O THR G 101 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA3 Processing sheet with id=AA4, first strand: chain 'C' and resid 42 through 43 removed outlier: 7.291A pdb=" N ARG C 42 " --> pdb=" O ILE D 89 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA4 Processing sheet with id=AA5, first strand: chain 'C' and resid 77 through 78 Processing sheet with id=AA6, first strand: chain 'C' and resid 100 through 102 Processing sheet with id=AA7, first strand: chain 'E' and resid 83 through 84 removed outlier: 7.013A pdb=" N ARG E 83 " --> pdb=" O VAL F 81 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA7 Processing sheet with id=AA8, first strand: chain 'E' and resid 118 through 119 Processing sheet with id=AA9, first strand: chain 'G' and resid 42 through 43 removed outlier: 7.220A pdb=" N ARG G 42 " --> pdb=" O ILE H 89 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA9 Processing sheet with id=AB1, first strand: chain 'G' and resid 77 through 78 343 hydrogen bonds defined for protein. 1017 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 317 hydrogen bonds 630 hydrogen bond angles 0 basepair planarities 129 basepair parallelities 276 stacking parallelities Total time for adding SS restraints: 4.64 Time building geometry restraints manager: 8.55 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: 2152 1.33 - 1.45: 4325 1.45 - 1.57: 5772 1.57 - 1.69: 588 1.69 - 1.81: 24 Bond restraints: 12861 Sorted by residual: bond pdb=" CB GLU D 113 " pdb=" CG GLU D 113 " ideal model delta sigma weight residual 1.520 1.598 -0.078 3.00e-02 1.11e+03 6.78e+00 bond pdb=" CB GLU H 113 " pdb=" CG GLU H 113 " ideal model delta sigma weight residual 1.520 1.597 -0.077 3.00e-02 1.11e+03 6.66e+00 bond pdb=" C3' DT I -59 " pdb=" O3' DT I -59 " ideal model delta sigma weight residual 1.422 1.492 -0.070 3.00e-02 1.11e+03 5.42e+00 bond pdb=" C3' DG I 20 " pdb=" O3' DG I 20 " ideal model delta sigma weight residual 1.422 1.488 -0.066 3.00e-02 1.11e+03 4.90e+00 bond pdb=" CG1 ILE B 50 " pdb=" CD1 ILE B 50 " ideal model delta sigma weight residual 1.513 1.435 0.078 3.90e-02 6.57e+02 4.04e+00 ... (remaining 12856 not shown) Histogram of bond angle deviations from ideal: 98.28 - 105.18: 1307 105.18 - 112.07: 6852 112.07 - 118.96: 3540 118.96 - 125.86: 5968 125.86 - 132.75: 962 Bond angle restraints: 18629 Sorted by residual: angle pdb=" CB GLU D 113 " pdb=" CG GLU D 113 " pdb=" CD GLU D 113 " ideal model delta sigma weight residual 112.60 123.17 -10.57 1.70e+00 3.46e-01 3.87e+01 angle pdb=" CB GLU H 113 " pdb=" CG GLU H 113 " pdb=" CD GLU H 113 " ideal model delta sigma weight residual 112.60 123.16 -10.56 1.70e+00 3.46e-01 3.86e+01 angle pdb=" CA GLU H 113 " pdb=" CB GLU H 113 " pdb=" CG GLU H 113 " ideal model delta sigma weight residual 114.10 125.26 -11.16 2.00e+00 2.50e-01 3.11e+01 angle pdb=" CA GLU D 113 " pdb=" CB GLU D 113 " pdb=" CG GLU D 113 " ideal model delta sigma weight residual 114.10 125.19 -11.09 2.00e+00 2.50e-01 3.08e+01 angle pdb=" C THR E 80 " pdb=" N ASP E 81 " pdb=" CA ASP E 81 " ideal model delta sigma weight residual 121.54 131.03 -9.49 1.91e+00 2.74e-01 2.47e+01 ... (remaining 18624 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 35.23: 5367 35.23 - 70.45: 1342 70.45 - 105.68: 8 105.68 - 140.90: 2 140.90 - 176.13: 2 Dihedral angle restraints: 6721 sinusoidal: 4519 harmonic: 2202 Sorted by residual: dihedral pdb=" CA GLU A 50 " pdb=" C GLU A 50 " pdb=" N ILE A 51 " pdb=" CA ILE A 51 " ideal model delta harmonic sigma weight residual 180.00 153.63 26.37 0 5.00e+00 4.00e-02 2.78e+01 dihedral pdb=" CA PRO E 121 " pdb=" C PRO E 121 " pdb=" N LYS E 122 " pdb=" CA LYS E 122 " ideal model delta harmonic sigma weight residual 180.00 158.63 21.37 0 5.00e+00 4.00e-02 1.83e+01 dihedral pdb=" CA ILE F 29 " pdb=" C ILE F 29 " pdb=" N THR F 30 " pdb=" CA THR F 30 " ideal model delta harmonic sigma weight residual 180.00 -159.33 -20.67 0 5.00e+00 4.00e-02 1.71e+01 ... (remaining 6718 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.051: 1523 0.051 - 0.102: 513 0.102 - 0.153: 71 0.153 - 0.204: 9 0.204 - 0.255: 3 Chirality restraints: 2119 Sorted by residual: chirality pdb=" CA ASP E 81 " pdb=" N ASP E 81 " pdb=" C ASP E 81 " pdb=" CB ASP E 81 " both_signs ideal model delta sigma weight residual False 2.51 2.25 0.26 2.00e-01 2.50e+01 1.63e+00 chirality pdb=" CA ASP A 81 " pdb=" N ASP A 81 " pdb=" C ASP A 81 " pdb=" CB ASP A 81 " both_signs ideal model delta sigma weight residual False 2.51 2.27 0.24 2.00e-01 2.50e+01 1.49e+00 chirality pdb=" P DC I -58 " pdb=" OP1 DC I -58 " pdb=" OP2 DC I -58 " pdb=" O5' DC I -58 " both_signs ideal model delta sigma weight residual True 2.35 -2.56 -0.22 2.00e-01 2.50e+01 1.16e+00 ... (remaining 2116 not shown) Planarity restraints: 1340 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB TYR F 88 " 0.017 2.00e-02 2.50e+03 1.83e-02 6.68e+00 pdb=" CG TYR F 88 " -0.045 2.00e-02 2.50e+03 pdb=" CD1 TYR F 88 " 0.017 2.00e-02 2.50e+03 pdb=" CD2 TYR F 88 " 0.008 2.00e-02 2.50e+03 pdb=" CE1 TYR F 88 " -0.004 2.00e-02 2.50e+03 pdb=" CE2 TYR F 88 " 0.004 2.00e-02 2.50e+03 pdb=" CZ TYR F 88 " 0.002 2.00e-02 2.50e+03 pdb=" OH TYR F 88 " 0.000 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TYR B 51 " -0.016 2.00e-02 2.50e+03 1.54e-02 4.72e+00 pdb=" CG TYR B 51 " 0.038 2.00e-02 2.50e+03 pdb=" CD1 TYR B 51 " -0.008 2.00e-02 2.50e+03 pdb=" CD2 TYR B 51 " -0.004 2.00e-02 2.50e+03 pdb=" CE1 TYR B 51 " -0.003 2.00e-02 2.50e+03 pdb=" CE2 TYR B 51 " -0.007 2.00e-02 2.50e+03 pdb=" CZ TYR B 51 " -0.002 2.00e-02 2.50e+03 pdb=" OH TYR B 51 " 0.003 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE C 25 " -0.012 2.00e-02 2.50e+03 1.46e-02 3.71e+00 pdb=" CG PHE C 25 " 0.033 2.00e-02 2.50e+03 pdb=" CD1 PHE C 25 " -0.010 2.00e-02 2.50e+03 pdb=" CD2 PHE C 25 " -0.011 2.00e-02 2.50e+03 pdb=" CE1 PHE C 25 " 0.001 2.00e-02 2.50e+03 pdb=" CE2 PHE C 25 " 0.001 2.00e-02 2.50e+03 pdb=" CZ PHE C 25 " -0.003 2.00e-02 2.50e+03 ... (remaining 1337 not shown) Histogram of nonbonded interaction distances: 2.34 - 2.85: 3453 2.85 - 3.37: 9714 3.37 - 3.88: 22898 3.88 - 4.39: 28712 4.39 - 4.90: 39523 Nonbonded interactions: 104300 Sorted by model distance: nonbonded pdb=" OG1 THR D 90 " pdb=" OE1 GLU D 93 " model vdw 2.343 2.440 nonbonded pdb=" OG1 THR H 90 " pdb=" OE1 GLU H 93 " model vdw 2.344 2.440 nonbonded pdb=" NE2 GLN A 93 " pdb=" OE2 GLU A 97 " model vdw 2.345 2.520 nonbonded pdb=" O6 DG I -21 " pdb=" N4 DC J 21 " model vdw 2.381 2.520 nonbonded pdb=" O6 DG I -10 " pdb=" N4 DC J 10 " model vdw 2.391 2.520 ... (remaining 104295 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = (chain 'A' and resid 38 through 134) selection = chain 'E' } ncs_group { reference = (chain 'B' and ((resid 23 and (name N or name CA or name C or name O or name CB \ )) or resid 24 through 101)) selection = chain 'F' } ncs_group { reference = chain 'C' selection = chain 'G' } ncs_group { reference = chain 'D' selection = chain 'H' } 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 P 294 5.49 5 S 14 5.16 5 C 6648 2.51 5 N 2268 2.21 5 O 2824 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.810 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: 4.550 Check model and map are aligned: 0.190 Convert atoms to be neutral: 0.110 Process input model: 46.000 Find NCS groups from input model: 0.460 Set up NCS constraints: 0.080 Set refine NCS operators: 0.000 Set scattering table: 0.010 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.010 Load rotamer database and sin/cos tables:3.570 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 55.800 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8253 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.009 0.080 12861 Z= 0.524 Angle : 0.871 11.158 18629 Z= 0.501 Chirality : 0.049 0.255 2119 Planarity : 0.005 0.044 1340 Dihedral : 26.969 176.131 5301 Min Nonbonded Distance : 2.343 Molprobity Statistics. All-atom Clashscore : 2.86 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.34 % Favored : 98.66 % Rotamer Outliers : 1.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: -5.83 (0.19), residues: 744 helix: -4.27 (0.11), residues: 524 sheet: None (None), residues: 0 loop : -2.04 (0.34), residues: 220 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 195 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 185 time to evaluate : 0.876 Fit side-chains outliers start: 10 outliers final: 0 residues processed: 189 average time/residue: 0.3845 time to fit residues: 91.6750 Evaluate side-chains 128 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 128 time to evaluate : 0.874 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.0185 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 83 optimal weight: 0.9990 chunk 75 optimal weight: 2.9990 chunk 41 optimal weight: 0.7980 chunk 25 optimal weight: 0.9980 chunk 50 optimal weight: 0.7980 chunk 40 optimal weight: 0.8980 chunk 77 optimal weight: 3.9990 chunk 30 optimal weight: 0.9990 chunk 47 optimal weight: 0.5980 chunk 57 optimal weight: 3.9990 chunk 89 optimal weight: 3.9990 overall best weight: 0.8180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 68 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 25 ASN ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 84 GLN C 104 GLN E 39 HIS E 108 ASN G 24 GLN G 82 HIS D 47 GLN H 95 GLN Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8214 moved from start: 0.2505 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.042 12861 Z= 0.216 Angle : 0.610 5.503 18629 Z= 0.367 Chirality : 0.037 0.154 2119 Planarity : 0.005 0.034 1340 Dihedral : 31.282 170.201 3803 Min Nonbonded Distance : 2.151 Molprobity Statistics. All-atom Clashscore : 7.52 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.21 % Favored : 98.79 % Rotamer Outliers : 1.91 % 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.26 (0.26), residues: 744 helix: -1.38 (0.19), residues: 532 sheet: None (None), residues: 0 loop : -1.31 (0.37), residues: 212 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 162 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 150 time to evaluate : 0.925 Fit side-chains outliers start: 12 outliers final: 7 residues processed: 159 average time/residue: 0.3434 time to fit residues: 70.7101 Evaluate side-chains 144 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 137 time to evaluate : 0.976 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.1153 time to fit residues: 2.5504 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 49 optimal weight: 2.9990 chunk 27 optimal weight: 0.9990 chunk 74 optimal weight: 3.9990 chunk 61 optimal weight: 3.9990 chunk 24 optimal weight: 0.8980 chunk 90 optimal weight: 1.9990 chunk 97 optimal weight: 0.8980 chunk 80 optimal weight: 4.9990 chunk 89 optimal weight: 1.9990 chunk 30 optimal weight: 0.8980 chunk 72 optimal weight: 4.9990 overall best weight: 1.1384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 104 GLN D 109 HIS H 95 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8238 moved from start: 0.2900 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.048 12861 Z= 0.246 Angle : 0.592 5.702 18629 Z= 0.353 Chirality : 0.037 0.161 2119 Planarity : 0.004 0.036 1340 Dihedral : 31.234 169.205 3803 Min Nonbonded Distance : 2.161 Molprobity Statistics. All-atom Clashscore : 7.20 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.48 % Favored : 98.52 % Rotamer Outliers : 0.79 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.45 (0.29), residues: 744 helix: 0.11 (0.22), residues: 533 sheet: None (None), residues: 0 loop : -1.10 (0.38), residues: 211 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 139 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 134 time to evaluate : 1.018 Fit side-chains outliers start: 5 outliers final: 2 residues processed: 135 average time/residue: 0.3465 time to fit residues: 60.3716 Evaluate side-chains 129 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 127 time to evaluate : 1.031 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.1353 time to fit residues: 1.6851 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 88 optimal weight: 1.9990 chunk 67 optimal weight: 7.9990 chunk 46 optimal weight: 0.8980 chunk 9 optimal weight: 0.6980 chunk 42 optimal weight: 0.7980 chunk 60 optimal weight: 4.9990 chunk 90 optimal weight: 1.9990 chunk 95 optimal weight: 0.9990 chunk 47 optimal weight: 0.5980 chunk 85 optimal weight: 3.9990 chunk 25 optimal weight: 0.9990 overall best weight: 0.7982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 27 GLN ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 104 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8213 moved from start: 0.3165 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.045 12861 Z= 0.200 Angle : 0.557 5.523 18629 Z= 0.336 Chirality : 0.035 0.148 2119 Planarity : 0.004 0.036 1340 Dihedral : 31.025 167.556 3803 Min Nonbonded Distance : 2.157 Molprobity Statistics. All-atom Clashscore : 7.38 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.02 % Favored : 97.98 % Rotamer Outliers : 0.95 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.67 (0.30), residues: 744 helix: 1.01 (0.23), residues: 530 sheet: None (None), residues: 0 loop : -0.86 (0.39), residues: 214 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 135 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 129 time to evaluate : 0.957 Fit side-chains outliers start: 6 outliers final: 3 residues processed: 131 average time/residue: 0.3558 time to fit residues: 60.0565 Evaluate side-chains 129 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 126 time to evaluate : 0.961 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1997 time to fit residues: 1.9992 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 79 optimal weight: 4.9990 chunk 54 optimal weight: 5.9990 chunk 1 optimal weight: 0.9990 chunk 71 optimal weight: 4.9990 chunk 39 optimal weight: 0.6980 chunk 81 optimal weight: 7.9990 chunk 66 optimal weight: 6.9990 chunk 0 optimal weight: 5.9990 chunk 48 optimal weight: 0.9990 chunk 85 optimal weight: 3.9990 chunk 24 optimal weight: 1.9990 overall best weight: 1.7388 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 104 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8270 moved from start: 0.3172 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.006 0.054 12861 Z= 0.316 Angle : 0.627 5.856 18629 Z= 0.369 Chirality : 0.039 0.168 2119 Planarity : 0.005 0.042 1340 Dihedral : 31.479 170.696 3803 Min Nonbonded Distance : 2.149 Molprobity Statistics. All-atom Clashscore : 6.46 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.55 % Favored : 97.45 % Rotamer Outliers : 1.75 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.65 (0.30), residues: 744 helix: 1.01 (0.23), residues: 532 sheet: None (None), residues: 0 loop : -0.95 (0.39), residues: 212 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 143 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 132 time to evaluate : 1.014 Fit side-chains outliers start: 11 outliers final: 8 residues processed: 141 average time/residue: 0.3286 time to fit residues: 60.8243 Evaluate side-chains 134 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 126 time to evaluate : 0.968 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 0 residues processed: 8 average time/residue: 0.1056 time to fit residues: 2.7650 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 32 optimal weight: 0.9990 chunk 86 optimal weight: 1.9990 chunk 18 optimal weight: 0.0670 chunk 56 optimal weight: 3.9990 chunk 23 optimal weight: 0.7980 chunk 95 optimal weight: 2.9990 chunk 79 optimal weight: 4.9990 chunk 44 optimal weight: 0.5980 chunk 7 optimal weight: 0.7980 chunk 31 optimal weight: 1.9990 chunk 50 optimal weight: 1.9990 overall best weight: 0.6520 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 104 GLN H 67 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8244 moved from start: 0.3262 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.043 12861 Z= 0.193 Angle : 0.577 5.580 18629 Z= 0.348 Chirality : 0.036 0.155 2119 Planarity : 0.004 0.038 1340 Dihedral : 31.395 169.777 3803 Min Nonbonded Distance : 2.144 Molprobity Statistics. All-atom Clashscore : 6.83 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.02 % Favored : 97.98 % Rotamer Outliers : 0.79 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.03 (0.30), residues: 744 helix: 1.32 (0.23), residues: 528 sheet: None (None), residues: 0 loop : -0.81 (0.38), residues: 216 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 132 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 127 time to evaluate : 1.021 Fit side-chains outliers start: 5 outliers final: 3 residues processed: 130 average time/residue: 0.3250 time to fit residues: 55.7868 Evaluate side-chains 126 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 123 time to evaluate : 0.964 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1027 time to fit residues: 1.7006 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 92 optimal weight: 1.9990 chunk 10 optimal weight: 0.6980 chunk 54 optimal weight: 5.9990 chunk 69 optimal weight: 5.9990 chunk 80 optimal weight: 4.9990 chunk 53 optimal weight: 5.9990 chunk 95 optimal weight: 1.9990 chunk 59 optimal weight: 4.9990 chunk 58 optimal weight: 3.9990 chunk 44 optimal weight: 2.9990 chunk 38 optimal weight: 0.7980 overall best weight: 1.6986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 104 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8270 moved from start: 0.3271 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.006 0.052 12861 Z= 0.312 Angle : 0.621 5.839 18629 Z= 0.366 Chirality : 0.039 0.162 2119 Planarity : 0.005 0.043 1340 Dihedral : 31.411 170.078 3803 Min Nonbonded Distance : 2.133 Molprobity Statistics. All-atom Clashscore : 6.23 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.28 % Favored : 97.72 % Rotamer Outliers : 0.79 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.88 (0.30), residues: 744 helix: 1.23 (0.23), residues: 530 sheet: None (None), residues: 0 loop : -0.97 (0.38), residues: 214 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 133 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 128 time to evaluate : 1.027 Fit side-chains outliers start: 5 outliers final: 4 residues processed: 129 average time/residue: 0.3463 time to fit residues: 57.9977 Evaluate side-chains 131 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 127 time to evaluate : 1.025 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.1702 time to fit residues: 2.1410 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 56 optimal weight: 3.9990 chunk 28 optimal weight: 0.2980 chunk 18 optimal weight: 0.0000 chunk 60 optimal weight: 4.9990 chunk 65 optimal weight: 4.9990 chunk 47 optimal weight: 0.6980 chunk 8 optimal weight: 0.1980 chunk 75 optimal weight: 3.9990 chunk 86 optimal weight: 1.9990 chunk 91 optimal weight: 3.9990 chunk 83 optimal weight: 0.8980 overall best weight: 0.4184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 104 GLN D 84 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8240 moved from start: 0.3363 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.044 12861 Z= 0.177 Angle : 0.574 5.500 18629 Z= 0.347 Chirality : 0.036 0.226 2119 Planarity : 0.004 0.036 1340 Dihedral : 31.371 169.331 3803 Min Nonbonded Distance : 2.113 Molprobity Statistics. All-atom Clashscore : 6.69 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.02 % Favored : 97.98 % Rotamer Outliers : 0.64 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.28 (0.30), residues: 744 helix: 1.54 (0.23), residues: 530 sheet: None (None), residues: 0 loop : -0.85 (0.38), residues: 214 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 129 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 125 time to evaluate : 0.997 Fit side-chains outliers start: 4 outliers final: 1 residues processed: 128 average time/residue: 0.3519 time to fit residues: 58.4338 Evaluate side-chains 123 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 122 time to evaluate : 1.078 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1033 time to fit residues: 1.3654 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 88 optimal weight: 0.1980 chunk 91 optimal weight: 2.9990 chunk 53 optimal weight: 6.9990 chunk 38 optimal weight: 0.0000 chunk 69 optimal weight: 5.9990 chunk 27 optimal weight: 1.9990 chunk 80 optimal weight: 4.9990 chunk 84 optimal weight: 0.7980 chunk 58 optimal weight: 3.9990 chunk 94 optimal weight: 3.9990 chunk 57 optimal weight: 3.9990 overall best weight: 1.1988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 104 GLN D 84 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8252 moved from start: 0.3403 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.047 12861 Z= 0.248 Angle : 0.583 5.714 18629 Z= 0.348 Chirality : 0.036 0.166 2119 Planarity : 0.004 0.038 1340 Dihedral : 31.255 168.634 3803 Min Nonbonded Distance : 2.061 Molprobity Statistics. All-atom Clashscore : 6.73 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.02 % Favored : 97.98 % Rotamer Outliers : 0.32 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.39 (0.30), residues: 744 helix: 1.61 (0.23), residues: 530 sheet: None (None), residues: 0 loop : -0.77 (0.38), residues: 214 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 130 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 128 time to evaluate : 1.009 Fit side-chains outliers start: 2 outliers final: 2 residues processed: 129 average time/residue: 0.3473 time to fit residues: 58.6644 Evaluate side-chains 128 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 126 time to evaluate : 0.938 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.1028 time to fit residues: 1.4881 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 44 optimal weight: 1.9990 chunk 65 optimal weight: 4.9990 chunk 98 optimal weight: 0.7980 chunk 90 optimal weight: 0.0970 chunk 78 optimal weight: 4.9990 chunk 8 optimal weight: 0.0570 chunk 60 optimal weight: 3.9990 chunk 48 optimal weight: 3.9990 chunk 62 optimal weight: 3.9990 chunk 83 optimal weight: 0.9980 chunk 24 optimal weight: 0.5980 overall best weight: 0.5096 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 104 GLN D 84 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8239 moved from start: 0.3432 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.040 12861 Z= 0.181 Angle : 0.567 5.475 18629 Z= 0.342 Chirality : 0.035 0.175 2119 Planarity : 0.004 0.037 1340 Dihedral : 31.250 168.383 3803 Min Nonbonded Distance : 2.153 Molprobity Statistics. All-atom Clashscore : 6.83 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.15 % Favored : 97.85 % Rotamer Outliers : 0.16 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.54 (0.31), residues: 744 helix: 1.72 (0.23), residues: 530 sheet: None (None), residues: 0 loop : -0.71 (0.38), residues: 214 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 125 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 124 time to evaluate : 1.012 Fit side-chains outliers start: 1 outliers final: 1 residues processed: 125 average time/residue: 0.3442 time to fit residues: 56.4276 Evaluate side-chains 124 residues out of total 632 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 123 time to evaluate : 0.947 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1032 time to fit residues: 1.3029 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 72 optimal weight: 4.9990 chunk 11 optimal weight: 0.7980 chunk 21 optimal weight: 0.6980 chunk 78 optimal weight: 4.9990 chunk 32 optimal weight: 0.9990 chunk 80 optimal weight: 4.9990 chunk 9 optimal weight: 0.3980 chunk 14 optimal weight: 2.9990 chunk 69 optimal weight: 5.9990 chunk 4 optimal weight: 0.9990 chunk 56 optimal weight: 3.9990 overall best weight: 0.7784 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 104 GLN D 84 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3606 r_free = 0.3606 target = 0.140628 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 37)----------------| | r_work = 0.3108 r_free = 0.3108 target = 0.102612 restraints weight = 15870.081| |-----------------------------------------------------------------------------| r_work (start): 0.3153 rms_B_bonded: 0.95 r_work: 0.3124 rms_B_bonded: 0.90 restraints_weight: 0.5000 r_work: 0.3105 rms_B_bonded: 1.10 restraints_weight: 0.2500 r_work: 0.3080 rms_B_bonded: 1.44 restraints_weight: 0.1250 r_work: 0.3046 rms_B_bonded: 1.97 restraints_weight: 0.0625 r_work: 0.2996 rms_B_bonded: 2.79 restraints_weight: 0.0312 r_work: 0.2919 rms_B_bonded: 4.17 restraints_weight: 0.0156 r_work (final): 0.2919 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8424 moved from start: 0.3439 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.043 12861 Z= 0.207 Angle : 0.566 5.587 18629 Z= 0.341 Chirality : 0.035 0.167 2119 Planarity : 0.004 0.038 1340 Dihedral : 31.229 168.425 3803 Min Nonbonded Distance : 2.032 Molprobity Statistics. All-atom Clashscore : 6.87 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.28 % Favored : 97.72 % Rotamer Outliers : 0.16 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.55 (0.31), residues: 744 helix: 1.72 (0.23), residues: 530 sheet: None (None), residues: 0 loop : -0.69 (0.38), residues: 214 =============================================================================== Job complete usr+sys time: 1957.52 seconds wall clock time: 36 minutes 21.95 seconds (2181.95 seconds total)