Starting phenix.real_space_refine (version: 1.21rc1) on Fri Apr 21 14:17:40 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/5xlp_6731/04_2023/5xlp_6731.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/5xlp_6731/04_2023/5xlp_6731.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=4.2 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/5xlp_6731/04_2023/5xlp_6731.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/5xlp_6731/04_2023/5xlp_6731.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/5xlp_6731/04_2023/5xlp_6731.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/5xlp_6731/04_2023/5xlp_6731.pdb" } resolution = 4.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.001 sd= 0.010 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 5 Type Number sf(0) Gaussians P 24 5.49 5 S 11 5.16 5 C 6107 2.51 5 N 1737 2.21 5 O 2073 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "C ARG 161": "NH1" <-> "NH2" Residue "C ARG 218": "NH1" <-> "NH2" Residue "C PHE 336": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D ARG 161": "NH1" <-> "NH2" Residue "D ARG 218": "NH1" <-> "NH2" Residue "D PHE 336": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "E ARG 161": "NH1" <-> "NH2" Residue "E ARG 218": "NH1" <-> "NH2" Residue "E PHE 336": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "F ARG 161": "NH1" <-> "NH2" Residue "F ARG 218": "NH1" <-> "NH2" Residue "F PHE 336": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "M ARG 75": "NH1" <-> "NH2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-4938/modules/chem_data/mon_lib" Total number of atoms: 9952 Number of models: 1 Model: "" Number of chains: 6 Chain: "C" Number of atoms: 2004 Number of conformers: 1 Conformer: "" Number of residues, atoms: 287, 2004 Classifications: {'peptide': 287} Incomplete info: {'truncation_to_alanine': 43} Link IDs: {'PTRANS': 9, 'TRANS': 277} Chain breaks: 1 Unresolved chain link angles: 1 Unresolved non-hydrogen bonds: 245 Unresolved non-hydrogen angles: 307 Unresolved non-hydrogen dihedrals: 228 Unresolved non-hydrogen chiralities: 2 Planarities with less than four sites: {'GLN:plan1': 4, 'HIS:plan': 1, 'TYR:plan': 4, 'TRP:plan': 4, 'ASP:plan': 1, 'PHE:plan': 1, 'GLU:plan': 4, 'ARG:plan': 14} Unresolved non-hydrogen planarities: 180 Chain: "D" Number of atoms: 2273 Number of conformers: 1 Conformer: "" Number of residues, atoms: 326, 2273 Classifications: {'peptide': 326} Incomplete info: {'truncation_to_alanine': 50} Link IDs: {'PTRANS': 12, 'TRANS': 313} Unresolved chain link angles: 1 Unresolved non-hydrogen bonds: 272 Unresolved non-hydrogen angles: 339 Unresolved non-hydrogen dihedrals: 250 Unresolved non-hydrogen chiralities: 3 Planarities with less than four sites: {'GLN:plan1': 5, 'HIS:plan': 1, 'TYR:plan': 4, 'ASN:plan1': 1, 'TRP:plan': 4, 'ASP:plan': 2, 'PHE:plan': 1, 'GLU:plan': 4, 'ARG:plan': 15} Unresolved non-hydrogen planarities: 195 Chain: "E" Number of atoms: 2273 Number of conformers: 1 Conformer: "" Number of residues, atoms: 326, 2273 Classifications: {'peptide': 326} Incomplete info: {'truncation_to_alanine': 50} Link IDs: {'PTRANS': 12, 'TRANS': 313} Unresolved chain link angles: 1 Unresolved non-hydrogen bonds: 272 Unresolved non-hydrogen angles: 339 Unresolved non-hydrogen dihedrals: 250 Unresolved non-hydrogen chiralities: 3 Planarities with less than four sites: {'GLN:plan1': 5, 'HIS:plan': 1, 'TYR:plan': 4, 'ASN:plan1': 1, 'TRP:plan': 4, 'ASP:plan': 2, 'PHE:plan': 1, 'GLU:plan': 4, 'ARG:plan': 15} Unresolved non-hydrogen planarities: 195 Chain: "F" Number of atoms: 2273 Number of conformers: 1 Conformer: "" Number of residues, atoms: 326, 2273 Classifications: {'peptide': 326} Incomplete info: {'truncation_to_alanine': 50} Link IDs: {'PTRANS': 12, 'TRANS': 313} Unresolved chain link angles: 1 Unresolved non-hydrogen bonds: 272 Unresolved non-hydrogen angles: 339 Unresolved non-hydrogen dihedrals: 250 Unresolved non-hydrogen chiralities: 3 Planarities with less than four sites: {'GLN:plan1': 5, 'HIS:plan': 1, 'TYR:plan': 4, 'ASN:plan1': 1, 'TRP:plan': 4, 'ASP:plan': 2, 'PHE:plan': 1, 'GLU:plan': 4, 'ARG:plan': 15} Unresolved non-hydrogen planarities: 195 Chain: "K" Number of atoms: 514 Number of conformers: 1 Conformer: "" Number of residues, atoms: 24, 514 Classifications: {'RNA': 24} Modifications used: {'rna2p_pur': 7, 'rna2p_pyr': 10, 'rna3p_pur': 6, 'rna3p_pyr': 1} Link IDs: {'rna2p': 16, 'rna3p': 7} Chain: "M" Number of atoms: 615 Number of conformers: 1 Conformer: "" Number of residues, atoms: 78, 615 Classifications: {'peptide': 78} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 2, 'TRANS': 75} Unresolved non-hydrogen bonds: 6 Unresolved non-hydrogen angles: 8 Unresolved non-hydrogen dihedrals: 6 Planarities with less than four sites: {'HIS:plan': 1} Unresolved non-hydrogen planarities: 5 Time building chain proxies: 5.88, per 1000 atoms: 0.59 Number of scatterers: 9952 At special positions: 0 Unit cell: (110.04, 111.35, 127.07, 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 11 16.00 P 24 15.00 O 2073 8.00 N 1737 7.00 C 6107 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: 4.14 Conformation dependent library (CDL) restraints added in 1.4 seconds 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2480 Finding SS restraints... Secondary structure from input PDB file: 45 helices and 20 sheets defined 29.3% alpha, 6.0% beta 0 base pairs and 4 stacking pairs defined. Time for finding SS restraints: 2.27 Creating SS restraints... Processing helix chain 'C' and resid 121 through 138 Processing helix chain 'C' and resid 141 through 151 Processing helix chain 'C' and resid 159 through 163 removed outlier: 3.505A pdb=" N ILE C 162 " --> pdb=" O GLU C 159 " (cutoff:3.500A) removed outlier: 3.744A pdb=" N ASN C 163 " --> pdb=" O VAL C 160 " (cutoff:3.500A) No H-bonds generated for 'chain 'C' and resid 159 through 163' Processing helix chain 'C' and resid 187 through 191 Processing helix chain 'C' and resid 199 through 214 removed outlier: 4.197A pdb=" N VAL C 214 " --> pdb=" O LEU C 210 " (cutoff:3.500A) Processing helix chain 'C' and resid 266 through 271 Processing helix chain 'C' and resid 272 through 274 No H-bonds generated for 'chain 'C' and resid 272 through 274' Processing helix chain 'C' and resid 293 through 297 removed outlier: 4.391A pdb=" N GLN C 297 " --> pdb=" O ALA C 294 " (cutoff:3.500A) Processing helix chain 'C' and resid 316 through 325 removed outlier: 4.900A pdb=" N GLN C 323 " --> pdb=" O ALA C 319 " (cutoff:3.500A) Processing helix chain 'C' and resid 331 through 336 Processing helix chain 'D' and resid 74 through 83 removed outlier: 3.933A pdb=" N THR D 78 " --> pdb=" O PRO D 74 " (cutoff:3.500A) removed outlier: 3.547A pdb=" N VAL D 79 " --> pdb=" O ASN D 75 " (cutoff:3.500A) Processing helix chain 'D' and resid 121 through 138 Processing helix chain 'D' and resid 141 through 151 Processing helix chain 'D' and resid 159 through 163 removed outlier: 3.505A pdb=" N ILE D 162 " --> pdb=" O GLU D 159 " (cutoff:3.500A) removed outlier: 3.743A pdb=" N ASN D 163 " --> pdb=" O VAL D 160 " (cutoff:3.500A) No H-bonds generated for 'chain 'D' and resid 159 through 163' Processing helix chain 'D' and resid 187 through 191 Processing helix chain 'D' and resid 199 through 214 removed outlier: 4.198A pdb=" N VAL D 214 " --> pdb=" O LEU D 210 " (cutoff:3.500A) Processing helix chain 'D' and resid 266 through 271 Processing helix chain 'D' and resid 272 through 274 No H-bonds generated for 'chain 'D' and resid 272 through 274' Processing helix chain 'D' and resid 293 through 297 removed outlier: 4.390A pdb=" N GLN D 297 " --> pdb=" O ALA D 294 " (cutoff:3.500A) Processing helix chain 'D' and resid 316 through 325 removed outlier: 4.901A pdb=" N GLN D 323 " --> pdb=" O ALA D 319 " (cutoff:3.500A) Processing helix chain 'D' and resid 331 through 336 Processing helix chain 'E' and resid 75 through 83 removed outlier: 3.520A pdb=" N VAL E 79 " --> pdb=" O ASN E 75 " (cutoff:3.500A) Processing helix chain 'E' and resid 121 through 138 Processing helix chain 'E' and resid 141 through 151 Processing helix chain 'E' and resid 159 through 163 removed outlier: 3.504A pdb=" N ILE E 162 " --> pdb=" O GLU E 159 " (cutoff:3.500A) removed outlier: 3.743A pdb=" N ASN E 163 " --> pdb=" O VAL E 160 " (cutoff:3.500A) No H-bonds generated for 'chain 'E' and resid 159 through 163' Processing helix chain 'E' and resid 187 through 191 Processing helix chain 'E' and resid 199 through 214 removed outlier: 4.198A pdb=" N VAL E 214 " --> pdb=" O LEU E 210 " (cutoff:3.500A) Processing helix chain 'E' and resid 266 through 271 Processing helix chain 'E' and resid 272 through 274 No H-bonds generated for 'chain 'E' and resid 272 through 274' Processing helix chain 'E' and resid 293 through 297 removed outlier: 4.391A pdb=" N GLN E 297 " --> pdb=" O ALA E 294 " (cutoff:3.500A) Processing helix chain 'E' and resid 316 through 325 removed outlier: 4.901A pdb=" N GLN E 323 " --> pdb=" O ALA E 319 " (cutoff:3.500A) Processing helix chain 'E' and resid 331 through 336 Processing helix chain 'F' and resid 73 through 81 removed outlier: 4.084A pdb=" N GLN F 77 " --> pdb=" O SER F 73 " (cutoff:3.500A) removed outlier: 3.697A pdb=" N VAL F 79 " --> pdb=" O ASN F 75 " (cutoff:3.500A) Processing helix chain 'F' and resid 121 through 138 Processing helix chain 'F' and resid 141 through 151 Processing helix chain 'F' and resid 159 through 163 removed outlier: 3.505A pdb=" N ILE F 162 " --> pdb=" O GLU F 159 " (cutoff:3.500A) removed outlier: 3.743A pdb=" N ASN F 163 " --> pdb=" O VAL F 160 " (cutoff:3.500A) No H-bonds generated for 'chain 'F' and resid 159 through 163' Processing helix chain 'F' and resid 187 through 191 Processing helix chain 'F' and resid 199 through 214 removed outlier: 4.198A pdb=" N VAL F 214 " --> pdb=" O LEU F 210 " (cutoff:3.500A) Processing helix chain 'F' and resid 266 through 271 Processing helix chain 'F' and resid 272 through 274 No H-bonds generated for 'chain 'F' and resid 272 through 274' Processing helix chain 'F' and resid 293 through 297 removed outlier: 4.391A pdb=" N GLN F 297 " --> pdb=" O ALA F 294 " (cutoff:3.500A) Processing helix chain 'F' and resid 316 through 325 removed outlier: 4.902A pdb=" N GLN F 323 " --> pdb=" O ALA F 319 " (cutoff:3.500A) Processing helix chain 'F' and resid 331 through 336 Processing helix chain 'M' and resid 46 through 56 removed outlier: 4.256A pdb=" N SER M 52 " --> pdb=" O GLU M 48 " (cutoff:3.500A) Processing helix chain 'M' and resid 61 through 74 removed outlier: 3.596A pdb=" N ILE M 65 " --> pdb=" O GLY M 61 " (cutoff:3.500A) removed outlier: 3.665A pdb=" N ARG M 74 " --> pdb=" O ASN M 70 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'C' and resid 36 through 37 removed outlier: 3.745A pdb=" N ARG C 50 " --> pdb=" O GLN C 37 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'C' and resid 107 through 108 removed outlier: 3.659A pdb=" N SER C 107 " --> pdb=" O ASP C 221 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'C' and resid 168 through 172 Processing sheet with id=AA4, first strand: chain 'C' and resid 174 through 175 Processing sheet with id=AA5, first strand: chain 'D' and resid 36 through 37 removed outlier: 3.745A pdb=" N ARG D 50 " --> pdb=" O GLN D 37 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'D' and resid 56 through 59 Processing sheet with id=AA7, first strand: chain 'D' and resid 107 through 108 removed outlier: 3.658A pdb=" N SER D 107 " --> pdb=" O ASP D 221 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'D' and resid 168 through 172 Processing sheet with id=AA9, first strand: chain 'D' and resid 174 through 175 Processing sheet with id=AB1, first strand: chain 'E' and resid 36 through 37 removed outlier: 3.746A pdb=" N ARG E 50 " --> pdb=" O GLN E 37 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'E' and resid 57 through 59 Processing sheet with id=AB3, first strand: chain 'E' and resid 107 through 108 removed outlier: 3.658A pdb=" N SER E 107 " --> pdb=" O ASP E 221 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'E' and resid 168 through 172 Processing sheet with id=AB5, first strand: chain 'E' and resid 174 through 175 Processing sheet with id=AB6, first strand: chain 'F' and resid 36 through 37 removed outlier: 3.744A pdb=" N ARG F 50 " --> pdb=" O GLN F 37 " (cutoff:3.500A) Processing sheet with id=AB7, first strand: chain 'F' and resid 55 through 58 removed outlier: 3.626A pdb=" N ASN F 55 " --> pdb=" O LEU F 91 " (cutoff:3.500A) Processing sheet with id=AB8, first strand: chain 'F' and resid 107 through 108 removed outlier: 3.658A pdb=" N SER F 107 " --> pdb=" O ASP F 221 " (cutoff:3.500A) Processing sheet with id=AB9, first strand: chain 'F' and resid 168 through 172 Processing sheet with id=AC1, first strand: chain 'F' and resid 174 through 175 Processing sheet with id=AC2, first strand: chain 'M' and resid 19 through 21 248 hydrogen bonds defined for protein. 684 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 4 stacking parallelities Total time for adding SS restraints: 2.31 Time building geometry restraints manager: 4.44 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: 2658 1.33 - 1.45: 1827 1.45 - 1.57: 5592 1.57 - 1.69: 47 1.69 - 1.80: 18 Bond restraints: 10142 Sorted by residual: bond pdb=" C HIS E 164 " pdb=" N ILE E 165 " ideal model delta sigma weight residual 1.334 1.306 0.028 1.24e-02 6.50e+03 5.00e+00 bond pdb=" C HIS C 164 " pdb=" N ILE C 165 " ideal model delta sigma weight residual 1.334 1.306 0.028 1.24e-02 6.50e+03 4.94e+00 bond pdb=" C LEU M 12 " pdb=" N ASN M 13 " ideal model delta sigma weight residual 1.331 1.282 0.049 2.22e-02 2.03e+03 4.85e+00 bond pdb=" C HIS F 164 " pdb=" N ILE F 165 " ideal model delta sigma weight residual 1.334 1.307 0.027 1.24e-02 6.50e+03 4.77e+00 bond pdb=" C HIS D 164 " pdb=" N ILE D 165 " ideal model delta sigma weight residual 1.334 1.307 0.027 1.24e-02 6.50e+03 4.72e+00 ... (remaining 10137 not shown) Histogram of bond angle deviations from ideal: 98.72 - 106.15: 372 106.15 - 113.58: 5705 113.58 - 121.02: 4542 121.02 - 128.45: 3246 128.45 - 135.88: 98 Bond angle restraints: 13963 Sorted by residual: angle pdb=" N VAL E 42 " pdb=" CA VAL E 42 " pdb=" C VAL E 42 " ideal model delta sigma weight residual 111.88 101.65 10.23 1.06e+00 8.90e-01 9.32e+01 angle pdb=" N VAL D 42 " pdb=" CA VAL D 42 " pdb=" C VAL D 42 " ideal model delta sigma weight residual 111.88 101.66 10.22 1.06e+00 8.90e-01 9.29e+01 angle pdb=" N VAL C 42 " pdb=" CA VAL C 42 " pdb=" C VAL C 42 " ideal model delta sigma weight residual 111.88 101.67 10.21 1.06e+00 8.90e-01 9.28e+01 angle pdb=" N VAL F 42 " pdb=" CA VAL F 42 " pdb=" C VAL F 42 " ideal model delta sigma weight residual 111.88 101.67 10.21 1.06e+00 8.90e-01 9.27e+01 angle pdb=" N ILE D 280 " pdb=" CA ILE D 280 " pdb=" C ILE D 280 " ideal model delta sigma weight residual 111.88 102.37 9.51 1.06e+00 8.90e-01 8.06e+01 ... (remaining 13958 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 35.85: 5914 35.85 - 71.71: 65 71.71 - 107.56: 13 107.56 - 143.41: 1 143.41 - 179.27: 6 Dihedral angle restraints: 5999 sinusoidal: 2158 harmonic: 3841 Sorted by residual: dihedral pdb=" CA HIS M 11 " pdb=" C HIS M 11 " pdb=" N LEU M 12 " pdb=" CA LEU M 12 " ideal model delta harmonic sigma weight residual -180.00 -118.49 -61.51 0 5.00e+00 4.00e-02 1.51e+02 dihedral pdb=" C5' C K 12 " pdb=" C4' C K 12 " pdb=" C3' C K 12 " pdb=" O3' C K 12 " ideal model delta sinusoidal sigma weight residual 147.00 84.74 62.26 1 8.00e+00 1.56e-02 8.02e+01 dihedral pdb=" O4' C K 5 " pdb=" C1' C K 5 " pdb=" N1 C K 5 " pdb=" C2 C K 5 " ideal model delta sinusoidal sigma weight residual 200.00 48.73 151.27 1 1.50e+01 4.44e-03 8.01e+01 ... (remaining 5996 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.083: 1424 0.083 - 0.165: 241 0.165 - 0.248: 21 0.248 - 0.330: 2 0.330 - 0.413: 3 Chirality restraints: 1691 Sorted by residual: chirality pdb=" C3' G K 18 " pdb=" C4' G K 18 " pdb=" O3' G K 18 " pdb=" C2' G K 18 " both_signs ideal model delta sigma weight residual False -2.74 -2.33 -0.41 2.00e-01 2.50e+01 4.26e+00 chirality pdb=" C3' A K 16 " pdb=" C4' A K 16 " pdb=" O3' A K 16 " pdb=" C2' A K 16 " both_signs ideal model delta sigma weight residual False -2.74 -2.36 -0.38 2.00e-01 2.50e+01 3.66e+00 chirality pdb=" C3' G K 9 " pdb=" C4' G K 9 " pdb=" O3' G K 9 " pdb=" C2' G K 9 " both_signs ideal model delta sigma weight residual False -2.48 -2.13 -0.35 2.00e-01 2.50e+01 3.09e+00 ... (remaining 1688 not shown) Planarity restraints: 1783 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA HIS M 11 " 0.023 2.00e-02 2.50e+03 4.63e-02 2.14e+01 pdb=" C HIS M 11 " -0.080 2.00e-02 2.50e+03 pdb=" O HIS M 11 " 0.029 2.00e-02 2.50e+03 pdb=" N LEU M 12 " 0.028 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C TYR F 271 " -0.051 5.00e-02 4.00e+02 7.60e-02 9.24e+00 pdb=" N PRO F 272 " 0.131 5.00e-02 4.00e+02 pdb=" CA PRO F 272 " -0.038 5.00e-02 4.00e+02 pdb=" CD PRO F 272 " -0.042 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C TYR E 271 " -0.051 5.00e-02 4.00e+02 7.60e-02 9.24e+00 pdb=" N PRO E 272 " 0.131 5.00e-02 4.00e+02 pdb=" CA PRO E 272 " -0.038 5.00e-02 4.00e+02 pdb=" CD PRO E 272 " -0.042 5.00e-02 4.00e+02 ... (remaining 1780 not shown) Histogram of nonbonded interaction distances: 2.19 - 2.74: 1094 2.74 - 3.28: 9866 3.28 - 3.82: 15966 3.82 - 4.36: 18829 4.36 - 4.90: 30316 Nonbonded interactions: 76071 Sorted by model distance: nonbonded pdb=" OD1 ASP F 80 " pdb=" OH TYR M 6 " model vdw 2.194 2.440 nonbonded pdb=" O ASN E 75 " pdb=" OG1 THR E 78 " model vdw 2.232 2.440 nonbonded pdb=" OG1 THR D 90 " pdb=" O ALA D 254 " model vdw 2.258 2.440 nonbonded pdb=" OG1 THR C 105 " pdb=" OE1 GLN C 223 " model vdw 2.278 2.440 nonbonded pdb=" OG1 THR D 105 " pdb=" OE1 GLN D 223 " model vdw 2.279 2.440 ... (remaining 76066 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 'C' selection = (chain 'D' and (resid 15 through 53 or resid 93 through 340)) selection = (chain 'E' and (resid 15 through 53 or resid 93 through 340)) selection = (chain 'F' and (resid 15 through 53 or resid 93 through 340)) } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.600 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: 2.470 Check model and map are aligned: 0.150 Set scattering table: 0.090 Process input model: 28.720 Find NCS groups from input model: 0.640 Set up NCS constraints: 0.040 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:2.990 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 35.710 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5858 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.095 10142 Z= 0.586 Angle : 1.358 15.262 13963 Z= 0.748 Chirality : 0.064 0.413 1691 Planarity : 0.009 0.076 1783 Dihedral : 14.752 179.265 3519 Min Nonbonded Distance : 2.194 Molprobity Statistics. All-atom Clashscore : 12.48 Ramachandran Plot: Outliers : 0.45 % Allowed : 18.18 % Favored : 81.37 % Rotamer Outliers : 0.57 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.85 % 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.53 (0.17), residues: 1331 helix: -3.73 (0.19), residues: 253 sheet: -3.98 (0.28), residues: 206 loop : -3.80 (0.16), residues: 872 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 319 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 314 time to evaluate : 1.183 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 5 outliers final: 1 residues processed: 318 average time/residue: 0.2731 time to fit residues: 113.4793 Evaluate side-chains 181 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 180 time to evaluate : 1.121 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0851 time to fit residues: 1.6326 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 112 optimal weight: 3.9990 chunk 100 optimal weight: 6.9990 chunk 55 optimal weight: 1.9990 chunk 34 optimal weight: 2.9990 chunk 67 optimal weight: 7.9990 chunk 53 optimal weight: 0.7980 chunk 104 optimal weight: 9.9990 chunk 40 optimal weight: 9.9990 chunk 63 optimal weight: 5.9990 chunk 77 optimal weight: 5.9990 chunk 120 optimal weight: 0.5980 overall best weight: 2.0786 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 140 HIS C 141 ASN C 144 ASN C 151 ASN ** D 83 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 141 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 144 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** E 151 ASN ** E 229 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** E 300 GLN F 140 HIS F 151 ASN F 258 GLN F 300 GLN ** M 32 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 10 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6200 moved from start: 0.4208 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.056 10142 Z= 0.296 Angle : 0.905 12.018 13963 Z= 0.466 Chirality : 0.046 0.280 1691 Planarity : 0.007 0.077 1783 Dihedral : 16.682 177.797 1691 Min Nonbonded Distance : 2.047 Molprobity Statistics. All-atom Clashscore : 20.62 Ramachandran Plot: Outliers : 0.38 % Allowed : 13.67 % Favored : 85.95 % Rotamer Outliers : 0.23 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.16 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.81 (0.19), residues: 1331 helix: -2.34 (0.26), residues: 293 sheet: -3.83 (0.29), residues: 180 loop : -3.58 (0.17), residues: 858 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 237 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 235 time to evaluate : 1.072 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 2 outliers final: 0 residues processed: 236 average time/residue: 0.2810 time to fit residues: 86.8724 Evaluate side-chains 168 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 168 time to evaluate : 1.094 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.4581 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 67 optimal weight: 3.9990 chunk 37 optimal weight: 8.9990 chunk 100 optimal weight: 5.9990 chunk 82 optimal weight: 0.9980 chunk 33 optimal weight: 7.9990 chunk 120 optimal weight: 8.9990 chunk 130 optimal weight: 0.9990 chunk 107 optimal weight: 6.9990 chunk 119 optimal weight: 10.0000 chunk 41 optimal weight: 4.9990 chunk 97 optimal weight: 7.9990 overall best weight: 3.3988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 144 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 241 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 208 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 141 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 144 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** E 223 GLN ** E 229 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** E 241 GLN E 256 HIS E 300 GLN F 208 HIS ** M 32 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6489 moved from start: 0.6288 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.006 0.073 10142 Z= 0.409 Angle : 1.001 10.437 13963 Z= 0.523 Chirality : 0.050 0.302 1691 Planarity : 0.008 0.078 1783 Dihedral : 16.871 178.640 1691 Min Nonbonded Distance : 1.869 Molprobity Statistics. All-atom Clashscore : 26.81 Ramachandran Plot: Outliers : 0.38 % Allowed : 17.58 % Favored : 82.04 % Rotamer Outliers : 0.80 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.47 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.76 (0.19), residues: 1331 helix: -2.16 (0.27), residues: 309 sheet: -4.01 (0.32), residues: 147 loop : -3.62 (0.18), residues: 875 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 213 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 206 time to evaluate : 1.187 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 7 outliers final: 1 residues processed: 212 average time/residue: 0.2933 time to fit residues: 81.9697 Evaluate side-chains 152 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 151 time to evaluate : 1.096 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1129 time to fit residues: 1.6620 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 119 optimal weight: 9.9990 chunk 90 optimal weight: 0.4980 chunk 62 optimal weight: 5.9990 chunk 13 optimal weight: 0.8980 chunk 57 optimal weight: 0.8980 chunk 81 optimal weight: 4.9990 chunk 121 optimal weight: 10.0000 chunk 128 optimal weight: 5.9990 chunk 63 optimal weight: 0.0470 chunk 115 optimal weight: 6.9990 chunk 34 optimal weight: 2.9990 overall best weight: 1.0680 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 32 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 208 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** D 262 ASN ** E 140 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 141 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 144 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 229 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** F 140 HIS ** F 229 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** F 329 ASN ** M 32 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6318 moved from start: 0.6740 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.036 10142 Z= 0.210 Angle : 0.810 9.832 13963 Z= 0.419 Chirality : 0.044 0.248 1691 Planarity : 0.006 0.069 1783 Dihedral : 16.227 177.189 1691 Min Nonbonded Distance : 2.021 Molprobity Statistics. All-atom Clashscore : 19.73 Ramachandran Plot: Outliers : 0.15 % Allowed : 13.90 % Favored : 85.95 % Rotamer Outliers : 0.11 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.31 (0.20), residues: 1331 helix: -1.40 (0.30), residues: 309 sheet: -3.81 (0.33), residues: 164 loop : -3.44 (0.18), residues: 858 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 217 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 216 time to evaluate : 1.157 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 1 outliers final: 0 residues processed: 217 average time/residue: 0.2851 time to fit residues: 81.0288 Evaluate side-chains 158 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 158 time to evaluate : 1.092 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.4577 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 107 optimal weight: 7.9990 chunk 72 optimal weight: 2.9990 chunk 1 optimal weight: 7.9990 chunk 95 optimal weight: 9.9990 chunk 53 optimal weight: 2.9990 chunk 109 optimal weight: 0.8980 chunk 88 optimal weight: 4.9990 chunk 0 optimal weight: 5.9990 chunk 65 optimal weight: 3.9990 chunk 115 optimal weight: 7.9990 chunk 32 optimal weight: 6.9990 overall best weight: 3.1788 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** D 127 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 140 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** F 329 ASN ** M 32 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6576 moved from start: 0.7676 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.074 10142 Z= 0.370 Angle : 0.962 11.920 13963 Z= 0.499 Chirality : 0.049 0.271 1691 Planarity : 0.007 0.088 1783 Dihedral : 16.518 179.217 1691 Min Nonbonded Distance : 2.006 Molprobity Statistics. All-atom Clashscore : 28.45 Ramachandran Plot: Outliers : 0.38 % Allowed : 19.46 % Favored : 80.17 % Rotamer Outliers : 0.11 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.16 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.53 (0.20), residues: 1331 helix: -1.71 (0.28), residues: 325 sheet: -3.55 (0.33), residues: 170 loop : -3.64 (0.18), residues: 836 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 190 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 189 time to evaluate : 1.105 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 1 outliers final: 1 residues processed: 190 average time/residue: 0.2947 time to fit residues: 72.5329 Evaluate side-chains 152 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 151 time to evaluate : 1.114 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0974 time to fit residues: 1.7730 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 43 optimal weight: 0.9990 chunk 115 optimal weight: 9.9990 chunk 25 optimal weight: 9.9990 chunk 75 optimal weight: 0.9990 chunk 31 optimal weight: 2.9990 chunk 128 optimal weight: 8.9990 chunk 106 optimal weight: 0.8980 chunk 59 optimal weight: 5.9990 chunk 10 optimal weight: 0.9990 chunk 42 optimal weight: 0.1980 chunk 67 optimal weight: 0.2980 overall best weight: 0.6784 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 127 ASN C 324 HIS D 151 ASN ** D 229 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** E 110 ASN ** E 141 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 144 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** F 329 ASN M 32 ASN Total number of N/Q/H flips: 6 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6327 moved from start: 0.7943 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.041 10142 Z= 0.190 Angle : 0.790 12.734 13963 Z= 0.404 Chirality : 0.044 0.233 1691 Planarity : 0.005 0.075 1783 Dihedral : 15.897 176.331 1691 Min Nonbonded Distance : 2.079 Molprobity Statistics. All-atom Clashscore : 18.72 Ramachandran Plot: Outliers : 0.23 % Allowed : 13.37 % Favored : 86.40 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.08 (0.21), residues: 1331 helix: -1.20 (0.29), residues: 336 sheet: -3.40 (0.34), residues: 170 loop : -3.41 (0.18), residues: 825 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 213 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 213 time to evaluate : 1.149 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 213 average time/residue: 0.2959 time to fit residues: 82.3246 Evaluate side-chains 162 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 162 time to evaluate : 1.113 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.4530 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 124 optimal weight: 9.9990 chunk 14 optimal weight: 9.9990 chunk 73 optimal weight: 2.9990 chunk 93 optimal weight: 10.0000 chunk 72 optimal weight: 9.9990 chunk 108 optimal weight: 9.9990 chunk 71 optimal weight: 2.9990 chunk 128 optimal weight: 7.9990 chunk 80 optimal weight: 7.9990 chunk 78 optimal weight: 3.9990 chunk 59 optimal weight: 5.9990 overall best weight: 4.7990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 241 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 324 HIS D 140 HIS ** D 151 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** E 322 GLN F 324 HIS Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6700 moved from start: 0.8818 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.008 0.097 10142 Z= 0.510 Angle : 1.110 13.375 13963 Z= 0.573 Chirality : 0.053 0.290 1691 Planarity : 0.008 0.087 1783 Dihedral : 16.870 179.420 1691 Min Nonbonded Distance : 1.987 Molprobity Statistics. All-atom Clashscore : 32.63 Ramachandran Plot: Outliers : 0.53 % Allowed : 21.19 % Favored : 78.29 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.39 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.77 (0.19), residues: 1331 helix: -1.69 (0.27), residues: 325 sheet: -3.63 (0.32), residues: 176 loop : -3.92 (0.17), residues: 830 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 197 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 197 time to evaluate : 1.069 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 197 average time/residue: 0.3067 time to fit residues: 78.3314 Evaluate side-chains 152 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 152 time to evaluate : 1.093 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.4811 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 79 optimal weight: 2.9990 chunk 51 optimal weight: 2.9990 chunk 76 optimal weight: 6.9990 chunk 38 optimal weight: 2.9990 chunk 25 optimal weight: 5.9990 chunk 24 optimal weight: 1.9990 chunk 81 optimal weight: 0.9990 chunk 87 optimal weight: 0.9980 chunk 63 optimal weight: 0.9990 chunk 11 optimal weight: 6.9990 chunk 100 optimal weight: 3.9990 overall best weight: 1.5988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 127 ASN ** C 241 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 127 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 151 ASN D 324 HIS E 258 GLN E 310 ASN E 322 GLN F 329 ASN Total number of N/Q/H flips: 7 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6479 moved from start: 0.8966 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.049 10142 Z= 0.242 Angle : 0.865 11.367 13963 Z= 0.443 Chirality : 0.046 0.248 1691 Planarity : 0.006 0.085 1783 Dihedral : 16.260 178.529 1691 Min Nonbonded Distance : 2.053 Molprobity Statistics. All-atom Clashscore : 21.31 Ramachandran Plot: Outliers : 0.23 % Allowed : 14.95 % Favored : 84.82 % Rotamer Outliers : 0.11 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.23 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.38 (0.20), residues: 1331 helix: -1.52 (0.28), residues: 348 sheet: -3.51 (0.32), residues: 186 loop : -3.59 (0.19), residues: 797 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 216 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 215 time to evaluate : 1.099 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 1 outliers final: 1 residues processed: 216 average time/residue: 0.2694 time to fit residues: 77.0748 Evaluate side-chains 151 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 150 time to evaluate : 1.066 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0863 time to fit residues: 1.5809 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 116 optimal weight: 20.0000 chunk 122 optimal weight: 9.9990 chunk 112 optimal weight: 0.8980 chunk 119 optimal weight: 8.9990 chunk 71 optimal weight: 0.8980 chunk 52 optimal weight: 4.9990 chunk 93 optimal weight: 9.9990 chunk 36 optimal weight: 7.9990 chunk 107 optimal weight: 0.9990 chunk 78 optimal weight: 10.0000 chunk 126 optimal weight: 3.9990 overall best weight: 2.3586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 241 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 262 ASN C 324 HIS D 140 HIS ** D 324 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** E 127 ASN E 258 GLN E 322 GLN Total number of N/Q/H flips: 6 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6566 moved from start: 0.9245 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.062 10142 Z= 0.304 Angle : 0.910 11.205 13963 Z= 0.466 Chirality : 0.047 0.254 1691 Planarity : 0.006 0.082 1783 Dihedral : 16.339 176.807 1691 Min Nonbonded Distance : 2.054 Molprobity Statistics. All-atom Clashscore : 25.22 Ramachandran Plot: Outliers : 0.38 % Allowed : 19.53 % Favored : 80.09 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.23 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.42 (0.20), residues: 1331 helix: -1.57 (0.28), residues: 349 sheet: -3.44 (0.33), residues: 176 loop : -3.65 (0.18), residues: 806 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 196 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 196 time to evaluate : 1.149 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 196 average time/residue: 0.2923 time to fit residues: 75.3242 Evaluate side-chains 156 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 156 time to evaluate : 1.126 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.5379 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 77 optimal weight: 10.0000 chunk 59 optimal weight: 0.1980 chunk 87 optimal weight: 0.0770 chunk 132 optimal weight: 3.9990 chunk 121 optimal weight: 8.9990 chunk 105 optimal weight: 3.9990 chunk 10 optimal weight: 0.9980 chunk 81 optimal weight: 6.9990 chunk 64 optimal weight: 7.9990 chunk 83 optimal weight: 9.9990 chunk 112 optimal weight: 0.9980 overall best weight: 1.2540 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 229 GLN ** C 241 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** D 324 HIS E 127 ASN ** E 229 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** E 322 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6455 moved from start: 0.9456 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.052 10142 Z= 0.230 Angle : 0.854 11.165 13963 Z= 0.435 Chirality : 0.045 0.237 1691 Planarity : 0.006 0.108 1783 Dihedral : 16.118 178.215 1691 Min Nonbonded Distance : 1.932 Molprobity Statistics. All-atom Clashscore : 21.68 Ramachandran Plot: Outliers : 0.23 % Allowed : 16.15 % Favored : 83.62 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.16 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.14 (0.21), residues: 1331 helix: -1.29 (0.29), residues: 340 sheet: -3.47 (0.32), residues: 186 loop : -3.41 (0.19), residues: 805 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2662 Ramachandran restraints generated. 1331 Oldfield, 0 Emsley, 1331 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 202 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 202 time to evaluate : 0.993 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 202 average time/residue: 0.2915 time to fit residues: 77.1547 Evaluate side-chains 156 residues out of total 1067 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 156 time to evaluate : 1.096 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.5881 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 133 random chunks: chunk 32 optimal weight: 0.0970 chunk 97 optimal weight: 0.0870 chunk 15 optimal weight: 1.9990 chunk 29 optimal weight: 0.0770 chunk 105 optimal weight: 9.9990 chunk 44 optimal weight: 6.9990 chunk 108 optimal weight: 10.0000 chunk 13 optimal weight: 7.9990 chunk 19 optimal weight: 9.9990 chunk 92 optimal weight: 20.0000 chunk 6 optimal weight: 0.3980 overall best weight: 0.5316 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 241 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** D 324 HIS E 127 ASN E 300 GLN E 322 GLN F 329 ASN M 32 ASN ** M 49 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** M 70 ASN Total number of N/Q/H flips: 7 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4446 r_free = 0.4446 target = 0.171738 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 47)----------------| | r_work = 0.4277 r_free = 0.4277 target = 0.158456 restraints weight = 25844.554| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 29)----------------| | r_work = 0.4275 r_free = 0.4275 target = 0.158692 restraints weight = 30043.142| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 28)----------------| | r_work = 0.4292 r_free = 0.4292 target = 0.159966 restraints weight = 24485.356| |-----------------------------------------------------------------------------| r_work (final): 0.4294 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6317 moved from start: 0.9628 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.053 10142 Z= 0.187 Angle : 0.803 9.728 13963 Z= 0.409 Chirality : 0.044 0.221 1691 Planarity : 0.006 0.101 1783 Dihedral : 15.766 178.221 1691 Min Nonbonded Distance : 1.956 Molprobity Statistics. All-atom Clashscore : 18.56 Ramachandran Plot: Outliers : 0.15 % Allowed : 15.33 % Favored : 84.52 % Rotamer Outliers : 0.11 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.16 % 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.85 (0.21), residues: 1331 helix: -1.00 (0.30), residues: 346 sheet: -3.38 (0.32), residues: 186 loop : -3.23 (0.19), residues: 799 =============================================================================== Job complete usr+sys time: 2153.95 seconds wall clock time: 39 minutes 47.40 seconds (2387.40 seconds total)