Starting phenix.real_space_refine (version: dev) on Fri Mar 17 14:26:50 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7a4a_11630/03_2023/7a4a_11630.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7a4a_11630/03_2023/7a4a_11630.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.76 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7a4a_11630/03_2023/7a4a_11630.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7a4a_11630/03_2023/7a4a_11630.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7a4a_11630/03_2023/7a4a_11630.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7a4a_11630/03_2023/7a4a_11630.pdb" } resolution = 3.76 write_initial_geo_file = False refinement { macro_cycles = 10 } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= -0.000 sd= 0.004 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 Process input model Symmetric amino acids flipped Residue "A ARG 50": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A ARG 65": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A ARG 75": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A ARG 86": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A ARG 117": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A ARG 146": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A TYR 148": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "A ARG 177": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A ARG 187": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A ARG 192": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A ARG 223": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A PHE 255": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "A PHE 259": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "A PHE 301": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "A ARG 305": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A ARG 381": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "A TYR 433": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "B ARG 50": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B ARG 65": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B ARG 75": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B ARG 86": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B ARG 117": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B ARG 146": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B TYR 148": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "B ARG 177": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B ARG 187": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B ARG 192": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B ARG 223": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B PHE 255": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "B PHE 259": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "B PHE 301": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "B ARG 305": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B ARG 381": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "B TYR 433": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "C ARG 50": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C ARG 65": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C ARG 75": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C ARG 86": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C ARG 117": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C ARG 146": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C TYR 148": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "C ARG 177": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C ARG 187": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C ARG 192": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C ARG 223": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C PHE 255": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "C PHE 259": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "C PHE 301": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "C ARG 305": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C ARG 381": "NH1" <-> "NH2" "HH11" <-> "HH21" "HH12" <-> "HH22" Residue "C TYR 433": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Time to flip residues: 0.04s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4893/modules/chem_data/mon_lib" Total number of atoms: 19785 Number of models: 1 Model: "" Number of chains: 6 Chain: "A" Number of atoms: 6581 Number of conformers: 1 Conformer: "" Number of residues, atoms: 439, 6581 Classifications: {'peptide': 439} Modifications used: {'NH3': 1} Link IDs: {'PCIS': 1, 'PTRANS': 17, 'TRANS': 420} Chain: "B" Number of atoms: 6581 Number of conformers: 1 Conformer: "" Number of residues, atoms: 439, 6581 Classifications: {'peptide': 439} Modifications used: {'NH3': 1} Link IDs: {'PCIS': 1, 'PTRANS': 17, 'TRANS': 420} Chain: "C" Number of atoms: 6581 Number of conformers: 1 Conformer: "" Number of residues, atoms: 439, 6581 Classifications: {'peptide': 439} Modifications used: {'NH3': 1} Link IDs: {'PCIS': 1, 'PTRANS': 17, 'TRANS': 420} Chain: "A" Number of atoms: 14 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 14 Unusual residues: {'NAG': 1} Classifications: {'undetermined': 1} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen chiralities: 1 Chain: "B" Number of atoms: 14 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 14 Unusual residues: {'NAG': 1} Classifications: {'undetermined': 1} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen chiralities: 1 Chain: "C" Number of atoms: 14 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 14 Unusual residues: {'NAG': 1} Classifications: {'undetermined': 1} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen chiralities: 1 Time building chain proxies: 12.00, per 1000 atoms: 0.61 Number of scatterers: 19785 At special positions: 0 Unit cell: (86.901, 86.901, 134.016, 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 105 16.00 O 2025 8.00 N 1743 7.00 C 6294 6.00 H 9618 1.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=45, symmetry=0 Simple disulfide: pdb=" SG CYS A 1 " - pdb=" SG CYS A 41 " distance=2.03 Simple disulfide: pdb=" SG CYS A 14 " - pdb=" SG CYS A 23 " distance=2.03 Simple disulfide: pdb=" SG CYS A 66 " - pdb=" SG CYS A 162 " distance=2.02 Simple disulfide: pdb=" SG CYS A 87 " - pdb=" SG CYS A 135 " distance=2.04 Simple disulfide: pdb=" SG CYS A 93 " - pdb=" SG CYS A 142 " distance=2.03 Simple disulfide: pdb=" SG CYS A 98 " - pdb=" SG CYS A 123 " distance=2.03 Simple disulfide: pdb=" SG CYS A 127 " - pdb=" SG CYS A 132 " distance=2.03 Simple disulfide: pdb=" SG CYS A 129 " - pdb=" SG CYS A 138 " distance=2.04 Simple disulfide: pdb=" SG CYS A 246 " - pdb=" SG CYS A 257 " distance=2.03 Simple disulfide: pdb=" SG CYS A 264 " - pdb=" SG CYS A 277 " distance=2.03 Simple disulfide: pdb=" SG CYS A 266 " - pdb=" SG CYS A 275 " distance=2.03 Simple disulfide: pdb=" SG CYS A 337 " - pdb=" SG CYS A 408 " distance=2.03 Simple disulfide: pdb=" SG CYS A 347 " - pdb=" SG CYS A 350 " distance=2.03 Simple disulfide: pdb=" SG CYS A 360 " - pdb=" SG CYS A 382 " distance=2.03 Simple disulfide: pdb=" SG CYS A 373 " - pdb=" SG CYS A 404 " distance=2.03 Simple disulfide: pdb=" SG CYS B 1 " - pdb=" SG CYS B 41 " distance=2.03 Simple disulfide: pdb=" SG CYS B 14 " - pdb=" SG CYS B 23 " distance=2.04 Simple disulfide: pdb=" SG CYS B 66 " - pdb=" SG CYS B 162 " distance=2.02 Simple disulfide: pdb=" SG CYS B 87 " - pdb=" SG CYS B 135 " distance=2.04 Simple disulfide: pdb=" SG CYS B 93 " - pdb=" SG CYS B 142 " distance=2.03 Simple disulfide: pdb=" SG CYS B 98 " - pdb=" SG CYS B 123 " distance=2.03 Simple disulfide: pdb=" SG CYS B 127 " - pdb=" SG CYS B 132 " distance=2.03 Simple disulfide: pdb=" SG CYS B 129 " - pdb=" SG CYS B 138 " distance=2.04 Simple disulfide: pdb=" SG CYS B 246 " - pdb=" SG CYS B 257 " distance=2.03 Simple disulfide: pdb=" SG CYS B 264 " - pdb=" SG CYS B 277 " distance=2.03 Simple disulfide: pdb=" SG CYS B 266 " - pdb=" SG CYS B 275 " distance=2.03 Simple disulfide: pdb=" SG CYS B 337 " - pdb=" SG CYS B 408 " distance=2.03 Simple disulfide: pdb=" SG CYS B 347 " - pdb=" SG CYS B 350 " distance=2.03 Simple disulfide: pdb=" SG CYS B 360 " - pdb=" SG CYS B 382 " distance=2.03 Simple disulfide: pdb=" SG CYS B 373 " - pdb=" SG CYS B 404 " distance=2.03 Simple disulfide: pdb=" SG CYS C 1 " - pdb=" SG CYS C 41 " distance=2.03 Simple disulfide: pdb=" SG CYS C 14 " - pdb=" SG CYS C 23 " distance=2.04 Simple disulfide: pdb=" SG CYS C 66 " - pdb=" SG CYS C 162 " distance=2.02 Simple disulfide: pdb=" SG CYS C 87 " - pdb=" SG CYS C 135 " distance=2.04 Simple disulfide: pdb=" SG CYS C 93 " - pdb=" SG CYS C 142 " distance=2.03 Simple disulfide: pdb=" SG CYS C 98 " - pdb=" SG CYS C 123 " distance=2.03 Simple disulfide: pdb=" SG CYS C 127 " - pdb=" SG CYS C 132 " distance=2.03 Simple disulfide: pdb=" SG CYS C 129 " - pdb=" SG CYS C 138 " distance=2.04 Simple disulfide: pdb=" SG CYS C 246 " - pdb=" SG CYS C 257 " distance=2.03 Simple disulfide: pdb=" SG CYS C 264 " - pdb=" SG CYS C 277 " distance=2.03 Simple disulfide: pdb=" SG CYS C 266 " - pdb=" SG CYS C 275 " distance=2.03 Simple disulfide: pdb=" SG CYS C 337 " - pdb=" SG CYS C 408 " distance=2.03 Simple disulfide: pdb=" SG CYS C 347 " - pdb=" SG CYS C 350 " distance=2.03 Simple disulfide: pdb=" SG CYS C 360 " - pdb=" SG CYS C 382 " distance=2.03 Simple disulfide: pdb=" SG CYS C 373 " - pdb=" SG CYS C 404 " distance=2.03 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 501 " - " ASN A 414 " " NAG B 501 " - " ASN B 414 " " NAG C 501 " - " ASN C 414 " Time building additional restraints: 14.51 Conformation dependent library (CDL) restraints added in 1.7 seconds 2622 Ramachandran restraints generated. 1311 Oldfield, 0 Emsley, 1311 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2514 Finding SS restraints... Secondary structure from input PDB file: 9 helices and 33 sheets defined 4.8% alpha, 37.1% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.82 Creating SS restraints... Processing helix chain 'A' and resid 112 through 117 removed outlier: 4.094A pdb=" N SER A 116 " --> pdb=" O GLU A 113 " (cutoff:3.500A) Processing helix chain 'A' and resid 248 through 254 Processing helix chain 'A' and resid 282 through 289 removed outlier: 3.927A pdb=" N LEU A 289 " --> pdb=" O HIS A 285 " (cutoff:3.500A) Processing helix chain 'B' and resid 112 through 117 removed outlier: 4.094A pdb=" N SER B 116 " --> pdb=" O GLU B 113 " (cutoff:3.500A) Processing helix chain 'B' and resid 248 through 254 Processing helix chain 'B' and resid 282 through 289 removed outlier: 3.928A pdb=" N LEU B 289 " --> pdb=" O HIS B 285 " (cutoff:3.500A) Processing helix chain 'C' and resid 112 through 117 removed outlier: 4.095A pdb=" N SER C 116 " --> pdb=" O GLU C 113 " (cutoff:3.500A) Processing helix chain 'C' and resid 248 through 254 Processing helix chain 'C' and resid 282 through 289 removed outlier: 3.928A pdb=" N LEU C 289 " --> pdb=" O HIS C 285 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 5 through 9 removed outlier: 4.028A pdb=" N ILE A 5 " --> pdb=" O THR A 32 " (cutoff:3.500A) removed outlier: 7.048A pdb=" N HIS A 320 " --> pdb=" O LEU A 203 " (cutoff:3.500A) removed outlier: 5.813A pdb=" N LEU A 203 " --> pdb=" O HIS A 320 " (cutoff:3.500A) removed outlier: 7.609A pdb=" N VAL A 322 " --> pdb=" O PHE A 201 " (cutoff:3.500A) removed outlier: 6.013A pdb=" N PHE A 201 " --> pdb=" O VAL A 322 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 20 through 24 removed outlier: 3.645A pdb=" N CYS A 14 " --> pdb=" O THR C 330 " (cutoff:3.500A) removed outlier: 4.367A pdb=" N VAL C 331 " --> pdb=" O CYS C 41 " (cutoff:3.500A) removed outlier: 6.095A pdb=" N CYS C 41 " --> pdb=" O VAL C 331 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 20 through 24 removed outlier: 3.645A pdb=" N CYS A 14 " --> pdb=" O THR C 330 " (cutoff:3.500A) removed outlier: 4.367A pdb=" N VAL C 331 " --> pdb=" O CYS C 41 " (cutoff:3.500A) removed outlier: 6.095A pdb=" N CYS C 41 " --> pdb=" O VAL C 331 " (cutoff:3.500A) removed outlier: 3.568A pdb=" N VAL C 300 " --> pdb=" O GLN C 293 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'A' and resid 184 through 188 removed outlier: 3.568A pdb=" N VAL A 300 " --> pdb=" O GLN A 293 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'A' and resid 184 through 188 removed outlier: 6.095A pdb=" N CYS A 41 " --> pdb=" O VAL A 331 " (cutoff:3.500A) removed outlier: 4.367A pdb=" N VAL A 331 " --> pdb=" O CYS A 41 " (cutoff:3.500A) removed outlier: 8.669A pdb=" N CYS B 14 " --> pdb=" O ILE A 328 " (cutoff:3.500A) removed outlier: 6.464A pdb=" N THR A 330 " --> pdb=" O CYS B 14 " (cutoff:3.500A) removed outlier: 7.537A pdb=" N ILE B 16 " --> pdb=" O THR A 330 " (cutoff:3.500A) removed outlier: 6.607A pdb=" N LYS A 332 " --> pdb=" O ILE B 16 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'A' and resid 69 through 74 removed outlier: 5.793A pdb=" N ILE A 70 " --> pdb=" O TYR A 161 " (cutoff:3.500A) removed outlier: 7.465A pdb=" N TYR A 161 " --> pdb=" O ILE A 70 " (cutoff:3.500A) removed outlier: 3.694A pdb=" N PHE A 72 " --> pdb=" O ILE A 159 " (cutoff:3.500A) removed outlier: 3.721A pdb=" N ALA A 217 " --> pdb=" O VAL A 228 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'A' and resid 77 through 85 Processing sheet with id=AA8, first strand: chain 'A' and resid 244 through 245 Processing sheet with id=AA9, first strand: chain 'A' and resid 265 through 267 removed outlier: 3.575A pdb=" N THR A 276 " --> pdb=" O THR A 265 " (cutoff:3.500A) removed outlier: 3.696A pdb=" N TYR A 433 " --> pdb=" O CYS A 277 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'A' and resid 338 through 340 Processing sheet with id=AB2, first strand: chain 'A' and resid 369 through 372 Processing sheet with id=AB3, first strand: chain 'A' and resid 399 through 402 Processing sheet with id=AB4, first strand: chain 'A' and resid 405 through 407 Processing sheet with id=AB5, first strand: chain 'B' and resid 5 through 9 removed outlier: 4.028A pdb=" N ILE B 5 " --> pdb=" O THR B 32 " (cutoff:3.500A) removed outlier: 7.048A pdb=" N HIS B 320 " --> pdb=" O LEU B 203 " (cutoff:3.500A) removed outlier: 5.812A pdb=" N LEU B 203 " --> pdb=" O HIS B 320 " (cutoff:3.500A) removed outlier: 7.609A pdb=" N VAL B 322 " --> pdb=" O PHE B 201 " (cutoff:3.500A) removed outlier: 6.013A pdb=" N PHE B 201 " --> pdb=" O VAL B 322 " (cutoff:3.500A) Processing sheet with id=AB6, first strand: chain 'B' and resid 184 through 188 removed outlier: 3.568A pdb=" N VAL B 300 " --> pdb=" O GLN B 293 " (cutoff:3.500A) Processing sheet with id=AB7, first strand: chain 'B' and resid 184 through 188 removed outlier: 6.095A pdb=" N CYS B 41 " --> pdb=" O VAL B 331 " (cutoff:3.500A) removed outlier: 4.368A pdb=" N VAL B 331 " --> pdb=" O CYS B 41 " (cutoff:3.500A) removed outlier: 8.708A pdb=" N CYS C 14 " --> pdb=" O ILE B 328 " (cutoff:3.500A) removed outlier: 6.506A pdb=" N THR B 330 " --> pdb=" O CYS C 14 " (cutoff:3.500A) removed outlier: 7.582A pdb=" N ILE C 16 " --> pdb=" O THR B 330 " (cutoff:3.500A) removed outlier: 6.646A pdb=" N LYS B 332 " --> pdb=" O ILE C 16 " (cutoff:3.500A) Processing sheet with id=AB8, first strand: chain 'B' and resid 69 through 74 removed outlier: 5.793A pdb=" N ILE B 70 " --> pdb=" O TYR B 161 " (cutoff:3.500A) removed outlier: 7.464A pdb=" N TYR B 161 " --> pdb=" O ILE B 70 " (cutoff:3.500A) removed outlier: 3.694A pdb=" N PHE B 72 " --> pdb=" O ILE B 159 " (cutoff:3.500A) removed outlier: 3.721A pdb=" N ALA B 217 " --> pdb=" O VAL B 228 " (cutoff:3.500A) Processing sheet with id=AB9, first strand: chain 'B' and resid 77 through 85 Processing sheet with id=AC1, first strand: chain 'B' and resid 244 through 245 Processing sheet with id=AC2, first strand: chain 'B' and resid 265 through 267 removed outlier: 3.575A pdb=" N THR B 276 " --> pdb=" O THR B 265 " (cutoff:3.500A) removed outlier: 3.696A pdb=" N TYR B 433 " --> pdb=" O CYS B 277 " (cutoff:3.500A) Processing sheet with id=AC3, first strand: chain 'B' and resid 338 through 340 Processing sheet with id=AC4, first strand: chain 'B' and resid 369 through 372 Processing sheet with id=AC5, first strand: chain 'B' and resid 399 through 402 Processing sheet with id=AC6, first strand: chain 'B' and resid 405 through 407 Processing sheet with id=AC7, first strand: chain 'C' and resid 5 through 9 removed outlier: 4.027A pdb=" N ILE C 5 " --> pdb=" O THR C 32 " (cutoff:3.500A) removed outlier: 7.049A pdb=" N HIS C 320 " --> pdb=" O LEU C 203 " (cutoff:3.500A) removed outlier: 5.812A pdb=" N LEU C 203 " --> pdb=" O HIS C 320 " (cutoff:3.500A) removed outlier: 7.609A pdb=" N VAL C 322 " --> pdb=" O PHE C 201 " (cutoff:3.500A) removed outlier: 6.013A pdb=" N PHE C 201 " --> pdb=" O VAL C 322 " (cutoff:3.500A) Processing sheet with id=AC8, first strand: chain 'C' and resid 69 through 74 removed outlier: 5.793A pdb=" N ILE C 70 " --> pdb=" O TYR C 161 " (cutoff:3.500A) removed outlier: 7.465A pdb=" N TYR C 161 " --> pdb=" O ILE C 70 " (cutoff:3.500A) removed outlier: 3.694A pdb=" N PHE C 72 " --> pdb=" O ILE C 159 " (cutoff:3.500A) removed outlier: 3.721A pdb=" N ALA C 217 " --> pdb=" O VAL C 228 " (cutoff:3.500A) Processing sheet with id=AC9, first strand: chain 'C' and resid 77 through 85 Processing sheet with id=AD1, first strand: chain 'C' and resid 244 through 245 Processing sheet with id=AD2, first strand: chain 'C' and resid 265 through 267 removed outlier: 3.574A pdb=" N THR C 276 " --> pdb=" O THR C 265 " (cutoff:3.500A) removed outlier: 3.696A pdb=" N TYR C 433 " --> pdb=" O CYS C 277 " (cutoff:3.500A) Processing sheet with id=AD3, first strand: chain 'C' and resid 338 through 340 Processing sheet with id=AD4, first strand: chain 'C' and resid 369 through 372 Processing sheet with id=AD5, first strand: chain 'C' and resid 399 through 402 Processing sheet with id=AD6, first strand: chain 'C' and resid 405 through 407 366 hydrogen bonds defined for protein. 861 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: 7.58 Time building geometry restraints manager: 15.28 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 0.83 - 1.03: 9618 1.03 - 1.23: 170 1.23 - 1.43: 4112 1.43 - 1.63: 5960 1.63 - 1.82: 120 Bond restraints: 19980 Sorted by residual: bond pdb=" CA ASN A 67 " pdb=" CB ASN A 67 " ideal model delta sigma weight residual 1.527 1.496 0.031 1.75e-02 3.27e+03 3.15e+00 bond pdb=" CA ASN B 67 " pdb=" CB ASN B 67 " ideal model delta sigma weight residual 1.527 1.496 0.031 1.75e-02 3.27e+03 3.13e+00 bond pdb=" CA ASN C 67 " pdb=" CB ASN C 67 " ideal model delta sigma weight residual 1.527 1.496 0.030 1.75e-02 3.27e+03 3.00e+00 bond pdb=" CB ILE B 207 " pdb=" CG2 ILE B 207 " ideal model delta sigma weight residual 1.521 1.475 0.046 3.30e-02 9.18e+02 1.90e+00 bond pdb=" CB ILE C 207 " pdb=" CG2 ILE C 207 " ideal model delta sigma weight residual 1.521 1.476 0.045 3.30e-02 9.18e+02 1.88e+00 ... (remaining 19975 not shown) Histogram of bond angle deviations from ideal: 99.83 - 106.70: 399 106.70 - 113.56: 23868 113.56 - 120.42: 5720 120.42 - 127.28: 5899 127.28 - 134.15: 102 Bond angle restraints: 35988 Sorted by residual: angle pdb=" N LEU A 270 " pdb=" CA LEU A 270 " pdb=" CB LEU A 270 " ideal model delta sigma weight residual 114.17 110.47 3.70 1.14e+00 7.69e-01 1.05e+01 angle pdb=" N LEU C 270 " pdb=" CA LEU C 270 " pdb=" CB LEU C 270 " ideal model delta sigma weight residual 114.17 110.50 3.67 1.14e+00 7.69e-01 1.04e+01 angle pdb=" N LEU B 270 " pdb=" CA LEU B 270 " pdb=" CB LEU B 270 " ideal model delta sigma weight residual 114.17 110.50 3.67 1.14e+00 7.69e-01 1.04e+01 angle pdb=" C TRP B 196 " pdb=" N ASN B 197 " pdb=" CA ASN B 197 " ideal model delta sigma weight residual 122.17 117.39 4.78 1.54e+00 4.22e-01 9.64e+00 angle pdb=" C TRP A 196 " pdb=" N ASN A 197 " pdb=" CA ASN A 197 " ideal model delta sigma weight residual 122.17 117.39 4.78 1.54e+00 4.22e-01 9.62e+00 ... (remaining 35983 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 14.74: 7370 14.74 - 29.48: 633 29.48 - 44.22: 136 44.22 - 58.96: 69 58.96 - 73.70: 15 Dihedral angle restraints: 8223 sinusoidal: 3810 harmonic: 4413 Sorted by residual: dihedral pdb=" CA GLU C 71 " pdb=" C GLU C 71 " pdb=" N PHE C 72 " pdb=" CA PHE C 72 " ideal model delta harmonic sigma weight residual -180.00 -152.07 -27.93 0 5.00e+00 4.00e-02 3.12e+01 dihedral pdb=" CA GLU B 71 " pdb=" C GLU B 71 " pdb=" N PHE B 72 " pdb=" CA PHE B 72 " ideal model delta harmonic sigma weight residual -180.00 -152.11 -27.89 0 5.00e+00 4.00e-02 3.11e+01 dihedral pdb=" CA GLU A 71 " pdb=" C GLU A 71 " pdb=" N PHE A 72 " pdb=" CA PHE A 72 " ideal model delta harmonic sigma weight residual 180.00 -152.12 -27.88 0 5.00e+00 4.00e-02 3.11e+01 ... (remaining 8220 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.030: 889 0.030 - 0.059: 506 0.059 - 0.089: 187 0.089 - 0.119: 56 0.119 - 0.149: 21 Chirality restraints: 1659 Sorted by residual: chirality pdb=" CA PRO A 190 " pdb=" N PRO A 190 " pdb=" C PRO A 190 " pdb=" CB PRO A 190 " both_signs ideal model delta sigma weight residual False 2.72 2.57 0.15 2.00e-01 2.50e+01 5.52e-01 chirality pdb=" CA PRO C 190 " pdb=" N PRO C 190 " pdb=" C PRO C 190 " pdb=" CB PRO C 190 " both_signs ideal model delta sigma weight residual False 2.72 2.57 0.15 2.00e-01 2.50e+01 5.50e-01 chirality pdb=" CA PRO B 190 " pdb=" N PRO B 190 " pdb=" C PRO B 190 " pdb=" CB PRO B 190 " both_signs ideal model delta sigma weight residual False 2.72 2.57 0.15 2.00e-01 2.50e+01 5.47e-01 ... (remaining 1656 not shown) Planarity restraints: 3081 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C LEU A 211 " -0.048 5.00e-02 4.00e+02 7.25e-02 8.42e+00 pdb=" N PRO A 212 " 0.125 5.00e-02 4.00e+02 pdb=" CA PRO A 212 " -0.038 5.00e-02 4.00e+02 pdb=" CD PRO A 212 " -0.039 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C LEU B 211 " 0.048 5.00e-02 4.00e+02 7.24e-02 8.39e+00 pdb=" N PRO B 212 " -0.125 5.00e-02 4.00e+02 pdb=" CA PRO B 212 " 0.038 5.00e-02 4.00e+02 pdb=" CD PRO B 212 " 0.039 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C LEU C 211 " 0.048 5.00e-02 4.00e+02 7.22e-02 8.35e+00 pdb=" N PRO C 212 " -0.125 5.00e-02 4.00e+02 pdb=" CA PRO C 212 " 0.038 5.00e-02 4.00e+02 pdb=" CD PRO C 212 " 0.039 5.00e-02 4.00e+02 ... (remaining 3078 not shown) Histogram of nonbonded interaction distances: 1.56 - 2.17: 1021 2.17 - 2.78: 39337 2.78 - 3.38: 49624 3.38 - 3.99: 67698 3.99 - 4.60: 104566 Nonbonded interactions: 262246 Sorted by model distance: nonbonded pdb=" O LYS A 332 " pdb=" H ILE B 16 " model vdw 1.558 1.850 nonbonded pdb=" O LYS B 332 " pdb=" H ILE C 16 " model vdw 1.573 1.850 nonbonded pdb=" HH TYR C 120 " pdb=" O PHE C 259 " model vdw 1.593 1.850 nonbonded pdb=" HH TYR B 120 " pdb=" O PHE B 259 " model vdw 1.603 1.850 nonbonded pdb=" HH TYR A 120 " pdb=" O PHE A 259 " model vdw 1.610 1.850 ... (remaining 262241 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' selection = chain 'B' selection = chain 'C' } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Set scattering table Set to: electron Number of scattering types: 5 Type Number sf(0) Gaussians S 105 5.16 5 C 6294 2.51 5 N 1743 2.21 5 O 2025 1.98 5 H 9618 0.53 5 sf(0) = scattering factor at diffraction angle 0. Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as 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.940 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.260 Extract box with map and model: 4.480 Check model and map are aligned: 0.290 Process input model: 67.810 Find NCS groups from input model: 0.880 Set up NCS constraints: 0.090 Set refine NCS operators: 0.000 Set scattering table: 0.190 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.860 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 77.800 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7833 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.007 0.055 10362 Z= 0.497 Angle : 0.710 6.857 14073 Z= 0.401 Chirality : 0.043 0.149 1659 Planarity : 0.006 0.073 1821 Dihedral : 11.654 73.698 3747 Min Nonbonded Distance : 2.109 Molprobity Statistics. All-atom Clashscore : 13.28 Ramachandran Plot: Outliers : 0.00 % Allowed : 12.81 % Favored : 87.19 % Rotamer Outliers : 2.96 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.56 % 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: -6.00 (0.16), residues: 1311 helix: -4.78 (0.19), residues: 81 sheet: -3.30 (0.17), residues: 543 loop : -4.57 (0.16), residues: 687 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2622 Ramachandran restraints generated. 1311 Oldfield, 0 Emsley, 1311 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2622 Ramachandran restraints generated. 1311 Oldfield, 0 Emsley, 1311 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Residue ILE 199 is missing expected H atoms. Skipping. Residue SER 202 is missing expected H atoms. Skipping. Residue LEU 203 is missing expected H atoms. Skipping. Residue LYS 272 is missing expected H atoms. Skipping. Residue ILE 199 is missing expected H atoms. Skipping. Residue SER 202 is missing expected H atoms. Skipping. Residue LEU 203 is missing expected H atoms. Skipping. Residue LYS 272 is missing expected H atoms. Skipping. Residue ILE 199 is missing expected H atoms. Skipping. Residue SER 202 is missing expected H atoms. Skipping. Residue LEU 203 is missing expected H atoms. Skipping. Residue LYS 272 is missing expected H atoms. Skipping. Evaluate side-chains 298 residues out of total 1182 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 35 poor density : 263 time to evaluate : 1.624 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 outliers start: 35 outliers final: 11 residues processed: 295 average time/residue: 0.4474 time to fit residues: 188.0370 Evaluate side-chains 143 residues out of total 1182 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 132 time to evaluate : 1.469 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 0 residues processed: 11 average time/residue: 0.1769 time to fit residues: 5.2747 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 129 random chunks: chunk 108 optimal weight: 0.9990 chunk 97 optimal weight: 9.9990 chunk 54 optimal weight: 9.9990 chunk 33 optimal weight: 4.9990 chunk 65 optimal weight: 0.9990 chunk 52 optimal weight: 20.0000 chunk 101 optimal weight: 1.9990 chunk 39 optimal weight: 3.9990 chunk 61 optimal weight: 2.9990 chunk 75 optimal weight: 4.9990 chunk 117 optimal weight: 2.9990 overall best weight: 1.9990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Sorry: Reduce crashed with command 'molprobity.reduce -quiet -build -allalt -'. Dumping stdin to file 'reduce_failure.pdb'. Return code: -15 Dumping stderr: