Starting phenix.real_space_refine on Sat Dec 9 00:01:29 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wli_32584/12_2023/7wli_32584_neut.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wli_32584/12_2023/7wli_32584.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.3 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wli_32584/12_2023/7wli_32584.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wli_32584/12_2023/7wli_32584.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wli_32584/12_2023/7wli_32584_neut.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wli_32584/12_2023/7wli_32584_neut.pdb" } resolution = 3.3 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.068 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: 6 Type Number sf(0) Gaussians Ca 2 9.91 5 P 1 5.49 5 S 78 5.16 5 C 6321 2.51 5 N 1501 2.21 5 O 1590 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A ARG 238": "NH1" <-> "NH2" Residue "A GLU 261": "OE1" <-> "OE2" Residue "A GLU 403": "OE1" <-> "OE2" Residue "A GLU 663": "OE1" <-> "OE2" Residue "A GLU 666": "OE1" <-> "OE2" Residue "A GLU 863": "OE1" <-> "OE2" Residue "A GLU 1190": "OE1" <-> "OE2" Residue "A GLU 1297": "OE1" <-> "OE2" Residue "A ARG 1466": "NH1" <-> "NH2" 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: 9493 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 9493 Number of conformers: 2 Conformer: "A" Number of residues, atoms: 1156, 9482 Unusual residues: {' CA': 2, '3PE': 11, 'NAG': 2, 'Y01': 6} Classifications: {'peptide': 1135, 'undetermined': 21} Incomplete info: {'backbone_only': 1, 'truncation_to_alanine': 8} Link IDs: {'PTRANS': 31, 'TRANS': 1103, None: 21} Not linked: pdbres="LYS A1738 " pdbres=" CA A2301 " Not linked: pdbres=" CA A2301 " pdbres=" CA A2302 " Not linked: pdbres=" CA A2302 " pdbres="NAG A2303 " Not linked: pdbres="NAG A2303 " pdbres="NAG A2304 " Not linked: pdbres="NAG A2304 " pdbres="3PE A2305 " ... (remaining 16 not shown) Chain breaks: 3 Unresolved chain link angles: 3 Unresolved non-hydrogen bonds: 466 Unresolved non-hydrogen angles: 543 Unresolved non-hydrogen dihedrals: 471 Unresolved non-hydrogen chiralities: 15 Planarities with less than four sites: {'3PE:plan-1': 10, '3PE:plan-2': 10, 'TRP:plan': 1, 'Y01:plan-2': 1, 'Y01:plan-1': 3, 'ARG:plan': 3} Unresolved non-hydrogen planarities: 102 Conformer: "B" Number of residues, atoms: 1156, 9482 Unusual residues: {' CA': 2, '3PE': 11, 'NAG': 2, 'Y01': 6} Classifications: {'peptide': 1135, 'undetermined': 21} Incomplete info: {'backbone_only': 1, 'truncation_to_alanine': 8} Link IDs: {'PTRANS': 31, 'TRANS': 1103, None: 21} Not linked: pdbres="LYS A1738 " pdbres=" CA A2301 " Not linked: pdbres=" CA A2301 " pdbres=" CA A2302 " Not linked: pdbres=" CA A2302 " pdbres="NAG A2303 " Not linked: pdbres="NAG A2303 " pdbres="NAG A2304 " Not linked: pdbres="NAG A2304 " pdbres="3PE A2305 " ... (remaining 16 not shown) Chain breaks: 3 Unresolved chain link angles: 3 Unresolved non-hydrogen bonds: 466 Unresolved non-hydrogen angles: 543 Unresolved non-hydrogen dihedrals: 471 Unresolved non-hydrogen chiralities: 15 Planarities with less than four sites: {'3PE:plan-1': 10, '3PE:plan-2': 10, 'TRP:plan': 1, 'Y01:plan-2': 1, 'Y01:plan-1': 3, 'ARG:plan': 3} Unresolved non-hydrogen planarities: 102 bond proxies already assigned to first conformer: 9676 Residues with excluded nonbonded symmetry interactions: 1 residue: pdb=" N APHE A 854 " occ=0.50 ... (20 atoms not shown) pdb=" CZ BPHE A 854 " occ=0.50 Time building chain proxies: 9.58, per 1000 atoms: 1.01 Number of scatterers: 9493 At special positions: 0 Unit cell: (112.32, 123.76, 109.2, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 6 Type Number sf(0) Ca 2 19.99 S 78 16.00 P 1 15.00 O 1590 8.00 N 1501 7.00 C 6321 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=5, symmetry=0 Simple disulfide: pdb=" SG CYS A 239 " - pdb=" SG CYS A 280 " distance=2.03 Simple disulfide: pdb=" SG CYS A 269 " - pdb=" SG CYS A 326 " distance=2.03 Simple disulfide: pdb=" SG CYS A1335 " - pdb=" SG CYS A1349 " distance=2.03 Simple disulfide: pdb=" SG CYS A1646 " - pdb=" SG CYS A1652 " distance=2.03 Simple disulfide: pdb=" SG CYS A1691 " - pdb=" SG CYS A1698 " 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 A2303 " - " ASN A1345 " " NAG A2304 " - " ASN A1342 " Time building additional restraints: 4.38 Conformation dependent library (CDL) restraints added in 3.4 seconds 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2150 Finding SS restraints... Secondary structure from input PDB file: 54 helices and 2 sheets defined 76.7% alpha, 0.5% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.95 Creating SS restraints... Processing helix chain 'A' and resid 66 through 78 removed outlier: 4.192A pdb=" N CYS A 72 " --> pdb=" O PRO A 68 " (cutoff:3.500A) removed outlier: 4.583A pdb=" N LYS A 74 " --> pdb=" O ASN A 70 " (cutoff:3.500A) removed outlier: 5.263A pdb=" N MET A 75 " --> pdb=" O TRP A 71 " (cutoff:3.500A) Processing helix chain 'A' and resid 78 through 98 removed outlier: 3.796A pdb=" N GLU A 82 " --> pdb=" O ASN A 78 " (cutoff:3.500A) Processing helix chain 'A' and resid 109 through 138 Processing helix chain 'A' and resid 144 through 147 Processing helix chain 'A' and resid 148 through 170 removed outlier: 3.586A pdb=" N ARG A 152 " --> pdb=" O ASP A 148 " (cutoff:3.500A) Processing helix chain 'A' and resid 177 through 190 removed outlier: 4.313A pdb=" N VAL A 182 " --> pdb=" O THR A 179 " (cutoff:3.500A) Proline residue: A 185 - end of helix Processing helix chain 'A' and resid 192 through 206 Processing helix chain 'A' and resid 209 through 231 removed outlier: 4.308A pdb=" N PHE A 223 " --> pdb=" O VAL A 219 " (cutoff:3.500A) Processing helix chain 'A' and resid 233 through 237 Processing helix chain 'A' and resid 260 through 264 removed outlier: 3.561A pdb=" N ASP A 263 " --> pdb=" O GLU A 260 " (cutoff:3.500A) Processing helix chain 'A' and resid 318 through 323 Processing helix chain 'A' and resid 332 through 337 Processing helix chain 'A' and resid 342 through 355 removed outlier: 3.779A pdb=" N ALA A 346 " --> pdb=" O ASN A 342 " (cutoff:3.500A) Processing helix chain 'A' and resid 358 through 370 removed outlier: 3.586A pdb=" N ILE A 362 " --> pdb=" O GLY A 358 " (cutoff:3.500A) Processing helix chain 'A' and resid 371 through 373 No H-bonds generated for 'chain 'A' and resid 371 through 373' Processing helix chain 'A' and resid 374 through 383 Processing helix chain 'A' and resid 389 through 413 removed outlier: 3.742A pdb=" N MET A 412 " --> pdb=" O GLU A 408 " (cutoff:3.500A) Processing helix chain 'A' and resid 628 through 638 removed outlier: 4.111A pdb=" N LYS A 632 " --> pdb=" O GLU A 628 " (cutoff:3.500A) removed outlier: 3.657A pdb=" N GLY A 635 " --> pdb=" O ALA A 631 " (cutoff:3.500A) removed outlier: 3.963A pdb=" N ILE A 636 " --> pdb=" O LYS A 632 " (cutoff:3.500A) removed outlier: 3.787A pdb=" N VAL A 637 " --> pdb=" O LEU A 633 " (cutoff:3.500A) Processing helix chain 'A' and resid 641 through 660 removed outlier: 4.200A pdb=" N GLU A 660 " --> pdb=" O SER A 656 " (cutoff:3.500A) Processing helix chain 'A' and resid 665 through 695 removed outlier: 3.768A pdb=" N ASN A 670 " --> pdb=" O GLU A 666 " (cutoff:3.500A) Processing helix chain 'A' and resid 695 through 702 removed outlier: 3.895A pdb=" N ASN A 702 " --> pdb=" O ASP A 698 " (cutoff:3.500A) Processing helix chain 'A' and resid 702 through 721 removed outlier: 3.662A pdb=" N ILE A 706 " --> pdb=" O ASN A 702 " (cutoff:3.500A) Processing helix chain 'A' and resid 725 through 730 Processing helix chain 'A' and resid 731 through 744 removed outlier: 3.861A pdb=" N LEU A 735 " --> pdb=" O THR A 732 " (cutoff:3.500A) Processing helix chain 'A' and resid 746 through 785 removed outlier: 3.643A pdb=" N ARG A 750 " --> pdb=" O PRO A 746 " (cutoff:3.500A) removed outlier: 3.707A pdb=" N LEU A 752 " --> pdb=" O LEU A 748 " (cutoff:3.500A) removed outlier: 3.971A pdb=" N ASN A 761 " --> pdb=" O LYS A 757 " (cutoff:3.500A) removed outlier: 5.635A pdb=" N THR A 764 " --> pdb=" O ASP A 760 " (cutoff:3.500A) removed outlier: 4.207A pdb=" N PHE A 765 " --> pdb=" O ASN A 761 " (cutoff:3.500A) WARNING: missing atoms! Processing helix chain 'A' and resid 806 through 820 Processing helix chain 'A' and resid 822 through 832 removed outlier: 3.599A pdb=" N VAL A 826 " --> pdb=" O ASP A 822 " (cutoff:3.500A) Processing helix chain 'A' and resid 835 through 837 No H-bonds generated for 'chain 'A' and resid 835 through 837' Processing helix chain 'A' and resid 838 through 850 removed outlier: 3.990A pdb=" N PHE A 842 " --> pdb=" O ALA A 838 " (cutoff:3.500A) Processing helix chain 'A' and resid 852 through 866 Processing helix chain 'A' and resid 1171 through 1188 Processing helix chain 'A' and resid 1198 through 1227 Processing helix chain 'A' and resid 1234 through 1237 removed outlier: 3.520A pdb=" N ARG A1237 " --> pdb=" O ALA A1234 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 1234 through 1237' Processing helix chain 'A' and resid 1238 through 1260 Processing helix chain 'A' and resid 1267 through 1281 Processing helix chain 'A' and resid 1289 through 1302 Processing helix chain 'A' and resid 1306 through 1328 Processing helix chain 'A' and resid 1345 through 1353 Processing helix chain 'A' and resid 1365 through 1378 Processing helix chain 'A' and resid 1381 through 1392 removed outlier: 3.695A pdb=" N ILE A1385 " --> pdb=" O GLY A1381 " (cutoff:3.500A) Processing helix chain 'A' and resid 1404 through 1406 No H-bonds generated for 'chain 'A' and resid 1404 through 1406' Processing helix chain 'A' and resid 1407 through 1463 removed outlier: 3.918A pdb=" N LEU A1415 " --> pdb=" O PHE A1411 " (cutoff:3.500A) removed outlier: 3.519A pdb=" N LYS A1436 " --> pdb=" O GLU A1432 " (cutoff:3.500A) removed outlier: 3.674A pdb=" N GLU A1451 " --> pdb=" O ARG A1447 " (cutoff:3.500A) Processing helix chain 'A' and resid 1486 through 1505 removed outlier: 3.564A pdb=" N ASP A1490 " --> pdb=" O SER A1486 " (cutoff:3.500A) removed outlier: 4.466A pdb=" N PHE A1492 " --> pdb=" O TYR A1488 " (cutoff:3.500A) removed outlier: 3.744A pdb=" N ILE A1493 " --> pdb=" O LEU A1489 " (cutoff:3.500A) Processing helix chain 'A' and resid 1512 through 1541 Processing helix chain 'A' and resid 1542 through 1549 removed outlier: 3.575A pdb=" N PHE A1546 " --> pdb=" O GLY A1542 " (cutoff:3.500A) removed outlier: 3.643A pdb=" N PHE A1547 " --> pdb=" O LEU A1543 " (cutoff:3.500A) Processing helix chain 'A' and resid 1551 through 1575 removed outlier: 3.577A pdb=" N LEU A1568 " --> pdb=" O MET A1564 " (cutoff:3.500A) removed outlier: 3.623A pdb=" N GLU A1569 " --> pdb=" O GLY A1565 " (cutoff:3.500A) removed outlier: 3.575A pdb=" N GLU A1570 " --> pdb=" O ILE A1566 " (cutoff:3.500A) removed outlier: 3.572A pdb=" N ILE A1571 " --> pdb=" O THR A1567 " (cutoff:3.500A) removed outlier: 3.791A pdb=" N GLU A1572 " --> pdb=" O LEU A1568 " (cutoff:3.500A) removed outlier: 3.665A pdb=" N ILE A1573 " --> pdb=" O GLU A1569 " (cutoff:3.500A) Processing helix chain 'A' and resid 1580 through 1590 removed outlier: 3.528A pdb=" N ILE A1584 " --> pdb=" O ASN A1580 " (cutoff:3.500A) Processing helix chain 'A' and resid 1591 through 1600 removed outlier: 3.660A pdb=" N LYS A1597 " --> pdb=" O ARG A1594 " (cutoff:3.500A) Processing helix chain 'A' and resid 1602 through 1614 Processing helix chain 'A' and resid 1615 through 1642 Processing helix chain 'A' and resid 1664 through 1676 removed outlier: 3.863A pdb=" N THR A1670 " --> pdb=" O MET A1666 " (cutoff:3.500A) Processing helix chain 'A' and resid 1679 through 1688 Processing helix chain 'A' and resid 1695 through 1700 removed outlier: 4.173A pdb=" N LEU A1699 " --> pdb=" O GLU A1695 " (cutoff:3.500A) Processing helix chain 'A' and resid 1701 through 1738 removed outlier: 3.576A pdb=" N VAL A1705 " --> pdb=" O SER A1701 " (cutoff:3.500A) Proline residue: A1707 - end of helix Processing sheet with id=AA1, first strand: chain 'A' and resid 256 through 257 removed outlier: 4.753A pdb=" N CYS A 239 " --> pdb=" O TYR A 257 " (cutoff:3.500A) removed outlier: 3.912A pdb=" N ARG A 238 " --> pdb=" O ARG A 327 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 1333 through 1334 616 hydrogen bonds defined for protein. 1755 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.99 Time building geometry restraints manager: 4.50 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.18 - 1.31: 1426 1.31 - 1.44: 2571 1.44 - 1.57: 5441 1.57 - 1.69: 137 1.69 - 1.82: 127 Bond restraints: 9702 Sorted by residual: bond pdb=" C ALA A1443 " pdb=" O ALA A1443 " ideal model delta sigma weight residual 1.236 1.356 -0.120 1.15e-02 7.56e+03 1.08e+02 bond pdb=" C GLN A1439 " pdb=" O GLN A1439 " ideal model delta sigma weight residual 1.236 1.328 -0.092 1.15e-02 7.56e+03 6.37e+01 bond pdb=" C LYS A1459 " pdb=" O LYS A1459 " ideal model delta sigma weight residual 1.237 1.330 -0.093 1.17e-02 7.31e+03 6.35e+01 bond pdb=" C21 3PE A2305 " pdb=" O21 3PE A2305 " ideal model delta sigma weight residual 1.316 1.468 -0.152 2.00e-02 2.50e+03 5.79e+01 bond pdb=" C GLU A1450 " pdb=" O GLU A1450 " ideal model delta sigma weight residual 1.237 1.316 -0.079 1.16e-02 7.43e+03 4.66e+01 ... (remaining 9697 not shown) Histogram of bond angle deviations from ideal: 98.40 - 105.56: 184 105.56 - 112.72: 5303 112.72 - 119.87: 3240 119.87 - 127.03: 4316 127.03 - 134.19: 87 Bond angle restraints: 13130 Sorted by residual: angle pdb=" N ALA A1443 " pdb=" CA ALA A1443 " pdb=" C ALA A1443 " ideal model delta sigma weight residual 110.97 123.96 -12.99 1.09e+00 8.42e-01 1.42e+02 angle pdb=" N LYS A1463 " pdb=" CA LYS A1463 " pdb=" C LYS A1463 " ideal model delta sigma weight residual 108.31 123.19 -14.88 1.52e+00 4.33e-01 9.58e+01 angle pdb=" C ALA A1446 " pdb=" N ARG A1447 " pdb=" CA ARG A1447 " ideal model delta sigma weight residual 120.28 109.12 11.16 1.44e+00 4.82e-01 6.00e+01 angle pdb=" N ALA A1446 " pdb=" CA ALA A1446 " pdb=" C ALA A1446 " ideal model delta sigma weight residual 111.07 103.08 7.99 1.07e+00 8.73e-01 5.57e+01 angle pdb=" N ARG A1453 " pdb=" CA ARG A1453 " pdb=" C ARG A1453 " ideal model delta sigma weight residual 112.23 103.94 8.29 1.26e+00 6.30e-01 4.33e+01 ... (remaining 13125 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 18.12: 5328 18.12 - 36.25: 606 36.25 - 54.37: 132 54.37 - 72.49: 21 72.49 - 90.61: 5 Dihedral angle restraints: 6092 sinusoidal: 2773 harmonic: 3319 Sorted by residual: dihedral pdb=" CB CYS A1646 " pdb=" SG CYS A1646 " pdb=" SG CYS A1652 " pdb=" CB CYS A1652 " ideal model delta sinusoidal sigma weight residual 93.00 174.76 -81.76 1 1.00e+01 1.00e-02 8.22e+01 dihedral pdb=" CB CYS A 239 " pdb=" SG CYS A 239 " pdb=" SG CYS A 280 " pdb=" CB CYS A 280 " ideal model delta sinusoidal sigma weight residual 93.00 145.97 -52.97 1 1.00e+01 1.00e-02 3.82e+01 dihedral pdb=" CA ASN A 332 " pdb=" C ASN A 332 " pdb=" N PRO A 333 " pdb=" CA PRO A 333 " ideal model delta harmonic sigma weight residual -180.00 -151.66 -28.34 0 5.00e+00 4.00e-02 3.21e+01 ... (remaining 6089 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.101: 1444 0.101 - 0.202: 64 0.202 - 0.303: 15 0.303 - 0.404: 6 0.404 - 0.505: 6 Chirality restraints: 1535 Sorted by residual: chirality pdb=" CBG Y01 A2312 " pdb=" CAQ Y01 A2312 " pdb=" CBD Y01 A2312 " pdb=" CBI Y01 A2312 " both_signs ideal model delta sigma weight residual False -2.33 -2.84 0.51 2.00e-01 2.50e+01 6.39e+00 chirality pdb=" CBG Y01 A2315 " pdb=" CAQ Y01 A2315 " pdb=" CBD Y01 A2315 " pdb=" CBI Y01 A2315 " both_signs ideal model delta sigma weight residual False -2.33 -2.82 0.49 2.00e-01 2.50e+01 6.04e+00 chirality pdb=" CBG Y01 A2309 " pdb=" CAQ Y01 A2309 " pdb=" CBD Y01 A2309 " pdb=" CBI Y01 A2309 " both_signs ideal model delta sigma weight residual False -2.33 -2.82 0.49 2.00e-01 2.50e+01 6.03e+00 ... (remaining 1532 not shown) Planarity restraints: 1579 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA ARG A1449 " -0.024 2.00e-02 2.50e+03 4.64e-02 2.16e+01 pdb=" C ARG A1449 " 0.080 2.00e-02 2.50e+03 pdb=" O ARG A1449 " -0.029 2.00e-02 2.50e+03 pdb=" N GLU A1450 " -0.028 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C ASN A 702 " 0.054 5.00e-02 4.00e+02 8.13e-02 1.06e+01 pdb=" N PRO A 703 " -0.141 5.00e-02 4.00e+02 pdb=" CA PRO A 703 " 0.041 5.00e-02 4.00e+02 pdb=" CD PRO A 703 " 0.045 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CG GLU A 403 " -0.016 2.00e-02 2.50e+03 3.24e-02 1.05e+01 pdb=" CD GLU A 403 " 0.056 2.00e-02 2.50e+03 pdb=" OE1 GLU A 403 " -0.020 2.00e-02 2.50e+03 pdb=" OE2 GLU A 403 " -0.021 2.00e-02 2.50e+03 ... (remaining 1576 not shown) Histogram of nonbonded interaction distances: 1.89 - 2.49: 47 2.49 - 3.09: 7232 3.09 - 3.70: 14554 3.70 - 4.30: 21308 4.30 - 4.90: 34299 Nonbonded interactions: 77440 Sorted by model distance: nonbonded pdb=" OE1 GLN A1441 " pdb=" OD2 ASP A1734 " model vdw 1.888 3.040 nonbonded pdb=" OD1 ASP A1380 " pdb=" ND2 ASN A1681 " model vdw 2.139 2.520 nonbonded pdb=" O HIS A1481 " pdb=" OG1 THR A1485 " model vdw 2.144 2.440 nonbonded pdb=" OG SER A 270 " pdb=" O ASP A 274 " model vdw 2.149 2.440 nonbonded pdb=" OH TYR A1524 " pdb=" OE1 GLU A1569 " model vdw 2.151 2.440 ... (remaining 77435 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. 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=0.50 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 1.010 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: 4.140 Check model and map are aligned: 0.140 Set scattering table: 0.100 Process input model: 35.630 Find NCS groups from input model: 0.260 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.860 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 44.210 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7600 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.012 0.152 9702 Z= 0.739 Angle : 1.012 14.878 13130 Z= 0.525 Chirality : 0.064 0.505 1535 Planarity : 0.005 0.081 1577 Dihedral : 15.399 90.614 3927 Min Nonbonded Distance : 1.888 Molprobity Statistics. All-atom Clashscore : 15.00 Ramachandran Plot: Outliers : 0.18 % Allowed : 7.28 % Favored : 92.55 % Rotamer: Outliers : 0.20 % Allowed : 0.60 % Favored : 99.20 % Cbeta Deviations : 0.09 % 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.17 (0.26), residues: 1130 helix: 0.53 (0.19), residues: 759 sheet: None (None), residues: 0 loop : -1.36 (0.36), residues: 371 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.015 0.001 TRP A 823 HIS 0.005 0.001 HIS A1358 PHE 0.043 0.002 PHE A1626 TYR 0.019 0.002 TYR A1361 ARG 0.004 0.000 ARG A 749 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 220 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 219 time to evaluate : 0.965 Fit side-chains revert: symmetry clash TARDY: cannot create tardy model for: "CYS A 766 " (corrupted residue). Skipping it. revert: symmetry clash revert: symmetry clash outliers start: 1 outliers final: 0 residues processed: 220 average time/residue: 0.1961 time to fit residues: 62.4000 Evaluate side-chains 185 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 185 time to evaluate : 1.126 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.4296 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 115 random chunks: chunk 97 optimal weight: 9.9990 chunk 87 optimal weight: 3.9990 chunk 48 optimal weight: 0.8980 chunk 29 optimal weight: 0.9990 chunk 58 optimal weight: 6.9990 chunk 46 optimal weight: 0.8980 chunk 90 optimal weight: 0.6980 chunk 34 optimal weight: 4.9990 chunk 54 optimal weight: 0.7980 chunk 67 optimal weight: 0.3980 chunk 104 optimal weight: 6.9990 overall best weight: 0.7380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 339 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 643 ASN A 653 ASN A 761 ASN A1241 ASN A1404 ASN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7583 moved from start: 0.1272 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.049 9702 Z= 0.226 Angle : 0.687 14.021 13130 Z= 0.343 Chirality : 0.042 0.243 1535 Planarity : 0.005 0.056 1577 Dihedral : 9.711 58.556 1807 Min Nonbonded Distance : 2.151 Molprobity Statistics. All-atom Clashscore : 14.01 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.83 % Favored : 93.17 % Rotamer: Outliers : 1.09 % Allowed : 9.74 % Favored : 89.17 % 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.18 (0.26), residues: 1130 helix: 0.93 (0.19), residues: 776 sheet: None (None), residues: 0 loop : -1.69 (0.35), residues: 354 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.001 TRP A1551 HIS 0.003 0.001 HIS A1358 PHE 0.033 0.002 PHE A1626 TYR 0.014 0.002 TYR A1709 ARG 0.006 0.000 ARG A 731 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 208 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 197 time to evaluate : 1.256 Fit side-chains TARDY: cannot create tardy model for: "CYS A 766 " (corrupted residue). Skipping it. revert: symmetry clash outliers start: 11 outliers final: 7 residues processed: 201 average time/residue: 0.1878 time to fit residues: 55.8502 Evaluate side-chains 188 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 181 time to evaluate : 1.077 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.0973 time to fit residues: 2.7832 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 115 random chunks: chunk 58 optimal weight: 1.9990 chunk 32 optimal weight: 7.9990 chunk 86 optimal weight: 10.0000 chunk 71 optimal weight: 0.8980 chunk 28 optimal weight: 0.8980 chunk 104 optimal weight: 5.9990 chunk 113 optimal weight: 5.9990 chunk 93 optimal weight: 9.9990 chunk 103 optimal weight: 0.7980 chunk 35 optimal weight: 0.9980 chunk 83 optimal weight: 0.9980 overall best weight: 0.9180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 339 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 761 ASN ** A1326 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7589 moved from start: 0.1619 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.043 9702 Z= 0.237 Angle : 0.645 13.558 13130 Z= 0.325 Chirality : 0.042 0.235 1535 Planarity : 0.004 0.049 1577 Dihedral : 9.208 59.184 1807 Min Nonbonded Distance : 2.157 Molprobity Statistics. All-atom Clashscore : 13.96 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.30 % Favored : 93.70 % Rotamer: Outliers : 0.70 % Allowed : 12.92 % Favored : 86.38 % 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.32 (0.26), residues: 1130 helix: 1.07 (0.19), residues: 779 sheet: None (None), residues: 0 loop : -1.71 (0.35), residues: 351 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.022 0.001 TRP A1551 HIS 0.004 0.001 HIS A1358 PHE 0.027 0.002 PHE A1626 TYR 0.012 0.001 TYR A1387 ARG 0.006 0.000 ARG A 731 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 199 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 192 time to evaluate : 1.374 Fit side-chains TARDY: cannot create tardy model for: "CYS A 766 " (corrupted residue). Skipping it. revert: symmetry clash outliers start: 7 outliers final: 3 residues processed: 195 average time/residue: 0.1922 time to fit residues: 55.3664 Evaluate side-chains 179 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 176 time to evaluate : 1.100 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1005 time to fit residues: 2.0491 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 115 random chunks: chunk 103 optimal weight: 0.5980 chunk 78 optimal weight: 1.9990 chunk 54 optimal weight: 0.7980 chunk 11 optimal weight: 3.9990 chunk 49 optimal weight: 0.6980 chunk 70 optimal weight: 0.8980 chunk 104 optimal weight: 7.9990 chunk 111 optimal weight: 0.6980 chunk 99 optimal weight: 5.9990 chunk 29 optimal weight: 0.0170 chunk 92 optimal weight: 1.9990 overall best weight: 0.5618 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 339 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 761 ASN A 866 GLN ** A1326 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1353 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7547 moved from start: 0.1937 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.040 9702 Z= 0.191 Angle : 0.628 19.043 13130 Z= 0.314 Chirality : 0.041 0.237 1535 Planarity : 0.004 0.048 1577 Dihedral : 8.729 59.499 1807 Min Nonbonded Distance : 2.193 Molprobity Statistics. All-atom Clashscore : 13.55 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.94 % Favored : 94.06 % Rotamer: Outliers : 0.70 % Allowed : 14.41 % Favored : 84.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.56 (0.26), residues: 1130 helix: 1.25 (0.19), residues: 779 sheet: None (None), residues: 0 loop : -1.64 (0.35), residues: 351 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.016 0.001 TRP A1551 HIS 0.003 0.001 HIS A1435 PHE 0.023 0.002 PHE A1434 TYR 0.013 0.001 TYR A1387 ARG 0.007 0.000 ARG A 731 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 216 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 209 time to evaluate : 1.147 Fit side-chains TARDY: cannot create tardy model for: "CYS A 766 " (corrupted residue). Skipping it. revert: symmetry clash outliers start: 7 outliers final: 3 residues processed: 210 average time/residue: 0.1874 time to fit residues: 58.0311 Evaluate side-chains 190 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 187 time to evaluate : 1.096 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1012 time to fit residues: 2.0731 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 115 random chunks: chunk 63 optimal weight: 4.9990 chunk 1 optimal weight: 1.9990 chunk 82 optimal weight: 0.5980 chunk 45 optimal weight: 0.8980 chunk 94 optimal weight: 0.0470 chunk 76 optimal weight: 2.9990 chunk 0 optimal weight: 10.0000 chunk 56 optimal weight: 10.0000 chunk 99 optimal weight: 2.9990 chunk 28 optimal weight: 0.4980 chunk 37 optimal weight: 4.9990 overall best weight: 0.8080 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 339 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1353 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7561 moved from start: 0.2151 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.042 9702 Z= 0.217 Angle : 0.621 13.118 13130 Z= 0.312 Chirality : 0.041 0.234 1535 Planarity : 0.004 0.048 1577 Dihedral : 8.638 59.589 1807 Min Nonbonded Distance : 2.152 Molprobity Statistics. All-atom Clashscore : 14.12 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.03 % Favored : 93.97 % Rotamer: Outliers : 0.99 % Allowed : 15.31 % Favored : 83.70 % 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.62 (0.26), residues: 1130 helix: 1.29 (0.19), residues: 778 sheet: None (None), residues: 0 loop : -1.60 (0.35), residues: 352 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP A1551 HIS 0.003 0.001 HIS A1358 PHE 0.021 0.002 PHE A1434 TYR 0.013 0.001 TYR A1387 ARG 0.007 0.000 ARG A 731 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 208 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 198 time to evaluate : 1.064 Fit side-chains revert: symmetry clash TARDY: cannot create tardy model for: "CYS A 766 " (corrupted residue). Skipping it. revert: symmetry clash outliers start: 10 outliers final: 4 residues processed: 200 average time/residue: 0.1837 time to fit residues: 54.4365 Evaluate side-chains 190 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 186 time to evaluate : 1.049 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.0983 time to fit residues: 2.1651 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 115 random chunks: chunk 100 optimal weight: 0.7980 chunk 21 optimal weight: 3.9990 chunk 65 optimal weight: 2.9990 chunk 27 optimal weight: 2.9990 chunk 111 optimal weight: 1.9990 chunk 92 optimal weight: 0.9980 chunk 51 optimal weight: 0.6980 chunk 9 optimal weight: 3.9990 chunk 36 optimal weight: 3.9990 chunk 58 optimal weight: 7.9990 chunk 107 optimal weight: 0.9980 overall best weight: 1.0982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A1353 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7584 moved from start: 0.2264 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.045 9702 Z= 0.253 Angle : 0.637 9.513 13130 Z= 0.322 Chirality : 0.042 0.233 1535 Planarity : 0.004 0.047 1577 Dihedral : 8.687 59.857 1807 Min Nonbonded Distance : 2.201 Molprobity Statistics. All-atom Clashscore : 13.55 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.30 % Favored : 93.70 % Rotamer: Outliers : 0.99 % Allowed : 16.60 % Favored : 82.41 % 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.57 (0.26), residues: 1130 helix: 1.24 (0.19), residues: 778 sheet: None (None), residues: 0 loop : -1.58 (0.35), residues: 352 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP A1551 HIS 0.003 0.001 HIS A1358 PHE 0.021 0.002 PHE A 848 TYR 0.016 0.001 TYR A1387 ARG 0.007 0.000 ARG A 731 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 212 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 202 time to evaluate : 1.010 Fit side-chains revert: symmetry clash TARDY: cannot create tardy model for: "CYS A 766 " (corrupted residue). Skipping it. revert: symmetry clash revert: symmetry clash outliers start: 10 outliers final: 7 residues processed: 206 average time/residue: 0.2054 time to fit residues: 62.7933 Evaluate side-chains 190 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 183 time to evaluate : 1.083 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.1586 time to fit residues: 3.7508 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 115 random chunks: chunk 12 optimal weight: 0.8980 chunk 63 optimal weight: 0.6980 chunk 81 optimal weight: 1.9990 chunk 62 optimal weight: 0.8980 chunk 93 optimal weight: 5.9990 chunk 110 optimal weight: 0.6980 chunk 69 optimal weight: 0.5980 chunk 67 optimal weight: 0.3980 chunk 51 optimal weight: 0.5980 chunk 68 optimal weight: 0.6980 chunk 44 optimal weight: 0.7980 overall best weight: 0.5980 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 339 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1353 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A1433 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7538 moved from start: 0.2452 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.035 9702 Z= 0.193 Angle : 0.613 9.565 13130 Z= 0.308 Chirality : 0.041 0.237 1535 Planarity : 0.004 0.048 1577 Dihedral : 8.357 59.816 1807 Min Nonbonded Distance : 2.210 Molprobity Statistics. All-atom Clashscore : 13.39 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.41 % Favored : 94.59 % Rotamer: Outliers : 0.30 % Allowed : 17.10 % Favored : 82.60 % 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.65 (0.26), residues: 1130 helix: 1.31 (0.19), residues: 779 sheet: None (None), residues: 0 loop : -1.58 (0.35), residues: 351 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP A1551 HIS 0.004 0.001 HIS A1358 PHE 0.028 0.002 PHE A1434 TYR 0.017 0.001 TYR A1387 ARG 0.007 0.000 ARG A 731 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 199 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 196 time to evaluate : 1.026 Fit side-chains TARDY: cannot create tardy model for: "CYS A 766 " (corrupted residue). Skipping it. revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 3 outliers final: 1 residues processed: 196 average time/residue: 0.1856 time to fit residues: 54.1941 Evaluate side-chains 185 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 184 time to evaluate : 1.105 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1122 time to fit residues: 1.7659 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 115 random chunks: chunk 66 optimal weight: 3.9990 chunk 33 optimal weight: 8.9990 chunk 21 optimal weight: 0.9980 chunk 70 optimal weight: 1.9990 chunk 75 optimal weight: 0.0270 chunk 54 optimal weight: 2.9990 chunk 10 optimal weight: 9.9990 chunk 87 optimal weight: 3.9990 chunk 100 optimal weight: 1.9990 chunk 106 optimal weight: 0.6980 chunk 96 optimal weight: 4.9990 overall best weight: 1.1442 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 339 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1353 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7582 moved from start: 0.2512 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.043 9702 Z= 0.263 Angle : 0.657 9.487 13130 Z= 0.333 Chirality : 0.042 0.233 1535 Planarity : 0.004 0.046 1577 Dihedral : 8.491 59.154 1807 Min Nonbonded Distance : 2.173 Molprobity Statistics. All-atom Clashscore : 14.12 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.39 % Favored : 93.61 % Rotamer: Outliers : 0.20 % Allowed : 18.99 % Favored : 80.82 % 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.67 (0.26), residues: 1130 helix: 1.33 (0.19), residues: 779 sheet: None (None), residues: 0 loop : -1.61 (0.35), residues: 351 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.015 0.001 TRP A 231 HIS 0.004 0.001 HIS A1358 PHE 0.027 0.002 PHE A1246 TYR 0.013 0.001 TYR A1387 ARG 0.008 0.000 ARG A 731 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 186 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 184 time to evaluate : 1.079 Fit side-chains TARDY: cannot create tardy model for: "CYS A 766 " (corrupted residue). Skipping it. revert: symmetry clash outliers start: 2 outliers final: 1 residues processed: 184 average time/residue: 0.1833 time to fit residues: 50.7922 Evaluate side-chains 174 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 173 time to evaluate : 1.128 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1209 time to fit residues: 1.7067 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 115 random chunks: chunk 103 optimal weight: 0.5980 chunk 106 optimal weight: 0.0270 chunk 62 optimal weight: 0.1980 chunk 44 optimal weight: 2.9990 chunk 81 optimal weight: 1.9990 chunk 31 optimal weight: 1.9990 chunk 93 optimal weight: 6.9990 chunk 97 optimal weight: 0.0370 chunk 102 optimal weight: 0.6980 chunk 67 optimal weight: 0.7980 chunk 109 optimal weight: 8.9990 overall best weight: 0.3116 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 339 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 816 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1353 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7505 moved from start: 0.2664 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.040 9702 Z= 0.177 Angle : 0.636 9.593 13130 Z= 0.322 Chirality : 0.041 0.238 1535 Planarity : 0.004 0.048 1577 Dihedral : 8.073 59.139 1807 Min Nonbonded Distance : 2.166 Molprobity Statistics. All-atom Clashscore : 13.24 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.68 % Favored : 94.32 % Rotamer: Outliers : 0.30 % Allowed : 19.68 % Favored : 80.02 % 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: 1130 helix: 1.41 (0.19), residues: 779 sheet: None (None), residues: 0 loop : -1.61 (0.35), residues: 351 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A1551 HIS 0.004 0.001 HIS A1358 PHE 0.029 0.001 PHE A1179 TYR 0.014 0.001 TYR A1387 ARG 0.008 0.000 ARG A 731 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 194 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 191 time to evaluate : 1.073 Fit side-chains TARDY: cannot create tardy model for: "CYS A 766 " (corrupted residue). Skipping it. revert: symmetry clash outliers start: 3 outliers final: 0 residues processed: 192 average time/residue: 0.1824 time to fit residues: 52.3118 Evaluate side-chains 184 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 184 time to evaluate : 1.032 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.4675 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 115 random chunks: chunk 66 optimal weight: 3.9990 chunk 51 optimal weight: 0.3980 chunk 75 optimal weight: 2.9990 chunk 114 optimal weight: 6.9990 chunk 105 optimal weight: 0.8980 chunk 91 optimal weight: 3.9990 chunk 9 optimal weight: 3.9990 chunk 70 optimal weight: 0.8980 chunk 55 optimal weight: 2.9990 chunk 72 optimal weight: 0.8980 chunk 97 optimal weight: 4.9990 overall best weight: 1.2182 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 339 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 352 GLN ** A1353 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1553 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7576 moved from start: 0.2679 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.041 9702 Z= 0.278 Angle : 0.679 9.458 13130 Z= 0.344 Chirality : 0.044 0.279 1535 Planarity : 0.004 0.045 1577 Dihedral : 8.312 58.057 1807 Min Nonbonded Distance : 2.142 Molprobity Statistics. All-atom Clashscore : 14.27 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.30 % Favored : 93.70 % Rotamer: Outliers : 0.10 % Allowed : 20.58 % Favored : 79.32 % 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.68 (0.26), residues: 1130 helix: 1.35 (0.19), residues: 779 sheet: None (None), residues: 0 loop : -1.62 (0.35), residues: 351 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP A 231 HIS 0.005 0.001 HIS A1358 PHE 0.047 0.002 PHE A 125 TYR 0.021 0.002 TYR A1361 ARG 0.009 0.000 ARG A 731 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2260 Ramachandran restraints generated. 1130 Oldfield, 0 Emsley, 1130 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 187 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 186 time to evaluate : 1.124 Fit side-chains TARDY: cannot create tardy model for: "CYS A 766 " (corrupted residue). Skipping it. revert: symmetry clash outliers start: 1 outliers final: 1 residues processed: 186 average time/residue: 0.1813 time to fit residues: 50.6922 Evaluate side-chains 179 residues out of total 1012 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 178 time to evaluate : 0.961 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0862 time to fit residues: 1.5748 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 115 random chunks: chunk 27 optimal weight: 0.9980 chunk 84 optimal weight: 5.9990 chunk 13 optimal weight: 2.9990 chunk 25 optimal weight: 0.6980 chunk 91 optimal weight: 6.9990 chunk 38 optimal weight: 6.9990 chunk 93 optimal weight: 0.4980 chunk 11 optimal weight: 0.3980 chunk 16 optimal weight: 0.8980 chunk 80 optimal weight: 0.6980 chunk 5 optimal weight: 6.9990 overall best weight: 0.6380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 190 ASN ** A 339 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 352 GLN ** A1353 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1398 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1553 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4426 r_free = 0.4426 target = 0.210020 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 54)----------------| | r_work = 0.3815 r_free = 0.3815 target = 0.151523 restraints weight = 17568.102| |-----------------------------------------------------------------------------| r_work (start): 0.3754 rms_B_bonded: 3.02 r_work: 0.3588 rms_B_bonded: 3.85 restraints_weight: 0.5000 r_work (final): 0.3588 ------------------------------------------------------------------------------- Occupancy refinement ******************** r_start: 0.3590 |-occupancy refinement: start-------------------------------------------------| | r_work = 0.3590 r_free = 0.3590 target_work(ls_wunit_k1) = 0.131 | | occupancies: max = 1.00 min = 0.50 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| |-occupancy refinement: end---------------------------------------------------| | r_work = 0.3590 r_free = 0.3590 target_work(ls_wunit_k1) = 0.131 | | occupancies: max = 1.00 min = 0.47 number of occupancies < 0.1: 0 | |-----------------------------------------------------------------------------| r_final: 0.3590 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7684 moved from start: 0.2726 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.044 9702 Z= 0.208 Angle : 0.662 9.544 13130 Z= 0.336 Chirality : 0.043 0.263 1535 Planarity : 0.004 0.047 1577 Dihedral : 8.139 59.164 1807 Min Nonbonded Distance : 2.156 Molprobity Statistics. All-atom Clashscore : 14.48 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.50 % Favored : 94.50 % Rotamer: Outliers : 0.00 % Allowed : 20.58 % Favored : 79.42 % 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.71 (0.26), residues: 1130 helix: 1.36 (0.19), residues: 779 sheet: None (None), residues: 0 loop : -1.58 (0.36), residues: 351 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP A1551 HIS 0.004 0.001 HIS A1358 PHE 0.028 0.002 PHE A 125 TYR 0.017 0.001 TYR A 364 ARG 0.009 0.000 ARG A 731 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 2348.97 seconds wall clock time: 43 minutes 15.14 seconds (2595.14 seconds total)