Starting phenix.real_space_refine (version: 1.21rc1) on Fri Apr 21 10:15:27 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7vr1_32096/04_2023/7vr1_32096.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7vr1_32096/04_2023/7vr1_32096.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.4 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7vr1_32096/04_2023/7vr1_32096.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7vr1_32096/04_2023/7vr1_32096.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7vr1_32096/04_2023/7vr1_32096.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7vr1_32096/04_2023/7vr1_32096.pdb" } resolution = 3.4 write_initial_geo_file = False refinement { macro_cycles = 10 } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= 0.001 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: 4 Type Number sf(0) Gaussians S 42 5.16 5 C 6110 2.51 5 N 1668 2.21 5 O 1702 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "B PHE 175": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 281": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B GLU 380": "OE1" <-> "OE2" Residue "B PHE 517": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 532": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 106": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 427": "OE1" <-> "OE2" Residue "A PHE 517": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-4938/modules/chem_data/mon_lib" Total number of atoms: 9522 Number of models: 1 Model: "" Number of chains: 2 Chain: "B" Number of atoms: 4761 Number of conformers: 1 Conformer: "" Number of residues, atoms: 598, 4761 Classifications: {'peptide': 598} Link IDs: {'PTRANS': 21, 'TRANS': 576} Chain breaks: 2 Chain: "A" Number of atoms: 4761 Number of conformers: 1 Conformer: "" Number of residues, atoms: 598, 4761 Classifications: {'peptide': 598} Link IDs: {'PTRANS': 21, 'TRANS': 576} Chain breaks: 2 Time building chain proxies: 4.88, per 1000 atoms: 0.51 Number of scatterers: 9522 At special positions: 0 Unit cell: (114.182, 73.5593, 146.021, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 4 Type Number sf(0) S 42 16.00 O 1702 8.00 N 1668 7.00 C 6110 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: 3.25 Conformation dependent library (CDL) restraints added in 1.3 seconds 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2244 Finding SS restraints... Secondary structure from input PDB file: 46 helices and 5 sheets defined 67.2% alpha, 2.0% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.42 Creating SS restraints... Processing helix chain 'B' and resid 70 through 83 removed outlier: 3.784A pdb=" N LEU B 81 " --> pdb=" O TRP B 77 " (cutoff:3.500A) Processing helix chain 'B' and resid 88 through 100 removed outlier: 3.516A pdb=" N HIS B 97 " --> pdb=" O LEU B 93 " (cutoff:3.500A) Processing helix chain 'B' and resid 101 through 123 removed outlier: 3.641A pdb=" N PHE B 106 " --> pdb=" O VAL B 102 " (cutoff:3.500A) removed outlier: 4.003A pdb=" N LEU B 107 " --> pdb=" O SER B 103 " (cutoff:3.500A) removed outlier: 3.611A pdb=" N TYR B 110 " --> pdb=" O PHE B 106 " (cutoff:3.500A) removed outlier: 3.544A pdb=" N ARG B 117 " --> pdb=" O ARG B 113 " (cutoff:3.500A) removed outlier: 3.589A pdb=" N VAL B 123 " --> pdb=" O ALA B 119 " (cutoff:3.500A) Processing helix chain 'B' and resid 126 through 140 Processing helix chain 'B' and resid 140 through 175 removed outlier: 3.546A pdb=" N SER B 149 " --> pdb=" O THR B 145 " (cutoff:3.500A) removed outlier: 3.562A pdb=" N ARG B 152 " --> pdb=" O ASN B 148 " (cutoff:3.500A) removed outlier: 3.673A pdb=" N ALA B 168 " --> pdb=" O SER B 164 " (cutoff:3.500A) removed outlier: 3.726A pdb=" N ARG B 172 " --> pdb=" O ALA B 168 " (cutoff:3.500A) removed outlier: 4.001A pdb=" N LEU B 173 " --> pdb=" O HIS B 169 " (cutoff:3.500A) Processing helix chain 'B' and resid 178 through 186 removed outlier: 3.887A pdb=" N SER B 184 " --> pdb=" O TYR B 180 " (cutoff:3.500A) Processing helix chain 'B' and resid 192 through 216 removed outlier: 3.632A pdb=" N SER B 196 " --> pdb=" O ASN B 192 " (cutoff:3.500A) removed outlier: 4.096A pdb=" N ASP B 200 " --> pdb=" O SER B 196 " (cutoff:3.500A) removed outlier: 4.965A pdb=" N VAL B 201 " --> pdb=" O LEU B 197 " (cutoff:3.500A) removed outlier: 3.912A pdb=" N VAL B 202 " --> pdb=" O THR B 198 " (cutoff:3.500A) removed outlier: 3.510A pdb=" N ALA B 203 " --> pdb=" O GLU B 199 " (cutoff:3.500A) removed outlier: 3.609A pdb=" N ALA B 209 " --> pdb=" O ALA B 205 " (cutoff:3.500A) removed outlier: 3.715A pdb=" N HIS B 210 " --> pdb=" O ALA B 206 " (cutoff:3.500A) Processing helix chain 'B' and resid 216 through 235 removed outlier: 3.782A pdb=" N LEU B 230 " --> pdb=" O SER B 226 " (cutoff:3.500A) Processing helix chain 'B' and resid 241 through 261 removed outlier: 3.859A pdb=" N ALA B 245 " --> pdb=" O ALA B 241 " (cutoff:3.500A) removed outlier: 3.576A pdb=" N ILE B 246 " --> pdb=" O TRP B 242 " (cutoff:3.500A) removed outlier: 3.596A pdb=" N LEU B 258 " --> pdb=" O THR B 254 " (cutoff:3.500A) removed outlier: 3.815A pdb=" N PHE B 261 " --> pdb=" O VAL B 257 " (cutoff:3.500A) Processing helix chain 'B' and resid 261 through 289 removed outlier: 3.505A pdb=" N PHE B 265 " --> pdb=" O PHE B 261 " (cutoff:3.500A) removed outlier: 4.080A pdb=" N VAL B 286 " --> pdb=" O MET B 282 " (cutoff:3.500A) removed outlier: 3.930A pdb=" N VAL B 287 " --> pdb=" O HIS B 283 " (cutoff:3.500A) Processing helix chain 'B' and resid 292 through 297 Processing helix chain 'B' and resid 298 through 337 removed outlier: 3.561A pdb=" N SER B 309 " --> pdb=" O LEU B 305 " (cutoff:3.500A) removed outlier: 3.809A pdb=" N LEU B 313 " --> pdb=" O SER B 309 " (cutoff:3.500A) removed outlier: 3.570A pdb=" N LEU B 319 " --> pdb=" O SER B 315 " (cutoff:3.500A) removed outlier: 4.179A pdb=" N TRP B 326 " --> pdb=" O LEU B 322 " (cutoff:3.500A) removed outlier: 3.606A pdb=" N TYR B 327 " --> pdb=" O GLU B 323 " (cutoff:3.500A) removed outlier: 4.032A pdb=" N MET B 329 " --> pdb=" O LEU B 325 " (cutoff:3.500A) removed outlier: 3.558A pdb=" N LEU B 334 " --> pdb=" O LEU B 330 " (cutoff:3.500A) Processing helix chain 'B' and resid 338 through 353 removed outlier: 3.654A pdb=" N MET B 346 " --> pdb=" O SER B 342 " (cutoff:3.500A) removed outlier: 4.140A pdb=" N VAL B 347 " --> pdb=" O GLY B 343 " (cutoff:3.500A) removed outlier: 3.629A pdb=" N VAL B 349 " --> pdb=" O LEU B 345 " (cutoff:3.500A) Proline residue: B 350 - end of helix Processing helix chain 'B' and resid 361 through 371 removed outlier: 3.719A pdb=" N LYS B 367 " --> pdb=" O GLU B 363 " (cutoff:3.500A) Processing helix chain 'B' and resid 377 through 407 removed outlier: 3.623A pdb=" N THR B 382 " --> pdb=" O VAL B 378 " (cutoff:3.500A) removed outlier: 3.694A pdb=" N GLU B 383 " --> pdb=" O SER B 379 " (cutoff:3.500A) removed outlier: 3.655A pdb=" N THR B 386 " --> pdb=" O THR B 382 " (cutoff:3.500A) removed outlier: 3.532A pdb=" N ALA B 388 " --> pdb=" O ALA B 384 " (cutoff:3.500A) removed outlier: 3.894A pdb=" N LEU B 392 " --> pdb=" O ALA B 388 " (cutoff:3.500A) Processing helix chain 'B' and resid 408 through 431 removed outlier: 3.781A pdb=" N GLY B 414 " --> pdb=" O THR B 410 " (cutoff:3.500A) removed outlier: 4.101A pdb=" N MET B 422 " --> pdb=" O ARG B 418 " (cutoff:3.500A) removed outlier: 3.684A pdb=" N GLU B 427 " --> pdb=" O PHE B 423 " (cutoff:3.500A) Processing helix chain 'B' and resid 535 through 539 Processing helix chain 'B' and resid 553 through 558 Processing helix chain 'B' and resid 575 through 584 removed outlier: 3.747A pdb=" N ILE B 579 " --> pdb=" O ASP B 575 " (cutoff:3.500A) Processing helix chain 'B' and resid 585 through 592 removed outlier: 3.922A pdb=" N LEU B 589 " --> pdb=" O LEU B 585 " (cutoff:3.500A) removed outlier: 3.830A pdb=" N GLN B 590 " --> pdb=" O HIS B 586 " (cutoff:3.500A) removed outlier: 3.613A pdb=" N ARG B 591 " --> pdb=" O HIS B 587 " (cutoff:3.500A) Processing helix chain 'B' and resid 600 through 605 removed outlier: 3.710A pdb=" N VAL B 604 " --> pdb=" O ASP B 600 " (cutoff:3.500A) removed outlier: 3.596A pdb=" N LEU B 605 " --> pdb=" O TRP B 601 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 600 through 605' Processing helix chain 'B' and resid 606 through 622 removed outlier: 3.577A pdb=" N GLN B 611 " --> pdb=" O GLY B 607 " (cutoff:3.500A) removed outlier: 3.618A pdb=" N MET B 618 " --> pdb=" O GLY B 614 " (cutoff:3.500A) removed outlier: 3.555A pdb=" N PHE B 619 " --> pdb=" O MET B 615 " (cutoff:3.500A) removed outlier: 3.787A pdb=" N TYR B 620 " --> pdb=" O ALA B 616 " (cutoff:3.500A) removed outlier: 3.599A pdb=" N HIS B 621 " --> pdb=" O ARG B 617 " (cutoff:3.500A) removed outlier: 3.542A pdb=" N ARG B 622 " --> pdb=" O MET B 618 " (cutoff:3.500A) Processing helix chain 'B' and resid 636 through 651 removed outlier: 3.739A pdb=" N LYS B 642 " --> pdb=" O ASP B 638 " (cutoff:3.500A) Processing helix chain 'B' and resid 660 through 666 removed outlier: 3.787A pdb=" N LYS B 665 " --> pdb=" O SER B 662 " (cutoff:3.500A) Processing helix chain 'B' and resid 685 through 688 removed outlier: 4.328A pdb=" N ALA B 688 " --> pdb=" O ASP B 685 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 685 through 688' Processing helix chain 'B' and resid 689 through 695 removed outlier: 3.572A pdb=" N THR B 693 " --> pdb=" O ARG B 689 " (cutoff:3.500A) Processing helix chain 'A' and resid 70 through 83 removed outlier: 3.514A pdb=" N LEU A 81 " --> pdb=" O TRP A 77 " (cutoff:3.500A) Processing helix chain 'A' and resid 88 through 123 removed outlier: 3.622A pdb=" N GLY A 92 " --> pdb=" O CYS A 88 " (cutoff:3.500A) removed outlier: 3.658A pdb=" N VAL A 123 " --> pdb=" O ALA A 119 " (cutoff:3.500A) Processing helix chain 'A' and resid 126 through 175 removed outlier: 3.612A pdb=" N GLN A 136 " --> pdb=" O TRP A 132 " (cutoff:3.500A) removed outlier: 3.986A pdb=" N LEU A 142 " --> pdb=" O LEU A 138 " (cutoff:3.500A) Proline residue: A 143 - end of helix removed outlier: 3.523A pdb=" N ARG A 152 " --> pdb=" O ASN A 148 " (cutoff:3.500A) removed outlier: 3.840A pdb=" N ARG A 172 " --> pdb=" O ALA A 168 " (cutoff:3.500A) removed outlier: 3.952A pdb=" N LEU A 173 " --> pdb=" O HIS A 169 " (cutoff:3.500A) Processing helix chain 'A' and resid 178 through 186 removed outlier: 3.815A pdb=" N SER A 184 " --> pdb=" O TYR A 180 " (cutoff:3.500A) Processing helix chain 'A' and resid 192 through 199 removed outlier: 3.603A pdb=" N SER A 196 " --> pdb=" O ASN A 192 " (cutoff:3.500A) Processing helix chain 'A' and resid 199 through 235 removed outlier: 3.605A pdb=" N ALA A 209 " --> pdb=" O ALA A 205 " (cutoff:3.500A) removed outlier: 3.600A pdb=" N HIS A 210 " --> pdb=" O ALA A 206 " (cutoff:3.500A) removed outlier: 3.786A pdb=" N LYS A 217 " --> pdb=" O SER A 213 " (cutoff:3.500A) Proline residue: A 218 - end of helix Processing helix chain 'A' and resid 241 through 260 removed outlier: 4.195A pdb=" N ALA A 245 " --> pdb=" O ALA A 241 " (cutoff:3.500A) removed outlier: 3.812A pdb=" N ILE A 246 " --> pdb=" O TRP A 242 " (cutoff:3.500A) Processing helix chain 'A' and resid 261 through 289 removed outlier: 3.705A pdb=" N PHE A 265 " --> pdb=" O PHE A 261 " (cutoff:3.500A) removed outlier: 3.730A pdb=" N VAL A 286 " --> pdb=" O MET A 282 " (cutoff:3.500A) removed outlier: 3.762A pdb=" N VAL A 287 " --> pdb=" O HIS A 283 " (cutoff:3.500A) Processing helix chain 'A' and resid 292 through 297 Processing helix chain 'A' and resid 298 through 337 removed outlier: 3.551A pdb=" N GLU A 302 " --> pdb=" O GLY A 298 " (cutoff:3.500A) removed outlier: 4.216A pdb=" N TRP A 326 " --> pdb=" O LEU A 322 " (cutoff:3.500A) removed outlier: 3.666A pdb=" N TYR A 327 " --> pdb=" O GLU A 323 " (cutoff:3.500A) removed outlier: 3.851A pdb=" N MET A 329 " --> pdb=" O LEU A 325 " (cutoff:3.500A) Processing helix chain 'A' and resid 338 through 353 removed outlier: 3.784A pdb=" N VAL A 347 " --> pdb=" O GLY A 343 " (cutoff:3.500A) removed outlier: 3.664A pdb=" N VAL A 349 " --> pdb=" O LEU A 345 " (cutoff:3.500A) Proline residue: A 350 - end of helix Processing helix chain 'A' and resid 363 through 372 removed outlier: 3.809A pdb=" N LYS A 367 " --> pdb=" O GLU A 363 " (cutoff:3.500A) removed outlier: 3.666A pdb=" N LEU A 370 " --> pdb=" O LYS A 366 " (cutoff:3.500A) removed outlier: 3.967A pdb=" N LYS A 372 " --> pdb=" O ALA A 368 " (cutoff:3.500A) Processing helix chain 'A' and resid 377 through 406 removed outlier: 3.576A pdb=" N ARG A 381 " --> pdb=" O LEU A 377 " (cutoff:3.500A) removed outlier: 3.580A pdb=" N THR A 382 " --> pdb=" O VAL A 378 " (cutoff:3.500A) removed outlier: 4.049A pdb=" N GLU A 383 " --> pdb=" O SER A 379 " (cutoff:3.500A) removed outlier: 3.751A pdb=" N THR A 386 " --> pdb=" O THR A 382 " (cutoff:3.500A) removed outlier: 3.627A pdb=" N ALA A 388 " --> pdb=" O ALA A 384 " (cutoff:3.500A) Processing helix chain 'A' and resid 408 through 431 removed outlier: 3.522A pdb=" N ARG A 431 " --> pdb=" O GLU A 427 " (cutoff:3.500A) Processing helix chain 'A' and resid 535 through 537 No H-bonds generated for 'chain 'A' and resid 535 through 537' Processing helix chain 'A' and resid 553 through 558 Processing helix chain 'A' and resid 575 through 584 removed outlier: 3.874A pdb=" N ILE A 579 " --> pdb=" O ASP A 575 " (cutoff:3.500A) Processing helix chain 'A' and resid 606 through 622 removed outlier: 3.581A pdb=" N GLN A 611 " --> pdb=" O GLY A 607 " (cutoff:3.500A) removed outlier: 3.727A pdb=" N MET A 615 " --> pdb=" O GLN A 611 " (cutoff:3.500A) Processing helix chain 'A' and resid 636 through 650 removed outlier: 3.898A pdb=" N LYS A 642 " --> pdb=" O ASP A 638 " (cutoff:3.500A) Processing helix chain 'A' and resid 660 through 664 Processing sheet with id=AA1, first strand: chain 'B' and resid 474 through 475 Processing sheet with id=AA2, first strand: chain 'B' and resid 504 through 505 removed outlier: 6.331A pdb=" N ALA B 626 " --> pdb=" O LEU B 655 " (cutoff:3.500A) removed outlier: 6.874A pdb=" N ILE B 657 " --> pdb=" O ALA B 626 " (cutoff:3.500A) removed outlier: 6.059A pdb=" N LEU B 628 " --> pdb=" O ILE B 657 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'B' and resid 672 through 673 Processing sheet with id=AA4, first strand: chain 'A' and resid 474 through 475 Processing sheet with id=AA5, first strand: chain 'A' and resid 539 through 542 removed outlier: 6.262A pdb=" N PHE A 540 " --> pdb=" O LEU A 627 " (cutoff:3.500A) 524 hydrogen bonds defined for protein. 1542 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.38 Time building geometry restraints manager: 3.65 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: 3110 1.34 - 1.46: 1706 1.46 - 1.58: 4840 1.58 - 1.70: 0 1.70 - 1.81: 70 Bond restraints: 9726 Sorted by residual: bond pdb=" C ARG B 622 " pdb=" O ARG B 622 " ideal model delta sigma weight residual 1.235 1.240 -0.005 4.70e-03 4.53e+04 1.04e+00 bond pdb=" CA ILE A 558 " pdb=" CB ILE A 558 " ideal model delta sigma weight residual 1.540 1.527 0.013 1.36e-02 5.41e+03 9.55e-01 bond pdb=" C GLN B 472 " pdb=" N GLY B 473 " ideal model delta sigma weight residual 1.322 1.333 -0.011 1.20e-02 6.94e+03 8.51e-01 bond pdb=" CG GLU B 471 " pdb=" CD GLU B 471 " ideal model delta sigma weight residual 1.516 1.496 0.020 2.50e-02 1.60e+03 6.15e-01 bond pdb=" CB MET A 548 " pdb=" CG MET A 548 " ideal model delta sigma weight residual 1.520 1.542 -0.022 3.00e-02 1.11e+03 5.49e-01 ... (remaining 9721 not shown) Histogram of bond angle deviations from ideal: 99.81 - 106.66: 246 106.66 - 113.51: 5335 113.51 - 120.36: 3709 120.36 - 127.21: 3772 127.21 - 134.06: 112 Bond angle restraints: 13174 Sorted by residual: angle pdb=" C GLN B 177 " pdb=" CA GLN B 177 " pdb=" CB GLN B 177 " ideal model delta sigma weight residual 110.42 116.01 -5.59 1.99e+00 2.53e-01 7.89e+00 angle pdb=" CB MET A 548 " pdb=" CG MET A 548 " pdb=" SD MET A 548 " ideal model delta sigma weight residual 112.70 120.51 -7.81 3.00e+00 1.11e-01 6.78e+00 angle pdb=" CB MET B 548 " pdb=" CG MET B 548 " pdb=" SD MET B 548 " ideal model delta sigma weight residual 112.70 119.83 -7.13 3.00e+00 1.11e-01 5.66e+00 angle pdb=" N GLU B 291 " pdb=" CA GLU B 291 " pdb=" C GLU B 291 " ideal model delta sigma weight residual 113.72 110.23 3.49 1.52e+00 4.33e-01 5.28e+00 angle pdb=" N GLY A 500 " pdb=" CA GLY A 500 " pdb=" C GLY A 500 " ideal model delta sigma weight residual 113.18 118.54 -5.36 2.37e+00 1.78e-01 5.12e+00 ... (remaining 13169 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 16.27: 5139 16.27 - 32.55: 526 32.55 - 48.82: 136 48.82 - 65.09: 10 65.09 - 81.36: 9 Dihedral angle restraints: 5820 sinusoidal: 2338 harmonic: 3482 Sorted by residual: dihedral pdb=" CA GLU B 471 " pdb=" C GLU B 471 " pdb=" N GLN B 472 " pdb=" CA GLN B 472 " ideal model delta harmonic sigma weight residual 180.00 150.21 29.79 0 5.00e+00 4.00e-02 3.55e+01 dihedral pdb=" CA TRP A 679 " pdb=" C TRP A 679 " pdb=" N LYS A 680 " pdb=" CA LYS A 680 " ideal model delta harmonic sigma weight residual 180.00 162.20 17.80 0 5.00e+00 4.00e-02 1.27e+01 dihedral pdb=" CA GLU A 499 " pdb=" C GLU A 499 " pdb=" N GLY A 500 " pdb=" CA GLY A 500 " ideal model delta harmonic sigma weight residual 180.00 -162.27 -17.73 0 5.00e+00 4.00e-02 1.26e+01 ... (remaining 5817 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.057: 1279 0.057 - 0.113: 191 0.113 - 0.170: 11 0.170 - 0.226: 0 0.226 - 0.283: 1 Chirality restraints: 1482 Sorted by residual: chirality pdb=" CB ILE B 558 " pdb=" CA ILE B 558 " pdb=" CG1 ILE B 558 " pdb=" CG2 ILE B 558 " both_signs ideal model delta sigma weight residual False 2.64 2.36 0.28 2.00e-01 2.50e+01 2.00e+00 chirality pdb=" CA ILE B 657 " pdb=" N ILE B 657 " pdb=" C ILE B 657 " pdb=" CB ILE B 657 " both_signs ideal model delta sigma weight residual False 2.43 2.58 -0.14 2.00e-01 2.50e+01 5.23e-01 chirality pdb=" CA PHE A 517 " pdb=" N PHE A 517 " pdb=" C PHE A 517 " pdb=" CB PHE A 517 " both_signs ideal model delta sigma weight residual False 2.51 2.37 0.14 2.00e-01 2.50e+01 4.98e-01 ... (remaining 1479 not shown) Planarity restraints: 1670 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA LEU B 523 " -0.006 2.00e-02 2.50e+03 1.23e-02 1.51e+00 pdb=" C LEU B 523 " 0.021 2.00e-02 2.50e+03 pdb=" O LEU B 523 " -0.008 2.00e-02 2.50e+03 pdb=" N TRP B 524 " -0.007 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA PHE A 333 " 0.006 2.00e-02 2.50e+03 1.20e-02 1.45e+00 pdb=" C PHE A 333 " -0.021 2.00e-02 2.50e+03 pdb=" O PHE A 333 " 0.008 2.00e-02 2.50e+03 pdb=" N LEU A 334 " 0.007 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C LYS B 533 " -0.020 5.00e-02 4.00e+02 2.96e-02 1.40e+00 pdb=" N PRO B 534 " 0.051 5.00e-02 4.00e+02 pdb=" CA PRO B 534 " -0.015 5.00e-02 4.00e+02 pdb=" CD PRO B 534 " -0.016 5.00e-02 4.00e+02 ... (remaining 1667 not shown) Histogram of nonbonded interaction distances: 2.25 - 2.78: 2114 2.78 - 3.31: 9280 3.31 - 3.84: 15127 3.84 - 4.37: 17111 4.37 - 4.90: 30172 Nonbonded interactions: 73804 Sorted by model distance: nonbonded pdb=" OE2 GLU B 90 " pdb=" OG SER B 161 " model vdw 2.253 2.440 nonbonded pdb=" OD1 ASP B 200 " pdb=" NH1 ARG B 418 " model vdw 2.271 2.520 nonbonded pdb=" NH2 ARG B 285 " pdb=" O MET A 548 " model vdw 2.303 2.520 nonbonded pdb=" O THR A 632 " pdb=" NH1 ARG A 660 " model vdw 2.307 2.520 nonbonded pdb=" OE2 GLU A 90 " pdb=" OG SER A 161 " model vdw 2.311 2.440 ... (remaining 73799 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = chain 'A' selection = chain 'B' } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.650 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.550 Check model and map are aligned: 0.140 Set scattering table: 0.080 Process input model: 24.820 Find NCS groups from input model: 0.470 Set up NCS constraints: 0.050 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.610 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 31.380 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8123 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.003 0.023 9726 Z= 0.173 Angle : 0.549 7.814 13174 Z= 0.276 Chirality : 0.039 0.283 1482 Planarity : 0.003 0.030 1670 Dihedral : 13.992 81.363 3576 Min Nonbonded Distance : 2.253 Molprobity Statistics. All-atom Clashscore : 7.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.83 % Favored : 94.17 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.76 (0.26), residues: 1184 helix: 0.49 (0.20), residues: 742 sheet: -2.64 (0.56), residues: 86 loop : -1.95 (0.32), residues: 356 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 108 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 108 time to evaluate : 1.160 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 108 average time/residue: 0.1915 time to fit residues: 30.9170 Evaluate side-chains 83 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 83 time to evaluate : 1.143 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.5647 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 99 optimal weight: 20.0000 chunk 89 optimal weight: 0.9980 chunk 49 optimal weight: 7.9990 chunk 30 optimal weight: 5.9990 chunk 60 optimal weight: 3.9990 chunk 47 optimal weight: 0.3980 chunk 92 optimal weight: 0.8980 chunk 35 optimal weight: 5.9990 chunk 56 optimal weight: 10.0000 chunk 68 optimal weight: 3.9990 chunk 107 optimal weight: 8.9990 overall best weight: 2.0584 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 587 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8140 moved from start: 0.0987 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.029 9726 Z= 0.260 Angle : 0.585 7.354 13174 Z= 0.297 Chirality : 0.041 0.190 1482 Planarity : 0.004 0.030 1670 Dihedral : 4.096 18.994 1330 Min Nonbonded Distance : 2.182 Molprobity Statistics. All-atom Clashscore : 8.77 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.17 % Favored : 93.83 % Rotamer Outliers : 0.40 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.40 (0.25), residues: 1184 helix: 0.80 (0.20), residues: 760 sheet: -2.57 (0.56), residues: 86 loop : -2.00 (0.33), residues: 338 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 100 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 96 time to evaluate : 1.144 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 4 outliers final: 2 residues processed: 98 average time/residue: 0.1959 time to fit residues: 28.6051 Evaluate side-chains 85 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 83 time to evaluate : 1.155 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.2347 time to fit residues: 2.1883 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 59 optimal weight: 7.9990 chunk 33 optimal weight: 2.9990 chunk 89 optimal weight: 1.9990 chunk 72 optimal weight: 0.9980 chunk 29 optimal weight: 7.9990 chunk 107 optimal weight: 6.9990 chunk 115 optimal weight: 9.9990 chunk 95 optimal weight: 9.9990 chunk 106 optimal weight: 0.5980 chunk 36 optimal weight: 7.9990 chunk 86 optimal weight: 0.9990 overall best weight: 1.5186 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 556 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 587 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8123 moved from start: 0.1265 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.027 9726 Z= 0.205 Angle : 0.542 7.868 13174 Z= 0.270 Chirality : 0.039 0.184 1482 Planarity : 0.003 0.028 1670 Dihedral : 3.969 17.383 1330 Min Nonbonded Distance : 2.195 Molprobity Statistics. All-atom Clashscore : 7.88 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.25 % Favored : 93.75 % Rotamer Outliers : 0.89 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.07 (0.26), residues: 1184 helix: 1.09 (0.20), residues: 760 sheet: -2.50 (0.64), residues: 74 loop : -1.98 (0.32), residues: 350 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 105 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 96 time to evaluate : 1.057 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 9 outliers final: 3 residues processed: 101 average time/residue: 0.1892 time to fit residues: 28.6542 Evaluate side-chains 88 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 85 time to evaluate : 1.211 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0884 time to fit residues: 2.1146 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 106 optimal weight: 0.4980 chunk 80 optimal weight: 2.9990 chunk 55 optimal weight: 10.0000 chunk 11 optimal weight: 0.7980 chunk 51 optimal weight: 5.9990 chunk 72 optimal weight: 5.9990 chunk 107 optimal weight: 7.9990 chunk 114 optimal weight: 6.9990 chunk 56 optimal weight: 6.9990 chunk 102 optimal weight: 8.9990 chunk 30 optimal weight: 10.0000 overall best weight: 3.2586 ------------------------------------------------------------------------------- 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.8182 moved from start: 0.1570 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.006 0.052 9726 Z= 0.378 Angle : 0.618 8.650 13174 Z= 0.311 Chirality : 0.044 0.479 1482 Planarity : 0.004 0.031 1670 Dihedral : 4.199 18.064 1330 Min Nonbonded Distance : 2.146 Molprobity Statistics. All-atom Clashscore : 9.35 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.76 % Favored : 93.24 % Rotamer Outliers : 0.99 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.24 (0.25), residues: 1184 helix: 1.00 (0.20), residues: 760 sheet: -2.59 (0.57), residues: 86 loop : -2.11 (0.32), residues: 338 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 97 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 87 time to evaluate : 1.088 Fit side-chains revert: symmetry clash outliers start: 10 outliers final: 6 residues processed: 94 average time/residue: 0.1830 time to fit residues: 26.7814 Evaluate side-chains 84 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 78 time to evaluate : 1.185 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 0 residues processed: 6 average time/residue: 0.1033 time to fit residues: 2.6723 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 94 optimal weight: 0.8980 chunk 64 optimal weight: 0.7980 chunk 1 optimal weight: 3.9990 chunk 84 optimal weight: 0.9990 chunk 47 optimal weight: 6.9990 chunk 97 optimal weight: 3.9990 chunk 78 optimal weight: 0.9980 chunk 0 optimal weight: 6.9990 chunk 58 optimal weight: 20.0000 chunk 102 optimal weight: 8.9990 chunk 28 optimal weight: 0.0470 overall best weight: 0.7480 ------------------------------------------------------------------------------- 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.8107 moved from start: 0.1602 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.043 9726 Z= 0.155 Angle : 0.540 8.289 13174 Z= 0.266 Chirality : 0.039 0.366 1482 Planarity : 0.003 0.027 1670 Dihedral : 4.018 18.366 1330 Min Nonbonded Distance : 2.236 Molprobity Statistics. All-atom Clashscore : 7.94 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.00 % Favored : 94.00 % Rotamer Outliers : 1.59 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.13 (0.26), residues: 1184 helix: 1.31 (0.20), residues: 758 sheet: -2.45 (0.58), residues: 86 loop : -2.02 (0.33), residues: 340 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 105 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 16 poor density : 89 time to evaluate : 1.213 Fit side-chains outliers start: 16 outliers final: 5 residues processed: 98 average time/residue: 0.1764 time to fit residues: 27.1731 Evaluate side-chains 85 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 80 time to evaluate : 1.033 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.0900 time to fit residues: 2.2449 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 38 optimal weight: 10.0000 chunk 102 optimal weight: 6.9990 chunk 22 optimal weight: 1.9990 chunk 66 optimal weight: 2.9990 chunk 28 optimal weight: 2.9990 chunk 114 optimal weight: 10.0000 chunk 94 optimal weight: 5.9990 chunk 52 optimal weight: 5.9990 chunk 9 optimal weight: 4.9990 chunk 37 optimal weight: 4.9990 chunk 59 optimal weight: 7.9990 overall best weight: 3.5990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 430 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 587 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8191 moved from start: 0.1794 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.006 0.050 9726 Z= 0.410 Angle : 0.640 8.150 13174 Z= 0.322 Chirality : 0.044 0.337 1482 Planarity : 0.004 0.031 1670 Dihedral : 4.289 20.551 1330 Min Nonbonded Distance : 2.145 Molprobity Statistics. All-atom Clashscore : 10.13 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.26 % Favored : 92.74 % Rotamer Outliers : 0.99 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.11 (0.25), residues: 1184 helix: 1.04 (0.19), residues: 760 sheet: -1.48 (0.78), residues: 54 loop : -2.23 (0.30), residues: 370 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 91 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 81 time to evaluate : 1.152 Fit side-chains outliers start: 10 outliers final: 4 residues processed: 90 average time/residue: 0.1812 time to fit residues: 25.0530 Evaluate side-chains 80 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 76 time to evaluate : 1.140 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.0882 time to fit residues: 2.3046 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 110 optimal weight: 1.9990 chunk 12 optimal weight: 3.9990 chunk 65 optimal weight: 0.9980 chunk 83 optimal weight: 0.9980 chunk 64 optimal weight: 0.6980 chunk 96 optimal weight: 5.9990 chunk 63 optimal weight: 0.7980 chunk 113 optimal weight: 3.9990 chunk 71 optimal weight: 1.9990 chunk 69 optimal weight: 0.7980 chunk 52 optimal weight: 0.9980 overall best weight: 0.8580 ------------------------------------------------------------------------------- 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.8117 moved from start: 0.1791 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.036 9726 Z= 0.162 Angle : 0.553 9.193 13174 Z= 0.271 Chirality : 0.040 0.301 1482 Planarity : 0.003 0.027 1670 Dihedral : 4.063 18.450 1330 Min Nonbonded Distance : 2.242 Molprobity Statistics. All-atom Clashscore : 8.35 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.66 % Favored : 94.34 % Rotamer Outliers : 0.50 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.15 (0.26), residues: 1184 helix: 1.34 (0.20), residues: 760 sheet: -2.49 (0.57), residues: 88 loop : -2.03 (0.32), residues: 336 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 91 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 86 time to evaluate : 1.261 Fit side-chains outliers start: 5 outliers final: 4 residues processed: 89 average time/residue: 0.1806 time to fit residues: 24.9082 Evaluate side-chains 85 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 81 time to evaluate : 1.008 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.0853 time to fit residues: 2.0649 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 70 optimal weight: 0.9980 chunk 45 optimal weight: 3.9990 chunk 67 optimal weight: 0.7980 chunk 34 optimal weight: 0.9990 chunk 22 optimal weight: 0.9980 chunk 72 optimal weight: 0.8980 chunk 77 optimal weight: 1.9990 chunk 56 optimal weight: 0.9990 chunk 10 optimal weight: 2.9990 chunk 89 optimal weight: 0.9980 chunk 103 optimal weight: 0.7980 overall best weight: 0.8980 ------------------------------------------------------------------------------- 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.8109 moved from start: 0.1844 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.033 9726 Z= 0.165 Angle : 0.559 12.171 13174 Z= 0.270 Chirality : 0.040 0.285 1482 Planarity : 0.003 0.028 1670 Dihedral : 3.952 18.469 1330 Min Nonbonded Distance : 2.234 Molprobity Statistics. All-atom Clashscore : 8.56 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.74 % Favored : 94.26 % Rotamer Outliers : 0.30 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.42 (0.26), residues: 1184 helix: 1.47 (0.20), residues: 762 sheet: -1.46 (0.77), residues: 54 loop : -2.06 (0.31), residues: 368 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 91 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 88 time to evaluate : 1.159 Fit side-chains revert: symmetry clash outliers start: 3 outliers final: 2 residues processed: 90 average time/residue: 0.1828 time to fit residues: 25.5702 Evaluate side-chains 88 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 86 time to evaluate : 1.156 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.1077 time to fit residues: 1.9898 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 108 optimal weight: 0.3980 chunk 99 optimal weight: 5.9990 chunk 105 optimal weight: 0.3980 chunk 63 optimal weight: 0.9990 chunk 46 optimal weight: 0.0170 chunk 83 optimal weight: 1.9990 chunk 32 optimal weight: 0.6980 chunk 95 optimal weight: 9.9990 chunk 100 optimal weight: 7.9990 chunk 69 optimal weight: 4.9990 chunk 112 optimal weight: 5.9990 overall best weight: 0.5020 ------------------------------------------------------------------------------- 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.8075 moved from start: 0.1952 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.036 9726 Z= 0.139 Angle : 0.549 12.138 13174 Z= 0.265 Chirality : 0.039 0.274 1482 Planarity : 0.003 0.027 1670 Dihedral : 3.870 24.284 1330 Min Nonbonded Distance : 2.216 Molprobity Statistics. All-atom Clashscore : 7.88 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.66 % Favored : 94.34 % Rotamer Outliers : 0.10 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.66 (0.26), residues: 1184 helix: 1.63 (0.20), residues: 764 sheet: -0.57 (1.03), residues: 34 loop : -2.00 (0.30), residues: 386 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 94 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 93 time to evaluate : 1.212 Fit side-chains outliers start: 1 outliers final: 1 residues processed: 93 average time/residue: 0.1855 time to fit residues: 26.3686 Evaluate side-chains 86 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 85 time to evaluate : 1.128 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1185 time to fit residues: 1.7550 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 68 optimal weight: 4.9990 chunk 53 optimal weight: 0.8980 chunk 77 optimal weight: 0.9980 chunk 117 optimal weight: 5.9990 chunk 108 optimal weight: 2.9990 chunk 93 optimal weight: 0.7980 chunk 9 optimal weight: 3.9990 chunk 72 optimal weight: 0.8980 chunk 57 optimal weight: 6.9990 chunk 74 optimal weight: 0.9990 chunk 99 optimal weight: 5.9990 overall best weight: 0.9182 ------------------------------------------------------------------------------- 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.8088 moved from start: 0.1993 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.032 9726 Z= 0.160 Angle : 0.558 11.933 13174 Z= 0.267 Chirality : 0.039 0.276 1482 Planarity : 0.003 0.027 1670 Dihedral : 3.843 21.354 1330 Min Nonbonded Distance : 2.227 Molprobity Statistics. All-atom Clashscore : 8.14 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.49 % Favored : 94.51 % Rotamer Outliers : 0.30 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.72 (0.26), residues: 1184 helix: 1.69 (0.20), residues: 764 sheet: -0.56 (1.03), residues: 34 loop : -1.99 (0.31), residues: 386 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2368 Ramachandran restraints generated. 1184 Oldfield, 0 Emsley, 1184 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 90 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 87 time to evaluate : 1.163 Fit side-chains outliers start: 3 outliers final: 1 residues processed: 89 average time/residue: 0.1744 time to fit residues: 24.1940 Evaluate side-chains 85 residues out of total 1006 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 84 time to evaluate : 1.150 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.3742 time to fit residues: 2.0252 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 118 random chunks: chunk 28 optimal weight: 0.9980 chunk 86 optimal weight: 0.7980 chunk 13 optimal weight: 4.9990 chunk 26 optimal weight: 0.7980 chunk 93 optimal weight: 3.9990 chunk 39 optimal weight: 4.9990 chunk 96 optimal weight: 9.9990 chunk 11 optimal weight: 0.9980 chunk 17 optimal weight: 0.9980 chunk 82 optimal weight: 6.9990 chunk 5 optimal weight: 0.9980 overall best weight: 0.9180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3894 r_free = 0.3894 target = 0.147599 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 40)----------------| | r_work = 0.3327 r_free = 0.3327 target = 0.103123 restraints weight = 13643.430| |-----------------------------------------------------------------------------| r_work (start): 0.3261 rms_B_bonded: 2.10 r_work: 0.3104 rms_B_bonded: 2.67 restraints_weight: 0.5000 r_work: 0.2958 rms_B_bonded: 4.27 restraints_weight: 0.2500 r_work (final): 0.2958 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8290 moved from start: 0.2034 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.030 9726 Z= 0.161 Angle : 0.554 11.550 13174 Z= 0.267 Chirality : 0.039 0.274 1482 Planarity : 0.003 0.027 1670 Dihedral : 3.820 21.386 1330 Min Nonbonded Distance : 2.229 Molprobity Statistics. All-atom Clashscore : 8.30 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.91 % Favored : 94.09 % Rotamer Outliers : 0.30 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.75 (0.26), residues: 1184 helix: 1.71 (0.20), residues: 764 sheet: -0.52 (1.03), residues: 34 loop : -2.00 (0.30), residues: 386 =============================================================================== Job complete usr+sys time: 1741.95 seconds wall clock time: 32 minutes 37.77 seconds (1957.77 seconds total)