Starting phenix.real_space_refine on Sun Jul 21 17:16:29 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/07_2024/7m6q_23704.cif Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/07_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/07_2024/7m6q_23704.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/07_2024/7m6q_23704.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/07_2024/7m6q_23704.cif" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7m6q_23704/07_2024/7m6q_23704.cif" } 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-5388/modules/chem_data/mon_lib" Total number of atoms: 14725 Number of models: 1 Model: "" Number of chains: 14 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: 5 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 5 Classifications: {'peptide': 1} Modifications used: {'COO': 1} Chain: "A" 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: "B" Number of atoms: 5 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 5 Classifications: {'peptide': 1} Modifications used: {'COO': 1} Chain: "B" 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 Chain: "C" Number of atoms: 5 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 5 Classifications: {'peptide': 1} Modifications used: {'COO': 1} Chain: "C" 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 Chain: "D" Number of atoms: 5 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 5 Classifications: {'peptide': 1} Modifications used: {'COO': 1} Chain: "D" 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 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: 8.30, per 1000 atoms: 0.56 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.30 Conformation dependent library (CDL) restraints added in 2.6 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... Secondary structure from input PDB file: 30 helices and 25 sheets defined 40.0% alpha, 33.3% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.97 Creating SS restraints... Processing helix chain 'A' and resid 33 through 43 removed outlier: 3.545A pdb=" N PHE A 37 " --> pdb=" O PRO A 33 " (cutoff:3.500A) removed outlier: 3.727A pdb=" N MET A 42 " --> pdb=" O LEU A 38 " (cutoff:3.500A) Processing helix chain 'A' and resid 245 through 251 removed outlier: 3.893A pdb=" N ILE A 249 " --> pdb=" O GLY A 245 " (cutoff:3.500A) Processing helix chain 'A' and resid 251 through 269 removed outlier: 3.510A pdb=" N SER A 255 " --> pdb=" O MET A 251 " (cutoff:3.500A) 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 274 through 295 removed outlier: 3.812A pdb=" N GLY A 278 " --> pdb=" O PRO A 274 " (cutoff:3.500A) removed outlier: 3.771A pdb=" N ARG A 295 " --> pdb=" O SER A 291 " (cutoff:3.500A) Processing helix chain 'A' and resid 305 through 341 Processing helix chain 'A' and resid 395 through 442 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) removed outlier: 3.886A pdb=" N HIS A 442 " --> pdb=" O SER A 438 " (cutoff:3.500A) Processing helix chain 'B' and resid 33 through 43 removed outlier: 3.544A pdb=" N PHE B 37 " --> pdb=" O PRO B 33 " (cutoff:3.500A) removed outlier: 3.728A pdb=" N MET B 42 " --> pdb=" O LEU B 38 " (cutoff:3.500A) Processing helix chain 'B' and resid 245 through 251 removed outlier: 3.893A pdb=" N ILE B 249 " --> pdb=" O GLY B 245 " (cutoff:3.500A) Processing helix chain 'B' and resid 251 through 269 removed outlier: 3.510A pdb=" N SER B 255 " --> pdb=" O MET B 251 " (cutoff:3.500A) 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 274 through 295 removed outlier: 3.812A pdb=" N GLY B 278 " --> pdb=" O PRO B 274 " (cutoff:3.500A) removed outlier: 3.771A pdb=" N ARG B 295 " --> pdb=" O SER B 291 " (cutoff:3.500A) Processing helix chain 'B' and resid 305 through 341 Processing helix chain 'B' and resid 395 through 442 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) removed outlier: 3.886A pdb=" N HIS B 442 " --> pdb=" O SER B 438 " (cutoff:3.500A) Processing helix chain 'C' and resid 33 through 43 removed outlier: 3.544A pdb=" N PHE C 37 " --> pdb=" O PRO C 33 " (cutoff:3.500A) removed outlier: 3.727A pdb=" N MET C 42 " --> pdb=" O LEU C 38 " (cutoff:3.500A) Processing helix chain 'C' and resid 245 through 251 removed outlier: 3.894A pdb=" N ILE C 249 " --> pdb=" O GLY C 245 " (cutoff:3.500A) Processing helix chain 'C' and resid 251 through 269 removed outlier: 3.511A pdb=" N SER C 255 " --> pdb=" O MET C 251 " (cutoff:3.500A) 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 274 through 295 removed outlier: 3.812A pdb=" N GLY C 278 " --> pdb=" O PRO C 274 " (cutoff:3.500A) removed outlier: 3.771A pdb=" N ARG C 295 " --> pdb=" O SER C 291 " (cutoff:3.500A) Processing helix chain 'C' and resid 305 through 341 Processing helix chain 'C' and resid 395 through 442 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) removed outlier: 3.886A pdb=" N HIS C 442 " --> pdb=" O SER C 438 " (cutoff:3.500A) Processing helix chain 'D' and resid 33 through 43 removed outlier: 3.545A pdb=" N PHE D 37 " --> pdb=" O PRO D 33 " (cutoff:3.500A) removed outlier: 3.727A pdb=" N MET D 42 " --> pdb=" O LEU D 38 " (cutoff:3.500A) Processing helix chain 'D' and resid 245 through 251 removed outlier: 3.894A pdb=" N ILE D 249 " --> pdb=" O GLY D 245 " (cutoff:3.500A) Processing helix chain 'D' and resid 251 through 269 removed outlier: 3.511A pdb=" N SER D 255 " --> pdb=" O MET D 251 " (cutoff:3.500A) 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 274 through 295 removed outlier: 3.812A pdb=" N GLY D 278 " --> pdb=" O PRO D 274 " (cutoff:3.500A) removed outlier: 3.772A pdb=" N ARG D 295 " --> pdb=" O SER D 291 " (cutoff:3.500A) Processing helix chain 'D' and resid 305 through 341 Processing helix chain 'D' and resid 395 through 442 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) removed outlier: 3.887A pdb=" N HIS D 442 " --> pdb=" O SER D 438 " (cutoff:3.500A) Processing helix chain 'E' and resid 33 through 43 removed outlier: 3.545A pdb=" N PHE E 37 " --> pdb=" O PRO E 33 " (cutoff:3.500A) removed outlier: 3.726A pdb=" N MET E 42 " --> pdb=" O LEU E 38 " (cutoff:3.500A) Processing helix chain 'E' and resid 245 through 251 removed outlier: 3.894A pdb=" N ILE E 249 " --> pdb=" O GLY E 245 " (cutoff:3.500A) Processing helix chain 'E' and resid 251 through 269 removed outlier: 3.511A pdb=" N SER E 255 " --> pdb=" O MET E 251 " (cutoff:3.500A) 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 274 through 295 removed outlier: 3.812A pdb=" N GLY E 278 " --> pdb=" O PRO E 274 " (cutoff:3.500A) removed outlier: 3.771A pdb=" N ARG E 295 " --> pdb=" O SER E 291 " (cutoff:3.500A) Processing helix chain 'E' and resid 305 through 341 Processing helix chain 'E' and resid 395 through 442 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) removed outlier: 3.886A pdb=" N HIS E 442 " --> pdb=" O SER E 438 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 106 through 108 removed outlier: 3.959A pdb=" N TRP A 92 " --> pdb=" O VAL A 150 " (cutoff:3.500A) removed outlier: 7.233A pdb=" N ARG A 83 " --> pdb=" O PHE A 72 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE A 72 " --> pdb=" O ARG A 83 " (cutoff:3.500A) removed outlier: 6.743A pdb=" N ASN A 85 " --> pdb=" O ASN A 70 " (cutoff:3.500A) removed outlier: 5.895A pdb=" N VAL A 63 " --> pdb=" O GLU A 193 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 127 through 131 removed outlier: 3.661A pdb=" N VAL A 159 " --> pdb=" O LYS A 128 " (cutoff:3.500A) removed outlier: 3.921A pdb=" N ALA A 130 " --> pdb=" O THR A 157 " (cutoff:3.500A) removed outlier: 6.096A pdb=" N THR A 157 " --> pdb=" O ALA A 130 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N TRP A 92 " --> pdb=" O VAL A 150 " (cutoff:3.500A) removed outlier: 7.233A pdb=" N ARG A 83 " --> pdb=" O PHE A 72 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE A 72 " --> pdb=" O ARG A 83 " (cutoff:3.500A) removed outlier: 6.743A pdb=" N ASN A 85 " --> pdb=" O ASN A 70 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 123 through 124 Processing sheet with id=AA4, first strand: chain 'A' and resid 123 through 124 removed outlier: 3.565A pdb=" N ARG A 220 " --> pdb=" O GLU A 235 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'A' and resid 224 through 225 Processing sheet with id=AA6, first strand: chain 'B' and resid 106 through 108 removed outlier: 3.959A pdb=" N TRP B 92 " --> pdb=" O VAL B 150 " (cutoff:3.500A) removed outlier: 7.233A pdb=" N ARG B 83 " --> pdb=" O PHE B 72 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE B 72 " --> pdb=" O ARG B 83 " (cutoff:3.500A) removed outlier: 6.742A pdb=" N ASN B 85 " --> pdb=" O ASN B 70 " (cutoff:3.500A) removed outlier: 5.894A pdb=" N VAL B 63 " --> pdb=" O GLU B 193 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 127 through 131 removed outlier: 3.662A pdb=" N VAL B 159 " --> pdb=" O LYS B 128 " (cutoff:3.500A) removed outlier: 3.921A pdb=" N ALA B 130 " --> pdb=" O THR B 157 " (cutoff:3.500A) removed outlier: 6.095A pdb=" N THR B 157 " --> pdb=" O ALA B 130 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N TRP B 92 " --> pdb=" O VAL B 150 " (cutoff:3.500A) removed outlier: 7.233A pdb=" N ARG B 83 " --> pdb=" O PHE B 72 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE B 72 " --> pdb=" O ARG B 83 " (cutoff:3.500A) removed outlier: 6.742A pdb=" N ASN B 85 " --> pdb=" O ASN B 70 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'B' and resid 123 through 124 Processing sheet with id=AA9, first strand: chain 'B' and resid 123 through 124 removed outlier: 3.565A pdb=" N ARG B 220 " --> pdb=" O GLU B 235 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'B' and resid 224 through 225 Processing sheet with id=AB2, first strand: chain 'C' and resid 106 through 108 removed outlier: 3.959A pdb=" N TRP C 92 " --> pdb=" O VAL C 150 " (cutoff:3.500A) removed outlier: 7.233A pdb=" N ARG C 83 " --> pdb=" O PHE C 72 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE C 72 " --> pdb=" O ARG C 83 " (cutoff:3.500A) removed outlier: 6.742A pdb=" N ASN C 85 " --> pdb=" O ASN C 70 " (cutoff:3.500A) removed outlier: 5.895A pdb=" N VAL C 63 " --> pdb=" O GLU C 193 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'C' and resid 127 through 131 removed outlier: 3.663A pdb=" N VAL C 159 " --> pdb=" O LYS C 128 " (cutoff:3.500A) removed outlier: 3.922A pdb=" N ALA C 130 " --> pdb=" O THR C 157 " (cutoff:3.500A) removed outlier: 6.095A pdb=" N THR C 157 " --> pdb=" O ALA C 130 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N TRP C 92 " --> pdb=" O VAL C 150 " (cutoff:3.500A) removed outlier: 7.233A pdb=" N ARG C 83 " --> pdb=" O PHE C 72 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE C 72 " --> pdb=" O ARG C 83 " (cutoff:3.500A) removed outlier: 6.742A pdb=" N ASN C 85 " --> pdb=" O ASN C 70 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'C' and resid 123 through 124 Processing sheet with id=AB5, first strand: chain 'C' and resid 123 through 124 removed outlier: 3.565A pdb=" N ARG C 220 " --> pdb=" O GLU C 235 " (cutoff:3.500A) Processing sheet with id=AB6, first strand: chain 'C' and resid 224 through 225 Processing sheet with id=AB7, first strand: chain 'D' and resid 106 through 108 removed outlier: 3.958A pdb=" N TRP D 92 " --> pdb=" O VAL D 150 " (cutoff:3.500A) removed outlier: 7.232A pdb=" N ARG D 83 " --> pdb=" O PHE D 72 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE D 72 " --> pdb=" O ARG D 83 " (cutoff:3.500A) removed outlier: 6.744A pdb=" N ASN D 85 " --> pdb=" O ASN D 70 " (cutoff:3.500A) removed outlier: 5.894A pdb=" N VAL D 63 " --> pdb=" O GLU D 193 " (cutoff:3.500A) Processing sheet with id=AB8, first strand: chain 'D' and resid 127 through 131 removed outlier: 3.662A pdb=" N VAL D 159 " --> pdb=" O LYS D 128 " (cutoff:3.500A) removed outlier: 3.921A pdb=" N ALA D 130 " --> pdb=" O THR D 157 " (cutoff:3.500A) removed outlier: 6.095A pdb=" N THR D 157 " --> pdb=" O ALA D 130 " (cutoff:3.500A) removed outlier: 3.958A pdb=" N TRP D 92 " --> pdb=" O VAL D 150 " (cutoff:3.500A) removed outlier: 7.232A pdb=" N ARG D 83 " --> pdb=" O PHE D 72 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE D 72 " --> pdb=" O ARG D 83 " (cutoff:3.500A) removed outlier: 6.744A pdb=" N ASN D 85 " --> pdb=" O ASN D 70 " (cutoff:3.500A) Processing sheet with id=AB9, first strand: chain 'D' and resid 123 through 124 Processing sheet with id=AC1, first strand: chain 'D' and resid 123 through 124 removed outlier: 3.565A pdb=" N ARG D 220 " --> pdb=" O GLU D 235 " (cutoff:3.500A) Processing sheet with id=AC2, first strand: chain 'D' and resid 224 through 225 Processing sheet with id=AC3, first strand: chain 'E' and resid 106 through 108 removed outlier: 3.959A pdb=" N TRP E 92 " --> pdb=" O VAL E 150 " (cutoff:3.500A) removed outlier: 7.232A pdb=" N ARG E 83 " --> pdb=" O PHE E 72 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE E 72 " --> pdb=" O ARG E 83 " (cutoff:3.500A) removed outlier: 6.742A pdb=" N ASN E 85 " --> pdb=" O ASN E 70 " (cutoff:3.500A) removed outlier: 5.894A pdb=" N VAL E 63 " --> pdb=" O GLU E 193 " (cutoff:3.500A) Processing sheet with id=AC4, first strand: chain 'E' and resid 127 through 131 removed outlier: 3.662A pdb=" N VAL E 159 " --> pdb=" O LYS E 128 " (cutoff:3.500A) removed outlier: 3.921A pdb=" N ALA E 130 " --> pdb=" O THR E 157 " (cutoff:3.500A) removed outlier: 6.095A pdb=" N THR E 157 " --> pdb=" O ALA E 130 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N TRP E 92 " --> pdb=" O VAL E 150 " (cutoff:3.500A) removed outlier: 7.232A pdb=" N ARG E 83 " --> pdb=" O PHE E 72 " (cutoff:3.500A) removed outlier: 4.897A pdb=" N PHE E 72 " --> pdb=" O ARG E 83 " (cutoff:3.500A) removed outlier: 6.742A pdb=" N ASN E 85 " --> pdb=" O ASN E 70 " (cutoff:3.500A) Processing sheet with id=AC5, first strand: chain 'E' and resid 123 through 124 Processing sheet with id=AC6, first strand: chain 'E' and resid 123 through 124 removed outlier: 3.566A pdb=" N ARG E 220 " --> pdb=" O GLU E 235 " (cutoff:3.500A) Processing sheet with id=AC7, first strand: chain 'E' and resid 224 through 225 860 hydrogen bonds defined for protein. 2505 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: 5.02 Time building geometry restraints manager: 5.83 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: 13838 3.33 - 3.86: 23538 3.86 - 4.38: 25372 4.38 - 4.90: 44769 Nonbonded interactions: 111349 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 111344 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.810 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.010 Extract box with map and model: 0.580 Check model and map are aligned: 0.110 Set scattering table: 0.130 Process input model: 40.110 Find NCS groups from input model: 1.240 Set up NCS constraints: 0.070 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:13.120 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 57.180 ------------------------------------------------------------------------------- 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.561 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.779 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.3468 time to fit residues: 339.6019 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.711 ------------------------------------------------------------------------------- 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-5388/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-5388/modules/cctbx_project/iotbx/cli_parser.py", line 946, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/modules/phenix/phenix/programs/real_space_refine.py", line 200, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/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-5388/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 702, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 766, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/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-5388/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-5388/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-5388/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-5388/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-5388/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-5388/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/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-5388/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-5388/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1477, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1379, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1260, 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