Starting phenix.real_space_refine on Sat Dec 9 03:51:39 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7cfs_30346/12_2023/7cfs_30346.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7cfs_30346/12_2023/7cfs_30346.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.56 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7cfs_30346/12_2023/7cfs_30346.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7cfs_30346/12_2023/7cfs_30346.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7cfs_30346/12_2023/7cfs_30346.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7cfs_30346/12_2023/7cfs_30346.pdb" } resolution = 3.56 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 S 81 5.16 5 Na 1 4.78 5 C 6402 2.51 5 N 1631 2.21 5 O 1859 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A GLU 219": "OE1" <-> "OE2" Residue "A GLU 254": "OE1" <-> "OE2" Residue "A GLU 340": "OE1" <-> "OE2" Residue "A GLU 375": "OE1" <-> "OE2" Residue "B GLU 219": "OE1" <-> "OE2" Residue "B GLU 254": "OE1" <-> "OE2" Residue "B GLU 315": "OE1" <-> "OE2" Residue "B GLU 340": "OE1" <-> "OE2" Residue "B GLU 375": "OE1" <-> "OE2" Residue "C GLU 219": "OE1" <-> "OE2" Residue "C GLU 254": "OE1" <-> "OE2" Residue "C GLU 315": "OE1" <-> "OE2" Residue "C GLU 340": "OE1" <-> "OE2" Residue "C GLU 375": "OE1" <-> "OE2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/chem_data/mon_lib" Total number of atoms: 9974 Number of models: 1 Model: "" Number of chains: 6 Chain: "A" Number of atoms: 3249 Number of conformers: 1 Conformer: "" Number of residues, atoms: 426, 3249 Classifications: {'peptide': 426} Incomplete info: {'truncation_to_alanine': 44} Link IDs: {'PCIS': 1, 'PTRANS': 17, 'TRANS': 407} Unresolved non-hydrogen bonds: 177 Unresolved non-hydrogen angles: 217 Unresolved non-hydrogen dihedrals: 149 Unresolved non-hydrogen chiralities: 7 Planarities with less than four sites: {'GLN:plan1': 1, 'HIS:plan': 1, 'TYR:plan': 2, 'ASP:plan': 6, 'PHE:plan': 2, 'GLU:plan': 11, 'ARG:plan': 3} Unresolved non-hydrogen planarities: 112 Chain: "B" Number of atoms: 3267 Number of conformers: 1 Conformer: "" Number of residues, atoms: 426, 3267 Classifications: {'peptide': 426} Incomplete info: {'truncation_to_alanine': 42} Link IDs: {'PCIS': 1, 'PTRANS': 17, 'TRANS': 407} Unresolved chain link angles: 1 Unresolved non-hydrogen bonds: 159 Unresolved non-hydrogen angles: 195 Unresolved non-hydrogen dihedrals: 134 Unresolved non-hydrogen chiralities: 8 Planarities with less than four sites: {'GLN:plan1': 1, 'HIS:plan': 1, 'TYR:plan': 2, 'ASP:plan': 6, 'PHE:plan': 1, 'GLU:plan': 8, 'ARG:plan': 2} Unresolved non-hydrogen planarities: 89 Chain: "C" Number of atoms: 3268 Number of conformers: 1 Conformer: "" Number of residues, atoms: 426, 3268 Classifications: {'peptide': 426} Incomplete info: {'truncation_to_alanine': 41} Link IDs: {'PCIS': 1, 'PTRANS': 17, 'TRANS': 407} Unresolved non-hydrogen bonds: 156 Unresolved non-hydrogen angles: 189 Unresolved non-hydrogen dihedrals: 130 Unresolved non-hydrogen chiralities: 9 Planarities with less than four sites: {'GLN:plan1': 1, 'HIS:plan': 1, 'TYR:plan': 1, 'ASP:plan': 6, 'PHE:plan': 1, 'GLU:plan': 7, 'ARG:plan': 3} Unresolved non-hydrogen planarities: 83 Chain: "A" Number of atoms: 64 Number of conformers: 1 Conformer: "" Number of residues, atoms: 4, 64 Unusual residues: {' NA': 1, 'NAG': 2, 'Y01': 1} Classifications: {'undetermined': 4} Link IDs: {None: 3} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Chain: "B" Number of atoms: 63 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 63 Unusual residues: {'NAG': 2, 'Y01': 1} Classifications: {'undetermined': 3} Link IDs: {None: 2} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Chain: "C" Number of atoms: 63 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 63 Unusual residues: {'NAG': 2, 'Y01': 1} Classifications: {'undetermined': 3} Link IDs: {None: 2} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Time building chain proxies: 5.89, per 1000 atoms: 0.59 Number of scatterers: 9974 At special positions: 0 Unit cell: (97.344, 96.33, 129.792, 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 81 16.00 Na 1 11.00 O 1859 8.00 N 1631 7.00 C 6402 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=21, symmetry=0 Simple disulfide: pdb=" SG CYS A 93 " - pdb=" SG CYS A 194 " distance=2.03 Simple disulfide: pdb=" SG CYS A 172 " - pdb=" SG CYS A 179 " distance=2.03 Simple disulfide: pdb=" SG CYS A 290 " - pdb=" SG CYS A 367 " distance=2.03 Simple disulfide: pdb=" SG CYS A 310 " - pdb=" SG CYS A 363 " distance=2.03 Simple disulfide: pdb=" SG CYS A 314 " - pdb=" SG CYS A 361 " distance=2.33 Simple disulfide: pdb=" SG CYS A 323 " - pdb=" SG CYS A 345 " distance=2.03 Simple disulfide: pdb=" SG CYS A 325 " - pdb=" SG CYS A 337 " distance=2.03 Simple disulfide: pdb=" SG CYS B 93 " - pdb=" SG CYS B 194 " distance=2.03 Simple disulfide: pdb=" SG CYS B 172 " - pdb=" SG CYS B 179 " distance=2.03 Simple disulfide: pdb=" SG CYS B 290 " - pdb=" SG CYS B 367 " distance=2.03 Simple disulfide: pdb=" SG CYS B 310 " - pdb=" SG CYS B 363 " distance=2.03 Simple disulfide: pdb=" SG CYS B 314 " - pdb=" SG CYS B 361 " distance=2.33 Simple disulfide: pdb=" SG CYS B 323 " - pdb=" SG CYS B 345 " distance=2.03 Simple disulfide: pdb=" SG CYS B 325 " - pdb=" SG CYS B 337 " distance=2.03 Simple disulfide: pdb=" SG CYS C 93 " - pdb=" SG CYS C 194 " distance=2.03 Simple disulfide: pdb=" SG CYS C 172 " - pdb=" SG CYS C 179 " distance=2.03 Simple disulfide: pdb=" SG CYS C 290 " - pdb=" SG CYS C 367 " distance=2.02 Simple disulfide: pdb=" SG CYS C 310 " - pdb=" SG CYS C 363 " distance=2.03 Simple disulfide: pdb=" SG CYS C 314 " - pdb=" SG CYS C 361 " distance=2.33 Simple disulfide: pdb=" SG CYS C 323 " - pdb=" SG CYS C 345 " distance=2.03 Simple disulfide: pdb=" SG CYS C 325 " - pdb=" SG CYS C 337 " 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 601 " - " ASN A 395 " " NAG A 602 " - " ASN A 368 " " NAG B 601 " - " ASN B 395 " " NAG B 602 " - " ASN B 368 " " NAG C 601 " - " ASN C 395 " " NAG C 602 " - " ASN C 368 " Time building additional restraints: 3.55 Conformation dependent library (CDL) restraints added in 1.9 seconds 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2400 Finding SS restraints... Secondary structure from input PDB file: 51 helices and 15 sheets defined 35.2% alpha, 23.5% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.24 Creating SS restraints... Processing helix chain 'A' and resid 45 through 61 removed outlier: 4.264A pdb=" N LEU A 51 " --> pdb=" O ALA A 47 " (cutoff:3.500A) removed outlier: 3.777A pdb=" N CYS A 59 " --> pdb=" O ALA A 55 " (cutoff:3.500A) Processing helix chain 'A' and resid 67 through 71 removed outlier: 3.511A pdb=" N HIS A 70 " --> pdb=" O TYR A 67 " (cutoff:3.500A) Processing helix chain 'A' and resid 99 through 103 removed outlier: 3.830A pdb=" N VAL A 103 " --> pdb=" O PHE A 100 " (cutoff:3.500A) Processing helix chain 'A' and resid 106 through 111 removed outlier: 3.543A pdb=" N HIS A 110 " --> pdb=" O ASN A 106 " (cutoff:3.500A) Processing helix chain 'A' and resid 133 through 140 removed outlier: 3.581A pdb=" N ILE A 137 " --> pdb=" O LYS A 133 " (cutoff:3.500A) removed outlier: 4.049A pdb=" N LEU A 138 " --> pdb=" O GLN A 134 " (cutoff:3.500A) removed outlier: 3.694A pdb=" N ASP A 140 " --> pdb=" O GLU A 136 " (cutoff:3.500A) Processing helix chain 'A' and resid 154 through 162 removed outlier: 3.706A pdb=" N ARG A 160 " --> pdb=" O GLU A 156 " (cutoff:3.500A) Processing helix chain 'A' and resid 164 through 169 Processing helix chain 'A' and resid 180 through 184 removed outlier: 3.583A pdb=" N ASP A 183 " --> pdb=" O SER A 180 " (cutoff:3.500A) removed outlier: 3.585A pdb=" N PHE A 184 " --> pdb=" O ALA A 181 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 180 through 184' Processing helix chain 'A' and resid 225 through 229 removed outlier: 3.810A pdb=" N TYR A 229 " --> pdb=" O GLN A 226 " (cutoff:3.500A) Processing helix chain 'A' and resid 257 through 262 Processing helix chain 'A' and resid 310 through 323 removed outlier: 3.559A pdb=" N TYR A 318 " --> pdb=" O CYS A 314 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N LEU A 319 " --> pdb=" O GLU A 315 " (cutoff:3.500A) removed outlier: 3.571A pdb=" N CYS A 323 " --> pdb=" O LEU A 319 " (cutoff:3.500A) Processing helix chain 'A' and resid 338 through 343 removed outlier: 3.784A pdb=" N LYS A 343 " --> pdb=" O PRO A 339 " (cutoff:3.500A) Processing helix chain 'A' and resid 346 through 352 removed outlier: 3.548A pdb=" N PHE A 352 " --> pdb=" O PRO A 348 " (cutoff:3.500A) Processing helix chain 'A' and resid 386 through 394 removed outlier: 3.620A pdb=" N ALA A 391 " --> pdb=" O ALA A 387 " (cutoff:3.500A) removed outlier: 3.836A pdb=" N LYS A 392 " --> pdb=" O LYS A 388 " (cutoff:3.500A) Processing helix chain 'A' and resid 397 through 404 removed outlier: 3.513A pdb=" N GLY A 402 " --> pdb=" O GLU A 398 " (cutoff:3.500A) removed outlier: 3.516A pdb=" N GLU A 403 " --> pdb=" O GLN A 399 " (cutoff:3.500A) Processing helix chain 'A' and resid 427 through 443 removed outlier: 4.012A pdb=" N GLY A 433 " --> pdb=" O ALA A 429 " (cutoff:3.500A) removed outlier: 3.803A pdb=" N GLN A 438 " --> pdb=" O ASP A 434 " (cutoff:3.500A) removed outlier: 3.853A pdb=" N MET A 439 " --> pdb=" O ILE A 435 " (cutoff:3.500A) Processing helix chain 'A' and resid 446 through 462 removed outlier: 3.822A pdb=" N LEU A 451 " --> pdb=" O ILE A 447 " (cutoff:3.500A) removed outlier: 3.742A pdb=" N TYR A 458 " --> pdb=" O PHE A 454 " (cutoff:3.500A) removed outlier: 4.054A pdb=" N LYS A 462 " --> pdb=" O TYR A 458 " (cutoff:3.500A) Processing helix chain 'B' and resid 45 through 61 removed outlier: 4.263A pdb=" N LEU B 51 " --> pdb=" O ALA B 47 " (cutoff:3.500A) removed outlier: 3.777A pdb=" N CYS B 59 " --> pdb=" O ALA B 55 " (cutoff:3.500A) Processing helix chain 'B' and resid 67 through 71 removed outlier: 3.509A pdb=" N HIS B 70 " --> pdb=" O TYR B 67 " (cutoff:3.500A) Processing helix chain 'B' and resid 99 through 103 removed outlier: 3.830A pdb=" N VAL B 103 " --> pdb=" O PHE B 100 " (cutoff:3.500A) Processing helix chain 'B' and resid 106 through 111 removed outlier: 3.542A pdb=" N HIS B 110 " --> pdb=" O ASN B 106 " (cutoff:3.500A) Processing helix chain 'B' and resid 133 through 140 removed outlier: 3.582A pdb=" N ILE B 137 " --> pdb=" O LYS B 133 " (cutoff:3.500A) removed outlier: 4.049A pdb=" N LEU B 138 " --> pdb=" O GLN B 134 " (cutoff:3.500A) removed outlier: 3.695A pdb=" N ASP B 140 " --> pdb=" O GLU B 136 " (cutoff:3.500A) Processing helix chain 'B' and resid 154 through 162 removed outlier: 3.707A pdb=" N ARG B 160 " --> pdb=" O GLU B 156 " (cutoff:3.500A) Processing helix chain 'B' and resid 164 through 169 Processing helix chain 'B' and resid 180 through 184 removed outlier: 3.582A pdb=" N ASP B 183 " --> pdb=" O SER B 180 " (cutoff:3.500A) removed outlier: 3.584A pdb=" N PHE B 184 " --> pdb=" O ALA B 181 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 180 through 184' Processing helix chain 'B' and resid 225 through 229 removed outlier: 3.810A pdb=" N TYR B 229 " --> pdb=" O GLN B 226 " (cutoff:3.500A) Processing helix chain 'B' and resid 257 through 262 Processing helix chain 'B' and resid 310 through 323 removed outlier: 3.559A pdb=" N TYR B 318 " --> pdb=" O CYS B 314 " (cutoff:3.500A) removed outlier: 3.959A pdb=" N LEU B 319 " --> pdb=" O GLU B 315 " (cutoff:3.500A) removed outlier: 3.571A pdb=" N CYS B 323 " --> pdb=" O LEU B 319 " (cutoff:3.500A) Processing helix chain 'B' and resid 338 through 343 removed outlier: 3.784A pdb=" N LYS B 343 " --> pdb=" O PRO B 339 " (cutoff:3.500A) Processing helix chain 'B' and resid 346 through 352 removed outlier: 3.548A pdb=" N PHE B 352 " --> pdb=" O PRO B 348 " (cutoff:3.500A) Processing helix chain 'B' and resid 386 through 394 removed outlier: 3.620A pdb=" N ALA B 391 " --> pdb=" O ALA B 387 " (cutoff:3.500A) removed outlier: 3.835A pdb=" N LYS B 392 " --> pdb=" O LYS B 388 " (cutoff:3.500A) Processing helix chain 'B' and resid 397 through 404 removed outlier: 3.513A pdb=" N GLY B 402 " --> pdb=" O GLU B 398 " (cutoff:3.500A) removed outlier: 3.515A pdb=" N GLU B 403 " --> pdb=" O GLN B 399 " (cutoff:3.500A) Processing helix chain 'B' and resid 427 through 443 removed outlier: 4.012A pdb=" N GLY B 433 " --> pdb=" O ALA B 429 " (cutoff:3.500A) removed outlier: 3.803A pdb=" N GLN B 438 " --> pdb=" O ASP B 434 " (cutoff:3.500A) removed outlier: 3.851A pdb=" N MET B 439 " --> pdb=" O ILE B 435 " (cutoff:3.500A) Processing helix chain 'B' and resid 446 through 462 removed outlier: 3.822A pdb=" N LEU B 451 " --> pdb=" O ILE B 447 " (cutoff:3.500A) removed outlier: 3.743A pdb=" N TYR B 458 " --> pdb=" O PHE B 454 " (cutoff:3.500A) removed outlier: 4.053A pdb=" N LYS B 462 " --> pdb=" O TYR B 458 " (cutoff:3.500A) Processing helix chain 'C' and resid 45 through 61 removed outlier: 4.263A pdb=" N LEU C 51 " --> pdb=" O ALA C 47 " (cutoff:3.500A) removed outlier: 3.777A pdb=" N CYS C 59 " --> pdb=" O ALA C 55 " (cutoff:3.500A) Processing helix chain 'C' and resid 67 through 71 removed outlier: 3.510A pdb=" N HIS C 70 " --> pdb=" O TYR C 67 " (cutoff:3.500A) Processing helix chain 'C' and resid 99 through 103 removed outlier: 3.830A pdb=" N VAL C 103 " --> pdb=" O PHE C 100 " (cutoff:3.500A) Processing helix chain 'C' and resid 106 through 111 removed outlier: 3.543A pdb=" N HIS C 110 " --> pdb=" O ASN C 106 " (cutoff:3.500A) Processing helix chain 'C' and resid 133 through 140 removed outlier: 3.581A pdb=" N ILE C 137 " --> pdb=" O LYS C 133 " (cutoff:3.500A) removed outlier: 4.049A pdb=" N LEU C 138 " --> pdb=" O GLN C 134 " (cutoff:3.500A) removed outlier: 3.694A pdb=" N ASP C 140 " --> pdb=" O GLU C 136 " (cutoff:3.500A) Processing helix chain 'C' and resid 154 through 162 removed outlier: 3.707A pdb=" N ARG C 160 " --> pdb=" O GLU C 156 " (cutoff:3.500A) Processing helix chain 'C' and resid 164 through 169 Processing helix chain 'C' and resid 180 through 184 removed outlier: 3.583A pdb=" N ASP C 183 " --> pdb=" O SER C 180 " (cutoff:3.500A) removed outlier: 3.584A pdb=" N PHE C 184 " --> pdb=" O ALA C 181 " (cutoff:3.500A) No H-bonds generated for 'chain 'C' and resid 180 through 184' Processing helix chain 'C' and resid 225 through 229 removed outlier: 3.811A pdb=" N TYR C 229 " --> pdb=" O GLN C 226 " (cutoff:3.500A) Processing helix chain 'C' and resid 257 through 262 Processing helix chain 'C' and resid 310 through 323 removed outlier: 3.560A pdb=" N TYR C 318 " --> pdb=" O CYS C 314 " (cutoff:3.500A) removed outlier: 3.958A pdb=" N LEU C 319 " --> pdb=" O GLU C 315 " (cutoff:3.500A) removed outlier: 3.570A pdb=" N CYS C 323 " --> pdb=" O LEU C 319 " (cutoff:3.500A) Processing helix chain 'C' and resid 338 through 343 removed outlier: 3.783A pdb=" N LYS C 343 " --> pdb=" O PRO C 339 " (cutoff:3.500A) Processing helix chain 'C' and resid 346 through 352 removed outlier: 3.547A pdb=" N PHE C 352 " --> pdb=" O PRO C 348 " (cutoff:3.500A) Processing helix chain 'C' and resid 386 through 394 removed outlier: 3.620A pdb=" N ALA C 391 " --> pdb=" O ALA C 387 " (cutoff:3.500A) removed outlier: 3.835A pdb=" N LYS C 392 " --> pdb=" O LYS C 388 " (cutoff:3.500A) Processing helix chain 'C' and resid 397 through 404 removed outlier: 3.513A pdb=" N GLY C 402 " --> pdb=" O GLU C 398 " (cutoff:3.500A) removed outlier: 3.515A pdb=" N GLU C 403 " --> pdb=" O GLN C 399 " (cutoff:3.500A) Processing helix chain 'C' and resid 427 through 443 removed outlier: 4.013A pdb=" N GLY C 433 " --> pdb=" O ALA C 429 " (cutoff:3.500A) removed outlier: 3.802A pdb=" N GLN C 438 " --> pdb=" O ASP C 434 " (cutoff:3.500A) removed outlier: 3.852A pdb=" N MET C 439 " --> pdb=" O ILE C 435 " (cutoff:3.500A) Processing helix chain 'C' and resid 446 through 462 removed outlier: 3.822A pdb=" N LEU C 451 " --> pdb=" O ILE C 447 " (cutoff:3.500A) removed outlier: 3.742A pdb=" N TYR C 458 " --> pdb=" O PHE C 454 " (cutoff:3.500A) removed outlier: 4.053A pdb=" N LYS C 462 " --> pdb=" O TYR C 458 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 73 through 81 removed outlier: 10.264A pdb=" N PHE A 269 " --> pdb=" O LEU A 408 " (cutoff:3.500A) removed outlier: 10.837A pdb=" N ILE A 410 " --> pdb=" O PHE A 269 " (cutoff:3.500A) removed outlier: 10.376A pdb=" N THR A 271 " --> pdb=" O ILE A 410 " (cutoff:3.500A) removed outlier: 10.747A pdb=" N PHE A 412 " --> pdb=" O THR A 271 " (cutoff:3.500A) removed outlier: 9.267A pdb=" N VAL A 273 " --> pdb=" O PHE A 412 " (cutoff:3.500A) removed outlier: 9.727A pdb=" N VAL A 414 " --> pdb=" O VAL A 273 " (cutoff:3.500A) removed outlier: 7.729A pdb=" N CYS A 275 " --> pdb=" O VAL A 414 " (cutoff:3.500A) removed outlier: 7.740A pdb=" N ASN A 416 " --> pdb=" O CYS A 275 " (cutoff:3.500A) removed outlier: 4.663A pdb=" N GLU A 277 " --> pdb=" O ASN A 416 " (cutoff:3.500A) removed outlier: 3.635A pdb=" N ALA A 274 " --> pdb=" O GLU A 375 " (cutoff:3.500A) removed outlier: 3.668A pdb=" N GLN A 276 " --> pdb=" O GLY A 373 " (cutoff:3.500A) removed outlier: 3.504A pdb=" N GLN A 278 " --> pdb=" O ARG A 371 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 221 through 223 removed outlier: 10.264A pdb=" N PHE A 269 " --> pdb=" O LEU A 408 " (cutoff:3.500A) removed outlier: 10.837A pdb=" N ILE A 410 " --> pdb=" O PHE A 269 " (cutoff:3.500A) removed outlier: 10.376A pdb=" N THR A 271 " --> pdb=" O ILE A 410 " (cutoff:3.500A) removed outlier: 10.747A pdb=" N PHE A 412 " --> pdb=" O THR A 271 " (cutoff:3.500A) removed outlier: 9.267A pdb=" N VAL A 273 " --> pdb=" O PHE A 412 " (cutoff:3.500A) removed outlier: 9.727A pdb=" N VAL A 414 " --> pdb=" O VAL A 273 " (cutoff:3.500A) removed outlier: 7.729A pdb=" N CYS A 275 " --> pdb=" O VAL A 414 " (cutoff:3.500A) removed outlier: 7.740A pdb=" N ASN A 416 " --> pdb=" O CYS A 275 " (cutoff:3.500A) removed outlier: 4.663A pdb=" N GLU A 277 " --> pdb=" O ASN A 416 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 85 through 86 Processing sheet with id=AA4, first strand: chain 'A' and resid 185 through 189 removed outlier: 7.261A pdb=" N LYS A 193 " --> pdb=" O ASN A 94 " (cutoff:3.500A) removed outlier: 5.276A pdb=" N ASN A 94 " --> pdb=" O LYS A 193 " (cutoff:3.500A) removed outlier: 3.609A pdb=" N HIS A 250 " --> pdb=" O ALA A 89 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'A' and resid 173 through 174 Processing sheet with id=AA6, first strand: chain 'B' and resid 73 through 81 removed outlier: 10.264A pdb=" N PHE B 269 " --> pdb=" O LEU B 408 " (cutoff:3.500A) removed outlier: 10.838A pdb=" N ILE B 410 " --> pdb=" O PHE B 269 " (cutoff:3.500A) removed outlier: 10.376A pdb=" N THR B 271 " --> pdb=" O ILE B 410 " (cutoff:3.500A) removed outlier: 10.747A pdb=" N PHE B 412 " --> pdb=" O THR B 271 " (cutoff:3.500A) removed outlier: 9.267A pdb=" N VAL B 273 " --> pdb=" O PHE B 412 " (cutoff:3.500A) removed outlier: 9.726A pdb=" N VAL B 414 " --> pdb=" O VAL B 273 " (cutoff:3.500A) removed outlier: 7.729A pdb=" N CYS B 275 " --> pdb=" O VAL B 414 " (cutoff:3.500A) removed outlier: 7.740A pdb=" N ASN B 416 " --> pdb=" O CYS B 275 " (cutoff:3.500A) removed outlier: 4.662A pdb=" N GLU B 277 " --> pdb=" O ASN B 416 " (cutoff:3.500A) removed outlier: 3.634A pdb=" N ALA B 274 " --> pdb=" O GLU B 375 " (cutoff:3.500A) removed outlier: 3.668A pdb=" N GLN B 276 " --> pdb=" O GLY B 373 " (cutoff:3.500A) removed outlier: 3.503A pdb=" N GLN B 278 " --> pdb=" O ARG B 371 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 221 through 223 removed outlier: 10.264A pdb=" N PHE B 269 " --> pdb=" O LEU B 408 " (cutoff:3.500A) removed outlier: 10.838A pdb=" N ILE B 410 " --> pdb=" O PHE B 269 " (cutoff:3.500A) removed outlier: 10.376A pdb=" N THR B 271 " --> pdb=" O ILE B 410 " (cutoff:3.500A) removed outlier: 10.747A pdb=" N PHE B 412 " --> pdb=" O THR B 271 " (cutoff:3.500A) removed outlier: 9.267A pdb=" N VAL B 273 " --> pdb=" O PHE B 412 " (cutoff:3.500A) removed outlier: 9.726A pdb=" N VAL B 414 " --> pdb=" O VAL B 273 " (cutoff:3.500A) removed outlier: 7.729A pdb=" N CYS B 275 " --> pdb=" O VAL B 414 " (cutoff:3.500A) removed outlier: 7.740A pdb=" N ASN B 416 " --> pdb=" O CYS B 275 " (cutoff:3.500A) removed outlier: 4.662A pdb=" N GLU B 277 " --> pdb=" O ASN B 416 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'B' and resid 85 through 86 Processing sheet with id=AA9, first strand: chain 'B' and resid 185 through 189 removed outlier: 7.262A pdb=" N LYS B 193 " --> pdb=" O ASN B 94 " (cutoff:3.500A) removed outlier: 5.276A pdb=" N ASN B 94 " --> pdb=" O LYS B 193 " (cutoff:3.500A) removed outlier: 3.609A pdb=" N HIS B 250 " --> pdb=" O ALA B 89 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'B' and resid 173 through 174 Processing sheet with id=AB2, first strand: chain 'C' and resid 73 through 81 removed outlier: 10.264A pdb=" N PHE C 269 " --> pdb=" O LEU C 408 " (cutoff:3.500A) removed outlier: 10.838A pdb=" N ILE C 410 " --> pdb=" O PHE C 269 " (cutoff:3.500A) removed outlier: 10.376A pdb=" N THR C 271 " --> pdb=" O ILE C 410 " (cutoff:3.500A) removed outlier: 10.747A pdb=" N PHE C 412 " --> pdb=" O THR C 271 " (cutoff:3.500A) removed outlier: 9.268A pdb=" N VAL C 273 " --> pdb=" O PHE C 412 " (cutoff:3.500A) removed outlier: 9.726A pdb=" N VAL C 414 " --> pdb=" O VAL C 273 " (cutoff:3.500A) removed outlier: 7.729A pdb=" N CYS C 275 " --> pdb=" O VAL C 414 " (cutoff:3.500A) removed outlier: 7.740A pdb=" N ASN C 416 " --> pdb=" O CYS C 275 " (cutoff:3.500A) removed outlier: 4.663A pdb=" N GLU C 277 " --> pdb=" O ASN C 416 " (cutoff:3.500A) removed outlier: 3.635A pdb=" N ALA C 274 " --> pdb=" O GLU C 375 " (cutoff:3.500A) removed outlier: 3.669A pdb=" N GLN C 276 " --> pdb=" O GLY C 373 " (cutoff:3.500A) removed outlier: 3.503A pdb=" N GLN C 278 " --> pdb=" O ARG C 371 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'C' and resid 221 through 223 removed outlier: 10.264A pdb=" N PHE C 269 " --> pdb=" O LEU C 408 " (cutoff:3.500A) removed outlier: 10.838A pdb=" N ILE C 410 " --> pdb=" O PHE C 269 " (cutoff:3.500A) removed outlier: 10.376A pdb=" N THR C 271 " --> pdb=" O ILE C 410 " (cutoff:3.500A) removed outlier: 10.747A pdb=" N PHE C 412 " --> pdb=" O THR C 271 " (cutoff:3.500A) removed outlier: 9.268A pdb=" N VAL C 273 " --> pdb=" O PHE C 412 " (cutoff:3.500A) removed outlier: 9.726A pdb=" N VAL C 414 " --> pdb=" O VAL C 273 " (cutoff:3.500A) removed outlier: 7.729A pdb=" N CYS C 275 " --> pdb=" O VAL C 414 " (cutoff:3.500A) removed outlier: 7.740A pdb=" N ASN C 416 " --> pdb=" O CYS C 275 " (cutoff:3.500A) removed outlier: 4.663A pdb=" N GLU C 277 " --> pdb=" O ASN C 416 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'C' and resid 85 through 86 Processing sheet with id=AB5, first strand: chain 'C' and resid 185 through 189 removed outlier: 7.262A pdb=" N LYS C 193 " --> pdb=" O ASN C 94 " (cutoff:3.500A) removed outlier: 5.276A pdb=" N ASN C 94 " --> pdb=" O LYS C 193 " (cutoff:3.500A) removed outlier: 3.610A pdb=" N HIS C 250 " --> pdb=" O ALA C 89 " (cutoff:3.500A) Processing sheet with id=AB6, first strand: chain 'C' and resid 173 through 174 303 hydrogen bonds defined for protein. 810 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: 2.77 Time building geometry restraints manager: 3.92 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: 3154 1.34 - 1.45: 1745 1.45 - 1.57: 5202 1.57 - 1.69: 0 1.69 - 1.81: 108 Bond restraints: 10209 Sorted by residual: bond pdb=" CAI Y01 B 603 " pdb=" CAK Y01 B 603 " ideal model delta sigma weight residual 1.492 1.386 0.106 2.00e-02 2.50e+03 2.80e+01 bond pdb=" CAI Y01 C 603 " pdb=" CAK Y01 C 603 " ideal model delta sigma weight residual 1.492 1.386 0.106 2.00e-02 2.50e+03 2.79e+01 bond pdb=" CAI Y01 A 603 " pdb=" CAK Y01 A 603 " ideal model delta sigma weight residual 1.492 1.387 0.105 2.00e-02 2.50e+03 2.78e+01 bond pdb=" CAI Y01 C 603 " pdb=" CAZ Y01 C 603 " ideal model delta sigma weight residual 1.332 1.380 -0.048 2.00e-02 2.50e+03 5.79e+00 bond pdb=" CAI Y01 A 603 " pdb=" CAZ Y01 A 603 " ideal model delta sigma weight residual 1.332 1.380 -0.048 2.00e-02 2.50e+03 5.78e+00 ... (remaining 10204 not shown) Histogram of bond angle deviations from ideal: 96.80 - 104.23: 175 104.23 - 111.65: 4341 111.65 - 119.08: 3783 119.08 - 126.51: 5464 126.51 - 133.94: 119 Bond angle restraints: 13882 Sorted by residual: angle pdb=" N ILE A 381 " pdb=" CA ILE A 381 " pdb=" C ILE A 381 " ideal model delta sigma weight residual 108.88 117.02 -8.14 2.16e+00 2.14e-01 1.42e+01 angle pdb=" N ILE C 381 " pdb=" CA ILE C 381 " pdb=" C ILE C 381 " ideal model delta sigma weight residual 108.88 116.98 -8.10 2.16e+00 2.14e-01 1.41e+01 angle pdb=" N ILE B 381 " pdb=" CA ILE B 381 " pdb=" C ILE B 381 " ideal model delta sigma weight residual 108.88 116.97 -8.09 2.16e+00 2.14e-01 1.40e+01 angle pdb=" CA ILE C 381 " pdb=" C ILE C 381 " pdb=" N PRO C 382 " ideal model delta sigma weight residual 118.88 124.19 -5.31 1.54e+00 4.22e-01 1.19e+01 angle pdb=" CA ILE B 381 " pdb=" C ILE B 381 " pdb=" N PRO B 382 " ideal model delta sigma weight residual 118.88 124.17 -5.29 1.54e+00 4.22e-01 1.18e+01 ... (remaining 13877 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 16.89: 5486 16.89 - 33.78: 656 33.78 - 50.68: 141 50.68 - 67.57: 39 67.57 - 84.46: 21 Dihedral angle restraints: 6343 sinusoidal: 2607 harmonic: 3736 Sorted by residual: dihedral pdb=" CB CYS C 323 " pdb=" SG CYS C 323 " pdb=" SG CYS C 345 " pdb=" CB CYS C 345 " ideal model delta sinusoidal sigma weight residual -86.00 -170.46 84.46 1 1.00e+01 1.00e-02 8.67e+01 dihedral pdb=" CB CYS B 323 " pdb=" SG CYS B 323 " pdb=" SG CYS B 345 " pdb=" CB CYS B 345 " ideal model delta sinusoidal sigma weight residual -86.00 -170.45 84.45 1 1.00e+01 1.00e-02 8.67e+01 dihedral pdb=" CB CYS A 323 " pdb=" SG CYS A 323 " pdb=" SG CYS A 345 " pdb=" CB CYS A 345 " ideal model delta sinusoidal sigma weight residual -86.00 -170.45 84.45 1 1.00e+01 1.00e-02 8.67e+01 ... (remaining 6340 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.067: 1248 0.067 - 0.135: 270 0.135 - 0.202: 33 0.202 - 0.270: 0 0.270 - 0.337: 6 Chirality restraints: 1557 Sorted by residual: chirality pdb=" CBG Y01 C 603 " pdb=" CAQ Y01 C 603 " pdb=" CBD Y01 C 603 " pdb=" CBI Y01 C 603 " both_signs ideal model delta sigma weight residual False -2.33 -2.67 0.34 2.00e-01 2.50e+01 2.84e+00 chirality pdb=" CBG Y01 B 603 " pdb=" CAQ Y01 B 603 " pdb=" CBD Y01 B 603 " pdb=" CBI Y01 B 603 " both_signs ideal model delta sigma weight residual False -2.33 -2.67 0.33 2.00e-01 2.50e+01 2.79e+00 chirality pdb=" CBG Y01 A 603 " pdb=" CAQ Y01 A 603 " pdb=" CBD Y01 A 603 " pdb=" CBI Y01 A 603 " both_signs ideal model delta sigma weight residual False -2.33 -2.66 0.33 2.00e-01 2.50e+01 2.77e+00 ... (remaining 1554 not shown) Planarity restraints: 1795 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C THR B 338 " -0.031 5.00e-02 4.00e+02 4.73e-02 3.58e+00 pdb=" N PRO B 339 " 0.082 5.00e-02 4.00e+02 pdb=" CA PRO B 339 " -0.024 5.00e-02 4.00e+02 pdb=" CD PRO B 339 " -0.026 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C THR A 338 " 0.031 5.00e-02 4.00e+02 4.73e-02 3.58e+00 pdb=" N PRO A 339 " -0.082 5.00e-02 4.00e+02 pdb=" CA PRO A 339 " 0.024 5.00e-02 4.00e+02 pdb=" CD PRO A 339 " 0.026 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C THR C 338 " -0.031 5.00e-02 4.00e+02 4.72e-02 3.56e+00 pdb=" N PRO C 339 " 0.082 5.00e-02 4.00e+02 pdb=" CA PRO C 339 " -0.024 5.00e-02 4.00e+02 pdb=" CD PRO C 339 " -0.026 5.00e-02 4.00e+02 ... (remaining 1792 not shown) Histogram of nonbonded interaction distances: 1.83 - 2.45: 25 2.45 - 3.06: 5484 3.06 - 3.67: 14301 3.67 - 4.29: 20002 4.29 - 4.90: 33905 Nonbonded interactions: 73717 Sorted by model distance: nonbonded pdb=" CE1 HIS C 73 " pdb=" OE1 GLN C 422 " model vdw 1.835 3.260 nonbonded pdb=" ND1 HIS C 73 " pdb=" OE1 GLN C 422 " model vdw 1.971 2.520 nonbonded pdb=" CD2 HIS A 73 " pdb=" OE1 GLN A 422 " model vdw 1.990 3.260 nonbonded pdb=" CD2 HIS B 73 " pdb=" OE1 GLN B 422 " model vdw 2.116 3.260 nonbonded pdb=" OD1 ASP A 434 " pdb="NA NA A 604 " model vdw 2.251 2.470 ... (remaining 73712 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 40 through 87 or (resid 88 through 89 and (name N or name \ CA or name C or name O or name CB )) or resid 90 through 198 or (resid 199 and ( \ name N or name CA or name C or name O or name CB )) or resid 200 through 290 or \ (resid 291 through 300 and (name N or name CA or name C or name O or name CB )) \ or resid 301 through 452 or (resid 453 through 455 and (name N or name CA or nam \ e C or name O or name CB )) or resid 456 through 464 or (resid 465 and (name N o \ r name CA or name C or name O or name CB )) or resid 601 through 603)) selection = (chain 'B' and (resid 40 through 42 or (resid 43 through 44 and (name N or name \ CA or name C or name O or name CB )) or resid 45 through 146 or (resid 147 throu \ gh 148 and (name N or name CA or name C or name O or name CB )) or resid 149 thr \ ough 181 or (resid 182 and (name N or name CA or name C or name O or name CB )) \ or resid 183 through 290 or (resid 291 through 300 and (name N or name CA or nam \ e C or name O or name CB )) or resid 301 through 314 or (resid 315 and (name N o \ r name CA or name C or name O or name CB )) or resid 316 through 363 or (resid 3 \ 64 and (name N or name CA or name C or name O or name CB )) or resid 365 through \ 452 or (resid 453 through 455 and (name N or name CA or name C or name O or nam \ e CB )) or resid 456 through 465 or resid 601 through 603)) selection = (chain 'C' and (resid 40 through 87 or (resid 88 through 89 and (name N or name \ CA or name C or name O or name CB )) or resid 90 through 146 or (resid 147 throu \ gh 148 and (name N or name CA or name C or name O or name CB )) or resid 149 thr \ ough 181 or (resid 182 and (name N or name CA or name C or name O or name CB )) \ or resid 183 through 198 or (resid 199 and (name N or name CA or name C or name \ O or name CB )) or resid 200 through 234 or (resid 235 and (name N or name CA or \ name C or name O or name CB )) or resid 236 through 314 or (resid 315 and (name \ N or name CA or name C or name O or name CB )) or resid 316 through 363 or (res \ id 364 and (name N or name CA or name C or name O or name CB )) or resid 365 thr \ ough 457 or (resid 458 through 459 and (name N or name CA or name C or name O or \ name CB )) or resid 460 through 465 or resid 601 through 603)) } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.350 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.010 Extract box with map and model: 2.140 Check model and map are aligned: 0.140 Set scattering table: 0.090 Process input model: 28.920 Find NCS groups from input model: 0.740 Set up NCS constraints: 0.060 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:2.870 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 35.320 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8220 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.008 0.106 10209 Z= 0.495 Angle : 0.936 9.704 13882 Z= 0.481 Chirality : 0.058 0.337 1557 Planarity : 0.005 0.047 1789 Dihedral : 15.513 84.427 3880 Min Nonbonded Distance : 1.835 Molprobity Statistics. All-atom Clashscore : 5.64 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.90 % Favored : 94.10 % Rotamer: Outliers : 0.60 % Allowed : 8.97 % Favored : 90.42 % 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: -5.35 (0.18), residues: 1272 helix: -4.61 (0.10), residues: 342 sheet: -3.63 (0.23), residues: 300 loop : -2.77 (0.21), residues: 630 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.004 0.001 TRP C 46 HIS 0.004 0.001 HIS B 329 PHE 0.018 0.002 PHE B 263 TYR 0.016 0.001 TYR C 426 ARG 0.003 0.000 ARG A 371 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 188 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 182 time to evaluate : 1.143 Fit side-chains outliers start: 6 outliers final: 0 residues processed: 188 average time/residue: 0.2691 time to fit residues: 67.1227 Evaluate side-chains 127 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 127 time to evaluate : 1.180 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.6456 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 126 random chunks: chunk 106 optimal weight: 0.5980 chunk 95 optimal weight: 4.9990 chunk 52 optimal weight: 0.7980 chunk 32 optimal weight: 9.9990 chunk 64 optimal weight: 0.9990 chunk 51 optimal weight: 0.8980 chunk 98 optimal weight: 3.9990 chunk 38 optimal weight: 5.9990 chunk 60 optimal weight: 2.9990 chunk 73 optimal weight: 0.6980 chunk 114 optimal weight: 0.2980 overall best weight: 0.6580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 73 HIS A 276 GLN ** A 438 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 73 HIS B 276 GLN B 438 GLN C 276 GLN C 329 HIS C 438 GLN Total number of N/Q/H flips: 8 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8137 moved from start: 0.1932 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.037 10209 Z= 0.155 Angle : 0.565 7.474 13882 Z= 0.286 Chirality : 0.040 0.129 1557 Planarity : 0.004 0.037 1789 Dihedral : 6.243 57.419 1721 Min Nonbonded Distance : 2.330 Molprobity Statistics. All-atom Clashscore : 5.43 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.07 % Favored : 96.93 % Rotamer: Outliers : 1.81 % Allowed : 16.53 % Favored : 81.65 % 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: -4.08 (0.20), residues: 1272 helix: -2.83 (0.19), residues: 387 sheet: -3.38 (0.25), residues: 279 loop : -2.12 (0.23), residues: 606 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP B 46 HIS 0.005 0.001 HIS B 73 PHE 0.017 0.001 PHE A 301 TYR 0.020 0.001 TYR C 67 ARG 0.003 0.000 ARG A 326 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 174 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 156 time to evaluate : 1.184 Fit side-chains outliers start: 18 outliers final: 11 residues processed: 173 average time/residue: 0.1875 time to fit residues: 47.9775 Evaluate side-chains 138 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 127 time to evaluate : 1.141 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 0 residues processed: 11 average time/residue: 0.1180 time to fit residues: 3.7077 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 126 random chunks: chunk 63 optimal weight: 5.9990 chunk 35 optimal weight: 3.9990 chunk 95 optimal weight: 4.9990 chunk 77 optimal weight: 4.9990 chunk 31 optimal weight: 5.9990 chunk 114 optimal weight: 6.9990 chunk 123 optimal weight: 4.9990 chunk 102 optimal weight: 0.9980 chunk 113 optimal weight: 8.9990 chunk 39 optimal weight: 3.9990 chunk 91 optimal weight: 0.8980 overall best weight: 2.9786 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 72 HIS A 252 GLN A 438 GLN B 66 GLN B 72 HIS C 72 HIS C 252 GLN C 438 GLN Total number of N/Q/H flips: 8 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8223 moved from start: 0.2282 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.053 10209 Z= 0.390 Angle : 0.667 14.729 13882 Z= 0.331 Chirality : 0.043 0.145 1557 Planarity : 0.003 0.028 1789 Dihedral : 6.224 58.865 1721 Min Nonbonded Distance : 2.227 Molprobity Statistics. All-atom Clashscore : 6.73 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.32 % Favored : 95.68 % Rotamer: Outliers : 2.42 % Allowed : 19.96 % Favored : 77.62 % 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: -3.35 (0.22), residues: 1272 helix: -1.66 (0.25), residues: 372 sheet: -3.21 (0.25), residues: 303 loop : -1.89 (0.25), residues: 597 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.001 TRP B 46 HIS 0.003 0.001 HIS B 250 PHE 0.012 0.002 PHE C 301 TYR 0.013 0.001 TYR C 67 ARG 0.005 0.000 ARG B 160 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 163 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 24 poor density : 139 time to evaluate : 1.148 Fit side-chains revert: symmetry clash outliers start: 24 outliers final: 11 residues processed: 154 average time/residue: 0.1726 time to fit residues: 40.9117 Evaluate side-chains 139 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 128 time to evaluate : 1.190 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 0 residues processed: 11 average time/residue: 0.1229 time to fit residues: 3.8400 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 126 random chunks: chunk 113 optimal weight: 8.9990 chunk 86 optimal weight: 3.9990 chunk 59 optimal weight: 3.9990 chunk 12 optimal weight: 7.9990 chunk 54 optimal weight: 4.9990 chunk 76 optimal weight: 2.9990 chunk 115 optimal weight: 5.9990 chunk 121 optimal weight: 0.2980 chunk 60 optimal weight: 1.9990 chunk 109 optimal weight: 0.4980 chunk 32 optimal weight: 0.9980 overall best weight: 1.3584 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 438 GLN C 66 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8173 moved from start: 0.2564 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.028 10209 Z= 0.209 Angle : 0.577 13.909 13882 Z= 0.282 Chirality : 0.041 0.132 1557 Planarity : 0.003 0.026 1789 Dihedral : 5.614 56.950 1721 Min Nonbonded Distance : 2.276 Molprobity Statistics. All-atom Clashscore : 6.26 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.36 % Favored : 97.64 % Rotamer: Outliers : 2.42 % Allowed : 20.77 % Favored : 76.81 % 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: -2.79 (0.23), residues: 1272 helix: -1.07 (0.26), residues: 390 sheet: -3.28 (0.26), residues: 300 loop : -1.41 (0.26), residues: 582 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP B 46 HIS 0.003 0.001 HIS C 73 PHE 0.008 0.001 PHE C 98 TYR 0.014 0.001 TYR C 67 ARG 0.003 0.000 ARG A 201 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 168 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 24 poor density : 144 time to evaluate : 1.198 Fit side-chains revert: symmetry clash outliers start: 24 outliers final: 14 residues processed: 158 average time/residue: 0.1581 time to fit residues: 39.3139 Evaluate side-chains 143 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 129 time to evaluate : 1.149 Switching outliers to nearest non-outliers outliers start: 14 outliers final: 0 residues processed: 14 average time/residue: 0.0900 time to fit residues: 3.9272 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 126 random chunks: chunk 101 optimal weight: 4.9990 chunk 69 optimal weight: 0.5980 chunk 1 optimal weight: 8.9990 chunk 90 optimal weight: 7.9990 chunk 50 optimal weight: 0.0070 chunk 103 optimal weight: 5.9990 chunk 84 optimal weight: 20.0000 chunk 0 optimal weight: 10.0000 chunk 62 optimal weight: 3.9990 chunk 109 optimal weight: 2.9990 chunk 30 optimal weight: 7.9990 overall best weight: 2.5204 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8218 moved from start: 0.2705 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.044 10209 Z= 0.340 Angle : 0.650 16.126 13882 Z= 0.316 Chirality : 0.043 0.133 1557 Planarity : 0.003 0.038 1789 Dihedral : 5.702 57.687 1721 Min Nonbonded Distance : 2.214 Molprobity Statistics. All-atom Clashscore : 6.68 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.25 % Favored : 95.75 % Rotamer: Outliers : 1.81 % Allowed : 22.18 % Favored : 76.01 % 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: -2.48 (0.24), residues: 1272 helix: -0.58 (0.28), residues: 372 sheet: -3.01 (0.27), residues: 303 loop : -1.43 (0.27), residues: 597 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP B 46 HIS 0.004 0.001 HIS C 73 PHE 0.013 0.001 PHE C 442 TYR 0.012 0.001 TYR A 318 ARG 0.003 0.000 ARG A 371 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 149 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 131 time to evaluate : 1.089 Fit side-chains outliers start: 18 outliers final: 8 residues processed: 141 average time/residue: 0.1613 time to fit residues: 36.1900 Evaluate side-chains 134 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 126 time to evaluate : 1.175 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 0 residues processed: 8 average time/residue: 0.1067 time to fit residues: 3.0573 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 126 random chunks: chunk 40 optimal weight: 0.9980 chunk 109 optimal weight: 0.0070 chunk 24 optimal weight: 0.0270 chunk 71 optimal weight: 5.9990 chunk 30 optimal weight: 3.9990 chunk 121 optimal weight: 1.9990 chunk 101 optimal weight: 1.9990 chunk 56 optimal weight: 6.9990 chunk 10 optimal weight: 0.0270 chunk 64 optimal weight: 0.8980 chunk 117 optimal weight: 1.9990 overall best weight: 0.3914 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 66 GLN A 329 HIS A 438 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8119 moved from start: 0.3017 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.026 10209 Z= 0.125 Angle : 0.554 14.925 13882 Z= 0.261 Chirality : 0.039 0.129 1557 Planarity : 0.002 0.031 1789 Dihedral : 5.094 55.928 1721 Min Nonbonded Distance : 2.340 Molprobity Statistics. All-atom Clashscore : 4.50 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.04 % Favored : 97.96 % Rotamer: Outliers : 0.81 % Allowed : 22.78 % Favored : 76.41 % 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: -2.11 (0.25), residues: 1272 helix: -0.68 (0.27), residues: 408 sheet: -2.88 (0.30), residues: 258 loop : -1.05 (0.28), residues: 606 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A 46 HIS 0.002 0.000 HIS B 329 PHE 0.007 0.001 PHE C 98 TYR 0.012 0.001 TYR A 417 ARG 0.002 0.000 ARG A 201 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 167 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 159 time to evaluate : 1.169 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 8 outliers final: 3 residues processed: 163 average time/residue: 0.1650 time to fit residues: 42.1003 Evaluate side-chains 145 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 142 time to evaluate : 1.188 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1002 time to fit residues: 2.0151 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 126 random chunks: chunk 13 optimal weight: 0.0770 chunk 69 optimal weight: 5.9990 chunk 89 optimal weight: 0.7980 chunk 68 optimal weight: 0.6980 chunk 102 optimal weight: 3.9990 chunk 121 optimal weight: 0.9990 chunk 76 optimal weight: 3.9990 chunk 74 optimal weight: 4.9990 chunk 56 optimal weight: 2.9990 chunk 75 optimal weight: 1.9990 chunk 48 optimal weight: 6.9990 overall best weight: 0.9142 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 438 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8137 moved from start: 0.3104 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.026 10209 Z= 0.165 Angle : 0.559 14.262 13882 Z= 0.266 Chirality : 0.040 0.133 1557 Planarity : 0.003 0.031 1789 Dihedral : 5.082 56.289 1721 Min Nonbonded Distance : 2.253 Molprobity Statistics. All-atom Clashscore : 5.02 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.97 % Favored : 98.03 % Rotamer: Outliers : 1.01 % Allowed : 23.99 % Favored : 75.00 % 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: -1.98 (0.25), residues: 1272 helix: -0.53 (0.27), residues: 411 sheet: -2.96 (0.29), residues: 279 loop : -0.82 (0.28), residues: 582 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.004 0.001 TRP C 233 HIS 0.002 0.000 HIS C 73 PHE 0.008 0.001 PHE C 442 TYR 0.012 0.001 TYR C 417 ARG 0.003 0.000 ARG B 326 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 156 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 146 time to evaluate : 1.224 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 10 outliers final: 5 residues processed: 151 average time/residue: 0.1647 time to fit residues: 38.8401 Evaluate side-chains 140 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 135 time to evaluate : 1.166 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.1009 time to fit residues: 2.4013 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 126 random chunks: chunk 72 optimal weight: 4.9990 chunk 36 optimal weight: 0.6980 chunk 23 optimal weight: 5.9990 chunk 77 optimal weight: 4.9990 chunk 82 optimal weight: 9.9990 chunk 60 optimal weight: 1.9990 chunk 11 optimal weight: 0.0060 chunk 95 optimal weight: 0.2980 chunk 110 optimal weight: 1.9990 chunk 116 optimal weight: 2.9990 chunk 106 optimal weight: 0.6980 overall best weight: 0.7398 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 322 ASN A 438 GLN C 438 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8119 moved from start: 0.3281 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.032 10209 Z= 0.154 Angle : 0.567 16.338 13882 Z= 0.263 Chirality : 0.040 0.186 1557 Planarity : 0.002 0.032 1789 Dihedral : 4.959 55.956 1721 Min Nonbonded Distance : 2.289 Molprobity Statistics. All-atom Clashscore : 4.71 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.57 % Favored : 98.43 % Rotamer: Outliers : 0.81 % Allowed : 24.90 % Favored : 74.29 % 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: -1.82 (0.25), residues: 1272 helix: -0.41 (0.27), residues: 414 sheet: -2.67 (0.31), residues: 258 loop : -0.93 (0.28), residues: 600 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP C 233 HIS 0.002 0.000 HIS C 73 PHE 0.008 0.001 PHE C 442 TYR 0.012 0.001 TYR A 417 ARG 0.002 0.000 ARG B 326 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 153 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 145 time to evaluate : 1.186 Fit side-chains revert: symmetry clash outliers start: 8 outliers final: 4 residues processed: 149 average time/residue: 0.1652 time to fit residues: 38.5088 Evaluate side-chains 140 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 136 time to evaluate : 1.040 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.0889 time to fit residues: 2.0578 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 126 random chunks: chunk 113 optimal weight: 0.9980 chunk 116 optimal weight: 3.9990 chunk 68 optimal weight: 0.7980 chunk 49 optimal weight: 10.0000 chunk 88 optimal weight: 10.0000 chunk 34 optimal weight: 0.9980 chunk 102 optimal weight: 0.0570 chunk 107 optimal weight: 0.6980 chunk 112 optimal weight: 8.9990 chunk 74 optimal weight: 1.9990 chunk 119 optimal weight: 0.9980 overall best weight: 0.7098 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 322 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8111 moved from start: 0.3360 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.033 10209 Z= 0.151 Angle : 0.560 15.535 13882 Z= 0.263 Chirality : 0.040 0.243 1557 Planarity : 0.003 0.053 1789 Dihedral : 4.963 55.716 1721 Min Nonbonded Distance : 2.272 Molprobity Statistics. All-atom Clashscore : 5.07 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.97 % Favored : 98.03 % Rotamer: Outliers : 0.30 % Allowed : 25.10 % Favored : 74.60 % 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: -1.71 (0.25), residues: 1272 helix: -0.30 (0.27), residues: 414 sheet: -2.57 (0.31), residues: 258 loop : -0.91 (0.28), residues: 600 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.004 0.001 TRP C 233 HIS 0.002 0.000 HIS C 73 PHE 0.007 0.001 PHE C 98 TYR 0.025 0.001 TYR A 318 ARG 0.002 0.000 ARG B 326 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 147 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 144 time to evaluate : 1.130 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 3 outliers final: 1 residues processed: 147 average time/residue: 0.1616 time to fit residues: 37.3697 Evaluate side-chains 136 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 135 time to evaluate : 1.094 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1329 time to fit residues: 1.7304 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 126 random chunks: chunk 73 optimal weight: 0.5980 chunk 56 optimal weight: 3.9990 chunk 83 optimal weight: 0.9980 chunk 125 optimal weight: 0.5980 chunk 115 optimal weight: 10.0000 chunk 99 optimal weight: 3.9990 chunk 10 optimal weight: 0.0570 chunk 77 optimal weight: 0.9980 chunk 61 optimal weight: 4.9990 chunk 79 optimal weight: 5.9990 chunk 106 optimal weight: 0.8980 overall best weight: 0.6298 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 110 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8096 moved from start: 0.3507 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.035 10209 Z= 0.147 Angle : 0.562 16.416 13882 Z= 0.261 Chirality : 0.040 0.263 1557 Planarity : 0.003 0.048 1789 Dihedral : 4.843 55.374 1721 Min Nonbonded Distance : 2.270 Molprobity Statistics. All-atom Clashscore : 5.23 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.89 % Favored : 98.11 % Rotamer: Outliers : 0.00 % Allowed : 26.41 % Favored : 73.59 % 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: -1.66 (0.25), residues: 1272 helix: -0.17 (0.27), residues: 414 sheet: -2.67 (0.29), residues: 279 loop : -0.78 (0.28), residues: 579 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP C 233 HIS 0.002 0.000 HIS A 329 PHE 0.007 0.001 PHE C 98 TYR 0.011 0.001 TYR C 417 ARG 0.002 0.000 ARG B 326 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2544 Ramachandran restraints generated. 1272 Oldfield, 0 Emsley, 1272 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 139 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 139 time to evaluate : 1.216 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 139 average time/residue: 0.1672 time to fit residues: 36.1317 Evaluate side-chains 129 residues out of total 1119 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 129 time to evaluate : 1.203 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.5636 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 126 random chunks: chunk 30 optimal weight: 0.4980 chunk 92 optimal weight: 4.9990 chunk 14 optimal weight: 3.9990 chunk 27 optimal weight: 4.9990 chunk 100 optimal weight: 0.9980 chunk 41 optimal weight: 0.7980 chunk 102 optimal weight: 4.9990 chunk 12 optimal weight: 8.9990 chunk 18 optimal weight: 3.9990 chunk 87 optimal weight: 1.9990 chunk 5 optimal weight: 2.9990 overall best weight: 1.4584 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 329 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3865 r_free = 0.3865 target = 0.161944 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 57)----------------| | r_work = 0.3314 r_free = 0.3314 target = 0.116592 restraints weight = 13069.929| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 45)----------------| | r_work = 0.3359 r_free = 0.3359 target = 0.120230 restraints weight = 7404.795| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 41)----------------| | r_work = 0.3386 r_free = 0.3386 target = 0.122378 restraints weight = 5495.382| |-----------------------------------------------------------------------------| r_work (final): 0.3360 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8106 moved from start: 0.3424 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.038 10209 Z= 0.222 Angle : 0.591 15.946 13882 Z= 0.279 Chirality : 0.042 0.323 1557 Planarity : 0.003 0.056 1789 Dihedral : 4.973 55.018 1721 Min Nonbonded Distance : 2.236 Molprobity Statistics. All-atom Clashscore : 6.21 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.44 % Favored : 97.56 % Rotamer: Outliers : 0.60 % Allowed : 26.01 % Favored : 73.39 % 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: -1.62 (0.25), residues: 1272 helix: -0.09 (0.27), residues: 411 sheet: -2.72 (0.29), residues: 279 loop : -0.75 (0.28), residues: 582 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP C 233 HIS 0.006 0.001 HIS B 329 PHE 0.008 0.001 PHE C 442 TYR 0.013 0.001 TYR A 318 ARG 0.002 0.000 ARG B 326 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 1983.13 seconds wall clock time: 36 minutes 49.26 seconds (2209.26 seconds total)