Starting phenix.real_space_refine (version: dev) on Wed Mar 15 01:06:25 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wah_32385/03_2023/7wah_32385.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wah_32385/03_2023/7wah_32385.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.45 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wah_32385/03_2023/7wah_32385.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wah_32385/03_2023/7wah_32385.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wah_32385/03_2023/7wah_32385.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7wah_32385/03_2023/7wah_32385.pdb" } resolution = 2.45 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.006 sd= 0.105 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Process input model Symmetric amino acids flipped Residue "A PHE 161": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 507": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 638": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 776": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 804": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 862": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 934": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 1042": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 1081": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 1126": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 1208": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 1228": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 1276": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 1534": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.04s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4893/modules/chem_data/mon_lib" Total number of atoms: 13217 Number of models: 1 Model: "" Number of chains: 3 Chain: "R" Number of atoms: 826 Number of conformers: 1 Conformer: "" Number of residues, atoms: 39, 826 Classifications: {'RNA': 39} Modifications used: {'5*END': 1, 'rna2p_pur': 14, 'rna2p_pyr': 8, 'rna3p_pur': 5, 'rna3p_pyr': 12} Link IDs: {'rna2p': 21, 'rna3p': 17} Chain: "D" Number of atoms: 492 Number of conformers: 1 Conformer: "" Number of residues, atoms: 23, 492 Classifications: {'RNA': 23} Modifications used: {'rna2p_pur': 3, 'rna2p_pyr': 1, 'rna3p_pur': 11, 'rna3p_pyr': 8} Link IDs: {'rna2p': 4, 'rna3p': 18} Chain: "A" Number of atoms: 11899 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1477, 11899 Unusual residues: {' ZN': 4} Classifications: {'peptide': 1473, 'undetermined': 4} Incomplete info: {'truncation_to_alanine': 4} Link IDs: {'PTRANS': 80, 'TRANS': 1392, None: 4} Not linked: pdbres="SER A1603 " pdbres=" ZN A1701 " Not linked: pdbres=" ZN A1701 " pdbres=" ZN A1702 " Not linked: pdbres=" ZN A1702 " pdbres=" ZN A1703 " Not linked: pdbres=" ZN A1703 " pdbres=" ZN A1704 " Chain breaks: 9 Unresolved non-hydrogen bonds: 20 Unresolved non-hydrogen angles: 24 Unresolved non-hydrogen dihedrals: 16 Planarities with less than four sites: {'GLU:plan': 2, 'ARG:plan': 2} Unresolved non-hydrogen planarities: 18 List of CYS excluded from plausible disulfide bonds: (reason: may participate in coordination) ATOM 2022 SG CYS A 86 48.068 76.770 17.624 1.00 74.00 S ATOM 2268 SG CYS A 115 48.227 74.846 14.354 1.00 96.17 S ATOM 2330 SG CYS A 123 50.567 73.956 17.218 1.00 70.51 S ATOM 2354 SG CYS A 126 46.859 73.185 17.504 1.00 60.85 S ATOM 4438 SG CYS A 463 42.328 90.654 47.896 1.00 71.44 S ATOM 4500 SG CYS A 472 42.013 91.085 44.076 1.00 57.69 S ATOM 4514 SG CYS A 474 45.077 92.311 45.949 1.00 59.04 S ATOM 4536 SG CYS A 477 44.492 88.593 45.583 1.00 60.88 S ATOM 6350 SG CYS A 706 31.561 96.931 73.123 1.00 95.71 S ATOM 6365 SG CYS A 708 34.219 98.389 75.450 1.00 77.27 S ATOM 6387 SG CYS A 711 35.034 95.356 73.340 1.00 74.09 S ATOM 8250 SG CYS A 965 30.141 81.704 108.993 1.00 40.04 S ATOM 11040 SG CYS A1312 28.017 82.679 105.973 1.00 54.59 S ATOM 11103 SG CYS A1342 29.737 85.304 108.209 1.00 31.20 S ATOM 11121 SG CYS A1345 31.771 83.230 105.905 1.00 60.12 S Time building chain proxies: 7.70, per 1000 atoms: 0.58 Number of scatterers: 13217 At special positions: 0 Unit cell: (84.66, 123.67, 172.64, 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) Zn 4 29.99 S 56 16.00 P 61 15.00 O 2618 8.00 N 2348 7.00 C 8130 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: 5.63 Conformation dependent library (CDL) restraints added in 1.8 seconds Dynamic metal coordination Zn2+ tetrahedral coordination pdb=" ZN A1701 " pdb="ZN ZN A1701 " - pdb=" SG CYS A 86 " pdb="ZN ZN A1701 " - pdb=" SG CYS A 123 " pdb="ZN ZN A1701 " - pdb=" SG CYS A 126 " pdb="ZN ZN A1701 " - pdb=" SG CYS A 115 " pdb=" ZN A1702 " pdb="ZN ZN A1702 " - pdb=" SG CYS A 474 " pdb="ZN ZN A1702 " - pdb=" SG CYS A 463 " pdb="ZN ZN A1702 " - pdb=" SG CYS A 477 " pdb="ZN ZN A1702 " - pdb=" SG CYS A 472 " pdb=" ZN A1703 " pdb="ZN ZN A1703 " - pdb=" ND1 HIS A 703 " pdb="ZN ZN A1703 " - pdb=" SG CYS A 708 " pdb="ZN ZN A1703 " - pdb=" SG CYS A 711 " pdb="ZN ZN A1703 " - pdb=" SG CYS A 706 " pdb=" ZN A1704 " pdb="ZN ZN A1704 " - pdb=" SG CYS A1345 " pdb="ZN ZN A1704 " - pdb=" SG CYS A1342 " pdb="ZN ZN A1704 " - pdb=" SG CYS A1312 " pdb="ZN ZN A1704 " - pdb=" SG CYS A 965 " Number of angles added : 21 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2736 Finding SS restraints... Secondary structure from input PDB file: 44 helices and 23 sheets defined 28.5% alpha, 29.7% beta 20 base pairs and 23 stacking pairs defined. Time for finding SS restraints: 4.51 Creating SS restraints... Processing helix chain 'A' and resid 30 through 36 removed outlier: 3.547A pdb=" N GLY A 36 " --> pdb=" O GLU A 32 " (cutoff:3.500A) Processing helix chain 'A' and resid 37 through 39 No H-bonds generated for 'chain 'A' and resid 37 through 39' Processing helix chain 'A' and resid 57 through 76 removed outlier: 3.705A pdb=" N LEU A 61 " --> pdb=" O THR A 57 " (cutoff:3.500A) Processing helix chain 'A' and resid 123 through 130 Processing helix chain 'A' and resid 163 through 167 removed outlier: 3.928A pdb=" N ILE A 167 " --> pdb=" O GLY A 163 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 163 through 167' Processing helix chain 'A' and resid 211 through 225 Processing helix chain 'A' and resid 261 through 276 Processing helix chain 'A' and resid 278 through 280 No H-bonds generated for 'chain 'A' and resid 278 through 280' Processing helix chain 'A' and resid 281 through 292 Processing helix chain 'A' and resid 293 through 300 removed outlier: 5.706A pdb=" N SER A 296 " --> pdb=" O ARG A 293 " (cutoff:3.500A) removed outlier: 6.822A pdb=" N LYS A 297 " --> pdb=" O ARG A 294 " (cutoff:3.500A) Processing helix chain 'A' and resid 328 through 338 Processing helix chain 'A' and resid 340 through 343 removed outlier: 3.626A pdb=" N LYS A 343 " --> pdb=" O LYS A 340 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 340 through 343' Processing helix chain 'A' and resid 344 through 364 Processing helix chain 'A' and resid 443 through 460 Processing helix chain 'A' and resid 474 through 481 Processing helix chain 'A' and resid 542 through 554 Processing helix chain 'A' and resid 562 through 566 Processing helix chain 'A' and resid 584 through 594 Processing helix chain 'A' and resid 595 through 598 removed outlier: 3.654A pdb=" N ARG A 598 " --> pdb=" O TRP A 595 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 595 through 598' Processing helix chain 'A' and resid 599 through 610 removed outlier: 3.540A pdb=" N GLU A 604 " --> pdb=" O GLU A 600 " (cutoff:3.500A) removed outlier: 3.505A pdb=" N GLU A 605 " --> pdb=" O LYS A 601 " (cutoff:3.500A) removed outlier: 3.798A pdb=" N LYS A 608 " --> pdb=" O GLU A 604 " (cutoff:3.500A) removed outlier: 3.984A pdb=" N ARG A 609 " --> pdb=" O GLU A 605 " (cutoff:3.500A) Processing helix chain 'A' and resid 643 through 649 Processing helix chain 'A' and resid 675 through 691 removed outlier: 3.951A pdb=" N ILE A 690 " --> pdb=" O ALA A 686 " (cutoff:3.500A) removed outlier: 3.602A pdb=" N HIS A 691 " --> pdb=" O VAL A 687 " (cutoff:3.500A) Processing helix chain 'A' and resid 708 through 715 Processing helix chain 'A' and resid 780 through 785 Processing helix chain 'A' and resid 785 through 801 Processing helix chain 'A' and resid 809 through 813 Processing helix chain 'A' and resid 946 through 962 removed outlier: 3.579A pdb=" N LEU A 950 " --> pdb=" O PRO A 946 " (cutoff:3.500A) removed outlier: 3.593A pdb=" N THR A 960 " --> pdb=" O SER A 956 " (cutoff:3.500A) Processing helix chain 'A' and resid 1025 through 1030 Processing helix chain 'A' and resid 1163 through 1168 Processing helix chain 'A' and resid 1179 through 1183 Processing helix chain 'A' and resid 1231 through 1247 Processing helix chain 'A' and resid 1254 through 1258 Processing helix chain 'A' and resid 1259 through 1264 Processing helix chain 'A' and resid 1300 through 1305 Processing helix chain 'A' and resid 1306 through 1310 Processing helix chain 'A' and resid 1342 through 1349 Processing helix chain 'A' and resid 1395 through 1400 removed outlier: 4.196A pdb=" N SER A1399 " --> pdb=" O ARG A1395 " (cutoff:3.500A) Processing helix chain 'A' and resid 1419 through 1426 removed outlier: 4.172A pdb=" N THR A1423 " --> pdb=" O HIS A1419 " (cutoff:3.500A) removed outlier: 4.070A pdb=" N ASP A1426 " --> pdb=" O LYS A1422 " (cutoff:3.500A) Processing helix chain 'A' and resid 1464 through 1477 Processing helix chain 'A' and resid 1518 through 1534 Processing helix chain 'A' and resid 1538 through 1549 Processing helix chain 'A' and resid 1565 through 1569 Processing helix chain 'A' and resid 1575 through 1580 Processing helix chain 'A' and resid 1588 through 1594 Processing sheet with id=AA1, first strand: chain 'A' and resid 148 through 150 Processing sheet with id=AA2, first strand: chain 'A' and resid 15 through 17 removed outlier: 3.737A pdb=" N ASP A 185 " --> pdb=" O ARG A 175 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 41 through 44 Processing sheet with id=AA4, first strand: chain 'A' and resid 79 through 80 Processing sheet with id=AA5, first strand: chain 'A' and resid 483 through 485 removed outlier: 5.027A pdb=" N GLU A 409 " --> pdb=" O ASN A 574 " (cutoff:3.500A) removed outlier: 5.385A pdb=" N ASN A 574 " --> pdb=" O GLU A 409 " (cutoff:3.500A) removed outlier: 6.754A pdb=" N VAL A 411 " --> pdb=" O MET A 572 " (cutoff:3.500A) removed outlier: 5.549A pdb=" N MET A 572 " --> pdb=" O VAL A 411 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'A' and resid 483 through 485 removed outlier: 5.027A pdb=" N GLU A 409 " --> pdb=" O ASN A 574 " (cutoff:3.500A) removed outlier: 5.385A pdb=" N ASN A 574 " --> pdb=" O GLU A 409 " (cutoff:3.500A) removed outlier: 6.754A pdb=" N VAL A 411 " --> pdb=" O MET A 572 " (cutoff:3.500A) removed outlier: 5.549A pdb=" N MET A 572 " --> pdb=" O VAL A 411 " (cutoff:3.500A) removed outlier: 3.600A pdb=" N ASN A 611 " --> pdb=" O TYR A 577 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'A' and resid 441 through 442 Processing sheet with id=AA8, first strand: chain 'A' and resid 497 through 504 removed outlier: 10.417A pdb=" N GLU A 497 " --> pdb=" O ASN A 518 " (cutoff:3.500A) removed outlier: 9.255A pdb=" N ASN A 518 " --> pdb=" O GLU A 497 " (cutoff:3.500A) removed outlier: 11.024A pdb=" N ARG A 499 " --> pdb=" O LEU A 516 " (cutoff:3.500A) removed outlier: 9.665A pdb=" N LEU A 516 " --> pdb=" O ARG A 499 " (cutoff:3.500A) removed outlier: 9.810A pdb=" N ARG A 501 " --> pdb=" O GLY A 514 " (cutoff:3.500A) removed outlier: 9.321A pdb=" N GLY A 514 " --> pdb=" O ARG A 501 " (cutoff:3.500A) removed outlier: 6.133A pdb=" N ARG A 503 " --> pdb=" O ALA A 512 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'A' and resid 659 through 664 removed outlier: 7.009A pdb=" N CYS A 672 " --> pdb=" O PHE A 730 " (cutoff:3.500A) removed outlier: 4.081A pdb=" N PHE A 730 " --> pdb=" O CYS A 672 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'A' and resid 659 through 664 removed outlier: 7.009A pdb=" N CYS A 672 " --> pdb=" O PHE A 730 " (cutoff:3.500A) removed outlier: 4.081A pdb=" N PHE A 730 " --> pdb=" O CYS A 672 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'A' and resid 641 through 642 removed outlier: 3.535A pdb=" N VAL A 656 " --> pdb=" O ASP A 642 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'A' and resid 692 through 693 Processing sheet with id=AB4, first strand: chain 'A' and resid 739 through 744 removed outlier: 11.509A pdb=" N PHE A 739 " --> pdb=" O LYS A 756 " (cutoff:3.500A) removed outlier: 10.323A pdb=" N LYS A 756 " --> pdb=" O PHE A 739 " (cutoff:3.500A) removed outlier: 10.100A pdb=" N HIS A 741 " --> pdb=" O LYS A 754 " (cutoff:3.500A) removed outlier: 9.398A pdb=" N LYS A 754 " --> pdb=" O HIS A 741 " (cutoff:3.500A) removed outlier: 5.602A pdb=" N ALA A 743 " --> pdb=" O ALA A 752 " (cutoff:3.500A) Processing sheet with id=AB5, first strand: chain 'A' and resid 1558 through 1559 removed outlier: 6.676A pdb=" N VAL A 861 " --> pdb=" O TRP A1600 " (cutoff:3.500A) Processing sheet with id=AB6, first strand: chain 'A' and resid 865 through 868 Processing sheet with id=AB7, first strand: chain 'A' and resid 938 through 939 removed outlier: 6.884A pdb=" N ILE A 943 " --> pdb=" O LEU A1364 " (cutoff:3.500A) removed outlier: 6.492A pdb=" N SER A1453 " --> pdb=" O THR A 901 " (cutoff:3.500A) removed outlier: 7.343A pdb=" N THR A 901 " --> pdb=" O SER A1453 " (cutoff:3.500A) removed outlier: 5.355A pdb=" N SER A1455 " --> pdb=" O LEU A 899 " (cutoff:3.500A) removed outlier: 6.733A pdb=" N LEU A 899 " --> pdb=" O SER A1455 " (cutoff:3.500A) removed outlier: 5.458A pdb=" N GLU A1457 " --> pdb=" O CYS A 897 " (cutoff:3.500A) removed outlier: 7.143A pdb=" N CYS A 897 " --> pdb=" O GLU A1457 " (cutoff:3.500A) removed outlier: 5.564A pdb=" N ALA A1459 " --> pdb=" O ILE A 895 " (cutoff:3.500A) removed outlier: 6.909A pdb=" N ILE A 895 " --> pdb=" O ALA A1459 " (cutoff:3.500A) removed outlier: 5.220A pdb=" N GLU A1461 " --> pdb=" O GLY A 893 " (cutoff:3.500A) removed outlier: 7.077A pdb=" N GLY A 893 " --> pdb=" O GLU A1461 " (cutoff:3.500A) removed outlier: 6.762A pdb=" N THR A 892 " --> pdb=" O LEU A1505 " (cutoff:3.500A) removed outlier: 4.476A pdb=" N LEU A1505 " --> pdb=" O THR A 892 " (cutoff:3.500A) removed outlier: 6.639A pdb=" N LYS A 894 " --> pdb=" O VAL A1503 " (cutoff:3.500A) removed outlier: 4.457A pdb=" N VAL A1503 " --> pdb=" O LYS A 894 " (cutoff:3.500A) removed outlier: 6.707A pdb=" N ARG A 896 " --> pdb=" O GLU A1501 " (cutoff:3.500A) Processing sheet with id=AB8, first strand: chain 'A' and resid 1357 through 1359 removed outlier: 6.492A pdb=" N SER A1453 " --> pdb=" O THR A 901 " (cutoff:3.500A) removed outlier: 7.343A pdb=" N THR A 901 " --> pdb=" O SER A1453 " (cutoff:3.500A) removed outlier: 5.355A pdb=" N SER A1455 " --> pdb=" O LEU A 899 " (cutoff:3.500A) removed outlier: 6.733A pdb=" N LEU A 899 " --> pdb=" O SER A1455 " (cutoff:3.500A) removed outlier: 5.458A pdb=" N GLU A1457 " --> pdb=" O CYS A 897 " (cutoff:3.500A) removed outlier: 7.143A pdb=" N CYS A 897 " --> pdb=" O GLU A1457 " (cutoff:3.500A) removed outlier: 5.564A pdb=" N ALA A1459 " --> pdb=" O ILE A 895 " (cutoff:3.500A) removed outlier: 6.909A pdb=" N ILE A 895 " --> pdb=" O ALA A1459 " (cutoff:3.500A) removed outlier: 5.220A pdb=" N GLU A1461 " --> pdb=" O GLY A 893 " (cutoff:3.500A) removed outlier: 7.077A pdb=" N GLY A 893 " --> pdb=" O GLU A1461 " (cutoff:3.500A) removed outlier: 6.762A pdb=" N THR A 892 " --> pdb=" O LEU A1505 " (cutoff:3.500A) removed outlier: 4.476A pdb=" N LEU A1505 " --> pdb=" O THR A 892 " (cutoff:3.500A) removed outlier: 6.639A pdb=" N LYS A 894 " --> pdb=" O VAL A1503 " (cutoff:3.500A) removed outlier: 4.457A pdb=" N VAL A1503 " --> pdb=" O LYS A 894 " (cutoff:3.500A) removed outlier: 6.707A pdb=" N ARG A 896 " --> pdb=" O GLU A1501 " (cutoff:3.500A) Processing sheet with id=AB9, first strand: chain 'A' and resid 908 through 909 removed outlier: 4.105A pdb=" N ALA A 935 " --> pdb=" O ASP A 909 " (cutoff:3.500A) Processing sheet with id=AC1, first strand: chain 'A' and resid 932 through 933 removed outlier: 6.884A pdb=" N LYS A 932 " --> pdb=" O THR A1444 " (cutoff:3.500A) Processing sheet with id=AC2, first strand: chain 'A' and resid 976 through 978 Processing sheet with id=AC3, first strand: chain 'A' and resid 991 through 994 removed outlier: 6.826A pdb=" N TYR A1276 " --> pdb=" O ILE A1287 " (cutoff:3.500A) removed outlier: 4.512A pdb=" N ILE A1287 " --> pdb=" O TYR A1276 " (cutoff:3.500A) removed outlier: 6.635A pdb=" N LYS A1278 " --> pdb=" O GLU A1285 " (cutoff:3.500A) Processing sheet with id=AC4, first strand: chain 'A' and resid 1092 through 1098 removed outlier: 6.348A pdb=" N SER A1080 " --> pdb=" O GLN A1076 " (cutoff:3.500A) removed outlier: 4.881A pdb=" N GLN A1076 " --> pdb=" O SER A1080 " (cutoff:3.500A) removed outlier: 7.085A pdb=" N TYR A1082 " --> pdb=" O ILE A1074 " (cutoff:3.500A) removed outlier: 6.994A pdb=" N ARG A1122 " --> pdb=" O LEU A1046 " (cutoff:3.500A) removed outlier: 5.277A pdb=" N LEU A1046 " --> pdb=" O ARG A1122 " (cutoff:3.500A) removed outlier: 6.739A pdb=" N SER A1124 " --> pdb=" O LYS A1044 " (cutoff:3.500A) removed outlier: 5.378A pdb=" N LYS A1044 " --> pdb=" O SER A1124 " (cutoff:3.500A) removed outlier: 5.605A pdb=" N TYR A1126 " --> pdb=" O PHE A1042 " (cutoff:3.500A) removed outlier: 3.762A pdb=" N PHE A1042 " --> pdb=" O TYR A1126 " (cutoff:3.500A) removed outlier: 19.741A pdb=" N CYS A1143 " --> pdb=" O LYS A1223 " (cutoff:3.500A) removed outlier: 14.425A pdb=" N LYS A1223 " --> pdb=" O CYS A1143 " (cutoff:3.500A) removed outlier: 9.584A pdb=" N GLU A1145 " --> pdb=" O ASP A1221 " (cutoff:3.500A) removed outlier: 4.519A pdb=" N ASP A1221 " --> pdb=" O GLU A1145 " (cutoff:3.500A) removed outlier: 6.241A pdb=" N LEU A1222 " --> pdb=" O PHE A1005 " (cutoff:3.500A) removed outlier: 6.463A pdb=" N PHE A1005 " --> pdb=" O LEU A1222 " (cutoff:3.500A) removed outlier: 8.054A pdb=" N GLU A1224 " --> pdb=" O GLN A1003 " (cutoff:3.500A) removed outlier: 7.308A pdb=" N GLN A1003 " --> pdb=" O GLU A1224 " (cutoff:3.500A) removed outlier: 11.323A pdb=" N GLU A1226 " --> pdb=" O HIS A1001 " (cutoff:3.500A) removed outlier: 9.694A pdb=" N HIS A1001 " --> pdb=" O GLU A1226 " (cutoff:3.500A) removed outlier: 47.243A pdb=" N ILE A1002 " --> pdb=" O CYS A1119 " (cutoff:3.500A) removed outlier: 40.819A pdb=" N VAL A1121 " --> pdb=" O ILE A1002 " (cutoff:3.500A) removed outlier: 34.447A pdb=" N LYS A1004 " --> pdb=" O VAL A1121 " (cutoff:3.500A) removed outlier: 27.522A pdb=" N ASP A1123 " --> pdb=" O LYS A1004 " (cutoff:3.500A) removed outlier: 22.294A pdb=" N SER A1006 " --> pdb=" O ASP A1123 " (cutoff:3.500A) removed outlier: 16.384A pdb=" N ARG A1125 " --> pdb=" O SER A1006 " (cutoff:3.500A) removed outlier: 13.328A pdb=" N THR A1008 " --> pdb=" O ARG A1125 " (cutoff:3.500A) removed outlier: 7.340A pdb=" N GLN A1127 " --> pdb=" O THR A1008 " (cutoff:3.500A) Processing sheet with id=AC5, first strand: chain 'A' and resid 1156 through 1159 removed outlier: 7.316A pdb=" N GLU A1196 " --> pdb=" O LYS A1189 " (cutoff:3.500A) 456 hydrogen bonds defined for protein. 1218 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 49 hydrogen bonds 98 hydrogen bond angles 0 basepair planarities 20 basepair parallelities 23 stacking parallelities Total time for adding SS restraints: 7.21 Time building geometry restraints manager: 5.94 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: 4055 1.34 - 1.46: 2665 1.46 - 1.58: 6730 1.58 - 1.70: 121 1.70 - 1.82: 80 Bond restraints: 13651 Sorted by residual: bond pdb=" N PRO A1138 " pdb=" CA PRO A1138 " ideal model delta sigma weight residual 1.468 1.479 -0.011 7.90e-03 1.60e+04 2.04e+00 bond pdb=" CA GLU A1196 " pdb=" CB GLU A1196 " ideal model delta sigma weight residual 1.522 1.530 -0.008 7.00e-03 2.04e+04 1.17e+00 bond pdb=" N PRO A1153 " pdb=" CA PRO A1153 " ideal model delta sigma weight residual 1.469 1.481 -0.013 1.28e-02 6.10e+03 9.94e-01 bond pdb=" N ASP A1315 " pdb=" CA ASP A1315 " ideal model delta sigma weight residual 1.457 1.470 -0.013 1.29e-02 6.01e+03 9.72e-01 bond pdb=" CB PRO A1368 " pdb=" CG PRO A1368 " ideal model delta sigma weight residual 1.492 1.538 -0.046 5.00e-02 4.00e+02 8.51e-01 ... (remaining 13646 not shown) Histogram of bond angle deviations from ideal: 99.84 - 106.67: 776 106.67 - 113.50: 7316 113.50 - 120.32: 4881 120.32 - 127.15: 5426 127.15 - 133.98: 297 Bond angle restraints: 18696 Sorted by residual: angle pdb=" C2' U R 1 " pdb=" C1' U R 1 " pdb=" N1 U R 1 " ideal model delta sigma weight residual 112.00 117.12 -5.12 1.50e+00 4.44e-01 1.16e+01 angle pdb=" O4' G R 19 " pdb=" C1' G R 19 " pdb=" N9 G R 19 " ideal model delta sigma weight residual 108.20 113.04 -4.84 1.50e+00 4.44e-01 1.04e+01 angle pdb=" N ILE A1424 " pdb=" CA ILE A1424 " pdb=" C ILE A1424 " ideal model delta sigma weight residual 112.80 109.62 3.18 1.15e+00 7.56e-01 7.65e+00 angle pdb=" N GLY A 58 " pdb=" CA GLY A 58 " pdb=" C GLY A 58 " ideal model delta sigma weight residual 114.66 111.27 3.39 1.24e+00 6.50e-01 7.46e+00 angle pdb=" N GLY A 641 " pdb=" CA GLY A 641 " pdb=" C GLY A 641 " ideal model delta sigma weight residual 112.14 115.47 -3.33 1.29e+00 6.01e-01 6.66e+00 ... (remaining 18691 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 23.25: 7818 23.25 - 46.51: 248 46.51 - 69.76: 39 69.76 - 93.01: 6 93.01 - 116.27: 2 Dihedral angle restraints: 8113 sinusoidal: 3872 harmonic: 4241 Sorted by residual: dihedral pdb=" O4' A R 38 " pdb=" C1' A R 38 " pdb=" N9 A R 38 " pdb=" C4 A R 38 " ideal model delta sinusoidal sigma weight residual 68.00 137.16 -69.16 1 1.70e+01 3.46e-03 2.14e+01 dihedral pdb=" CA ASN A 640 " pdb=" C ASN A 640 " pdb=" N GLY A 641 " pdb=" CA GLY A 641 " ideal model delta harmonic sigma weight residual 180.00 -160.54 -19.46 0 5.00e+00 4.00e-02 1.51e+01 dihedral pdb=" O4' A R 12 " pdb=" C1' A R 12 " pdb=" N9 A R 12 " pdb=" C4 A R 12 " ideal model delta sinusoidal sigma weight residual -106.00 -156.28 50.28 1 1.70e+01 3.46e-03 1.20e+01 ... (remaining 8110 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.049: 1667 0.049 - 0.099: 278 0.099 - 0.148: 63 0.148 - 0.197: 0 0.197 - 0.247: 2 Chirality restraints: 2010 Sorted by residual: chirality pdb=" C1' G R 19 " pdb=" O4' G R 19 " pdb=" C2' G R 19 " pdb=" N9 G R 19 " both_signs ideal model delta sigma weight residual False 2.44 2.19 0.25 2.00e-01 2.50e+01 1.52e+00 chirality pdb=" C1' U R 1 " pdb=" O4' U R 1 " pdb=" C2' U R 1 " pdb=" N1 U R 1 " both_signs ideal model delta sigma weight residual False 2.47 2.26 0.20 2.00e-01 2.50e+01 1.04e+00 chirality pdb=" CA ILE A 628 " pdb=" N ILE A 628 " pdb=" C ILE A 628 " pdb=" CB ILE A 628 " both_signs ideal model delta sigma weight residual False 2.43 2.58 -0.14 2.00e-01 2.50e+01 5.08e-01 ... (remaining 2007 not shown) Planarity restraints: 2203 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C1' U R 1 " 0.038 2.00e-02 2.50e+03 1.98e-02 8.78e+00 pdb=" N1 U R 1 " -0.043 2.00e-02 2.50e+03 pdb=" C2 U R 1 " -0.003 2.00e-02 2.50e+03 pdb=" O2 U R 1 " -0.004 2.00e-02 2.50e+03 pdb=" N3 U R 1 " 0.004 2.00e-02 2.50e+03 pdb=" C4 U R 1 " 0.005 2.00e-02 2.50e+03 pdb=" O4 U R 1 " 0.010 2.00e-02 2.50e+03 pdb=" C5 U R 1 " 0.000 2.00e-02 2.50e+03 pdb=" C6 U R 1 " -0.007 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C THR A 903 " 0.046 5.00e-02 4.00e+02 6.99e-02 7.81e+00 pdb=" N PRO A 904 " -0.121 5.00e-02 4.00e+02 pdb=" CA PRO A 904 " 0.036 5.00e-02 4.00e+02 pdb=" CD PRO A 904 " 0.039 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C1' G R 19 " 0.028 2.00e-02 2.50e+03 1.33e-02 5.28e+00 pdb=" N9 G R 19 " -0.035 2.00e-02 2.50e+03 pdb=" C8 G R 19 " 0.000 2.00e-02 2.50e+03 pdb=" N7 G R 19 " 0.001 2.00e-02 2.50e+03 pdb=" C5 G R 19 " 0.000 2.00e-02 2.50e+03 pdb=" C6 G R 19 " 0.001 2.00e-02 2.50e+03 pdb=" O6 G R 19 " 0.008 2.00e-02 2.50e+03 pdb=" N1 G R 19 " 0.002 2.00e-02 2.50e+03 pdb=" C2 G R 19 " -0.002 2.00e-02 2.50e+03 pdb=" N2 G R 19 " -0.002 2.00e-02 2.50e+03 pdb=" N3 G R 19 " -0.001 2.00e-02 2.50e+03 pdb=" C4 G R 19 " -0.001 2.00e-02 2.50e+03 ... (remaining 2200 not shown) Histogram of nonbonded interaction distances: 2.24 - 2.77: 2356 2.77 - 3.30: 12038 3.30 - 3.84: 22979 3.84 - 4.37: 28609 4.37 - 4.90: 46429 Nonbonded interactions: 112411 Sorted by model distance: nonbonded pdb=" OP1 U R 16 " pdb=" OG SER A 564 " model vdw 2.239 2.440 nonbonded pdb=" OH TYR A1081 " pdb=" OE1 GLN A1105 " model vdw 2.262 2.440 nonbonded pdb=" NE2 HIS A 931 " pdb=" OG1 THR A1444 " model vdw 2.263 2.520 nonbonded pdb=" O ALA A 288 " pdb=" OG SER A 291 " model vdw 2.271 2.440 nonbonded pdb=" OP1 U R 25 " pdb=" OG SER A 948 " model vdw 2.278 2.440 ... (remaining 112406 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 Set scattering table Set to: electron Number of scattering types: 6 Type Number sf(0) Gaussians Zn 4 6.06 5 P 61 5.49 5 S 56 5.16 5 C 8130 2.51 5 N 2348 2.21 5 O 2618 1.98 5 sf(0) = scattering factor at diffraction angle 0. Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as 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: 10.930 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: 28.970 Check model and map are aligned: 0.180 Process input model: 44.570 Find NCS groups from input model: 0.320 Set up NCS constraints: 0.050 Set refine NCS operators: 0.000 Set scattering table: 0.120 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:11.300 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 96.450 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7962 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.046 13651 Z= 0.150 Angle : 0.509 5.117 18696 Z= 0.285 Chirality : 0.039 0.247 2010 Planarity : 0.004 0.070 2203 Dihedral : 11.560 116.269 5377 Min Nonbonded Distance : 2.239 Molprobity Statistics. All-atom Clashscore : 7.07 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.72 % Favored : 98.28 % Rotamer Outliers : 0.55 % 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.28 (0.22), residues: 1453 helix: 0.93 (0.28), residues: 353 sheet: 0.51 (0.29), residues: 290 loop : -0.20 (0.22), residues: 810 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 250 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 243 time to evaluate : 1.462 Fit side-chains revert: symmetry clash outliers start: 7 outliers final: 4 residues processed: 249 average time/residue: 1.6276 time to fit residues: 435.7171 Evaluate side-chains 132 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 128 time to evaluate : 1.511 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 4 residues processed: 0 time to fit residues: 2.1223 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 152 random chunks: chunk 128 optimal weight: 0.4980 chunk 115 optimal weight: 1.9990 chunk 63 optimal weight: 0.1980 chunk 39 optimal weight: 3.9990 chunk 77 optimal weight: 0.6980 chunk 61 optimal weight: 1.9990 chunk 119 optimal weight: 0.8980 chunk 46 optimal weight: 3.9990 chunk 72 optimal weight: 3.9990 chunk 88 optimal weight: 3.9990 chunk 138 optimal weight: 0.9990 overall best weight: 0.6582 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 21 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 30 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 170 GLN A 306 HIS A 626 HIS A1246 ASN A1476 GLN ** A1483 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7991 moved from start: 0.1343 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.044 13651 Z= 0.170 Angle : 0.519 6.000 18696 Z= 0.278 Chirality : 0.041 0.291 2010 Planarity : 0.004 0.066 2203 Dihedral : 10.919 114.350 2345 Min Nonbonded Distance : 2.155 Molprobity Statistics. All-atom Clashscore : 7.89 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.13 % Favored : 97.87 % Rotamer Outliers : 2.27 % 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.55 (0.22), residues: 1453 helix: 1.30 (0.28), residues: 351 sheet: 0.38 (0.28), residues: 304 loop : 0.03 (0.23), residues: 798 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 165 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 29 poor density : 136 time to evaluate : 1.597 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 29 outliers final: 13 residues processed: 154 average time/residue: 1.5019 time to fit residues: 251.6279 Evaluate side-chains 139 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 126 time to evaluate : 1.494 Switching outliers to nearest non-outliers outliers start: 13 outliers final: 8 residues processed: 6 average time/residue: 0.3361 time to fit residues: 4.5725 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 152 random chunks: chunk 76 optimal weight: 2.9990 chunk 42 optimal weight: 0.6980 chunk 114 optimal weight: 4.9990 chunk 93 optimal weight: 4.9990 chunk 38 optimal weight: 0.9980 chunk 138 optimal weight: 1.9990 chunk 149 optimal weight: 3.9990 chunk 123 optimal weight: 0.8980 chunk 137 optimal weight: 2.9990 chunk 47 optimal weight: 6.9990 chunk 110 optimal weight: 0.7980 overall best weight: 1.0782 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 21 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 184 HIS A1084 ASN ** A1483 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8027 moved from start: 0.1579 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.047 13651 Z= 0.221 Angle : 0.544 6.213 18696 Z= 0.291 Chirality : 0.042 0.304 2010 Planarity : 0.005 0.062 2203 Dihedral : 10.913 114.139 2345 Min Nonbonded Distance : 2.211 Molprobity Statistics. All-atom Clashscore : 8.28 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.13 % Favored : 97.87 % Rotamer Outliers : 3.28 % 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.47 (0.22), residues: 1453 helix: 1.23 (0.28), residues: 350 sheet: 0.32 (0.27), residues: 315 loop : -0.03 (0.23), residues: 788 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 176 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 42 poor density : 134 time to evaluate : 1.559 Fit side-chains revert: symmetry clash outliers start: 42 outliers final: 17 residues processed: 159 average time/residue: 1.4307 time to fit residues: 248.2198 Evaluate side-chains 137 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 17 poor density : 120 time to evaluate : 1.366 Switching outliers to nearest non-outliers outliers start: 17 outliers final: 11 residues processed: 6 average time/residue: 0.1752 time to fit residues: 3.6342 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 152 random chunks: chunk 136 optimal weight: 0.8980 chunk 103 optimal weight: 10.0000 chunk 71 optimal weight: 0.8980 chunk 15 optimal weight: 3.9990 chunk 65 optimal weight: 2.9990 chunk 92 optimal weight: 0.7980 chunk 138 optimal weight: 0.5980 chunk 146 optimal weight: 0.8980 chunk 72 optimal weight: 0.5980 chunk 131 optimal weight: 2.9990 chunk 39 optimal weight: 3.9990 overall best weight: 0.7580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 21 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1483 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8007 moved from start: 0.1715 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.043 13651 Z= 0.176 Angle : 0.513 5.933 18696 Z= 0.276 Chirality : 0.041 0.290 2010 Planarity : 0.004 0.061 2203 Dihedral : 10.862 114.372 2345 Min Nonbonded Distance : 2.208 Molprobity Statistics. All-atom Clashscore : 8.08 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.93 % Favored : 98.07 % Rotamer Outliers : 2.97 % 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.52 (0.22), residues: 1453 helix: 1.31 (0.28), residues: 351 sheet: 0.31 (0.27), residues: 325 loop : -0.00 (0.23), residues: 777 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 163 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 38 poor density : 125 time to evaluate : 1.671 Fit side-chains revert: symmetry clash outliers start: 38 outliers final: 18 residues processed: 151 average time/residue: 1.5154 time to fit residues: 249.3048 Evaluate side-chains 135 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 117 time to evaluate : 1.569 Switching outliers to nearest non-outliers outliers start: 18 outliers final: 15 residues processed: 3 average time/residue: 0.1584 time to fit residues: 2.8468 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 152 random chunks: chunk 122 optimal weight: 0.9990 chunk 83 optimal weight: 0.6980 chunk 2 optimal weight: 4.9990 chunk 109 optimal weight: 2.9990 chunk 60 optimal weight: 3.9990 chunk 125 optimal weight: 4.9990 chunk 101 optimal weight: 2.9990 chunk 0 optimal weight: 5.9990 chunk 75 optimal weight: 0.5980 chunk 131 optimal weight: 0.9980 chunk 37 optimal weight: 0.9980 overall best weight: 0.8582 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 21 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1483 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8017 moved from start: 0.1829 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.043 13651 Z= 0.189 Angle : 0.520 5.997 18696 Z= 0.278 Chirality : 0.041 0.293 2010 Planarity : 0.004 0.061 2203 Dihedral : 10.845 114.296 2345 Min Nonbonded Distance : 2.218 Molprobity Statistics. All-atom Clashscore : 7.81 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.20 % Favored : 97.80 % Rotamer Outliers : 3.28 % 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.49 (0.22), residues: 1453 helix: 1.32 (0.28), residues: 352 sheet: 0.30 (0.27), residues: 325 loop : -0.04 (0.23), residues: 776 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 164 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 42 poor density : 122 time to evaluate : 1.579 Fit side-chains revert: symmetry clash outliers start: 42 outliers final: 24 residues processed: 151 average time/residue: 1.4079 time to fit residues: 231.8825 Evaluate side-chains 140 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 24 poor density : 116 time to evaluate : 1.450 Switching outliers to nearest non-outliers outliers start: 24 outliers final: 20 residues processed: 4 average time/residue: 0.1862 time to fit residues: 3.2089 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 152 random chunks: chunk 49 optimal weight: 4.9990 chunk 132 optimal weight: 0.0000 chunk 29 optimal weight: 0.7980 chunk 86 optimal weight: 0.8980 chunk 36 optimal weight: 2.9990 chunk 147 optimal weight: 0.4980 chunk 122 optimal weight: 0.7980 chunk 68 optimal weight: 3.9990 chunk 12 optimal weight: 1.9990 chunk 48 optimal weight: 2.9990 chunk 77 optimal weight: 0.9990 overall best weight: 0.5984 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 21 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1483 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7995 moved from start: 0.1914 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.042 13651 Z= 0.156 Angle : 0.496 5.839 18696 Z= 0.266 Chirality : 0.040 0.284 2010 Planarity : 0.004 0.060 2203 Dihedral : 10.795 114.602 2345 Min Nonbonded Distance : 2.222 Molprobity Statistics. All-atom Clashscore : 7.73 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.20 % Favored : 97.80 % Rotamer Outliers : 2.81 % 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.22), residues: 1453 helix: 1.47 (0.28), residues: 352 sheet: 0.40 (0.28), residues: 309 loop : -0.05 (0.23), residues: 792 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 161 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 36 poor density : 125 time to evaluate : 1.475 Fit side-chains revert: symmetry clash outliers start: 36 outliers final: 23 residues processed: 153 average time/residue: 1.3198 time to fit residues: 219.2305 Evaluate side-chains 136 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 23 poor density : 113 time to evaluate : 0.954 Switching outliers to nearest non-outliers outliers start: 23 outliers final: 18 residues processed: 5 average time/residue: 0.1495 time to fit residues: 2.4615 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 152 random chunks: chunk 141 optimal weight: 2.9990 chunk 16 optimal weight: 5.9990 chunk 83 optimal weight: 0.8980 chunk 107 optimal weight: 4.9990 chunk 123 optimal weight: 0.7980 chunk 82 optimal weight: 0.9980 chunk 146 optimal weight: 0.8980 chunk 91 optimal weight: 3.9990 chunk 89 optimal weight: 2.9990 chunk 67 optimal weight: 5.9990 chunk 90 optimal weight: 0.6980 overall best weight: 0.8580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 21 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1483 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8016 moved from start: 0.2020 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.043 13651 Z= 0.190 Angle : 0.521 6.048 18696 Z= 0.279 Chirality : 0.041 0.296 2010 Planarity : 0.005 0.060 2203 Dihedral : 10.787 114.364 2345 Min Nonbonded Distance : 2.218 Molprobity Statistics. All-atom Clashscore : 8.20 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.27 % Favored : 97.73 % Rotamer Outliers : 2.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.51 (0.22), residues: 1453 helix: 1.39 (0.28), residues: 352 sheet: 0.35 (0.28), residues: 313 loop : -0.07 (0.23), residues: 788 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2906 Ramachandran restraints generated. 1453 Oldfield, 0 Emsley, 1453 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 153 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 37 poor density : 116 time to evaluate : 1.714 Fit side-chains revert: symmetry clash outliers start: 37 outliers final: 27 residues processed: 144 average time/residue: 1.4655 time to fit residues: 230.6429 Evaluate side-chains 143 residues out of total 1284 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 27 poor density : 116 time to evaluate : 1.781 Switching outliers to nearest non-outliers outliers start: 27 outliers final: 21 residues processed: 6 average time/residue: 0.4141 time to fit residues: 5.3257 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 152 random chunks: chunk 58 optimal weight: 4.9990 chunk 87 optimal weight: 0.0270 chunk 44 optimal weight: 0.0770 chunk 28 optimal weight: 0.6980 chunk 93 optimal weight: 2.9990 chunk 99 optimal weight: 5.9990 chunk 72 optimal weight: 2.9990 chunk 13 optimal weight: 0.5980 chunk 115 optimal weight: 1.9990 chunk 133 optimal weight: 2.9990 chunk 140 optimal weight: 0.8980 overall best weight: 0.4596 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Sorry: Reduce crashed with command 'molprobity.reduce -quiet -trim -'. Dumping stdin to file 'reduce_failure.pdb'. Return code: -15 Dumping stderr: