Starting phenix.real_space_refine on Mon Dec 11 05:50:47 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7v9a_31813/12_2023/7v9a_31813.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7v9a_31813/12_2023/7v9a_31813.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.94 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7v9a_31813/12_2023/7v9a_31813.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7v9a_31813/12_2023/7v9a_31813.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7v9a_31813/12_2023/7v9a_31813.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7v9a_31813/12_2023/7v9a_31813.pdb" } resolution = 3.94 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.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 Set scattering table Set to: electron Number of scattering types: 5 Type Number sf(0) Gaussians P 146 5.49 5 S 86 5.16 5 C 10128 2.51 5 N 2999 2.21 5 O 3517 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "B PHE 150": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 196": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 211": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 229": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B GLU 263": "OE1" <-> "OE2" Residue "B TYR 271": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ASP 288": "OD1" <-> "OD2" Residue "B PHE 293": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 303": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 318": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 340": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 345": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ASP 356": "OD1" <-> "OD2" Residue "B ASP 378": "OD1" <-> "OD2" Residue "B PHE 382": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 383": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 416": "NH1" <-> "NH2" Residue "B PHE 419": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ASP 420": "OD1" <-> "OD2" Residue "B ASP 448": "OD1" <-> "OD2" Residue "B PHE 457": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C ARG 65": "NH1" <-> "NH2" Residue "C TYR 86": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 91": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C ARG 138": "NH1" <-> "NH2" Residue "C GLU 152": "OE1" <-> "OE2" Residue "C TYR 153": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C ARG 195": "NH1" <-> "NH2" Residue "C ARG 197": "NH1" <-> "NH2" Residue "C GLU 201": "OE1" <-> "OE2" Residue "C GLU 210": "OE1" <-> "OE2" Residue "C ARG 212": "NH1" <-> "NH2" Residue "C TYR 225": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C ARG 248": "NH1" <-> "NH2" Residue "C ASP 258": "OD1" <-> "OD2" Residue "C TYR 287": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C ARG 322": "NH1" <-> "NH2" Residue "C GLU 324": "OE1" <-> "OE2" Residue "C ASP 325": "OD1" <-> "OD2" Residue "C GLU 341": "OE1" <-> "OE2" Residue "C ARG 373": "NH1" <-> "NH2" Residue "C TYR 376": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C ASP 399": "OD1" <-> "OD2" Residue "C TYR 419": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D ARG 46": "NH1" <-> "NH2" Residue "D PHE 49": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D PHE 58": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D ASP 114": "OD1" <-> "OD2" Residue "D PHE 117": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D ASP 122": "OD1" <-> "OD2" Residue "D PHE 123": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D ARG 156": "NH1" <-> "NH2" Residue "E PHE 70": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "E ASP 84": "OD1" <-> "OD2" Residue "E GLU 99": "OE1" <-> "OE2" Residue "E GLU 133": "OE1" <-> "OE2" Residue "E GLU 134": "OE1" <-> "OE2" Residue "E TYR 135": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "E GLU 144": "OE1" <-> "OE2" Residue "F TYR 5": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "F TYR 6": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "F TYR 41": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "F ARG 51": "NH1" <-> "NH2" Residue "F PHE 52": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "G ARG 65": "NH1" <-> "NH2" Residue "G TYR 86": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "G PHE 91": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "G ARG 114": "NH1" <-> "NH2" Residue "G ASP 125": "OD1" <-> "OD2" Residue "G GLU 137": "OE1" <-> "OE2" Residue "G ARG 138": "NH1" <-> "NH2" Residue "G ARG 141": "NH1" <-> "NH2" Residue "G GLU 152": "OE1" <-> "OE2" Residue "G TYR 153": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "G GLU 164": "OE1" <-> "OE2" Residue "G GLU 174": "OE1" <-> "OE2" Residue "G ARG 183": "NH1" <-> "NH2" Residue "G ARG 195": "NH1" <-> "NH2" Residue "G ARG 197": "NH1" <-> "NH2" Residue "G GLU 201": "OE1" <-> "OE2" Residue "G TYR 225": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "G ARG 248": "NH1" <-> "NH2" Residue "G ASP 258": "OD1" <-> "OD2" Residue "G GLU 279": "OE1" <-> "OE2" Residue "G ARG 298": "NH1" <-> "NH2" Residue "G TYR 311": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "G ARG 322": "NH1" <-> "NH2" Residue "G GLU 324": "OE1" <-> "OE2" Residue "G ASP 325": "OD1" <-> "OD2" Residue "G GLU 341": "OE1" <-> "OE2" Residue "G ARG 373": "NH1" <-> "NH2" Residue "G ASP 374": "OD1" <-> "OD2" Residue "G TYR 376": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "H GLU 81": "OE1" <-> "OE2" Residue "H GLU 92": "OE1" <-> "OE2" Residue "H TYR 97": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "H TYR 103": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "H GLU 105": "OE1" <-> "OE2" Residue "H GLU 131": "OE1" <-> "OE2" Residue "H PHE 138": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "I GLU 21": "OE1" <-> "OE2" Residue "I GLU 26": "OE1" <-> "OE2" Residue "I GLU 66": "OE1" <-> "OE2" Residue "I PHE 70": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "I TYR 105": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "I TYR 135": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "J PHE 2": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "J TYR 5": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "J TYR 6": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "J ASP 12": "OD1" <-> "OD2" Residue "J TYR 41": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "J PHE 52": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.05s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/chem_data/mon_lib" Total number of atoms: 16876 Number of models: 1 Model: "" Number of chains: 10 Chain: "B" Number of atoms: 2938 Number of conformers: 1 Conformer: "" Number of residues, atoms: 377, 2938 Classifications: {'peptide': 377} Link IDs: {'PTRANS': 25, 'TRANS': 351} Chain: "C" Number of atoms: 3118 Number of conformers: 1 Conformer: "" Number of residues, atoms: 393, 3118 Classifications: {'peptide': 393} Link IDs: {'PTRANS': 16, 'TRANS': 376} Chain: "D" Number of atoms: 940 Number of conformers: 1 Conformer: "" Number of residues, atoms: 117, 940 Classifications: {'peptide': 117} Link IDs: {'PTRANS': 9, 'TRANS': 107} Chain: "E" Number of atoms: 1031 Number of conformers: 1 Conformer: "" Number of residues, atoms: 129, 1031 Classifications: {'peptide': 129} Link IDs: {'PCIS': 1, 'PTRANS': 8, 'TRANS': 119} Chain: "F" Number of atoms: 526 Number of conformers: 1 Conformer: "" Number of residues, atoms: 62, 526 Classifications: {'peptide': 62} Link IDs: {'PTRANS': 5, 'TRANS': 56} Chain: "G" Number of atoms: 2762 Number of conformers: 1 Conformer: "" Number of residues, atoms: 349, 2762 Classifications: {'peptide': 349} Link IDs: {'PTRANS': 13, 'TRANS': 335} Chain: "H" Number of atoms: 811 Number of conformers: 1 Conformer: "" Number of residues, atoms: 98, 811 Classifications: {'peptide': 98} Link IDs: {'PTRANS': 9, 'TRANS': 88} Chain: "I" Number of atoms: 1065 Number of conformers: 1 Conformer: "" Number of residues, atoms: 134, 1065 Classifications: {'peptide': 134} Link IDs: {'PCIS': 1, 'PTRANS': 9, 'TRANS': 123} Chain: "J" Number of atoms: 541 Number of conformers: 1 Conformer: "" Number of residues, atoms: 64, 541 Classifications: {'peptide': 64} Link IDs: {'PTRANS': 5, 'TRANS': 58} Chain: "R" Number of atoms: 3144 Number of conformers: 1 Conformer: "" Number of residues, atoms: 146, 3144 Classifications: {'RNA': 146} Modifications used: {'rna2p_pur': 12, 'rna2p_pyr': 7, 'rna3p_pur': 70, 'rna3p_pyr': 57} Link IDs: {'rna2p': 19, 'rna3p': 126} Chain breaks: 2 Time building chain proxies: 9.28, per 1000 atoms: 0.55 Number of scatterers: 16876 At special positions: 0 Unit cell: (183.7, 96.8, 132, 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 86 16.00 P 146 15.00 O 3517 8.00 N 2999 7.00 C 10128 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 6.37 Conformation dependent library (CDL) restraints added in 2.7 seconds 3410 Ramachandran restraints generated. 1705 Oldfield, 0 Emsley, 1705 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 3200 Finding SS restraints... Secondary structure from input PDB file: 45 helices and 21 sheets defined 24.7% alpha, 23.0% beta 16 base pairs and 57 stacking pairs defined. Time for finding SS restraints: 5.50 Creating SS restraints... Processing helix chain 'B' and resid 200 through 206 removed outlier: 4.259A pdb=" N GLY B 206 " --> pdb=" O LEU B 202 " (cutoff:3.500A) Processing helix chain 'B' and resid 490 through 498 Processing helix chain 'C' and resid 53 through 57 removed outlier: 3.700A pdb=" N LYS C 57 " --> pdb=" O LEU C 54 " (cutoff:3.500A) Processing helix chain 'C' and resid 82 through 89 Processing helix chain 'C' and resid 101 through 114 Processing helix chain 'C' and resid 139 through 141 No H-bonds generated for 'chain 'C' and resid 139 through 141' Processing helix chain 'C' and resid 142 through 148 removed outlier: 3.785A pdb=" N SER C 148 " --> pdb=" O LYS C 144 " (cutoff:3.500A) Processing helix chain 'C' and resid 165 through 176 removed outlier: 3.910A pdb=" N LEU C 169 " --> pdb=" O GLY C 165 " (cutoff:3.500A) removed outlier: 3.522A pdb=" N LEU C 176 " --> pdb=" O ALA C 172 " (cutoff:3.500A) Processing helix chain 'C' and resid 225 through 238 removed outlier: 3.535A pdb=" N LEU C 229 " --> pdb=" O TYR C 225 " (cutoff:3.500A) Processing helix chain 'C' and resid 262 through 276 removed outlier: 3.525A pdb=" N VAL C 266 " --> pdb=" O THR C 262 " (cutoff:3.500A) removed outlier: 3.592A pdb=" N GLN C 270 " --> pdb=" O VAL C 266 " (cutoff:3.500A) removed outlier: 3.635A pdb=" N TRP C 271 " --> pdb=" O LEU C 267 " (cutoff:3.500A) Processing helix chain 'C' and resid 279 through 286 removed outlier: 3.746A pdb=" N VAL C 285 " --> pdb=" O TYR C 281 " (cutoff:3.500A) Processing helix chain 'C' and resid 288 through 295 removed outlier: 3.632A pdb=" N LEU C 293 " --> pdb=" O GLU C 290 " (cutoff:3.500A) removed outlier: 3.511A pdb=" N THR C 294 " --> pdb=" O LYS C 291 " (cutoff:3.500A) removed outlier: 4.575A pdb=" N SER C 295 " --> pdb=" O LEU C 292 " (cutoff:3.500A) Processing helix chain 'C' and resid 305 through 312 Processing helix chain 'C' and resid 351 through 356 Processing helix chain 'C' and resid 385 through 396 Processing helix chain 'C' and resid 409 through 414 removed outlier: 3.538A pdb=" N LYS C 413 " --> pdb=" O PRO C 409 " (cutoff:3.500A) Processing helix chain 'E' and resid 23 through 30 Processing helix chain 'E' and resid 42 through 58 removed outlier: 4.232A pdb=" N LYS E 56 " --> pdb=" O LYS E 52 " (cutoff:3.500A) removed outlier: 3.554A pdb=" N GLN E 57 " --> pdb=" O LYS E 53 " (cutoff:3.500A) removed outlier: 3.683A pdb=" N LYS E 58 " --> pdb=" O ALA E 54 " (cutoff:3.500A) Processing helix chain 'E' and resid 63 through 74 removed outlier: 3.736A pdb=" N LYS E 73 " --> pdb=" O LYS E 69 " (cutoff:3.500A) Processing helix chain 'E' and resid 93 through 101 Processing helix chain 'E' and resid 111 through 120 Processing helix chain 'E' and resid 134 through 137 removed outlier: 4.139A pdb=" N GLU E 137 " --> pdb=" O GLU E 134 " (cutoff:3.500A) No H-bonds generated for 'chain 'E' and resid 134 through 137' Processing helix chain 'E' and resid 138 through 148 removed outlier: 3.624A pdb=" N CYS E 142 " --> pdb=" O ALA E 138 " (cutoff:3.500A) Processing helix chain 'F' and resid 41 through 53 removed outlier: 3.713A pdb=" N ARG F 45 " --> pdb=" O TYR F 41 " (cutoff:3.500A) removed outlier: 3.604A pdb=" N PHE F 52 " --> pdb=" O ILE F 48 " (cutoff:3.500A) Processing helix chain 'G' and resid 53 through 57 removed outlier: 4.096A pdb=" N LYS G 57 " --> pdb=" O LEU G 54 " (cutoff:3.500A) Processing helix chain 'G' and resid 82 through 89 Processing helix chain 'G' and resid 101 through 114 removed outlier: 3.660A pdb=" N ALA G 107 " --> pdb=" O HIS G 103 " (cutoff:3.500A) Processing helix chain 'G' and resid 139 through 141 No H-bonds generated for 'chain 'G' and resid 139 through 141' Processing helix chain 'G' and resid 142 through 148 removed outlier: 3.735A pdb=" N SER G 148 " --> pdb=" O LYS G 144 " (cutoff:3.500A) Processing helix chain 'G' and resid 165 through 176 removed outlier: 4.085A pdb=" N LEU G 169 " --> pdb=" O GLY G 165 " (cutoff:3.500A) removed outlier: 3.505A pdb=" N LEU G 176 " --> pdb=" O ALA G 172 " (cutoff:3.500A) Processing helix chain 'G' and resid 226 through 237 Processing helix chain 'G' and resid 263 through 276 Processing helix chain 'G' and resid 279 through 284 Processing helix chain 'G' and resid 288 through 292 removed outlier: 3.702A pdb=" N LYS G 291 " --> pdb=" O PRO G 288 " (cutoff:3.500A) Processing helix chain 'G' and resid 303 through 312 removed outlier: 4.360A pdb=" N ASN G 307 " --> pdb=" O ASP G 303 " (cutoff:3.500A) removed outlier: 4.092A pdb=" N ALA G 308 " --> pdb=" O SER G 304 " (cutoff:3.500A) Processing helix chain 'G' and resid 351 through 356 Processing helix chain 'G' and resid 385 through 395 Processing helix chain 'H' and resid 148 through 151 removed outlier: 3.809A pdb=" N LEU H 151 " --> pdb=" O PRO H 148 " (cutoff:3.500A) No H-bonds generated for 'chain 'H' and resid 148 through 151' Processing helix chain 'I' and resid 20 through 30 Processing helix chain 'I' and resid 40 through 56 removed outlier: 4.128A pdb=" N LYS I 56 " --> pdb=" O LYS I 52 " (cutoff:3.500A) Processing helix chain 'I' and resid 63 through 74 removed outlier: 3.788A pdb=" N LYS I 73 " --> pdb=" O LYS I 69 " (cutoff:3.500A) Processing helix chain 'I' and resid 93 through 101 Processing helix chain 'I' and resid 111 through 119 removed outlier: 3.502A pdb=" N ALA I 118 " --> pdb=" O LEU I 114 " (cutoff:3.500A) removed outlier: 3.672A pdb=" N GLY I 119 " --> pdb=" O GLY I 115 " (cutoff:3.500A) Processing helix chain 'I' and resid 136 through 147 removed outlier: 3.999A pdb=" N GLU I 141 " --> pdb=" O GLU I 137 " (cutoff:3.500A) Processing helix chain 'J' and resid 42 through 53 removed outlier: 3.911A pdb=" N ILE J 46 " --> pdb=" O SER J 42 " (cutoff:3.500A) removed outlier: 3.669A pdb=" N THR J 47 " --> pdb=" O ARG J 43 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'B' and resid 155 through 160 removed outlier: 3.513A pdb=" N ARG B 155 " --> pdb=" O TRP B 518 " (cutoff:3.500A) removed outlier: 7.232A pdb=" N LEU B 516 " --> pdb=" O LEU B 157 " (cutoff:3.500A) removed outlier: 4.434A pdb=" N GLY B 159 " --> pdb=" O LEU B 514 " (cutoff:3.500A) removed outlier: 6.736A pdb=" N LEU B 514 " --> pdb=" O GLY B 159 " (cutoff:3.500A) removed outlier: 3.739A pdb=" N GLY B 466 " --> pdb=" O ALA B 479 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'B' and resid 173 through 177 removed outlier: 3.582A pdb=" N ASN B 187 " --> pdb=" O LYS B 173 " (cutoff:3.500A) removed outlier: 6.299A pdb=" N LEU B 185 " --> pdb=" O CYS B 175 " (cutoff:3.500A) removed outlier: 4.846A pdb=" N TRP B 177 " --> pdb=" O CYS B 183 " (cutoff:3.500A) removed outlier: 6.988A pdb=" N CYS B 183 " --> pdb=" O TRP B 177 " (cutoff:3.500A) removed outlier: 3.994A pdb=" N LEU B 218 " --> pdb=" O ILE B 195 " (cutoff:3.500A) removed outlier: 7.152A pdb=" N ASN B 197 " --> pdb=" O PRO B 216 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'B' and resid 226 through 231 removed outlier: 7.140A pdb=" N SER B 248 " --> pdb=" O TYR B 227 " (cutoff:3.500A) removed outlier: 4.838A pdb=" N TYR B 229 " --> pdb=" O ALA B 246 " (cutoff:3.500A) removed outlier: 7.163A pdb=" N ALA B 246 " --> pdb=" O TYR B 229 " (cutoff:3.500A) removed outlier: 5.946A pdb=" N TRP B 231 " --> pdb=" O TYR B 244 " (cutoff:3.500A) removed outlier: 8.041A pdb=" N TYR B 244 " --> pdb=" O TRP B 231 " (cutoff:3.500A) removed outlier: 6.582A pdb=" N ILE B 254 " --> pdb=" O SER B 267 " (cutoff:3.500A) removed outlier: 4.561A pdb=" N SER B 267 " --> pdb=" O ILE B 254 " (cutoff:3.500A) removed outlier: 7.198A pdb=" N ILE B 256 " --> pdb=" O ARG B 265 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'B' and resid 283 through 285 removed outlier: 3.646A pdb=" N CYS B 284 " --> pdb=" O PHE B 293 " (cutoff:3.500A) removed outlier: 3.518A pdb=" N PHE B 293 " --> pdb=" O CYS B 284 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'B' and resid 328 through 333 removed outlier: 3.521A pdb=" N LEU B 351 " --> pdb=" O ALA B 361 " (cutoff:3.500A) removed outlier: 4.128A pdb=" N ALA B 361 " --> pdb=" O LEU B 351 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'B' and resid 370 through 373 removed outlier: 6.711A pdb=" N GLY B 385 " --> pdb=" O THR B 371 " (cutoff:3.500A) removed outlier: 5.004A pdb=" N LEU B 373 " --> pdb=" O PHE B 383 " (cutoff:3.500A) removed outlier: 7.134A pdb=" N PHE B 383 " --> pdb=" O LEU B 373 " (cutoff:3.500A) removed outlier: 3.694A pdb=" N SER B 384 " --> pdb=" O LEU B 393 " (cutoff:3.500A) removed outlier: 4.335A pdb=" N LEU B 393 " --> pdb=" O SER B 384 " (cutoff:3.500A) removed outlier: 5.666A pdb=" N ALA B 386 " --> pdb=" O GLU B 391 " (cutoff:3.500A) removed outlier: 7.389A pdb=" N GLU B 391 " --> pdb=" O ALA B 386 " (cutoff:3.500A) removed outlier: 4.128A pdb=" N TRP B 405 " --> pdb=" O CYS B 394 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 419 through 421 removed outlier: 3.565A pdb=" N VAL B 429 " --> pdb=" O ASP B 420 " (cutoff:3.500A) removed outlier: 4.223A pdb=" N LEU B 455 " --> pdb=" O VAL B 439 " (cutoff:3.500A) removed outlier: 5.975A pdb=" N ASP B 441 " --> pdb=" O PRO B 453 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'C' and resid 63 through 64 removed outlier: 3.581A pdb=" N TYR C 323 " --> pdb=" O ASN C 63 " (cutoff:3.500A) removed outlier: 6.750A pdb=" N VAL C 300 " --> pdb=" O LEU C 321 " (cutoff:3.500A) removed outlier: 5.309A pdb=" N TYR C 323 " --> pdb=" O ARG C 298 " (cutoff:3.500A) removed outlier: 7.587A pdb=" N ARG C 298 " --> pdb=" O TYR C 323 " (cutoff:3.500A) removed outlier: 6.037A pdb=" N LYS C 297 " --> pdb=" O VAL C 334 " (cutoff:3.500A) removed outlier: 6.818A pdb=" N ILE C 336 " --> pdb=" O LYS C 297 " (cutoff:3.500A) removed outlier: 5.494A pdb=" N LEU C 299 " --> pdb=" O ILE C 336 " (cutoff:3.500A) removed outlier: 3.594A pdb=" N GLN C 331 " --> pdb=" O ALA C 348 " (cutoff:3.500A) removed outlier: 3.521A pdb=" N ARG C 368 " --> pdb=" O MET C 345 " (cutoff:3.500A) removed outlier: 6.523A pdb=" N ILE C 347 " --> pdb=" O ILE C 366 " (cutoff:3.500A) removed outlier: 6.487A pdb=" N ILE C 366 " --> pdb=" O ILE C 347 " (cutoff:3.500A) removed outlier: 4.262A pdb=" N ALA C 364 " --> pdb=" O ILE C 315 " (cutoff:3.500A) removed outlier: 3.584A pdb=" N ILE C 315 " --> pdb=" O ALA C 364 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'C' and resid 91 through 95 removed outlier: 3.587A pdb=" N CYS C 131 " --> pdb=" O ARG C 248 " (cutoff:3.500A) removed outlier: 6.955A pdb=" N ARG C 250 " --> pdb=" O CYS C 131 " (cutoff:3.500A) removed outlier: 9.054A pdb=" N ILE C 133 " --> pdb=" O ARG C 250 " (cutoff:3.500A) removed outlier: 6.534A pdb=" N GLN C 244 " --> pdb=" O ARG C 158 " (cutoff:3.500A) removed outlier: 7.692A pdb=" N ARG C 158 " --> pdb=" O GLN C 244 " (cutoff:3.500A) removed outlier: 6.103A pdb=" N LEU C 246 " --> pdb=" O ILE C 156 " (cutoff:3.500A) removed outlier: 7.092A pdb=" N ILE C 156 " --> pdb=" O LEU C 246 " (cutoff:3.500A) removed outlier: 4.910A pdb=" N ARG C 248 " --> pdb=" O VAL C 154 " (cutoff:3.500A) removed outlier: 6.383A pdb=" N VAL C 154 " --> pdb=" O ARG C 248 " (cutoff:3.500A) removed outlier: 3.967A pdb=" N ARG C 250 " --> pdb=" O GLU C 152 " (cutoff:3.500A) removed outlier: 3.579A pdb=" N GLY C 155 " --> pdb=" O PHE C 216 " (cutoff:3.500A) removed outlier: 3.577A pdb=" N GLY C 214 " --> pdb=" O VAL C 157 " (cutoff:3.500A) removed outlier: 6.887A pdb=" N LEU C 213 " --> pdb=" O TYR C 207 " (cutoff:3.500A) removed outlier: 4.837A pdb=" N TYR C 207 " --> pdb=" O LEU C 213 " (cutoff:3.500A) removed outlier: 7.207A pdb=" N ILE C 215 " --> pdb=" O ILE C 205 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'C' and resid 119 through 120 removed outlier: 3.587A pdb=" N CYS C 131 " --> pdb=" O ARG C 248 " (cutoff:3.500A) removed outlier: 6.955A pdb=" N ARG C 250 " --> pdb=" O CYS C 131 " (cutoff:3.500A) removed outlier: 9.054A pdb=" N ILE C 133 " --> pdb=" O ARG C 250 " (cutoff:3.500A) removed outlier: 3.815A pdb=" N MET C 254 " --> pdb=" O SER C 251 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'D' and resid 71 through 80 removed outlier: 5.767A pdb=" N PHE D 76 " --> pdb=" O LYS D 87 " (cutoff:3.500A) removed outlier: 6.489A pdb=" N LYS D 87 " --> pdb=" O PHE D 76 " (cutoff:3.500A) removed outlier: 3.972A pdb=" N HIS D 78 " --> pdb=" O VAL D 85 " (cutoff:3.500A) removed outlier: 3.578A pdb=" N VAL D 85 " --> pdb=" O HIS D 78 " (cutoff:3.500A) removed outlier: 7.141A pdb=" N TYR D 124 " --> pdb=" O ILE D 116 " (cutoff:3.500A) removed outlier: 5.025A pdb=" N ILE D 116 " --> pdb=" O TYR D 124 " (cutoff:3.500A) removed outlier: 6.568A pdb=" N SER D 126 " --> pdb=" O ASP D 114 " (cutoff:3.500A) removed outlier: 3.851A pdb=" N GLY D 111 " --> pdb=" O VAL D 102 " (cutoff:3.500A) removed outlier: 6.717A pdb=" N PHE D 144 " --> pdb=" O LEU D 73 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'D' and resid 94 through 96 removed outlier: 3.881A pdb=" N VAL D 95 " --> pdb=" O LEU D 152 " (cutoff:3.500A) removed outlier: 3.538A pdb=" N LEU D 152 " --> pdb=" O VAL D 95 " (cutoff:3.500A) No H-bonds generated for sheet with id=AB3 Processing sheet with id=AB4, first strand: chain 'E' and resid 60 through 61 removed outlier: 3.517A pdb=" N VAL E 129 " --> pdb=" O ILE E 78 " (cutoff:3.500A) removed outlier: 3.692A pdb=" N ILE E 78 " --> pdb=" O VAL E 129 " (cutoff:3.500A) removed outlier: 6.727A pdb=" N MET E 79 " --> pdb=" O VAL E 106 " (cutoff:3.500A) removed outlier: 7.742A pdb=" N ILE E 108 " --> pdb=" O MET E 79 " (cutoff:3.500A) removed outlier: 6.859A pdb=" N LEU E 81 " --> pdb=" O ILE E 108 " (cutoff:3.500A) Processing sheet with id=AB5, first strand: chain 'F' and resid 13 through 16 removed outlier: 3.611A pdb=" N VAL F 14 " --> pdb=" O TYR F 6 " (cutoff:3.500A) removed outlier: 3.532A pdb=" N TYR F 5 " --> pdb=" O CYS F 28 " (cutoff:3.500A) removed outlier: 3.768A pdb=" N CYS F 28 " --> pdb=" O TYR F 5 " (cutoff:3.500A) Processing sheet with id=AB6, first strand: chain 'G' and resid 63 through 64 removed outlier: 3.559A pdb=" N TYR G 323 " --> pdb=" O ASN G 63 " (cutoff:3.500A) Processing sheet with id=AB7, first strand: chain 'G' and resid 119 through 120 removed outlier: 3.733A pdb=" N GLY G 119 " --> pdb=" O CYS G 135 " (cutoff:3.500A) removed outlier: 3.635A pdb=" N LEU G 132 " --> pdb=" O LEU G 94 " (cutoff:3.500A) removed outlier: 3.542A pdb=" N LEU G 94 " --> pdb=" O LEU G 132 " (cutoff:3.500A) removed outlier: 9.801A pdb=" N VAL G 261 " --> pdb=" O PHE G 91 " (cutoff:3.500A) removed outlier: 6.637A pdb=" N ASN G 93 " --> pdb=" O VAL G 261 " (cutoff:3.500A) Processing sheet with id=AB8, first strand: chain 'G' and resid 178 through 182 removed outlier: 3.828A pdb=" N ARG G 197 " --> pdb=" O LEU G 180 " (cutoff:3.500A) removed outlier: 3.521A pdb=" N GLU G 201 " --> pdb=" O SER G 219 " (cutoff:3.500A) removed outlier: 6.675A pdb=" N GLU G 152 " --> pdb=" O VAL G 249 " (cutoff:3.500A) removed outlier: 3.669A pdb=" N VAL G 249 " --> pdb=" O GLU G 152 " (cutoff:3.500A) removed outlier: 3.798A pdb=" N GLU G 245 " --> pdb=" O ILE G 156 " (cutoff:3.500A) removed outlier: 6.022A pdb=" N ARG G 158 " --> pdb=" O MET G 243 " (cutoff:3.500A) removed outlier: 5.562A pdb=" N MET G 243 " --> pdb=" O ARG G 158 " (cutoff:3.500A) removed outlier: 3.609A pdb=" N MET G 254 " --> pdb=" O SER G 251 " (cutoff:3.500A) Processing sheet with id=AB9, first strand: chain 'G' and resid 298 through 300 removed outlier: 5.841A pdb=" N LEU G 299 " --> pdb=" O ILE G 336 " (cutoff:3.500A) removed outlier: 3.670A pdb=" N GLN G 331 " --> pdb=" O ALA G 348 " (cutoff:3.500A) removed outlier: 3.666A pdb=" N CYS G 344 " --> pdb=" O VAL G 335 " (cutoff:3.500A) removed outlier: 3.594A pdb=" N ARG G 368 " --> pdb=" O MET G 345 " (cutoff:3.500A) removed outlier: 6.187A pdb=" N ILE G 347 " --> pdb=" O ILE G 366 " (cutoff:3.500A) removed outlier: 6.041A pdb=" N ILE G 366 " --> pdb=" O ILE G 347 " (cutoff:3.500A) Processing sheet with id=AC1, first strand: chain 'H' and resid 71 through 80 removed outlier: 5.278A pdb=" N PHE H 76 " --> pdb=" O LYS H 87 " (cutoff:3.500A) removed outlier: 5.791A pdb=" N LYS H 87 " --> pdb=" O PHE H 76 " (cutoff:3.500A) removed outlier: 7.401A pdb=" N TYR H 124 " --> pdb=" O ILE H 116 " (cutoff:3.500A) removed outlier: 5.526A pdb=" N ILE H 116 " --> pdb=" O TYR H 124 " (cutoff:3.500A) removed outlier: 4.915A pdb=" N SER H 126 " --> pdb=" O ASP H 114 " (cutoff:3.500A) removed outlier: 3.658A pdb=" N GLY H 111 " --> pdb=" O VAL H 102 " (cutoff:3.500A) removed outlier: 6.634A pdb=" N PHE H 144 " --> pdb=" O LEU H 73 " (cutoff:3.500A) Processing sheet with id=AC2, first strand: chain 'I' and resid 60 through 61 removed outlier: 3.891A pdb=" N ILE I 78 " --> pdb=" O VAL I 129 " (cutoff:3.500A) Processing sheet with id=AC3, first strand: chain 'J' and resid 13 through 16 removed outlier: 3.629A pdb=" N VAL J 14 " --> pdb=" O TYR J 6 " (cutoff:3.500A) removed outlier: 3.580A pdb=" N GLN J 4 " --> pdb=" O THR J 16 " (cutoff:3.500A) removed outlier: 3.639A pdb=" N CYS J 28 " --> pdb=" O TYR J 5 " (cutoff:3.500A) 385 hydrogen bonds defined for protein. 1071 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 42 hydrogen bonds 84 hydrogen bond angles 0 basepair planarities 16 basepair parallelities 57 stacking parallelities Total time for adding SS restraints: 5.21 Time building geometry restraints manager: 6.93 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: 4860 1.34 - 1.46: 3643 1.46 - 1.57: 8636 1.57 - 1.69: 289 1.69 - 1.81: 125 Bond restraints: 17553 Sorted by residual: bond pdb=" C LEU H 158 " pdb=" N PRO H 159 " ideal model delta sigma weight residual 1.334 1.354 -0.020 1.11e-02 8.12e+03 3.33e+00 bond pdb=" CA TYR E 91 " pdb=" CB TYR E 91 " ideal model delta sigma weight residual 1.528 1.568 -0.039 2.61e-02 1.47e+03 2.27e+00 bond pdb=" C GLY B 444 " pdb=" N PRO B 445 " ideal model delta sigma weight residual 1.335 1.354 -0.018 1.36e-02 5.41e+03 1.79e+00 bond pdb=" N LEU G 54 " pdb=" CA LEU G 54 " ideal model delta sigma weight residual 1.457 1.474 -0.017 1.29e-02 6.01e+03 1.70e+00 bond pdb=" N LEU C 54 " pdb=" CA LEU C 54 " ideal model delta sigma weight residual 1.459 1.474 -0.015 1.16e-02 7.43e+03 1.65e+00 ... (remaining 17548 not shown) Histogram of bond angle deviations from ideal: 100.09 - 106.87: 1362 106.87 - 113.65: 9799 113.65 - 120.43: 6430 120.43 - 127.20: 6403 127.20 - 133.98: 464 Bond angle restraints: 24458 Sorted by residual: angle pdb=" CA ASP D 91 " pdb=" C ASP D 91 " pdb=" N GLU D 92 " ideal model delta sigma weight residual 118.81 113.40 5.41 1.15e+00 7.56e-01 2.21e+01 angle pdb=" N VAL C 239 " pdb=" CA VAL C 239 " pdb=" C VAL C 239 " ideal model delta sigma weight residual 113.71 109.66 4.05 9.50e-01 1.11e+00 1.82e+01 angle pdb=" N VAL G 239 " pdb=" CA VAL G 239 " pdb=" C VAL G 239 " ideal model delta sigma weight residual 113.71 109.85 3.86 9.50e-01 1.11e+00 1.65e+01 angle pdb=" C PRO C 53 " pdb=" N LEU C 54 " pdb=" CA LEU C 54 " ideal model delta sigma weight residual 120.79 126.37 -5.58 1.39e+00 5.18e-01 1.61e+01 angle pdb=" N VAL C 363 " pdb=" CA VAL C 363 " pdb=" C VAL C 363 " ideal model delta sigma weight residual 113.53 109.70 3.83 9.80e-01 1.04e+00 1.53e+01 ... (remaining 24453 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 34.19: 10238 34.19 - 68.38: 511 68.38 - 102.57: 77 102.57 - 136.76: 1 136.76 - 170.96: 6 Dihedral angle restraints: 10833 sinusoidal: 5847 harmonic: 4986 Sorted by residual: dihedral pdb=" C4' C R 398 " pdb=" C3' C R 398 " pdb=" C2' C R 398 " pdb=" C1' C R 398 " ideal model delta sinusoidal sigma weight residual -35.00 33.12 -68.12 1 8.00e+00 1.56e-02 9.41e+01 dihedral pdb=" C5' C R 398 " pdb=" C4' C R 398 " pdb=" C3' C R 398 " pdb=" O3' C R 398 " ideal model delta sinusoidal sigma weight residual 147.00 81.29 65.71 1 8.00e+00 1.56e-02 8.83e+01 dihedral pdb=" O4' C R 398 " pdb=" C4' C R 398 " pdb=" C3' C R 398 " pdb=" C2' C R 398 " ideal model delta sinusoidal sigma weight residual 24.00 -34.36 58.36 1 8.00e+00 1.56e-02 7.13e+01 ... (remaining 10830 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.055: 2519 0.055 - 0.111: 262 0.111 - 0.166: 23 0.166 - 0.222: 0 0.222 - 0.277: 2 Chirality restraints: 2806 Sorted by residual: chirality pdb=" C3' C R 398 " pdb=" C4' C R 398 " pdb=" O3' C R 398 " pdb=" C2' C R 398 " both_signs ideal model delta sigma weight residual False -2.74 -2.47 -0.28 2.00e-01 2.50e+01 1.92e+00 chirality pdb=" C2' C R 398 " pdb=" C3' C R 398 " pdb=" O2' C R 398 " pdb=" C1' C R 398 " both_signs ideal model delta sigma weight residual False -2.52 -2.75 0.23 2.00e-01 2.50e+01 1.28e+00 chirality pdb=" P C R 394 " pdb=" OP1 C R 394 " pdb=" OP2 C R 394 " pdb=" O5' C R 394 " both_signs ideal model delta sigma weight residual True 2.41 -2.56 -0.15 2.00e-01 2.50e+01 5.90e-01 ... (remaining 2803 not shown) Planarity restraints: 2586 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C LEU H 158 " -0.027 5.00e-02 4.00e+02 4.06e-02 2.64e+00 pdb=" N PRO H 159 " 0.070 5.00e-02 4.00e+02 pdb=" CA PRO H 159 " -0.020 5.00e-02 4.00e+02 pdb=" CD PRO H 159 " -0.023 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C TRP G 52 " 0.025 5.00e-02 4.00e+02 3.85e-02 2.37e+00 pdb=" N PRO G 53 " -0.067 5.00e-02 4.00e+02 pdb=" CA PRO G 53 " 0.020 5.00e-02 4.00e+02 pdb=" CD PRO G 53 " 0.022 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C LEU B 500 " 0.025 5.00e-02 4.00e+02 3.71e-02 2.20e+00 pdb=" N PRO B 501 " -0.064 5.00e-02 4.00e+02 pdb=" CA PRO B 501 " 0.019 5.00e-02 4.00e+02 pdb=" CD PRO B 501 " 0.021 5.00e-02 4.00e+02 ... (remaining 2583 not shown) Histogram of nonbonded interaction distances: 2.26 - 2.79: 3643 2.79 - 3.32: 13897 3.32 - 3.84: 26781 3.84 - 4.37: 30435 4.37 - 4.90: 51103 Nonbonded interactions: 125859 Sorted by model distance: nonbonded pdb=" O LYS G 291 " pdb=" OG1 THR G 294 " model vdw 2.260 2.440 nonbonded pdb=" OD2 ASP B 422 " pdb=" OG1 THR B 424 " model vdw 2.285 2.440 nonbonded pdb=" O GLY C 383 " pdb=" OG SER C 387 " model vdw 2.290 2.440 nonbonded pdb=" OD1 ASP B 180 " pdb=" OG SER B 182 " model vdw 2.292 2.440 nonbonded pdb=" O ASP G 85 " pdb=" OG1 THR G 89 " model vdw 2.295 2.440 ... (remaining 125854 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = (chain 'C' and resid 47 through 395) selection = chain 'G' } ncs_group { reference = chain 'E' selection = (chain 'I' and resid 22 through 150) } ncs_group { reference = chain 'F' selection = (chain 'J' and resid 1 through 62) } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 2.990 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: 16.450 Check model and map are aligned: 0.260 Set scattering table: 0.170 Process input model: 47.890 Find NCS groups from input model: 0.770 Set up NCS constraints: 0.080 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:3.210 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 71.840 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.5619 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.002 0.039 17553 Z= 0.133 Angle : 0.579 6.823 24458 Z= 0.353 Chirality : 0.035 0.277 2806 Planarity : 0.004 0.041 2586 Dihedral : 18.684 170.955 7633 Min Nonbonded Distance : 2.260 Molprobity Statistics. All-atom Clashscore : 6.87 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.22 % Favored : 94.78 % Rotamer: Outliers : 0.00 % Allowed : 3.45 % Favored : 96.55 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.98 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.05 (0.16), residues: 1705 helix: -1.49 (0.23), residues: 379 sheet: -2.27 (0.27), residues: 301 loop : -3.49 (0.14), residues: 1025 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP C 52 HIS 0.002 0.000 HIS H 78 PHE 0.009 0.001 PHE B 419 TYR 0.008 0.001 TYR E 91 ARG 0.002 0.000 ARG C 141 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3410 Ramachandran restraints generated. 1705 Oldfield, 0 Emsley, 1705 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3410 Ramachandran restraints generated. 1705 Oldfield, 0 Emsley, 1705 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 319 residues out of total 1509 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 319 time to evaluate : 1.723 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 319 average time/residue: 0.4255 time to fit residues: 181.6547 Evaluate side-chains 139 residues out of total 1509 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 139 time to evaluate : 1.732 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 2.2580 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 181 random chunks: chunk 152 optimal weight: 10.0000 chunk 137 optimal weight: 0.4980 chunk 76 optimal weight: 20.0000 chunk 46 optimal weight: 8.9990 chunk 92 optimal weight: 10.0000 chunk 73 optimal weight: 6.9990 chunk 141 optimal weight: 0.8980 chunk 54 optimal weight: 9.9990 chunk 86 optimal weight: 9.9990 chunk 105 optimal weight: 9.9990 chunk 164 optimal weight: 3.9990 overall best weight: 4.2786 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 252 ASN ** B 515 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 93 ASN C 120 HIS C 147 GLN C 160 HIS C 244 GLN ** E 30 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 31 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** E 68 GLN F 4 GLN G 63 ASN G 160 HIS ** G 242 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** G 244 GLN G 307 ASN G 331 GLN H 64 GLN ** H 93 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** H 132 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** I 31 GLN I 32 ASN I 59 GLN I 68 GLN J 59 GLN Total number of N/Q/H flips: 19 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6351 moved from start: 0.4094 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.057 17553 Z= 0.364 Angle : 0.804 8.989 24458 Z= 0.408 Chirality : 0.045 0.231 2806 Planarity : 0.006 0.062 2586 Dihedral : 19.912 174.330 4178 Min Nonbonded Distance : 2.004 Molprobity Statistics. All-atom Clashscore : 15.16 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.92 % Favored : 92.08 % Rotamer: Outliers : 3.18 % Allowed : 17.43 % Favored : 79.39 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.98 % 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: -3.47 (0.17), residues: 1705 helix: -0.96 (0.24), residues: 382 sheet: -2.08 (0.26), residues: 346 loop : -3.07 (0.16), residues: 977 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.002 TRP G 380 HIS 0.009 0.002 HIS C 120 PHE 0.024 0.003 PHE C 59 TYR 0.020 0.003 TYR H 145 ARG 0.015 0.001 ARG G 227 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3410 Ramachandran restraints generated. 1705 Oldfield, 0 Emsley, 1705 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3410 Ramachandran restraints generated. 1705 Oldfield, 0 Emsley, 1705 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 190 residues out of total 1509 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 48 poor density : 142 time to evaluate : 1.732 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 48 outliers final: 31 residues processed: 177 average time/residue: 0.3228 time to fit residues: 84.2665 Evaluate side-chains 140 residues out of total 1509 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 31 poor density : 109 time to evaluate : 1.680 Switching outliers to nearest non-outliers outliers start: 31 outliers final: 0 residues processed: 31 average time/residue: 0.1756 time to fit residues: 11.8373 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 181 random chunks: chunk 91 optimal weight: 7.9990 chunk 51 optimal weight: 8.9990 chunk 136 optimal weight: 8.9990 chunk 111 optimal weight: 2.9990 chunk 45 optimal weight: 1.9990 chunk 164 optimal weight: 5.9990 chunk 177 optimal weight: 10.0000 chunk 146 optimal weight: 0.7980 chunk 163 optimal weight: 9.9990 chunk 56 optimal weight: 10.0000 chunk 132 optimal weight: 0.9980 overall best weight: 2.5586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 460 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 515 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 93 ASN ** D 59 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 30 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 31 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** E 68 GLN G 242 GLN H 106 ASN J 44 HIS Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6359 moved from start: 0.4468 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.050 17553 Z= 0.210 Angle : 0.604 10.098 24458 Z= 0.297 Chirality : 0.039 0.157 2806 Planarity : 0.004 0.037 2586 Dihedral : 19.820 170.831 4178 Min Nonbonded Distance : 2.129 Molprobity Statistics. All-atom Clashscore : 12.22 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.45 % Favored : 93.55 % Rotamer: Outliers : 2.25 % Allowed : 19.02 % Favored : 78.73 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.98 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.89 (0.18), residues: 1705 helix: -0.28 (0.26), residues: 382 sheet: -1.77 (0.26), residues: 334 loop : -2.78 (0.17), residues: 989 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP G 52 HIS 0.008 0.001 HIS E 132 PHE 0.013 0.001 PHE C 59 TYR 0.014 0.002 TYR F 5 ARG 0.005 0.001 ARG C 212 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3410 Ramachandran restraints generated. 1705 Oldfield, 0 Emsley, 1705 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3410 Ramachandran restraints generated. 1705 Oldfield, 0 Emsley, 1705 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 161 residues out of total 1509 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 34 poor density : 127 time to evaluate : 1.739 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 34 outliers final: 16 residues processed: 152 average time/residue: 0.2999 time to fit residues: 69.9168 Evaluate side-chains 121 residues out of total 1509 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 16 poor density : 105 time to evaluate : 1.481 Switching outliers to nearest non-outliers outliers start: 16 outliers final: 0 residues processed: 16 average time/residue: 0.1670 time to fit residues: 6.9043 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 181 random chunks: chunk 162 optimal weight: 2.9990 chunk 123 optimal weight: 3.9990 chunk 85 optimal weight: 7.9990 chunk 18 optimal weight: 20.0000 chunk 78 optimal weight: 10.0000 chunk 110 optimal weight: 6.9990 chunk 165 optimal weight: 9.9990 chunk 174 optimal weight: 9.9990 chunk 86 optimal weight: 10.0000 chunk 156 optimal weight: 4.9990 chunk 47 optimal weight: 5.9990 overall best weight: 4.9990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 272 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 460 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 515 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 264 HIS C 270 GLN ** D 59 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 78 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 142 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 30 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 31 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** E 132 HIS G 276 HIS ** H 142 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6599 moved from start: 0.6011 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.083 17553 Z= 0.361 Angle : 0.779 10.159 24458 Z= 0.391 Chirality : 0.045 0.165 2806 Planarity : 0.006 0.040 2586 Dihedral : 20.112 169.138 4178 Min Nonbonded Distance : 2.030 Molprobity Statistics. All-atom Clashscore : 17.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.56 % Favored : 90.44 % Rotamer: Outliers : 3.98 % Allowed : 21.87 % Favored : 74.16 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 1.98 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.97 (0.18), residues: 1705 helix: -0.73 (0.25), residues: 386 sheet: -1.55 (0.27), residues: 325 loop : -2.75 (0.17), residues: 994 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.017 0.002 TRP B 257 HIS 0.011 0.002 HIS B 376 PHE 0.021 0.002 PHE C 59 TYR 0.032 0.002 TYR F 5 ARG 0.012 0.001 ARG B 298 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3410 Ramachandran restraints generated. 1705 Oldfield, 0 Emsley, 1705 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3410 Ramachandran restraints generated. 1705 Oldfield, 0 Emsley, 1705 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 169 residues out of total 1509 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 60 poor density : 109 time to evaluate : 1.827 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 60 outliers final: 35 residues processed: 156 average time/residue: 0.2994 time to fit residues: 71.8688 Evaluate side-chains 134 residues out of total 1509 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 35 poor density : 99 time to evaluate : 1.740 Switching outliers to nearest non-outliers outliers start: 35 outliers final: 0 residues processed: 35 average time/residue: 0.1707 time to fit residues: 13.4742 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/modules/cctbx_project/iotbx/cli_parser.py", line 931, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/phenix/phenix/programs/real_space_refine.py", line 191, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 297, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 312, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 673, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 737, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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-1.21rc1-5172/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: 70.3286 > 50: distance: 2 - 26: 14.318 distance: 8 - 34: 30.793 distance: 15 - 41: 26.131 distance: 22 - 26: 16.763 distance: 26 - 27: 19.882 distance: 27 - 28: 10.068 distance: 27 - 30: 8.530 distance: 28 - 29: 20.395 distance: 28 - 34: 15.776 distance: 30 - 31: 24.684 distance: 31 - 32: 22.944 distance: 31 - 33: 3.982 distance: 34 - 35: 10.841 distance: 35 - 36: 21.731 distance: 35 - 38: 4.762 distance: 36 - 37: 4.349 distance: 36 - 41: 48.868 distance: 38 - 39: 41.407 distance: 38 - 40: 28.041 distance: 41 - 42: 22.813 distance: 42 - 43: 25.791 distance: 42 - 45: 11.096 distance: 43 - 44: 43.729 distance: 43 - 49: 43.382 distance: 45 - 46: 5.593 distance: 46 - 48: 12.839 distance: 49 - 50: 39.089 distance: 50 - 51: 32.061 distance: 50 - 53: 22.587 distance: 51 - 52: 52.307 distance: 51 - 58: 23.168 distance: 53 - 54: 4.655 distance: 54 - 55: 63.602 distance: 55 - 56: 7.509 distance: 55 - 57: 6.002 distance: 58 - 59: 18.865 distance: 59 - 60: 16.336 distance: 59 - 62: 7.717 distance: 60 - 61: 10.287 distance: 60 - 66: 25.455 distance: 62 - 63: 15.702 distance: 63 - 64: 26.713 distance: 66 - 67: 23.948 distance: 66 - 72: 15.390 distance: 67 - 68: 6.153 distance: 67 - 70: 12.124 distance: 68 - 69: 19.985 distance: 68 - 73: 15.613 distance: 70 - 71: 9.387 distance: 71 - 72: 27.916 distance: 73 - 74: 8.261 distance: 74 - 75: 17.965 distance: 74 - 77: 22.157 distance: 75 - 76: 12.054 distance: 75 - 81: 5.497 distance: 77 - 78: 40.517 distance: 77 - 79: 21.891 distance: 78 - 80: 24.928 distance: 81 - 82: 31.577 distance: 82 - 83: 9.687 distance: 83 - 84: 6.035 distance: 83 - 86: 27.053 distance: 86 - 87: 9.288 distance: 87 - 88: 20.790 distance: 87 - 90: 30.679 distance: 88 - 89: 3.243 distance: 88 - 95: 40.428 distance: 90 - 91: 17.517 distance: 91 - 92: 13.379 distance: 95 - 96: 5.550 distance: 95 - 101: 3.884 distance: 96 - 97: 13.290 distance: 96 - 99: 5.365 distance: 97 - 98: 14.239 distance: 97 - 102: 23.408 distance: 100 - 101: 6.008