Starting phenix.real_space_refine (version: dev) on Mon Feb 20 15:21:59 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6nma_0445/02_2023/6nma_0445.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6nma_0445/02_2023/6nma_0445.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.38 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6nma_0445/02_2023/6nma_0445.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6nma_0445/02_2023/6nma_0445.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6nma_0445/02_2023/6nma_0445.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6nma_0445/02_2023/6nma_0445.pdb" } resolution = 3.38 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.000 sd= 0.004 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Process input model Symmetric amino acids flipped Residue "A ARG 217": "NH1" <-> "NH2" Residue "A GLU 227": "OE1" <-> "OE2" Residue "B GLU 45": "OE1" <-> "OE2" Residue "B GLU 125": "OE1" <-> "OE2" Residue "B GLU 159": "OE1" <-> "OE2" Residue "B GLU 165": "OE1" <-> "OE2" Residue "B ARG 174": "NH1" <-> "NH2" Residue "B GLU 178": "OE1" <-> "OE2" Residue "B GLU 204": "OE1" <-> "OE2" Residue "B GLU 223": "OE1" <-> "OE2" Residue "B GLU 247": "OE1" <-> "OE2" Residue "B GLU 250": "OE1" <-> "OE2" Residue "B TYR 294": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B GLU 365": "OE1" <-> "OE2" Residue "B TYR 366": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B GLU 379": "OE1" <-> "OE2" Residue "B GLU 402": "OE1" <-> "OE2" Residue "B GLU 412": "OE1" <-> "OE2" Residue "B GLU 484": "OE1" <-> "OE2" Residue "B TYR 504": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 548": "NH1" <-> "NH2" Residue "B ARG 737": "NH1" <-> "NH2" Residue "B GLU 743": "OE1" <-> "OE2" Residue "B ARG 747": "NH1" <-> "NH2" Residue "B TYR 872": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 883": "NH1" <-> "NH2" Residue "B PHE 884": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 887": "NH1" <-> "NH2" Residue "B GLU 981": "OE1" <-> "OE2" Residue "B GLU 1097": "OE1" <-> "OE2" Residue "B GLU 1117": "OE1" <-> "OE2" Residue "B GLU 1171": "OE1" <-> "OE2" Residue "B GLU 1202": "OE1" <-> "OE2" Residue "C GLU 174": "OE1" <-> "OE2" Residue "C GLU 176": "OE1" <-> "OE2" Residue "C GLU 184": "OE1" <-> "OE2" Residue "C ARG 217": "NH1" <-> "NH2" Time to flip residues: 0.04s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4877/modules/chem_data/mon_lib" Total number of atoms: 12404 Number of models: 1 Model: "" Number of chains: 6 Chain: "A" Number of atoms: 951 Number of conformers: 1 Conformer: "" Number of residues, atoms: 115, 951 Classifications: {'peptide': 115} Link IDs: {'PTRANS': 7, 'TRANS': 107} Chain: "B" Number of atoms: 9974 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1206, 9974 Classifications: {'peptide': 1206} Link IDs: {'PTRANS': 25, 'TRANS': 1180} Chain breaks: 4 Chain: "G" Number of atoms: 535 Number of conformers: 1 Conformer: "" Number of residues, atoms: 25, 535 Classifications: {'RNA': 25} Modifications used: {'rna2p_pur': 6, 'rna2p_pyr': 3, 'rna3p_pur': 9, 'rna3p_pyr': 7} Link IDs: {'rna2p': 9, 'rna3p': 15} Chain: "C" Number of atoms: 942 Number of conformers: 1 Conformer: "" Number of residues, atoms: 114, 942 Classifications: {'peptide': 114} Link IDs: {'PTRANS': 7, 'TRANS': 106} Chain: "A" Number of atoms: 1 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 1 Unusual residues: {' MG': 1} Classifications: {'undetermined': 1} Chain: "G" Number of atoms: 1 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 1 Unusual residues: {' MG': 1} Classifications: {'undetermined': 1} Time building chain proxies: 7.45, per 1000 atoms: 0.60 Number of scatterers: 12404 At special positions: 0 Unit cell: (126.854, 108.732, 127.92, 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) S 41 16.00 P 25 15.00 Mg 2 11.99 O 2417 8.00 N 2049 7.00 C 7870 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.30 Conformation dependent library (CDL) restraints added in 1.6 seconds 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2752 Finding SS restraints... Secondary structure from input PDB file: 61 helices and 12 sheets defined 45.4% alpha, 9.7% beta 6 base pairs and 7 stacking pairs defined. Time for finding SS restraints: 3.41 Creating SS restraints... Processing helix chain 'A' and resid 120 through 125 removed outlier: 4.091A pdb=" N VAL A 124 " --> pdb=" O MET A 120 " (cutoff:3.500A) removed outlier: 3.998A pdb=" N GLN A 125 " --> pdb=" O TYR A 121 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 120 through 125' Processing helix chain 'A' and resid 126 through 129 removed outlier: 5.018A pdb=" N ARG A 129 " --> pdb=" O ALA A 126 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 126 through 129' Processing helix chain 'A' and resid 130 through 148 removed outlier: 3.715A pdb=" N LYS A 134 " --> pdb=" O ALA A 130 " (cutoff:3.500A) removed outlier: 3.631A pdb=" N LEU A 138 " --> pdb=" O LYS A 134 " (cutoff:3.500A) removed outlier: 3.658A pdb=" N ASN A 145 " --> pdb=" O GLU A 141 " (cutoff:3.500A) removed outlier: 3.809A pdb=" N ILE A 146 " --> pdb=" O GLN A 142 " (cutoff:3.500A) removed outlier: 3.551A pdb=" N LYS A 148 " --> pdb=" O TYR A 144 " (cutoff:3.500A) Processing helix chain 'A' and resid 151 through 155 removed outlier: 3.636A pdb=" N ARG A 155 " --> pdb=" O TYR A 152 " (cutoff:3.500A) Processing helix chain 'A' and resid 200 through 204 Processing helix chain 'A' and resid 225 through 231 Processing helix chain 'B' and resid 26 through 35 removed outlier: 3.688A pdb=" N ASN B 30 " --> pdb=" O LYS B 26 " (cutoff:3.500A) Processing helix chain 'B' and resid 37 through 66 removed outlier: 3.638A pdb=" N ALA B 44 " --> pdb=" O ASP B 40 " (cutoff:3.500A) removed outlier: 5.686A pdb=" N LYS B 48 " --> pdb=" O ALA B 44 " (cutoff:3.500A) removed outlier: 5.160A pdb=" N GLY B 49 " --> pdb=" O GLU B 45 " (cutoff:3.500A) Processing helix chain 'B' and resid 73 through 82 removed outlier: 4.385A pdb=" N SER B 79 " --> pdb=" O ASN B 75 " (cutoff:3.500A) removed outlier: 4.080A pdb=" N ARG B 82 " --> pdb=" O ILE B 78 " (cutoff:3.500A) Processing helix chain 'B' and resid 90 through 104 removed outlier: 3.984A pdb=" N LEU B 94 " --> pdb=" O GLU B 90 " (cutoff:3.500A) Processing helix chain 'B' and resid 106 through 112 Processing helix chain 'B' and resid 121 through 126 removed outlier: 3.644A pdb=" N GLU B 125 " --> pdb=" O LYS B 121 " (cutoff:3.500A) Processing helix chain 'B' and resid 126 through 131 removed outlier: 3.889A pdb=" N GLU B 130 " --> pdb=" O THR B 126 " (cutoff:3.500A) Processing helix chain 'B' and resid 135 through 144 removed outlier: 4.038A pdb=" N VAL B 141 " --> pdb=" O GLU B 137 " (cutoff:3.500A) removed outlier: 3.651A pdb=" N PHE B 144 " --> pdb=" O LEU B 140 " (cutoff:3.500A) Processing helix chain 'B' and resid 147 through 150 Processing helix chain 'B' and resid 151 through 158 Processing helix chain 'B' and resid 173 through 178 Processing helix chain 'B' and resid 179 through 194 removed outlier: 3.645A pdb=" N VAL B 193 " --> pdb=" O ILE B 189 " (cutoff:3.500A) Processing helix chain 'B' and resid 198 through 209 removed outlier: 3.779A pdb=" N VAL B 202 " --> pdb=" O ASP B 198 " (cutoff:3.500A) removed outlier: 5.397A pdb=" N GLU B 204 " --> pdb=" O HIS B 200 " (cutoff:3.500A) removed outlier: 5.379A pdb=" N ILE B 205 " --> pdb=" O GLU B 201 " (cutoff:3.500A) removed outlier: 3.960A pdb=" N LYS B 208 " --> pdb=" O GLU B 204 " (cutoff:3.500A) Processing helix chain 'B' and resid 223 through 227 Processing helix chain 'B' and resid 230 through 242 removed outlier: 3.522A pdb=" N ILE B 234 " --> pdb=" O THR B 230 " (cutoff:3.500A) Processing helix chain 'B' and resid 254 through 262 Processing helix chain 'B' and resid 263 through 266 Processing helix chain 'B' and resid 296 through 309 Processing helix chain 'B' and resid 312 through 327 Processing helix chain 'B' and resid 341 through 351 removed outlier: 4.134A pdb=" N SER B 345 " --> pdb=" O GLY B 341 " (cutoff:3.500A) removed outlier: 3.891A pdb=" N THR B 346 " --> pdb=" O PRO B 342 " (cutoff:3.500A) removed outlier: 3.891A pdb=" N LYS B 349 " --> pdb=" O SER B 345 " (cutoff:3.500A) Processing helix chain 'B' and resid 356 through 370 removed outlier: 3.561A pdb=" N ASP B 360 " --> pdb=" O ASN B 356 " (cutoff:3.500A) Processing helix chain 'B' and resid 378 through 392 Processing helix chain 'B' and resid 397 through 403 Processing helix chain 'B' and resid 410 through 431 Processing helix chain 'B' and resid 432 through 437 removed outlier: 3.809A pdb=" N PHE B 436 " --> pdb=" O SER B 432 " (cutoff:3.500A) removed outlier: 3.683A pdb=" N ASP B 437 " --> pdb=" O GLU B 433 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 432 through 437' Processing helix chain 'B' and resid 452 through 471 removed outlier: 3.721A pdb=" N SER B 462 " --> pdb=" O ASP B 458 " (cutoff:3.500A) removed outlier: 3.908A pdb=" N TYR B 469 " --> pdb=" O SER B 465 " (cutoff:3.500A) Processing helix chain 'B' and resid 483 through 498 removed outlier: 4.738A pdb=" N ASP B 489 " --> pdb=" O SER B 485 " (cutoff:3.500A) removed outlier: 3.970A pdb=" N ILE B 496 " --> pdb=" O LEU B 492 " (cutoff:3.500A) Processing helix chain 'B' and resid 499 through 512 removed outlier: 3.597A pdb=" N ARG B 508 " --> pdb=" O TYR B 504 " (cutoff:3.500A) removed outlier: 3.764A pdb=" N ASN B 509 " --> pdb=" O ASP B 505 " (cutoff:3.500A) removed outlier: 4.467A pdb=" N TYR B 510 " --> pdb=" O ALA B 506 " (cutoff:3.500A) Processing helix chain 'B' and resid 538 through 543 Processing helix chain 'B' and resid 588 through 598 Proline residue: B 594 - end of helix Processing helix chain 'B' and resid 610 through 619 Processing helix chain 'B' and resid 629 through 644 removed outlier: 4.005A pdb=" N HIS B 633 " --> pdb=" O LEU B 629 " (cutoff:3.500A) Processing helix chain 'B' and resid 646 through 651 Processing helix chain 'B' and resid 666 through 677 Processing helix chain 'B' and resid 687 through 696 Processing helix chain 'B' and resid 718 through 727 removed outlier: 3.516A pdb=" N MET B 722 " --> pdb=" O ASN B 718 " (cutoff:3.500A) Processing helix chain 'B' and resid 786 through 790 Processing helix chain 'B' and resid 811 through 822 removed outlier: 3.941A pdb=" N ARG B 817 " --> pdb=" O ASN B 813 " (cutoff:3.500A) Processing helix chain 'B' and resid 873 through 885 removed outlier: 3.808A pdb=" N GLU B 885 " --> pdb=" O LYS B 881 " (cutoff:3.500A) Processing helix chain 'B' and resid 894 through 900 Processing helix chain 'B' and resid 907 through 919 removed outlier: 3.813A pdb=" N LEU B 914 " --> pdb=" O LYS B 910 " (cutoff:3.500A) removed outlier: 4.152A pdb=" N GLU B 916 " --> pdb=" O CYS B 912 " (cutoff:3.500A) Processing helix chain 'B' and resid 927 through 930 removed outlier: 4.265A pdb=" N GLY B 930 " --> pdb=" O LEU B 927 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 927 through 930' Processing helix chain 'B' and resid 939 through 948 removed outlier: 3.819A pdb=" N LYS B 945 " --> pdb=" O GLN B 941 " (cutoff:3.500A) Processing helix chain 'B' and resid 950 through 955 Processing helix chain 'B' and resid 1020 through 1031 removed outlier: 3.645A pdb=" N SER B1024 " --> pdb=" O SER B1020 " (cutoff:3.500A) removed outlier: 3.658A pdb=" N PHE B1031 " --> pdb=" O PHE B1027 " (cutoff:3.500A) Processing helix chain 'B' and resid 1091 through 1102 removed outlier: 3.597A pdb=" N TYR B1102 " --> pdb=" O LEU B1098 " (cutoff:3.500A) Processing helix chain 'B' and resid 1111 through 1116 Processing helix chain 'B' and resid 1120 through 1135 removed outlier: 4.468A pdb=" N TYR B1124 " --> pdb=" O ASP B1120 " (cutoff:3.500A) removed outlier: 3.950A pdb=" N SER B1126 " --> pdb=" O ALA B1122 " (cutoff:3.500A) removed outlier: 3.593A pdb=" N MET B1128 " --> pdb=" O TYR B1124 " (cutoff:3.500A) removed outlier: 4.164A pdb=" N ALA B1129 " --> pdb=" O SER B1125 " (cutoff:3.500A) removed outlier: 3.556A pdb=" N LEU B1135 " --> pdb=" O MET B1131 " (cutoff:3.500A) Processing helix chain 'B' and resid 1163 through 1167 Processing helix chain 'B' and resid 1180 through 1199 removed outlier: 3.671A pdb=" N LYS B1190 " --> pdb=" O ASN B1186 " (cutoff:3.500A) removed outlier: 3.847A pdb=" N VAL B1191 " --> pdb=" O ILE B1187 " (cutoff:3.500A) Processing helix chain 'B' and resid 1205 through 1209 Processing helix chain 'B' and resid 1214 through 1225 Processing helix chain 'C' and resid 123 through 149 removed outlier: 3.733A pdb=" N ALA C 130 " --> pdb=" O ALA C 126 " (cutoff:3.500A) removed outlier: 3.712A pdb=" N ILE C 133 " --> pdb=" O ARG C 129 " (cutoff:3.500A) removed outlier: 4.046A pdb=" N ASN C 145 " --> pdb=" O GLU C 141 " (cutoff:3.500A) removed outlier: 3.962A pdb=" N ILE C 146 " --> pdb=" O GLN C 142 " (cutoff:3.500A) Processing helix chain 'C' and resid 151 through 155 removed outlier: 3.980A pdb=" N ARG C 155 " --> pdb=" O TYR C 152 " (cutoff:3.500A) Processing helix chain 'C' and resid 225 through 231 Processing sheet with id=AA1, first strand: chain 'A' and resid 158 through 161 Processing sheet with id=AA2, first strand: chain 'A' and resid 196 through 197 removed outlier: 3.589A pdb=" N VAL A 207 " --> pdb=" O PHE A 196 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'B' and resid 520 through 522 removed outlier: 4.168A pdb=" N TYR B 794 " --> pdb=" O ALA B 21 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'B' and resid 520 through 522 removed outlier: 3.666A pdb=" N ARG B 737 " --> pdb=" O ALA B 801 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'B' and resid 337 through 339 Processing sheet with id=AA6, first strand: chain 'B' and resid 552 through 556 removed outlier: 4.923A pdb=" N ILE B 546 " --> pdb=" O ILE B 581 " (cutoff:3.500A) removed outlier: 3.506A pdb=" N ILE B 581 " --> pdb=" O ILE B 546 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 758 through 759 removed outlier: 3.551A pdb=" N VAL B 783 " --> pdb=" O HIS B 759 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'B' and resid 850 through 855 removed outlier: 3.539A pdb=" N GLU B 852 " --> pdb=" O VAL B 843 " (cutoff:3.500A) removed outlier: 4.162A pdb=" N GLU B 925 " --> pdb=" O ILE B 831 " (cutoff:3.500A) removed outlier: 6.230A pdb=" N ILE B 922 " --> pdb=" O PHE B 994 " (cutoff:3.500A) removed outlier: 6.921A pdb=" N ILE B 996 " --> pdb=" O ILE B 922 " (cutoff:3.500A) removed outlier: 5.918A pdb=" N LEU B 924 " --> pdb=" O ILE B 996 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'B' and resid 858 through 863 Processing sheet with id=AB1, first strand: chain 'B' and resid 1033 through 1037 removed outlier: 3.544A pdb=" N LEU B1042 " --> pdb=" O VAL B1037 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'B' and resid 1070 through 1073 Processing sheet with id=AB3, first strand: chain 'C' and resid 206 through 207 removed outlier: 6.467A pdb=" N CYS C 195 " --> pdb=" O CYS C 185 " (cutoff:3.500A) removed outlier: 4.818A pdb=" N CYS C 185 " --> pdb=" O CYS C 195 " (cutoff:3.500A) removed outlier: 5.747A pdb=" N TYR C 197 " --> pdb=" O GLY C 183 " (cutoff:3.500A) removed outlier: 3.764A pdb=" N ASN C 157 " --> pdb=" O VAL C 213 " (cutoff:3.500A) 442 hydrogen bonds defined for protein. 1236 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 13 hydrogen bonds 22 hydrogen bond angles 0 basepair planarities 6 basepair parallelities 7 stacking parallelities Total time for adding SS restraints: 4.64 Time building geometry restraints manager: 5.53 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.20 - 1.32: 2341 1.32 - 1.44: 3250 1.44 - 1.57: 7003 1.57 - 1.69: 49 1.69 - 1.81: 67 Bond restraints: 12710 Sorted by residual: bond pdb=" N SER B 855 " pdb=" CA SER B 855 " ideal model delta sigma weight residual 1.454 1.488 -0.034 1.19e-02 7.06e+03 8.18e+00 bond pdb=" N GLU B 882 " pdb=" CA GLU B 882 " ideal model delta sigma weight residual 1.459 1.492 -0.033 1.20e-02 6.94e+03 7.65e+00 bond pdb=" CA GLN B 975 " pdb=" CB GLN B 975 " ideal model delta sigma weight residual 1.527 1.459 0.068 2.48e-02 1.63e+03 7.48e+00 bond pdb=" N LEU B 856 " pdb=" CA LEU B 856 " ideal model delta sigma weight residual 1.456 1.489 -0.034 1.26e-02 6.30e+03 7.23e+00 bond pdb=" N ARG B1112 " pdb=" CA ARG B1112 " ideal model delta sigma weight residual 1.459 1.488 -0.030 1.16e-02 7.43e+03 6.64e+00 ... (remaining 12705 not shown) Histogram of bond angle deviations from ideal: 98.56 - 105.65: 281 105.65 - 112.74: 6561 112.74 - 119.84: 4624 119.84 - 126.93: 5603 126.93 - 134.03: 155 Bond angle restraints: 17224 Sorted by residual: angle pdb=" N ILE B 904 " pdb=" CA ILE B 904 " pdb=" C ILE B 904 " ideal model delta sigma weight residual 111.91 107.11 4.80 8.90e-01 1.26e+00 2.91e+01 angle pdb=" N ILE B1141 " pdb=" CA ILE B1141 " pdb=" C ILE B1141 " ideal model delta sigma weight residual 113.53 109.79 3.74 9.80e-01 1.04e+00 1.45e+01 angle pdb=" C ASN C 169 " pdb=" N LEU C 170 " pdb=" CA LEU C 170 " ideal model delta sigma weight residual 121.54 127.86 -6.32 1.91e+00 2.74e-01 1.09e+01 angle pdb=" N GLY B 902 " pdb=" CA GLY B 902 " pdb=" C GLY B 902 " ideal model delta sigma weight residual 113.18 105.36 7.82 2.37e+00 1.78e-01 1.09e+01 angle pdb=" N ILE B1111 " pdb=" CA ILE B1111 " pdb=" C ILE B1111 " ideal model delta sigma weight residual 111.91 109.04 2.87 8.90e-01 1.26e+00 1.04e+01 ... (remaining 17219 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 32.08: 7410 32.08 - 64.17: 232 64.17 - 96.25: 23 96.25 - 128.33: 0 128.33 - 160.42: 1 Dihedral angle restraints: 7666 sinusoidal: 3399 harmonic: 4267 Sorted by residual: dihedral pdb=" O4' U G 16 " pdb=" C1' U G 16 " pdb=" N1 U G 16 " pdb=" C2 U G 16 " ideal model delta sinusoidal sigma weight residual 232.00 71.58 160.42 1 1.70e+01 3.46e-03 6.45e+01 dihedral pdb=" CA ASN B 933 " pdb=" C ASN B 933 " pdb=" N SER B 934 " pdb=" CA SER B 934 " ideal model delta harmonic sigma weight residual 180.00 150.55 29.45 0 5.00e+00 4.00e-02 3.47e+01 dihedral pdb=" CA ARG B 86 " pdb=" C ARG B 86 " pdb=" N THR B 87 " pdb=" CA THR B 87 " ideal model delta harmonic sigma weight residual -180.00 -151.68 -28.32 0 5.00e+00 4.00e-02 3.21e+01 ... (remaining 7663 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.048: 1271 0.048 - 0.096: 509 0.096 - 0.144: 77 0.144 - 0.192: 7 0.192 - 0.240: 1 Chirality restraints: 1865 Sorted by residual: chirality pdb=" C1' C G 11 " pdb=" O4' C G 11 " pdb=" C2' C G 11 " pdb=" N1 C G 11 " both_signs ideal model delta sigma weight residual False 2.47 2.23 0.24 2.00e-01 2.50e+01 1.44e+00 chirality pdb=" CA ILE B1111 " pdb=" N ILE B1111 " pdb=" C ILE B1111 " pdb=" CB ILE B1111 " both_signs ideal model delta sigma weight residual False 2.43 2.61 -0.17 2.00e-01 2.50e+01 7.38e-01 chirality pdb=" CA ILE B1034 " pdb=" N ILE B1034 " pdb=" C ILE B1034 " pdb=" CB ILE B1034 " both_signs ideal model delta sigma weight residual False 2.43 2.60 -0.17 2.00e-01 2.50e+01 7.09e-01 ... (remaining 1862 not shown) Planarity restraints: 2110 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C1' C G 11 " 0.027 2.00e-02 2.50e+03 1.60e-02 5.75e+00 pdb=" N1 C G 11 " -0.037 2.00e-02 2.50e+03 pdb=" C2 C G 11 " 0.003 2.00e-02 2.50e+03 pdb=" O2 C G 11 " 0.001 2.00e-02 2.50e+03 pdb=" N3 C G 11 " -0.004 2.00e-02 2.50e+03 pdb=" C4 C G 11 " 0.001 2.00e-02 2.50e+03 pdb=" N4 C G 11 " 0.011 2.00e-02 2.50e+03 pdb=" C5 C G 11 " 0.002 2.00e-02 2.50e+03 pdb=" C6 C G 11 " -0.003 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C1' A G 21 " -0.029 2.00e-02 2.50e+03 1.23e-02 4.17e+00 pdb=" N9 A G 21 " 0.021 2.00e-02 2.50e+03 pdb=" C8 A G 21 " 0.010 2.00e-02 2.50e+03 pdb=" N7 A G 21 " 0.004 2.00e-02 2.50e+03 pdb=" C5 A G 21 " -0.005 2.00e-02 2.50e+03 pdb=" C6 A G 21 " -0.007 2.00e-02 2.50e+03 pdb=" N6 A G 21 " -0.012 2.00e-02 2.50e+03 pdb=" N1 A G 21 " 0.006 2.00e-02 2.50e+03 pdb=" C2 A G 21 " 0.003 2.00e-02 2.50e+03 pdb=" N3 A G 21 " 0.005 2.00e-02 2.50e+03 pdb=" C4 A G 21 " 0.004 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C GLN A 193 " 0.033 5.00e-02 4.00e+02 4.91e-02 3.85e+00 pdb=" N PRO A 194 " -0.085 5.00e-02 4.00e+02 pdb=" CA PRO A 194 " 0.025 5.00e-02 4.00e+02 pdb=" CD PRO A 194 " 0.027 5.00e-02 4.00e+02 ... (remaining 2107 not shown) Histogram of nonbonded interaction distances: 1.97 - 2.56: 141 2.56 - 3.14: 10455 3.14 - 3.73: 18256 3.73 - 4.31: 26135 4.31 - 4.90: 42863 Nonbonded interactions: 97850 Sorted by model distance: nonbonded pdb=" NZ LYS A 202 " pdb="MG MG A 301 " model vdw 1.973 2.250 nonbonded pdb=" OP2 A G 19 " pdb="MG MG G1301 " model vdw 2.020 2.170 nonbonded pdb=" O THR B1055 " pdb=" OG1 THR B1055 " model vdw 2.104 2.440 nonbonded pdb=" O THR B 246 " pdb=" OG1 THR B 246 " model vdw 2.197 2.440 nonbonded pdb=" O LYS B 595 " pdb=" OG SER B 599 " model vdw 2.200 2.440 ... (remaining 97845 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = (chain 'A' and resid 120 through 233) selection = chain 'C' } 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 P 25 5.49 5 Mg 2 5.21 5 S 41 5.16 5 C 7870 2.51 5 N 2049 2.21 5 O 2417 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: 0.670 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: 3.930 Check model and map are aligned: 0.180 Process input model: 36.590 Find NCS groups from input model: 0.350 Set up NCS constraints: 0.050 Set refine NCS operators: 0.000 Set scattering table: 0.110 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.750 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 44.640 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7771 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.009 0.078 12710 Z= 0.600 Angle : 0.780 7.816 17224 Z= 0.455 Chirality : 0.048 0.240 1865 Planarity : 0.004 0.049 2110 Dihedral : 15.303 160.416 4914 Min Nonbonded Distance : 1.973 Molprobity Statistics. All-atom Clashscore : 10.70 Ramachandran Plot: Outliers : 0.07 % Allowed : 10.98 % Favored : 88.95 % Rotamer Outliers : 8.09 % 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: -4.94 (0.17), residues: 1421 helix: -2.79 (0.16), residues: 586 sheet: -3.82 (0.33), residues: 163 loop : -3.58 (0.19), residues: 672 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 318 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 106 poor density : 212 time to evaluate : 1.322 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 106 outliers final: 61 residues processed: 298 average time/residue: 0.2615 time to fit residues: 110.2416 Evaluate side-chains 219 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 61 poor density : 158 time to evaluate : 1.476 Switching outliers to nearest non-outliers outliers start: 61 outliers final: 0 residues processed: 61 average time/residue: 0.1379 time to fit residues: 16.9188 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 144 random chunks: chunk 121 optimal weight: 1.9990 chunk 109 optimal weight: 0.5980 chunk 60 optimal weight: 0.9990 chunk 37 optimal weight: 9.9990 chunk 73 optimal weight: 6.9990 chunk 58 optimal weight: 0.7980 chunk 112 optimal weight: 2.9990 chunk 43 optimal weight: 0.8980 chunk 68 optimal weight: 0.2980 chunk 84 optimal weight: 0.7980 chunk 130 optimal weight: 0.8980 overall best weight: 0.6780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 157 ASN A 168 ASN A 208 ASN B 33 ASN B 63 ASN B 157 ASN B 263 ASN B 567 GLN B 651 ASN B 676 GLN ** B 793 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 797 HIS B 808 ASN B 861 ASN B 864 ASN B 906 GLN ** C 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** C 181 GLN Total number of N/Q/H flips: 16 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7680 moved from start: 0.1743 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.030 12710 Z= 0.170 Angle : 0.552 6.533 17224 Z= 0.295 Chirality : 0.042 0.157 1865 Planarity : 0.004 0.045 2110 Dihedral : 8.876 167.231 1868 Min Nonbonded Distance : 2.170 Molprobity Statistics. All-atom Clashscore : 9.56 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.67 % Favored : 92.33 % Rotamer Outliers : 2.21 % 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: -3.34 (0.21), residues: 1421 helix: -1.22 (0.20), residues: 593 sheet: -3.57 (0.32), residues: 174 loop : -2.74 (0.22), residues: 654 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 224 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 29 poor density : 195 time to evaluate : 1.505 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 29 outliers final: 14 residues processed: 216 average time/residue: 0.2328 time to fit residues: 74.6975 Evaluate side-chains 175 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 161 time to evaluate : 1.418 Switching outliers to nearest non-outliers outliers start: 14 outliers final: 0 residues processed: 14 average time/residue: 0.1495 time to fit residues: 5.4316 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 144 random chunks: chunk 72 optimal weight: 0.5980 chunk 40 optimal weight: 0.2980 chunk 108 optimal weight: 3.9990 chunk 89 optimal weight: 0.0970 chunk 36 optimal weight: 3.9990 chunk 131 optimal weight: 2.9990 chunk 141 optimal weight: 6.9990 chunk 116 optimal weight: 0.8980 chunk 129 optimal weight: 3.9990 chunk 44 optimal weight: 0.7980 chunk 105 optimal weight: 0.5980 overall best weight: 0.4778 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 793 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7651 moved from start: 0.2332 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.026 12710 Z= 0.139 Angle : 0.494 5.638 17224 Z= 0.266 Chirality : 0.040 0.154 1865 Planarity : 0.003 0.044 2110 Dihedral : 8.693 173.670 1868 Min Nonbonded Distance : 2.202 Molprobity Statistics. All-atom Clashscore : 7.93 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.60 % Favored : 92.40 % Rotamer Outliers : 2.06 % 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: -2.46 (0.22), residues: 1421 helix: -0.31 (0.22), residues: 592 sheet: -3.24 (0.35), residues: 160 loop : -2.46 (0.23), residues: 669 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 222 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 27 poor density : 195 time to evaluate : 1.863 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 27 outliers final: 10 residues processed: 215 average time/residue: 0.2331 time to fit residues: 74.3767 Evaluate side-chains 177 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 167 time to evaluate : 1.428 Switching outliers to nearest non-outliers outliers start: 10 outliers final: 0 residues processed: 10 average time/residue: 0.1191 time to fit residues: 4.3018 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 144 random chunks: chunk 129 optimal weight: 0.9980 chunk 98 optimal weight: 0.6980 chunk 67 optimal weight: 5.9990 chunk 14 optimal weight: 0.8980 chunk 62 optimal weight: 0.7980 chunk 87 optimal weight: 3.9990 chunk 131 optimal weight: 2.9990 chunk 139 optimal weight: 4.9990 chunk 68 optimal weight: 8.9990 chunk 124 optimal weight: 3.9990 chunk 37 optimal weight: 10.0000 overall best weight: 1.2782 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 28 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7692 moved from start: 0.2348 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.033 12710 Z= 0.241 Angle : 0.539 7.979 17224 Z= 0.283 Chirality : 0.042 0.166 1865 Planarity : 0.003 0.045 2110 Dihedral : 8.677 173.235 1868 Min Nonbonded Distance : 2.201 Molprobity Statistics. All-atom Clashscore : 8.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.09 % Favored : 91.91 % Rotamer Outliers : 2.14 % 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: -2.18 (0.22), residues: 1421 helix: -0.09 (0.22), residues: 595 sheet: -3.09 (0.36), residues: 156 loop : -2.31 (0.23), residues: 670 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 199 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 28 poor density : 171 time to evaluate : 1.570 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 28 outliers final: 18 residues processed: 190 average time/residue: 0.2439 time to fit residues: 69.1691 Evaluate side-chains 179 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 161 time to evaluate : 1.421 Switching outliers to nearest non-outliers outliers start: 18 outliers final: 0 residues processed: 18 average time/residue: 0.1270 time to fit residues: 6.3813 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 144 random chunks: chunk 115 optimal weight: 0.0870 chunk 79 optimal weight: 2.9990 chunk 2 optimal weight: 0.8980 chunk 103 optimal weight: 0.3980 chunk 57 optimal weight: 4.9990 chunk 118 optimal weight: 0.5980 chunk 96 optimal weight: 0.7980 chunk 0 optimal weight: 10.0000 chunk 71 optimal weight: 7.9990 chunk 124 optimal weight: 6.9990 chunk 35 optimal weight: 9.9990 overall best weight: 0.5558 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 157 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7645 moved from start: 0.2744 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.026 12710 Z= 0.143 Angle : 0.504 9.495 17224 Z= 0.264 Chirality : 0.040 0.167 1865 Planarity : 0.003 0.041 2110 Dihedral : 8.603 176.173 1868 Min Nonbonded Distance : 2.209 Molprobity Statistics. All-atom Clashscore : 8.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.39 % Favored : 92.61 % Rotamer Outliers : 0.99 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.83 (0.23), residues: 1421 helix: 0.23 (0.22), residues: 599 sheet: -2.98 (0.36), residues: 152 loop : -2.15 (0.23), residues: 670 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 188 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 175 time to evaluate : 1.520 Fit side-chains revert: symmetry clash outliers start: 13 outliers final: 6 residues processed: 184 average time/residue: 0.2425 time to fit residues: 66.5542 Evaluate side-chains 162 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 156 time to evaluate : 1.481 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 0 residues processed: 6 average time/residue: 0.1223 time to fit residues: 3.4100 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 144 random chunks: chunk 46 optimal weight: 1.9990 chunk 125 optimal weight: 0.8980 chunk 27 optimal weight: 0.0770 chunk 81 optimal weight: 3.9990 chunk 34 optimal weight: 0.7980 chunk 139 optimal weight: 5.9990 chunk 115 optimal weight: 0.9990 chunk 64 optimal weight: 4.9990 chunk 11 optimal weight: 4.9990 chunk 73 optimal weight: 10.0000 chunk 134 optimal weight: 2.9990 overall best weight: 0.9542 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7669 moved from start: 0.2771 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.029 12710 Z= 0.192 Angle : 0.522 8.858 17224 Z= 0.270 Chirality : 0.041 0.172 1865 Planarity : 0.003 0.040 2110 Dihedral : 8.590 175.746 1868 Min Nonbonded Distance : 2.218 Molprobity Statistics. All-atom Clashscore : 8.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.23 % Favored : 91.77 % Rotamer Outliers : 0.99 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.67 (0.23), residues: 1421 helix: 0.37 (0.23), residues: 602 sheet: -2.89 (0.37), residues: 146 loop : -2.09 (0.24), residues: 673 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 180 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 167 time to evaluate : 1.501 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 13 outliers final: 10 residues processed: 175 average time/residue: 0.2520 time to fit residues: 65.9652 Evaluate side-chains 168 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 158 time to evaluate : 1.574 Switching outliers to nearest non-outliers outliers start: 10 outliers final: 0 residues processed: 10 average time/residue: 0.1316 time to fit residues: 4.7114 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 144 random chunks: chunk 15 optimal weight: 0.9990 chunk 79 optimal weight: 1.9990 chunk 101 optimal weight: 0.9980 chunk 78 optimal weight: 0.9980 chunk 117 optimal weight: 5.9990 chunk 77 optimal weight: 5.9990 chunk 138 optimal weight: 9.9990 chunk 86 optimal weight: 2.9990 chunk 84 optimal weight: 1.9990 chunk 64 optimal weight: 6.9990 chunk 85 optimal weight: 1.9990 overall best weight: 1.3986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 157 ASN ** C 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7693 moved from start: 0.2699 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.035 12710 Z= 0.255 Angle : 0.539 6.789 17224 Z= 0.283 Chirality : 0.042 0.174 1865 Planarity : 0.003 0.044 2110 Dihedral : 8.621 173.698 1868 Min Nonbonded Distance : 2.173 Molprobity Statistics. All-atom Clashscore : 8.50 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.66 % Favored : 91.34 % Rotamer Outliers : 0.61 % 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: -1.66 (0.23), residues: 1421 helix: 0.42 (0.23), residues: 602 sheet: -2.86 (0.37), residues: 160 loop : -2.11 (0.24), residues: 659 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 167 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 159 time to evaluate : 1.480 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 8 outliers final: 5 residues processed: 163 average time/residue: 0.2491 time to fit residues: 60.7187 Evaluate side-chains 155 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 150 time to evaluate : 1.355 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.1420 time to fit residues: 3.3952 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 144 random chunks: chunk 55 optimal weight: 0.5980 chunk 82 optimal weight: 0.7980 chunk 41 optimal weight: 2.9990 chunk 27 optimal weight: 4.9990 chunk 26 optimal weight: 1.9990 chunk 88 optimal weight: 2.9990 chunk 94 optimal weight: 1.9990 chunk 68 optimal weight: 10.0000 chunk 12 optimal weight: 1.9990 chunk 109 optimal weight: 0.9990 chunk 126 optimal weight: 1.9990 overall best weight: 1.2786 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 732 ASN ** C 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7689 moved from start: 0.2763 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.036 12710 Z= 0.238 Angle : 0.535 6.732 17224 Z= 0.281 Chirality : 0.042 0.175 1865 Planarity : 0.003 0.038 2110 Dihedral : 8.614 173.315 1868 Min Nonbonded Distance : 2.179 Molprobity Statistics. All-atom Clashscore : 8.70 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.30 % Favored : 91.70 % Rotamer Outliers : 0.69 % 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: -1.58 (0.23), residues: 1421 helix: 0.51 (0.23), residues: 602 sheet: -2.85 (0.36), residues: 160 loop : -2.08 (0.24), residues: 659 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 172 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 163 time to evaluate : 1.504 Fit side-chains revert: symmetry clash outliers start: 9 outliers final: 7 residues processed: 168 average time/residue: 0.2499 time to fit residues: 62.8038 Evaluate side-chains 158 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 151 time to evaluate : 1.529 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.1303 time to fit residues: 3.7951 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 144 random chunks: chunk 132 optimal weight: 3.9990 chunk 121 optimal weight: 3.9990 chunk 129 optimal weight: 0.0570 chunk 77 optimal weight: 4.9990 chunk 56 optimal weight: 1.9990 chunk 101 optimal weight: 0.8980 chunk 39 optimal weight: 2.9990 chunk 116 optimal weight: 3.9990 chunk 122 optimal weight: 0.7980 chunk 128 optimal weight: 0.9990 chunk 84 optimal weight: 0.0000 overall best weight: 0.5504 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 732 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7648 moved from start: 0.3030 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.027 12710 Z= 0.144 Angle : 0.512 6.650 17224 Z= 0.268 Chirality : 0.040 0.170 1865 Planarity : 0.003 0.038 2110 Dihedral : 8.596 176.214 1868 Min Nonbonded Distance : 2.228 Molprobity Statistics. All-atom Clashscore : 8.13 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.53 % Favored : 92.47 % Rotamer Outliers : 0.69 % 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: -1.36 (0.23), residues: 1421 helix: 0.73 (0.23), residues: 599 sheet: -2.81 (0.37), residues: 146 loop : -2.00 (0.24), residues: 676 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 176 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 167 time to evaluate : 1.514 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 9 outliers final: 5 residues processed: 172 average time/residue: 0.2421 time to fit residues: 61.6520 Evaluate side-chains 159 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 154 time to evaluate : 1.315 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.1319 time to fit residues: 3.0398 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 144 random chunks: chunk 136 optimal weight: 0.6980 chunk 83 optimal weight: 3.9990 chunk 64 optimal weight: 6.9990 chunk 95 optimal weight: 2.9990 chunk 143 optimal weight: 3.9990 chunk 132 optimal weight: 3.9990 chunk 114 optimal weight: 3.9990 chunk 11 optimal weight: 4.9990 chunk 88 optimal weight: 2.9990 chunk 70 optimal weight: 0.1980 chunk 90 optimal weight: 0.6980 overall best weight: 1.5184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 975 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7694 moved from start: 0.2845 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.044 12710 Z= 0.276 Angle : 0.562 7.780 17224 Z= 0.293 Chirality : 0.042 0.167 1865 Planarity : 0.003 0.038 2110 Dihedral : 8.625 173.481 1868 Min Nonbonded Distance : 2.166 Molprobity Statistics. All-atom Clashscore : 9.07 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.59 % Favored : 91.41 % Rotamer Outliers : 0.31 % 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: -1.48 (0.23), residues: 1421 helix: 0.66 (0.23), residues: 602 sheet: -2.79 (0.36), residues: 160 loop : -2.11 (0.24), residues: 659 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2842 Ramachandran restraints generated. 1421 Oldfield, 0 Emsley, 1421 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 153 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 149 time to evaluate : 1.422 Fit side-chains revert: symmetry clash outliers start: 4 outliers final: 4 residues processed: 153 average time/residue: 0.2686 time to fit residues: 60.8357 Evaluate side-chains 151 residues out of total 1310 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 147 time to evaluate : 1.308 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.1284 time to fit residues: 2.8203 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 144 random chunks: chunk 121 optimal weight: 5.9990 chunk 34 optimal weight: 3.9990 chunk 105 optimal weight: 0.3980 chunk 16 optimal weight: 0.9990 chunk 31 optimal weight: 6.9990 chunk 114 optimal weight: 2.9990 chunk 47 optimal weight: 4.9990 chunk 117 optimal weight: 6.9990 chunk 14 optimal weight: 0.0470 chunk 21 optimal weight: 4.9990 chunk 100 optimal weight: 0.5980 overall best weight: 1.0082 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 145 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4004 r_free = 0.4004 target = 0.171203 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 62)----------------| | r_work = 0.3296 r_free = 0.3296 target = 0.110100 restraints weight = 15581.886| |-----------------------------------------------------------------------------| r_work (start): 0.3276 rms_B_bonded: 2.62 r_work: 0.3113 rms_B_bonded: 3.16 restraints_weight: 0.5000 r_work: 0.2956 rms_B_bonded: 4.96 restraints_weight: 0.2500 r_work (final): 0.2956 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7935 moved from start: 0.2969 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.036 12710 Z= 0.199 Angle : 0.534 7.768 17224 Z= 0.279 Chirality : 0.041 0.166 1865 Planarity : 0.003 0.038 2110 Dihedral : 8.594 174.326 1868 Min Nonbonded Distance : 2.188 Molprobity Statistics. All-atom Clashscore : 8.70 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.95 % Favored : 92.05 % Rotamer Outliers : 0.46 % 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: -1.40 (0.23), residues: 1421 helix: 0.74 (0.23), residues: 602 sheet: -2.81 (0.36), residues: 160 loop : -2.07 (0.24), residues: 659 =============================================================================== Job complete usr+sys time: 2802.15 seconds wall clock time: 52 minutes 1.93 seconds (3121.93 seconds total)