Starting phenix.real_space_refine on Sat Mar 16 13:01:01 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/03_2024/7m6q_23704_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/03_2024/7m6q_23704.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.91 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/03_2024/7m6q_23704.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/03_2024/7m6q_23704.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/03_2024/7m6q_23704_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/03_2024/7m6q_23704_updated.pdb" } resolution = 2.91 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= 0.000 sd= 0.001 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 4 Type Number sf(0) Gaussians S 100 5.16 5 C 9590 2.51 5 N 2390 2.21 5 O 2645 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped None Time to flip residues: 0.04s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/chem_data/mon_lib" Total number of atoms: 14725 Number of models: 1 Model: "" Number of chains: 10 Chain: "A" Number of atoms: 2903 Number of conformers: 1 Conformer: "" Number of residues, atoms: 359, 2903 Classifications: {'peptide': 359} Incomplete info: {'truncation_to_alanine': 5} Link IDs: {'PTRANS': 16, 'TRANS': 342} Chain breaks: 1 Unresolved non-hydrogen bonds: 22 Unresolved non-hydrogen angles: 27 Unresolved non-hydrogen dihedrals: 19 Unresolved non-hydrogen chiralities: 2 Planarities with less than four sites: {'PHE:plan': 1, 'ARG:plan': 1} Unresolved non-hydrogen planarities: 11 Chain: "B" Number of atoms: 2903 Number of conformers: 1 Conformer: "" Number of residues, atoms: 359, 2903 Classifications: {'peptide': 359} Incomplete info: {'truncation_to_alanine': 5} Link IDs: {'PTRANS': 16, 'TRANS': 342} Chain breaks: 1 Unresolved non-hydrogen bonds: 22 Unresolved non-hydrogen angles: 27 Unresolved non-hydrogen dihedrals: 19 Unresolved non-hydrogen chiralities: 2 Planarities with less than four sites: {'PHE:plan': 1, 'ARG:plan': 1} Unresolved non-hydrogen planarities: 11 Chain: "C" Number of atoms: 2903 Number of conformers: 1 Conformer: "" Number of residues, atoms: 359, 2903 Classifications: {'peptide': 359} Incomplete info: {'truncation_to_alanine': 5} Link IDs: {'PTRANS': 16, 'TRANS': 342} Chain breaks: 1 Unresolved non-hydrogen bonds: 22 Unresolved non-hydrogen angles: 27 Unresolved non-hydrogen dihedrals: 19 Unresolved non-hydrogen chiralities: 2 Planarities with less than four sites: {'PHE:plan': 1, 'ARG:plan': 1} Unresolved non-hydrogen planarities: 11 Chain: "D" Number of atoms: 2903 Number of conformers: 1 Conformer: "" Number of residues, atoms: 359, 2903 Classifications: {'peptide': 359} Incomplete info: {'truncation_to_alanine': 5} Link IDs: {'PTRANS': 16, 'TRANS': 342} Chain breaks: 1 Unresolved non-hydrogen bonds: 22 Unresolved non-hydrogen angles: 27 Unresolved non-hydrogen dihedrals: 19 Unresolved non-hydrogen chiralities: 2 Planarities with less than four sites: {'PHE:plan': 1, 'ARG:plan': 1} Unresolved non-hydrogen planarities: 11 Chain: "E" Number of atoms: 2903 Number of conformers: 1 Conformer: "" Number of residues, atoms: 359, 2903 Classifications: {'peptide': 359} Incomplete info: {'truncation_to_alanine': 5} Link IDs: {'PTRANS': 16, 'TRANS': 342} Chain breaks: 1 Unresolved non-hydrogen bonds: 22 Unresolved non-hydrogen angles: 27 Unresolved non-hydrogen dihedrals: 19 Unresolved non-hydrogen chiralities: 2 Planarities with less than four sites: {'PHE:plan': 1, 'ARG:plan': 1} Unresolved non-hydrogen planarities: 11 Chain: "A" Number of atoms: 47 Number of conformers: 1 Conformer: "" Number of residues, atoms: 4, 47 Unusual residues: {'NAG': 1, 'TCI': 1} Classifications: {'peptide': 2, 'undetermined': 2} Modifications used: {'COO': 2} Link IDs: {None: 3} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen dihedrals: 3 Unresolved non-hydrogen chiralities: 1 Chain: "B" Number of atoms: 42 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 42 Unusual residues: {'NAG': 1, 'TCI': 1} Classifications: {'peptide': 1, 'undetermined': 2} Modifications used: {'COO': 1} Link IDs: {None: 2} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen dihedrals: 3 Unresolved non-hydrogen chiralities: 1 Chain: "C" Number of atoms: 42 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 42 Unusual residues: {'NAG': 1, 'TCI': 1} Classifications: {'peptide': 1, 'undetermined': 2} Modifications used: {'COO': 1} Link IDs: {None: 2} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen dihedrals: 3 Unresolved non-hydrogen chiralities: 1 Chain: "D" Number of atoms: 42 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 42 Unusual residues: {'NAG': 1, 'TCI': 1} Classifications: {'peptide': 1, 'undetermined': 2} Modifications used: {'COO': 1} Link IDs: {None: 2} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen dihedrals: 3 Unresolved non-hydrogen chiralities: 1 Chain: "E" Number of atoms: 37 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 37 Unusual residues: {'NAG': 1, 'TCI': 1} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen dihedrals: 3 Unresolved non-hydrogen chiralities: 1 Time building chain proxies: 7.45, per 1000 atoms: 0.51 Number of scatterers: 14725 At special positions: 0 Unit cell: (94.6, 92.4, 139.7, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 4 Type Number sf(0) S 100 16.00 O 2645 8.00 N 2390 7.00 C 9590 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=5, symmetry=0 Simple disulfide: pdb=" SG CYS A 162 " - pdb=" SG CYS A 176 " distance=2.04 Simple disulfide: pdb=" SG CYS B 162 " - pdb=" SG CYS B 176 " distance=2.04 Simple disulfide: pdb=" SG CYS C 162 " - pdb=" SG CYS C 176 " distance=2.04 Simple disulfide: pdb=" SG CYS D 162 " - pdb=" SG CYS D 176 " distance=2.04 Simple disulfide: pdb=" SG CYS E 162 " - pdb=" SG CYS E 176 " distance=2.04 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 Links applied NAG-ASN " NAG A 802 " - " ASN A 62 " " NAG B 802 " - " ASN B 62 " " NAG C 802 " - " ASN C 62 " " NAG D 802 " - " ASN D 62 " " NAG E 501 " - " ASN E 62 " Time building additional restraints: 5.03 Conformation dependent library (CDL) restraints added in 2.4 seconds 3550 Ramachandran restraints generated. 1775 Oldfield, 0 Emsley, 1775 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 3450 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 30 helices and 15 sheets defined 36.9% alpha, 17.8% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 2.56 Creating SS restraints... Processing helix chain 'A' and resid 34 through 42 removed outlier: 3.727A pdb=" N MET A 42 " --> pdb=" O LEU A 38 " (cutoff:3.500A) Processing helix chain 'A' and resid 246 through 250 Processing helix chain 'A' and resid 252 through 268 removed outlier: 4.406A pdb=" N SER A 265 " --> pdb=" O LEU A 261 " (cutoff:3.500A) removed outlier: 4.791A pdb=" N PHE A 266 " --> pdb=" O SER A 262 " (cutoff:3.500A) removed outlier: 3.950A pdb=" N TRP A 267 " --> pdb=" O TRP A 263 " (cutoff:3.500A) removed outlier: 3.805A pdb=" N ILE A 268 " --> pdb=" O VAL A 264 " (cutoff:3.500A) Processing helix chain 'A' and resid 275 through 294 Processing helix chain 'A' and resid 306 through 340 Processing helix chain 'A' and resid 395 through 441 removed outlier: 3.962A pdb=" N MET A 399 " --> pdb=" O THR A 395 " (cutoff:3.500A) removed outlier: 3.911A pdb=" N LYS A 401 " --> pdb=" O GLU A 397 " (cutoff:3.500A) removed outlier: 4.645A pdb=" N LEU A 402 " --> pdb=" O GLU A 398 " (cutoff:3.500A) Proline residue: A 419 - end of helix removed outlier: 3.528A pdb=" N SER A 438 " --> pdb=" O LYS A 434 " (cutoff:3.500A) Processing helix chain 'B' and resid 34 through 42 removed outlier: 3.728A pdb=" N MET B 42 " --> pdb=" O LEU B 38 " (cutoff:3.500A) Processing helix chain 'B' and resid 246 through 250 Processing helix chain 'B' and resid 252 through 268 removed outlier: 4.408A pdb=" N SER B 265 " --> pdb=" O LEU B 261 " (cutoff:3.500A) removed outlier: 4.792A pdb=" N PHE B 266 " --> pdb=" O SER B 262 " (cutoff:3.500A) removed outlier: 3.949A pdb=" N TRP B 267 " --> pdb=" O TRP B 263 " (cutoff:3.500A) removed outlier: 3.805A pdb=" N ILE B 268 " --> pdb=" O VAL B 264 " (cutoff:3.500A) Processing helix chain 'B' and resid 275 through 294 Processing helix chain 'B' and resid 306 through 340 Processing helix chain 'B' and resid 395 through 441 removed outlier: 3.962A pdb=" N MET B 399 " --> pdb=" O THR B 395 " (cutoff:3.500A) removed outlier: 3.911A pdb=" N LYS B 401 " --> pdb=" O GLU B 397 " (cutoff:3.500A) removed outlier: 4.644A pdb=" N LEU B 402 " --> pdb=" O GLU B 398 " (cutoff:3.500A) Proline residue: B 419 - end of helix removed outlier: 3.528A pdb=" N SER B 438 " --> pdb=" O LYS B 434 " (cutoff:3.500A) Processing helix chain 'C' and resid 34 through 42 removed outlier: 3.727A pdb=" N MET C 42 " --> pdb=" O LEU C 38 " (cutoff:3.500A) Processing helix chain 'C' and resid 246 through 250 Processing helix chain 'C' and resid 252 through 268 removed outlier: 4.407A pdb=" N SER C 265 " --> pdb=" O LEU C 261 " (cutoff:3.500A) removed outlier: 4.792A pdb=" N PHE C 266 " --> pdb=" O SER C 262 " (cutoff:3.500A) removed outlier: 3.949A pdb=" N TRP C 267 " --> pdb=" O TRP C 263 " (cutoff:3.500A) removed outlier: 3.806A pdb=" N ILE C 268 " --> pdb=" O VAL C 264 " (cutoff:3.500A) Processing helix chain 'C' and resid 275 through 294 Processing helix chain 'C' and resid 306 through 340 Processing helix chain 'C' and resid 395 through 441 removed outlier: 3.962A pdb=" N MET C 399 " --> pdb=" O THR C 395 " (cutoff:3.500A) removed outlier: 3.910A pdb=" N LYS C 401 " --> pdb=" O GLU C 397 " (cutoff:3.500A) removed outlier: 4.645A pdb=" N LEU C 402 " --> pdb=" O GLU C 398 " (cutoff:3.500A) Proline residue: C 419 - end of helix removed outlier: 3.528A pdb=" N SER C 438 " --> pdb=" O LYS C 434 " (cutoff:3.500A) Processing helix chain 'D' and resid 34 through 42 removed outlier: 3.727A pdb=" N MET D 42 " --> pdb=" O LEU D 38 " (cutoff:3.500A) Processing helix chain 'D' and resid 246 through 250 Processing helix chain 'D' and resid 252 through 268 removed outlier: 4.407A pdb=" N SER D 265 " --> pdb=" O LEU D 261 " (cutoff:3.500A) removed outlier: 4.792A pdb=" N PHE D 266 " --> pdb=" O SER D 262 " (cutoff:3.500A) removed outlier: 3.950A pdb=" N TRP D 267 " --> pdb=" O TRP D 263 " (cutoff:3.500A) removed outlier: 3.806A pdb=" N ILE D 268 " --> pdb=" O VAL D 264 " (cutoff:3.500A) Processing helix chain 'D' and resid 275 through 294 Processing helix chain 'D' and resid 306 through 340 Processing helix chain 'D' and resid 395 through 441 removed outlier: 3.963A pdb=" N MET D 399 " --> pdb=" O THR D 395 " (cutoff:3.500A) removed outlier: 3.910A pdb=" N LYS D 401 " --> pdb=" O GLU D 397 " (cutoff:3.500A) removed outlier: 4.645A pdb=" N LEU D 402 " --> pdb=" O GLU D 398 " (cutoff:3.500A) Proline residue: D 419 - end of helix removed outlier: 3.528A pdb=" N SER D 438 " --> pdb=" O LYS D 434 " (cutoff:3.500A) Processing helix chain 'E' and resid 34 through 42 removed outlier: 3.726A pdb=" N MET E 42 " --> pdb=" O LEU E 38 " (cutoff:3.500A) Processing helix chain 'E' and resid 246 through 250 Processing helix chain 'E' and resid 252 through 268 removed outlier: 4.407A pdb=" N SER E 265 " --> pdb=" O LEU E 261 " (cutoff:3.500A) removed outlier: 4.792A pdb=" N PHE E 266 " --> pdb=" O SER E 262 " (cutoff:3.500A) removed outlier: 3.950A pdb=" N TRP E 267 " --> pdb=" O TRP E 263 " (cutoff:3.500A) removed outlier: 3.806A pdb=" N ILE E 268 " --> pdb=" O VAL E 264 " (cutoff:3.500A) Processing helix chain 'E' and resid 275 through 294 Processing helix chain 'E' and resid 306 through 340 Processing helix chain 'E' and resid 395 through 441 removed outlier: 3.962A pdb=" N MET E 399 " --> pdb=" O THR E 395 " (cutoff:3.500A) removed outlier: 3.911A pdb=" N LYS E 401 " --> pdb=" O GLU E 397 " (cutoff:3.500A) removed outlier: 4.644A pdb=" N LEU E 402 " --> pdb=" O GLU E 398 " (cutoff:3.500A) Proline residue: E 419 - end of helix removed outlier: 3.527A pdb=" N SER E 438 " --> pdb=" O LYS E 434 " (cutoff:3.500A) Processing sheet with id= A, first strand: chain 'A' and resid 192 through 194 removed outlier: 5.895A pdb=" N VAL A 63 " --> pdb=" O GLU A 193 " (cutoff:3.500A) removed outlier: 6.743A pdb=" N ASN A 85 " --> pdb=" O ASN A 70 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE A 72 " --> pdb=" O ARG A 83 " (cutoff:3.500A) removed outlier: 7.233A pdb=" N ARG A 83 " --> pdb=" O PHE A 72 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N TRP A 92 " --> pdb=" O VAL A 150 " (cutoff:3.500A) Processing sheet with id= B, first strand: chain 'A' and resid 67 through 69 Processing sheet with id= C, first strand: chain 'A' and resid 173 through 181 Processing sheet with id= D, first strand: chain 'B' and resid 192 through 194 removed outlier: 5.894A pdb=" N VAL B 63 " --> pdb=" O GLU B 193 " (cutoff:3.500A) removed outlier: 6.742A pdb=" N ASN B 85 " --> pdb=" O ASN B 70 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE B 72 " --> pdb=" O ARG B 83 " (cutoff:3.500A) removed outlier: 7.233A pdb=" N ARG B 83 " --> pdb=" O PHE B 72 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N TRP B 92 " --> pdb=" O VAL B 150 " (cutoff:3.500A) Processing sheet with id= E, first strand: chain 'B' and resid 67 through 69 Processing sheet with id= F, first strand: chain 'B' and resid 173 through 181 Processing sheet with id= G, first strand: chain 'C' and resid 192 through 194 removed outlier: 5.895A pdb=" N VAL C 63 " --> pdb=" O GLU C 193 " (cutoff:3.500A) removed outlier: 6.742A pdb=" N ASN C 85 " --> pdb=" O ASN C 70 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE C 72 " --> pdb=" O ARG C 83 " (cutoff:3.500A) removed outlier: 7.233A pdb=" N ARG C 83 " --> pdb=" O PHE C 72 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N TRP C 92 " --> pdb=" O VAL C 150 " (cutoff:3.500A) Processing sheet with id= H, first strand: chain 'C' and resid 67 through 69 Processing sheet with id= I, first strand: chain 'C' and resid 173 through 181 Processing sheet with id= J, first strand: chain 'D' and resid 192 through 194 removed outlier: 5.894A pdb=" N VAL D 63 " --> pdb=" O GLU D 193 " (cutoff:3.500A) removed outlier: 6.744A pdb=" N ASN D 85 " --> pdb=" O ASN D 70 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE D 72 " --> pdb=" O ARG D 83 " (cutoff:3.500A) removed outlier: 7.232A pdb=" N ARG D 83 " --> pdb=" O PHE D 72 " (cutoff:3.500A) removed outlier: 3.958A pdb=" N TRP D 92 " --> pdb=" O VAL D 150 " (cutoff:3.500A) Processing sheet with id= K, first strand: chain 'D' and resid 67 through 69 Processing sheet with id= L, first strand: chain 'D' and resid 173 through 181 Processing sheet with id= M, first strand: chain 'E' and resid 192 through 194 removed outlier: 5.894A pdb=" N VAL E 63 " --> pdb=" O GLU E 193 " (cutoff:3.500A) removed outlier: 6.742A pdb=" N ASN E 85 " --> pdb=" O ASN E 70 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE E 72 " --> pdb=" O ARG E 83 " (cutoff:3.500A) removed outlier: 7.232A pdb=" N ARG E 83 " --> pdb=" O PHE E 72 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N TRP E 92 " --> pdb=" O VAL E 150 " (cutoff:3.500A) Processing sheet with id= N, first strand: chain 'E' and resid 67 through 69 Processing sheet with id= O, first strand: chain 'E' and resid 173 through 181 665 hydrogen bonds defined for protein. 1980 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 0 stacking parallelities Total time for adding SS restraints: 4.40 Time building geometry restraints manager: 5.55 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.34: 4487 1.34 - 1.46: 3294 1.46 - 1.58: 7134 1.58 - 1.70: 0 1.70 - 1.82: 170 Bond restraints: 15085 Sorted by residual: bond pdb=" CB VAL C 259 " pdb=" CG2 VAL C 259 " ideal model delta sigma weight residual 1.521 1.460 0.061 3.30e-02 9.18e+02 3.46e+00 bond pdb=" CB VAL B 259 " pdb=" CG2 VAL B 259 " ideal model delta sigma weight residual 1.521 1.460 0.061 3.30e-02 9.18e+02 3.40e+00 bond pdb=" CB VAL E 259 " pdb=" CG2 VAL E 259 " ideal model delta sigma weight residual 1.521 1.460 0.061 3.30e-02 9.18e+02 3.40e+00 bond pdb=" CB VAL D 259 " pdb=" CG2 VAL D 259 " ideal model delta sigma weight residual 1.521 1.460 0.061 3.30e-02 9.18e+02 3.38e+00 bond pdb=" CB VAL A 259 " pdb=" CG2 VAL A 259 " ideal model delta sigma weight residual 1.521 1.461 0.060 3.30e-02 9.18e+02 3.32e+00 ... (remaining 15080 not shown) Histogram of bond angle deviations from ideal: 98.90 - 105.97: 362 105.97 - 113.03: 7998 113.03 - 120.09: 5513 120.09 - 127.15: 6417 127.15 - 134.21: 175 Bond angle restraints: 20465 Sorted by residual: angle pdb=" CA GLU B 215 " pdb=" CB GLU B 215 " pdb=" CG GLU B 215 " ideal model delta sigma weight residual 114.10 122.90 -8.80 2.00e+00 2.50e-01 1.94e+01 angle pdb=" CA GLU A 215 " pdb=" CB GLU A 215 " pdb=" CG GLU A 215 " ideal model delta sigma weight residual 114.10 122.86 -8.76 2.00e+00 2.50e-01 1.92e+01 angle pdb=" CA GLU C 215 " pdb=" CB GLU C 215 " pdb=" CG GLU C 215 " ideal model delta sigma weight residual 114.10 122.86 -8.76 2.00e+00 2.50e-01 1.92e+01 angle pdb=" CA GLU D 215 " pdb=" CB GLU D 215 " pdb=" CG GLU D 215 " ideal model delta sigma weight residual 114.10 122.85 -8.75 2.00e+00 2.50e-01 1.91e+01 angle pdb=" CA GLU E 215 " pdb=" CB GLU E 215 " pdb=" CG GLU E 215 " ideal model delta sigma weight residual 114.10 122.84 -8.74 2.00e+00 2.50e-01 1.91e+01 ... (remaining 20460 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 35.91: 8960 35.91 - 71.82: 170 71.82 - 107.74: 15 107.74 - 143.65: 0 143.65 - 179.56: 10 Dihedral angle restraints: 9155 sinusoidal: 3845 harmonic: 5310 Sorted by residual: dihedral pdb=" C11 TCI D 803 " pdb=" C12 TCI D 803 " pdb=" C13 TCI D 803 " pdb=" C15 TCI D 803 " ideal model delta sinusoidal sigma weight residual 64.70 -114.86 179.56 1 3.00e+01 1.11e-03 2.13e+01 dihedral pdb=" C11 TCI A 803 " pdb=" C12 TCI A 803 " pdb=" C13 TCI A 803 " pdb=" C15 TCI A 803 " ideal model delta sinusoidal sigma weight residual 64.70 -114.84 179.54 1 3.00e+01 1.11e-03 2.13e+01 dihedral pdb=" C11 TCI E 502 " pdb=" C12 TCI E 502 " pdb=" C13 TCI E 502 " pdb=" C15 TCI E 502 " ideal model delta sinusoidal sigma weight residual 64.70 -114.84 179.54 1 3.00e+01 1.11e-03 2.13e+01 ... (remaining 9152 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.088: 1868 0.088 - 0.177: 411 0.177 - 0.265: 11 0.265 - 0.354: 0 0.354 - 0.442: 10 Chirality restraints: 2300 Sorted by residual: chirality pdb=" CA PHE E 418 " pdb=" N PHE E 418 " pdb=" C PHE E 418 " pdb=" CB PHE E 418 " both_signs ideal model delta sigma weight residual False 2.51 2.07 0.44 2.00e-01 2.50e+01 4.88e+00 chirality pdb=" CA PHE A 418 " pdb=" N PHE A 418 " pdb=" C PHE A 418 " pdb=" CB PHE A 418 " both_signs ideal model delta sigma weight residual False 2.51 2.07 0.44 2.00e-01 2.50e+01 4.83e+00 chirality pdb=" CA PHE C 418 " pdb=" N PHE C 418 " pdb=" C PHE C 418 " pdb=" CB PHE C 418 " both_signs ideal model delta sigma weight residual False 2.51 2.07 0.44 2.00e-01 2.50e+01 4.79e+00 ... (remaining 2297 not shown) Planarity restraints: 2550 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA THR E 283 " 0.021 2.00e-02 2.50e+03 4.11e-02 1.69e+01 pdb=" C THR E 283 " -0.071 2.00e-02 2.50e+03 pdb=" O THR E 283 " 0.026 2.00e-02 2.50e+03 pdb=" N VAL E 284 " 0.024 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA THR A 283 " 0.021 2.00e-02 2.50e+03 4.10e-02 1.68e+01 pdb=" C THR A 283 " -0.071 2.00e-02 2.50e+03 pdb=" O THR A 283 " 0.026 2.00e-02 2.50e+03 pdb=" N VAL A 284 " 0.024 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA THR C 283 " -0.020 2.00e-02 2.50e+03 4.07e-02 1.66e+01 pdb=" C THR C 283 " 0.070 2.00e-02 2.50e+03 pdb=" O THR C 283 " -0.026 2.00e-02 2.50e+03 pdb=" N VAL C 284 " -0.024 2.00e-02 2.50e+03 ... (remaining 2547 not shown) Histogram of nonbonded interaction distances: 2.29 - 2.81: 3832 2.81 - 3.33: 13893 3.33 - 3.86: 23588 3.86 - 4.38: 25477 4.38 - 4.90: 44779 Nonbonded interactions: 111569 Sorted by model distance: nonbonded pdb=" OG1 THR C 186 " pdb=" OD1 ASN C 188 " model vdw 2.289 2.440 nonbonded pdb=" OG1 THR D 186 " pdb=" OD1 ASN D 188 " model vdw 2.289 2.440 nonbonded pdb=" OG1 THR A 186 " pdb=" OD1 ASN A 188 " model vdw 2.289 2.440 nonbonded pdb=" OG1 THR E 186 " pdb=" OD1 ASN E 188 " model vdw 2.290 2.440 nonbonded pdb=" OG1 THR B 186 " pdb=" OD1 ASN B 188 " model vdw 2.290 2.440 ... (remaining 111564 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 32 through 442 or resid 802 through 803)) selection = (chain 'B' and (resid 32 through 442 or resid 802 through 803)) selection = (chain 'C' and (resid 32 through 442 or resid 802 through 803)) selection = (chain 'D' and (resid 32 through 442 or resid 802 through 803)) selection = chain 'E' } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 1.300 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: 6.360 Check model and map are aligned: 0.200 Set scattering table: 0.120 Process input model: 38.220 Find NCS groups from input model: 1.180 Set up NCS constraints: 0.110 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.010 Load rotamer database and sin/cos tables:12.200 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 59.710 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8062 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.071 15085 Z= 0.563 Angle : 1.159 11.390 20465 Z= 0.632 Chirality : 0.070 0.442 2300 Planarity : 0.008 0.065 2545 Dihedral : 15.900 179.560 5690 Min Nonbonded Distance : 2.289 Molprobity Statistics. All-atom Clashscore : 5.71 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.10 % Favored : 96.90 % Rotamer: Outliers : 0.00 % Allowed : 3.45 % Favored : 96.55 % Cbeta Deviations : 0.29 % 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.39 (0.19), residues: 1775 helix: 1.54 (0.19), residues: 665 sheet: 0.19 (0.27), residues: 380 loop : -2.33 (0.19), residues: 730 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.002 TRP C 430 HIS 0.004 0.001 HIS C 442 PHE 0.033 0.004 PHE E 169 TYR 0.027 0.003 TYR A 433 ARG 0.010 0.001 ARG C 220 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3550 Ramachandran restraints generated. 1775 Oldfield, 0 Emsley, 1775 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3550 Ramachandran restraints generated. 1775 Oldfield, 0 Emsley, 1775 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 230 residues out of total 1620 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 230 time to evaluate : 1.641 Fit side-chains REVERT: A 267 TRP cc_start: 0.6770 (t-100) cc_final: 0.6330 (t-100) REVERT: A 428 PHE cc_start: 0.8463 (t80) cc_final: 0.8139 (m-80) REVERT: A 437 ARG cc_start: 0.7604 (mtp180) cc_final: 0.6383 (mtt180) REVERT: B 266 PHE cc_start: 0.6845 (m-10) cc_final: 0.6068 (m-80) REVERT: B 287 MET cc_start: 0.9094 (tpp) cc_final: 0.8809 (tpt) REVERT: B 428 PHE cc_start: 0.8589 (t80) cc_final: 0.8012 (t80) REVERT: B 437 ARG cc_start: 0.7359 (mtp180) cc_final: 0.6465 (mtt180) REVERT: C 257 LEU cc_start: 0.8513 (tp) cc_final: 0.8140 (mt) REVERT: C 287 MET cc_start: 0.8974 (tpp) cc_final: 0.8633 (tpt) REVERT: C 428 PHE cc_start: 0.8587 (t80) cc_final: 0.7996 (t80) REVERT: C 437 ARG cc_start: 0.7408 (mtp180) cc_final: 0.6335 (mtm180) REVERT: D 190 LEU cc_start: 0.8641 (tt) cc_final: 0.8415 (tm) REVERT: D 287 MET cc_start: 0.9063 (tpp) cc_final: 0.8850 (tpt) REVERT: D 428 PHE cc_start: 0.8651 (t80) cc_final: 0.8095 (t80) REVERT: D 437 ARG cc_start: 0.7423 (mtp180) cc_final: 0.6463 (mtm180) REVERT: E 113 MET cc_start: 0.8532 (pmm) cc_final: 0.8165 (ppp) REVERT: E 308 ASP cc_start: 0.8690 (m-30) cc_final: 0.8356 (m-30) REVERT: E 428 PHE cc_start: 0.8630 (t80) cc_final: 0.8110 (t80) REVERT: E 437 ARG cc_start: 0.7489 (mtp180) cc_final: 0.6586 (mtm180) outliers start: 0 outliers final: 6 residues processed: 230 average time/residue: 1.2149 time to fit residues: 307.2680 Evaluate side-chains 132 residues out of total 1620 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 126 time to evaluate : 1.766 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 115 ASP Chi-restraints excluded: chain A residue 154 ILE Chi-restraints excluded: chain B residue 204 ASP Chi-restraints excluded: chain D residue 115 ASP Chi-restraints excluded: chain D residue 204 ASP Chi-restraints excluded: chain E residue 115 ASP Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 8, in run_program(real_space_refine.Program) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/iotbx/cli_parser.py", line 944, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/phenix/phenix/programs/real_space_refine.py", line 189, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 306, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 700, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 764, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 51, in __init__ gradients_method = gradients_method) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/mmtbx/refinement/real_space/weight.py", line 87, in __init__ rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 368, in refine rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 112, in refine states_collector = self.states_accumulator) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 174, in __init__ exception_handling_params=lbfgs_exception_handling_params) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 131, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 245, in compute_functional_and_gradients compute_gradients=True) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1479, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1381, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5265/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1262, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 51.7431 > 50: distance: 49 - 53: 14.071 distance: 53 - 54: 10.838 distance: 54 - 55: 3.011 distance: 54 - 57: 5.023 distance: 55 - 62: 14.267 distance: 57 - 58: 11.824 distance: 58 - 59: 15.045 distance: 59 - 60: 12.340 distance: 59 - 61: 17.338 distance: 62 - 63: 17.980 distance: 63 - 64: 22.765 distance: 63 - 66: 10.283 distance: 64 - 65: 26.255 distance: 64 - 74: 21.856 distance: 66 - 67: 3.618 distance: 67 - 68: 6.370 distance: 67 - 69: 16.869 distance: 68 - 70: 7.200 distance: 69 - 71: 13.735 distance: 70 - 72: 15.460 distance: 71 - 72: 7.918 distance: 72 - 73: 11.487 distance: 74 - 75: 11.328 distance: 74 - 80: 12.573 distance: 75 - 76: 8.557 distance: 75 - 78: 8.234 distance: 76 - 77: 14.308 distance: 76 - 81: 16.213 distance: 78 - 79: 5.363 distance: 79 - 80: 13.604 distance: 81 - 82: 12.307 distance: 82 - 83: 6.540 distance: 83 - 84: 12.405 distance: 83 - 89: 14.556 distance: 85 - 86: 12.998 distance: 86 - 87: 14.095 distance: 86 - 88: 9.981 distance: 89 - 90: 9.884 distance: 90 - 91: 12.454 distance: 90 - 93: 17.291 distance: 91 - 92: 16.750 distance: 91 - 97: 15.951 distance: 93 - 94: 25.743 distance: 94 - 95: 5.132 distance: 94 - 96: 31.223 distance: 97 - 98: 11.802 distance: 98 - 99: 4.386 distance: 98 - 101: 15.963 distance: 99 - 100: 10.057 distance: 99 - 103: 15.928 distance: 101 - 102: 9.531 distance: 103 - 104: 14.998 distance: 103 - 172: 15.138 distance: 104 - 105: 5.905 distance: 104 - 107: 9.534 distance: 105 - 106: 11.596 distance: 105 - 111: 5.516 distance: 106 - 169: 9.051 distance: 107 - 108: 5.580 distance: 108 - 109: 6.008 distance: 108 - 110: 5.345 distance: 111 - 112: 5.335 distance: 112 - 113: 15.748 distance: 112 - 115: 7.393 distance: 113 - 114: 31.632 distance: 113 - 119: 18.756 distance: 115 - 116: 7.868 distance: 116 - 117: 13.070 distance: 116 - 118: 7.163 distance: 119 - 120: 8.535 distance: 120 - 121: 10.090 distance: 120 - 123: 11.284 distance: 121 - 122: 11.633 distance: 121 - 127: 5.148 distance: 123 - 124: 18.385 distance: 124 - 125: 11.076 distance: 124 - 126: 16.856 distance: 127 - 128: 8.460 distance: 128 - 129: 9.445 distance: 128 - 131: 7.261 distance: 129 - 130: 14.624 distance: 129 - 135: 5.654 distance: 131 - 132: 6.249 distance: 132 - 133: 6.925 distance: 132 - 134: 7.260