Starting phenix.real_space_refine (version: dev) on Sat May 14 09:13:59 2022 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7q2x_13783/05_2022/7q2x_13783_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7q2x_13783/05_2022/7q2x_13783.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.0 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7q2x_13783/05_2022/7q2x_13783.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7q2x_13783/05_2022/7q2x_13783.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7q2x_13783/05_2022/7q2x_13783_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7q2x_13783/05_2022/7q2x_13783_updated.pdb" } resolution = 3.0 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.007 sd= 0.059 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 Convert atoms to be neutral Process input model Symmetric amino acids flipped Residue "A GLU 1113": "OE1" <-> "OE2" Residue "C GLU 96": "OE1" <-> "OE2" Time to flip residues: 0.05s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4597/modules/chem_data/mon_lib" Total number of atoms: 18240 Number of models: 1 Model: "" Number of chains: 8 Chain: "A" Number of atoms: 3413 Number of conformers: 1 Conformer: "" Number of residues, atoms: 431, 3413 Classifications: {'peptide': 431} Link IDs: {'PTRANS': 12, 'TRANS': 418} Chain breaks: 1 Chain: "B" Number of atoms: 3102 Number of conformers: 1 Conformer: "" Number of residues, atoms: 388, 3102 Classifications: {'peptide': 388} Incomplete info: {'truncation_to_alanine': 2} Link IDs: {'PTRANS': 12, 'TRANS': 375} Chain breaks: 2 Unresolved non-hydrogen bonds: 8 Unresolved non-hydrogen angles: 9 Unresolved non-hydrogen dihedrals: 7 Planarities with less than four sites: {'GLU:plan': 1} Unresolved non-hydrogen planarities: 4 Chain: "C" Number of atoms: 2279 Number of conformers: 1 Conformer: "" Number of residues, atoms: 285, 2279 Classifications: {'peptide': 285} Incomplete info: {'truncation_to_alanine': 2} Link IDs: {'PTRANS': 1, 'TRANS': 283} Chain breaks: 5 Unresolved non-hydrogen bonds: 8 Unresolved non-hydrogen angles: 8 Unresolved non-hydrogen dihedrals: 8 Chain: "D" Number of atoms: 7906 Number of conformers: 1 Conformer: "" Number of residues, atoms: 994, 7906 Classifications: {'peptide': 994} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 27, 'TRANS': 966} Chain breaks: 10 Unresolved non-hydrogen bonds: 4 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 4 Chain: "F" Number of atoms: 756 Number of conformers: 1 Conformer: "" Number of residues, atoms: 36, 756 Classifications: {'DNA': 36} Link IDs: {'rna3p': 35} Chain: "G" Number of atoms: 720 Number of conformers: 1 Conformer: "" Number of residues, atoms: 36, 720 Classifications: {'DNA': 36} Link IDs: {'rna3p': 35} Chain: "A" Number of atoms: 32 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 32 Unusual residues: {' MG': 1, 'BEF': 1, 'ADP': 1} Classifications: {'undetermined': 3} Link IDs: {None: 2} Chain: "B" Number of atoms: 32 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 32 Unusual residues: {' MG': 1, 'BEF': 1, 'ADP': 1} Classifications: {'undetermined': 3} Link IDs: {None: 2} Time building chain proxies: 10.72, per 1000 atoms: 0.59 Number of scatterers: 18240 At special positions: 0 Unit cell: (121.98, 147.66, 161.57, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 8 Type Number sf(0) S 91 16.00 P 76 15.00 Mg 2 11.99 F 6 9.00 O 3633 8.00 N 3057 7.00 C 11373 6.00 Be 2 4.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=1, symmetry=0 Simple disulfide: pdb=" SG CYS D1059 " - pdb=" SG CYS D1071 " distance=2.88 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: 8.22 Conformation dependent library (CDL) restraints added in 2.9 seconds 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 4020 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 93 helices and 9 sheets defined 60.2% alpha, 6.7% beta 36 base pairs and 64 stacking pairs defined. Time for finding SS restraints: 6.56 Creating SS restraints... Processing helix chain 'A' and resid 38 through 48 removed outlier: 3.929A pdb=" N LEU A 42 " --> pdb=" O LYS A 38 " (cutoff:3.500A) Processing helix chain 'A' and resid 54 through 57 Processing helix chain 'A' and resid 62 through 65 Processing helix chain 'A' and resid 124 through 134 removed outlier: 3.870A pdb=" N VAL A 134 " --> pdb=" O LEU A 130 " (cutoff:3.500A) Processing helix chain 'A' and resid 153 through 155 No H-bonds generated for 'chain 'A' and resid 153 through 155' Processing helix chain 'A' and resid 158 through 168 removed outlier: 3.623A pdb=" N ALA A 168 " --> pdb=" O LEU A 164 " (cutoff:3.500A) Processing helix chain 'A' and resid 175 through 227 Proline residue: A 205 - end of helix removed outlier: 3.985A pdb=" N ARG A 215 " --> pdb=" O ARG A 211 " (cutoff:3.500A) removed outlier: 4.008A pdb=" N MET A 216 " --> pdb=" O ASN A 212 " (cutoff:3.500A) Processing helix chain 'A' and resid 969 through 977 removed outlier: 4.074A pdb=" N LYS A 974 " --> pdb=" O GLN A 970 " (cutoff:3.500A) Processing helix chain 'A' and resid 986 through 1047 removed outlier: 3.835A pdb=" N GLU A1027 " --> pdb=" O GLU A1023 " (cutoff:3.500A) Processing helix chain 'A' and resid 1063 through 1065 No H-bonds generated for 'chain 'A' and resid 1063 through 1065' Processing helix chain 'A' and resid 1086 through 1103 Processing helix chain 'A' and resid 1120 through 1131 Processing helix chain 'A' and resid 1146 through 1149 No H-bonds generated for 'chain 'A' and resid 1146 through 1149' Processing helix chain 'B' and resid 135 through 141 Processing helix chain 'B' and resid 191 through 201 removed outlier: 3.925A pdb=" N ILE B 195 " --> pdb=" O LYS B 191 " (cutoff:3.500A) Processing helix chain 'B' and resid 206 through 209 Processing helix chain 'B' and resid 279 through 288 Processing helix chain 'B' and resid 323 through 331 Processing helix chain 'B' and resid 335 through 349 removed outlier: 4.302A pdb=" N LEU B 339 " --> pdb=" O ASN B 335 " (cutoff:3.500A) Processing helix chain 'B' and resid 1247 through 1286 Processing helix chain 'B' and resid 1301 through 1304 No H-bonds generated for 'chain 'B' and resid 1301 through 1304' Processing helix chain 'B' and resid 1325 through 1342 Processing helix chain 'B' and resid 1359 through 1371 Processing helix chain 'B' and resid 1385 through 1389 Processing helix chain 'C' and resid 23 through 40 Processing helix chain 'C' and resid 53 through 58 Processing helix chain 'C' and resid 72 through 105 Processing helix chain 'C' and resid 170 through 173 Processing helix chain 'C' and resid 186 through 197 Processing helix chain 'C' and resid 276 through 287 Processing helix chain 'C' and resid 300 through 310 removed outlier: 4.004A pdb=" N ARG C 306 " --> pdb=" O GLU C 303 " (cutoff:3.500A) Processing helix chain 'C' and resid 313 through 322 removed outlier: 4.425A pdb=" N ASP C 317 " --> pdb=" O ASN C 313 " (cutoff:3.500A) Processing helix chain 'C' and resid 647 through 668 Processing helix chain 'C' and resid 692 through 700 Processing helix chain 'C' and resid 705 through 710 removed outlier: 3.718A pdb=" N ASP C 710 " --> pdb=" O ASP C 706 " (cutoff:3.500A) Processing helix chain 'C' and resid 713 through 726 Processing helix chain 'D' and resid 6 through 14 Processing helix chain 'D' and resid 27 through 40 removed outlier: 4.151A pdb=" N VAL D 40 " --> pdb=" O ASP D 36 " (cutoff:3.500A) Processing helix chain 'D' and resid 42 through 44 No H-bonds generated for 'chain 'D' and resid 42 through 44' Processing helix chain 'D' and resid 51 through 60 Processing helix chain 'D' and resid 62 through 64 No H-bonds generated for 'chain 'D' and resid 62 through 64' Processing helix chain 'D' and resid 67 through 91 Processing helix chain 'D' and resid 105 through 127 Processing helix chain 'D' and resid 154 through 176 Processing helix chain 'D' and resid 186 through 205 removed outlier: 3.545A pdb=" N PHE D 196 " --> pdb=" O PHE D 192 " (cutoff:3.500A) Proline residue: D 199 - end of helix removed outlier: 3.805A pdb=" N GLU D 205 " --> pdb=" O PHE D 201 " (cutoff:3.500A) Processing helix chain 'D' and resid 207 through 211 Processing helix chain 'D' and resid 213 through 229 Processing helix chain 'D' and resid 235 through 246 Processing helix chain 'D' and resid 251 through 263 Processing helix chain 'D' and resid 269 through 280 Processing helix chain 'D' and resid 289 through 304 Processing helix chain 'D' and resid 306 through 310 Processing helix chain 'D' and resid 313 through 321 removed outlier: 4.006A pdb=" N THR D 318 " --> pdb=" O ASN D 314 " (cutoff:3.500A) removed outlier: 3.979A pdb=" N LEU D 319 " --> pdb=" O LEU D 315 " (cutoff:3.500A) removed outlier: 3.852A pdb=" N LEU D 320 " --> pdb=" O VAL D 316 " (cutoff:3.500A) removed outlier: 4.335A pdb=" N ASN D 321 " --> pdb=" O ILE D 317 " (cutoff:3.500A) Processing helix chain 'D' and resid 325 through 344 removed outlier: 3.577A pdb=" N ALA D 343 " --> pdb=" O VAL D 339 " (cutoff:3.500A) removed outlier: 4.448A pdb=" N GLN D 344 " --> pdb=" O ALA D 340 " (cutoff:3.500A) Processing helix chain 'D' and resid 346 through 351 Processing helix chain 'D' and resid 353 through 367 removed outlier: 3.748A pdb=" N PHE D 367 " --> pdb=" O LEU D 363 " (cutoff:3.500A) Processing helix chain 'D' and resid 372 through 386 Processing helix chain 'D' and resid 394 through 408 removed outlier: 3.970A pdb=" N LYS D 398 " --> pdb=" O LYS D 394 " (cutoff:3.500A) removed outlier: 4.315A pdb=" N GLN D 408 " --> pdb=" O VAL D 404 " (cutoff:3.500A) Processing helix chain 'D' and resid 412 through 427 Processing helix chain 'D' and resid 441 through 456 Processing helix chain 'D' and resid 518 through 551 Processing helix chain 'D' and resid 556 through 571 Processing helix chain 'D' and resid 577 through 587 removed outlier: 3.851A pdb=" N HIS D 586 " --> pdb=" O LYS D 582 " (cutoff:3.500A) removed outlier: 4.270A pdb=" N LEU D 587 " --> pdb=" O LYS D 583 " (cutoff:3.500A) Processing helix chain 'D' and resid 600 through 611 Processing helix chain 'D' and resid 621 through 637 Processing helix chain 'D' and resid 641 through 656 Processing helix chain 'D' and resid 663 through 675 removed outlier: 4.007A pdb=" N ALA D 675 " --> pdb=" O ILE D 671 " (cutoff:3.500A) Processing helix chain 'D' and resid 699 through 713 removed outlier: 3.851A pdb=" N ALA D 713 " --> pdb=" O MET D 709 " (cutoff:3.500A) Processing helix chain 'D' and resid 715 through 719 Processing helix chain 'D' and resid 722 through 728 Processing helix chain 'D' and resid 732 through 736 removed outlier: 3.757A pdb=" N LYS D 736 " --> pdb=" O ALA D 732 " (cutoff:3.500A) No H-bonds generated for 'chain 'D' and resid 732 through 736' Processing helix chain 'D' and resid 738 through 749 Processing helix chain 'D' and resid 765 through 780 removed outlier: 4.450A pdb=" N VAL D 770 " --> pdb=" O GLU D 766 " (cutoff:3.500A) Processing helix chain 'D' and resid 786 through 802 removed outlier: 3.665A pdb=" N TYR D 789 " --> pdb=" O ASN D 786 " (cutoff:3.500A) Proline residue: D 790 - end of helix removed outlier: 4.544A pdb=" N GLU D 793 " --> pdb=" O PRO D 790 " (cutoff:3.500A) Processing helix chain 'D' and resid 806 through 820 Processing helix chain 'D' and resid 839 through 874 removed outlier: 3.554A pdb=" N GLU D 874 " --> pdb=" O LYS D 870 " (cutoff:3.500A) Processing helix chain 'D' and resid 909 through 922 Processing helix chain 'D' and resid 933 through 945 Proline residue: D 938 - end of helix Processing helix chain 'D' and resid 947 through 950 No H-bonds generated for 'chain 'D' and resid 947 through 950' Processing helix chain 'D' and resid 953 through 967 Processing helix chain 'D' and resid 971 through 986 Proline residue: D 979 - end of helix Processing helix chain 'D' and resid 991 through 1006 Processing helix chain 'D' and resid 1009 through 1012 No H-bonds generated for 'chain 'D' and resid 1009 through 1012' Processing helix chain 'D' and resid 1015 through 1023 Processing helix chain 'D' and resid 1028 through 1043 Processing helix chain 'D' and resid 1053 through 1061 removed outlier: 3.826A pdb=" N CYS D1059 " --> pdb=" O GLU D1055 " (cutoff:3.500A) removed outlier: 4.690A pdb=" N LEU D1060 " --> pdb=" O MET D1056 " (cutoff:3.500A) removed outlier: 5.016A pdb=" N ASP D1061 " --> pdb=" O ALA D1057 " (cutoff:3.500A) Processing helix chain 'D' and resid 1065 through 1077 removed outlier: 4.285A pdb=" N MET D1070 " --> pdb=" O GLY D1066 " (cutoff:3.500A) removed outlier: 4.545A pdb=" N CYS D1071 " --> pdb=" O ILE D1067 " (cutoff:3.500A) Processing helix chain 'D' and resid 1081 through 1099 removed outlier: 4.055A pdb=" N ALA D1084 " --> pdb=" O LYS D1081 " (cutoff:3.500A) removed outlier: 5.653A pdb=" N ILE D1085 " --> pdb=" O ASP D1082 " (cutoff:3.500A) removed outlier: 5.684A pdb=" N ILE D1090 " --> pdb=" O ASN D1087 " (cutoff:3.500A) Processing helix chain 'D' and resid 1104 through 1117 removed outlier: 3.853A pdb=" N LYS D1109 " --> pdb=" O GLU D1105 " (cutoff:3.500A) removed outlier: 4.167A pdb=" N PHE D1117 " --> pdb=" O PHE D1113 " (cutoff:3.500A) Processing helix chain 'D' and resid 1123 through 1136 removed outlier: 3.814A pdb=" N ARG D1136 " --> pdb=" O VAL D1132 " (cutoff:3.500A) Processing helix chain 'D' and resid 1141 through 1146 Processing helix chain 'D' and resid 1149 through 1153 Processing helix chain 'D' and resid 1159 through 1167 Processing sheet with id= A, first strand: chain 'A' and resid 17 through 20 removed outlier: 6.975A pdb=" N VAL A 82 " --> pdb=" O GLU A 4 " (cutoff:3.500A) removed outlier: 4.525A pdb=" N LEU A 6 " --> pdb=" O THR A 80 " (cutoff:3.500A) removed outlier: 6.610A pdb=" N THR A 80 " --> pdb=" O LEU A 6 " (cutoff:3.500A) removed outlier: 4.601A pdb=" N ILE A 8 " --> pdb=" O SER A 78 " (cutoff:3.500A) removed outlier: 6.864A pdb=" N SER A 78 " --> pdb=" O ILE A 8 " (cutoff:3.500A) removed outlier: 4.120A pdb=" N TYR A 115 " --> pdb=" O ALA A 122 " (cutoff:3.500A) Processing sheet with id= B, first strand: chain 'A' and resid 1165 through 1169 Processing sheet with id= C, first strand: chain 'A' and resid 1051 through 1057 Processing sheet with id= D, first strand: chain 'B' and resid 153 through 156 Processing sheet with id= E, first strand: chain 'B' and resid 170 through 173 Processing sheet with id= F, first strand: chain 'B' and resid 1403 through 1410 removed outlier: 7.572A pdb=" N PHE B 180 " --> pdb=" O GLN B1393 " (cutoff:3.500A) removed outlier: 8.663A pdb=" N VAL B1395 " --> pdb=" O PHE B 180 " (cutoff:3.500A) removed outlier: 6.683A pdb=" N ALA B 182 " --> pdb=" O VAL B1395 " (cutoff:3.500A) removed outlier: 8.324A pdb=" N VAL B1397 " --> pdb=" O ALA B 182 " (cutoff:3.500A) removed outlier: 7.343A pdb=" N VAL B 184 " --> pdb=" O VAL B1397 " (cutoff:3.500A) removed outlier: 9.381A pdb=" N LYS B1399 " --> pdb=" O VAL B 184 " (cutoff:3.500A) removed outlier: 6.031A pdb=" N PHE B1378 " --> pdb=" O SER B 181 " (cutoff:3.500A) removed outlier: 7.371A pdb=" N VAL B 183 " --> pdb=" O PHE B1378 " (cutoff:3.500A) removed outlier: 6.529A pdb=" N VAL B1380 " --> pdb=" O VAL B 183 " (cutoff:3.500A) removed outlier: 8.502A pdb=" N GLY B 185 " --> pdb=" O VAL B1380 " (cutoff:3.500A) Processing sheet with id= G, first strand: chain 'B' and resid 239 through 241 Processing sheet with id= H, first strand: chain 'B' and resid 1290 through 1295 Processing sheet with id= I, first strand: chain 'C' and resid 730 through 733 934 hydrogen bonds defined for protein. 2649 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 72 hydrogen bonds 144 hydrogen bond angles 0 basepair planarities 36 basepair parallelities 64 stacking parallelities Total time for adding SS restraints: 7.32 Time building geometry restraints manager: 8.70 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.20 - 1.33: 3160 1.33 - 1.45: 4673 1.45 - 1.57: 10517 1.57 - 1.69: 148 1.69 - 1.81: 153 Bond restraints: 18651 Sorted by residual: bond pdb=" F2 BEF A2002 " pdb="BE BEF A2002 " ideal model delta sigma weight residual 1.476 1.560 -0.084 2.00e-02 2.50e+03 1.74e+01 bond pdb=" F2 BEF B2002 " pdb="BE BEF B2002 " ideal model delta sigma weight residual 1.476 1.559 -0.083 2.00e-02 2.50e+03 1.71e+01 bond pdb=" F3 BEF A2002 " pdb="BE BEF A2002 " ideal model delta sigma weight residual 1.476 1.546 -0.070 2.00e-02 2.50e+03 1.24e+01 bond pdb=" F3 BEF B2002 " pdb="BE BEF B2002 " ideal model delta sigma weight residual 1.476 1.546 -0.070 2.00e-02 2.50e+03 1.22e+01 bond pdb=" F1 BEF B2002 " pdb="BE BEF B2002 " ideal model delta sigma weight residual 1.476 1.531 -0.055 2.00e-02 2.50e+03 7.59e+00 ... (remaining 18646 not shown) Histogram of bond angle deviations from ideal: 99.69 - 107.00: 794 107.00 - 114.31: 11230 114.31 - 121.61: 9472 121.61 - 128.92: 3797 128.92 - 136.23: 136 Bond angle restraints: 25429 Sorted by residual: angle pdb=" F2 BEF A2002 " pdb="BE BEF A2002 " pdb=" F3 BEF A2002 " ideal model delta sigma weight residual 119.96 107.49 12.47 3.00e+00 1.11e-01 1.73e+01 angle pdb=" F2 BEF B2002 " pdb="BE BEF B2002 " pdb=" F3 BEF B2002 " ideal model delta sigma weight residual 119.96 107.58 12.38 3.00e+00 1.11e-01 1.70e+01 angle pdb=" N PRO A 158 " pdb=" CA PRO A 158 " pdb=" CB PRO A 158 " ideal model delta sigma weight residual 103.30 99.74 3.56 1.12e+00 7.97e-01 1.01e+01 angle pdb=" F1 BEF A2002 " pdb="BE BEF A2002 " pdb=" F3 BEF A2002 " ideal model delta sigma weight residual 120.12 110.99 9.13 3.00e+00 1.11e-01 9.26e+00 angle pdb=" F1 BEF B2002 " pdb="BE BEF B2002 " pdb=" F3 BEF B2002 " ideal model delta sigma weight residual 120.12 110.99 9.13 3.00e+00 1.11e-01 9.25e+00 ... (remaining 25424 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 35.96: 10734 35.96 - 71.93: 442 71.93 - 107.89: 19 107.89 - 143.85: 1 143.85 - 179.82: 3 Dihedral angle restraints: 11199 sinusoidal: 5052 harmonic: 6147 Sorted by residual: dihedral pdb=" C5' ADP B2001 " pdb=" O5' ADP B2001 " pdb=" PA ADP B2001 " pdb=" O2A ADP B2001 " ideal model delta sinusoidal sigma weight residual -60.00 101.82 -161.81 1 2.00e+01 2.50e-03 4.68e+01 dihedral pdb=" O2A ADP B2001 " pdb=" O3A ADP B2001 " pdb=" PA ADP B2001 " pdb=" PB ADP B2001 " ideal model delta sinusoidal sigma weight residual -60.00 85.28 -145.28 1 2.00e+01 2.50e-03 4.37e+01 dihedral pdb=" O1B ADP B2001 " pdb=" O3A ADP B2001 " pdb=" PB ADP B2001 " pdb=" PA ADP B2001 " ideal model delta sinusoidal sigma weight residual -60.00 69.77 -129.78 1 2.00e+01 2.50e-03 3.94e+01 ... (remaining 11196 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.032: 2057 0.032 - 0.063: 625 0.063 - 0.095: 195 0.095 - 0.127: 80 0.127 - 0.159: 6 Chirality restraints: 2963 Sorted by residual: chirality pdb=" CA THR D 186 " pdb=" N THR D 186 " pdb=" C THR D 186 " pdb=" CB THR D 186 " both_signs ideal model delta sigma weight residual False 2.53 2.37 0.16 2.00e-01 2.50e+01 6.30e-01 chirality pdb=" CA ILE C 731 " pdb=" N ILE C 731 " pdb=" C ILE C 731 " pdb=" CB ILE C 731 " both_signs ideal model delta sigma weight residual False 2.43 2.57 -0.14 2.00e-01 2.50e+01 4.98e-01 chirality pdb=" CA ILE B 258 " pdb=" N ILE B 258 " pdb=" C ILE B 258 " pdb=" CB ILE B 258 " both_signs ideal model delta sigma weight residual False 2.43 2.57 -0.14 2.00e-01 2.50e+01 4.74e-01 ... (remaining 2960 not shown) Planarity restraints: 2951 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C SER D 305 " 0.065 5.00e-02 4.00e+02 9.79e-02 1.53e+01 pdb=" N PRO D 306 " -0.169 5.00e-02 4.00e+02 pdb=" CA PRO D 306 " 0.051 5.00e-02 4.00e+02 pdb=" CD PRO D 306 " 0.054 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C LYS B 312 " -0.039 5.00e-02 4.00e+02 5.90e-02 5.57e+00 pdb=" N PRO B 313 " 0.102 5.00e-02 4.00e+02 pdb=" CA PRO B 313 " -0.031 5.00e-02 4.00e+02 pdb=" CD PRO B 313 " -0.033 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C THR B1345 " -0.033 5.00e-02 4.00e+02 5.00e-02 4.00e+00 pdb=" N PRO B1346 " 0.087 5.00e-02 4.00e+02 pdb=" CA PRO B1346 " -0.026 5.00e-02 4.00e+02 pdb=" CD PRO B1346 " -0.028 5.00e-02 4.00e+02 ... (remaining 2948 not shown) Histogram of nonbonded interaction distances: 2.04 - 2.61: 154 2.61 - 3.18: 15992 3.18 - 3.75: 29288 3.75 - 4.33: 38824 4.33 - 4.90: 63004 Nonbonded interactions: 147262 Sorted by model distance: nonbonded pdb=" F1 BEF B2002 " pdb="MG MG B2003 " model vdw 2.035 2.120 nonbonded pdb=" O2B ADP A2001 " pdb="MG MG A2003 " model vdw 2.063 2.170 nonbonded pdb=" OG SER B 192 " pdb="MG MG B2003 " model vdw 2.085 2.170 nonbonded pdb=" OE1 GLN B 302 " pdb="MG MG B2003 " model vdw 2.133 2.170 nonbonded pdb=" OG SER A 39 " pdb="MG MG A2003 " model vdw 2.143 2.170 ... (remaining 147257 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: 8 Type Number sf(0) Gaussians P 76 5.49 5 Mg 2 5.21 5 S 91 5.16 5 Be 2 3.05 5 C 11373 2.51 5 N 3057 2.21 5 O 3633 1.98 5 F 6 1.80 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.270 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.010 Extract box with map and model: 2.730 Check model and map are aligned: 0.260 Convert atoms to be neutral: 0.160 Process input model: 54.400 Find NCS groups from input model: 0.490 Set up NCS constraints: 0.070 Set refine NCS operators: 0.000 Set scattering table: 0.010 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:10.670 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 69.070 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6973 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.084 18651 Z= 0.144 Angle : 0.546 12.467 25429 Z= 0.285 Chirality : 0.037 0.159 2963 Planarity : 0.004 0.098 2951 Dihedral : 17.413 179.817 7176 Min Nonbonded Distance : 2.035 Molprobity Statistics. All-atom Clashscore : 7.71 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.58 % Favored : 97.42 % Rotamer Outliers : 0.11 % 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.57 (0.19), residues: 2054 helix: 2.71 (0.15), residues: 1260 sheet: 1.48 (0.39), residues: 180 loop : -0.44 (0.24), residues: 614 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 440 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 438 time to evaluate : 2.010 Fit side-chains outliers start: 2 outliers final: 0 residues processed: 438 average time/residue: 1.1884 time to fit residues: 588.2021 Evaluate side-chains 385 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 385 time to evaluate : 2.136 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 3.1489 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 214 random chunks: chunk 180 optimal weight: 2.9990 chunk 162 optimal weight: 8.9990 chunk 90 optimal weight: 7.9990 chunk 55 optimal weight: 9.9990 chunk 109 optimal weight: 3.9990 chunk 86 optimal weight: 7.9990 chunk 167 optimal weight: 3.9990 chunk 64 optimal weight: 7.9990 chunk 101 optimal weight: 0.0970 chunk 124 optimal weight: 2.9990 chunk 194 optimal weight: 4.9990 overall best weight: 2.8186 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 71 GLN A 106 GLN B 171 GLN B 266 ASN B 307 ASN C 307 ASN C 653 ASN C 725 ASN C 738 ASN D 70 GLN D 83 ASN D 163 GLN ** D 371 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 408 GLN ** D 852 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** D1062 ASN Total number of N/Q/H flips: 14 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7534 moved from start: 0.3341 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.006 0.087 18651 Z= 0.357 Angle : 0.687 10.270 25429 Z= 0.358 Chirality : 0.044 0.176 2963 Planarity : 0.005 0.076 2951 Dihedral : 18.443 175.721 2941 Min Nonbonded Distance : 1.772 Molprobity Statistics. All-atom Clashscore : 9.23 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.92 % Favored : 97.08 % Rotamer Outliers : 3.15 % 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.97 (0.18), residues: 2054 helix: 2.19 (0.14), residues: 1266 sheet: 1.17 (0.39), residues: 183 loop : -0.57 (0.24), residues: 605 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 583 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 60 poor density : 523 time to evaluate : 2.294 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 60 outliers final: 16 residues processed: 547 average time/residue: 1.2372 time to fit residues: 760.1517 Evaluate side-chains 483 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 16 poor density : 467 time to evaluate : 1.971 Switching outliers to nearest non-outliers outliers start: 16 outliers final: 10 residues processed: 7 average time/residue: 0.5981 time to fit residues: 7.9647 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 214 random chunks: chunk 108 optimal weight: 6.9990 chunk 60 optimal weight: 0.0980 chunk 161 optimal weight: 1.9990 chunk 132 optimal weight: 0.8980 chunk 53 optimal weight: 0.9990 chunk 194 optimal weight: 3.9990 chunk 210 optimal weight: 4.9990 chunk 173 optimal weight: 4.9990 chunk 193 optimal weight: 4.9990 chunk 66 optimal weight: 0.9990 chunk 156 optimal weight: 4.9990 overall best weight: 0.9986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 28 ASN A 132 GLN A 983 ASN A1138 GLN B1283 GLN C 314 GLN C 660 ASN D 163 GLN ** D 254 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 322 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 610 GLN D 634 ASN D 781 ASN D 852 GLN Total number of N/Q/H flips: 12 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7510 moved from start: 0.3734 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.044 18651 Z= 0.178 Angle : 0.549 9.221 25429 Z= 0.283 Chirality : 0.038 0.153 2963 Planarity : 0.004 0.063 2951 Dihedral : 18.194 176.545 2941 Min Nonbonded Distance : 1.863 Molprobity Statistics. All-atom Clashscore : 8.16 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.82 % Favored : 97.18 % Rotamer Outliers : 3.36 % 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.19 (0.18), residues: 2054 helix: 2.41 (0.14), residues: 1260 sheet: 1.24 (0.39), residues: 181 loop : -0.55 (0.23), residues: 613 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 529 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 64 poor density : 465 time to evaluate : 2.121 Fit side-chains revert: symmetry clash outliers start: 64 outliers final: 24 residues processed: 491 average time/residue: 1.3549 time to fit residues: 746.5143 Evaluate side-chains 480 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 24 poor density : 456 time to evaluate : 2.019 Switching outliers to nearest non-outliers outliers start: 24 outliers final: 13 residues processed: 12 average time/residue: 0.9278 time to fit residues: 15.6774 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 214 random chunks: chunk 192 optimal weight: 6.9990 chunk 146 optimal weight: 2.9990 chunk 101 optimal weight: 0.1980 chunk 21 optimal weight: 10.0000 chunk 92 optimal weight: 5.9990 chunk 130 optimal weight: 5.9990 chunk 195 optimal weight: 4.9990 chunk 206 optimal weight: 0.4980 chunk 102 optimal weight: 0.7980 chunk 185 optimal weight: 2.9990 chunk 55 optimal weight: 9.9990 overall best weight: 1.4984 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 132 GLN ** A1138 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A1153 ASN B 211 GLN B1283 GLN B1393 GLN C 660 ASN ** D 254 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 312 GLN D 322 ASN ** D 634 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 781 ASN Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7561 moved from start: 0.4208 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.043 18651 Z= 0.220 Angle : 0.570 13.410 25429 Z= 0.293 Chirality : 0.039 0.157 2963 Planarity : 0.004 0.056 2951 Dihedral : 18.207 173.452 2941 Min Nonbonded Distance : 1.848 Molprobity Statistics. All-atom Clashscore : 8.60 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.87 % Favored : 97.13 % Rotamer Outliers : 3.84 % 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.17 (0.18), residues: 2054 helix: 2.39 (0.14), residues: 1257 sheet: 1.09 (0.39), residues: 181 loop : -0.47 (0.23), residues: 616 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 555 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 73 poor density : 482 time to evaluate : 2.058 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 73 outliers final: 28 residues processed: 518 average time/residue: 1.2491 time to fit residues: 727.6669 Evaluate side-chains 496 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 28 poor density : 468 time to evaluate : 2.048 Switching outliers to nearest non-outliers outliers start: 28 outliers final: 19 residues processed: 10 average time/residue: 0.5047 time to fit residues: 9.4026 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 214 random chunks: chunk 172 optimal weight: 8.9990 chunk 117 optimal weight: 6.9990 chunk 3 optimal weight: 3.9990 chunk 154 optimal weight: 3.9990 chunk 85 optimal weight: 6.9990 chunk 176 optimal weight: 0.9980 chunk 142 optimal weight: 0.6980 chunk 0 optimal weight: 8.9990 chunk 105 optimal weight: 8.9990 chunk 185 optimal weight: 1.9990 chunk 52 optimal weight: 0.9980 overall best weight: 1.7384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 106 GLN A 132 GLN A 142 ASN A 972 ASN A1138 GLN B1393 GLN C 653 ASN ** D 254 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 634 ASN D 781 ASN Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7595 moved from start: 0.4634 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.046 18651 Z= 0.240 Angle : 0.589 10.477 25429 Z= 0.302 Chirality : 0.040 0.156 2963 Planarity : 0.004 0.054 2951 Dihedral : 18.224 173.272 2941 Min Nonbonded Distance : 1.846 Molprobity Statistics. All-atom Clashscore : 8.79 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.26 % Favored : 96.74 % Rotamer Outliers : 4.05 % 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.07 (0.19), residues: 2054 helix: 2.32 (0.14), residues: 1253 sheet: 0.92 (0.39), residues: 181 loop : -0.50 (0.24), residues: 620 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 563 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 77 poor density : 486 time to evaluate : 2.256 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 77 outliers final: 33 residues processed: 523 average time/residue: 1.2337 time to fit residues: 726.6127 Evaluate side-chains 499 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 33 poor density : 466 time to evaluate : 2.149 Switching outliers to nearest non-outliers outliers start: 33 outliers final: 18 residues processed: 17 average time/residue: 1.1862 time to fit residues: 26.0554 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 214 random chunks: chunk 69 optimal weight: 9.9990 chunk 186 optimal weight: 0.4980 chunk 40 optimal weight: 3.9990 chunk 121 optimal weight: 0.6980 chunk 51 optimal weight: 0.0030 chunk 207 optimal weight: 20.0000 chunk 171 optimal weight: 1.9990 chunk 95 optimal weight: 10.0000 chunk 17 optimal weight: 4.9990 chunk 68 optimal weight: 1.9990 chunk 108 optimal weight: 1.9990 overall best weight: 1.0394 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 106 GLN A 132 GLN ** A 142 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 221 GLN A1022 ASN A1138 GLN A1153 ASN B1393 GLN C 212 ASN C 292 GLN C 725 ASN ** D 254 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 10 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7574 moved from start: 0.4774 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.037 18651 Z= 0.183 Angle : 0.580 13.298 25429 Z= 0.293 Chirality : 0.039 0.163 2963 Planarity : 0.004 0.054 2951 Dihedral : 18.185 171.793 2941 Min Nonbonded Distance : 1.887 Molprobity Statistics. All-atom Clashscore : 8.74 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.21 % Favored : 96.79 % Rotamer Outliers : 3.84 % 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.09 (0.19), residues: 2054 helix: 2.36 (0.15), residues: 1255 sheet: 0.90 (0.39), residues: 181 loop : -0.53 (0.24), residues: 618 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 547 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 73 poor density : 474 time to evaluate : 2.359 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 73 outliers final: 36 residues processed: 509 average time/residue: 1.1909 time to fit residues: 685.4822 Evaluate side-chains 496 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 36 poor density : 460 time to evaluate : 2.045 Switching outliers to nearest non-outliers outliers start: 36 outliers final: 19 residues processed: 19 average time/residue: 0.5944 time to fit residues: 16.9155 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 214 random chunks: chunk 199 optimal weight: 7.9990 chunk 23 optimal weight: 1.9990 chunk 117 optimal weight: 0.4980 chunk 151 optimal weight: 10.0000 chunk 174 optimal weight: 3.9990 chunk 115 optimal weight: 5.9990 chunk 206 optimal weight: 20.0000 chunk 129 optimal weight: 2.9990 chunk 125 optimal weight: 4.9990 chunk 95 optimal weight: 0.9990 chunk 127 optimal weight: 7.9990 overall best weight: 2.0988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 71 GLN A 106 GLN A 132 GLN ** A 142 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A1022 ASN A1138 GLN A1153 ASN B 132 ASN B 265 ASN B1393 GLN ** C 725 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** C 738 ASN ** D 254 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 10 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7624 moved from start: 0.5107 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.047 18651 Z= 0.277 Angle : 0.643 16.995 25429 Z= 0.321 Chirality : 0.041 0.223 2963 Planarity : 0.004 0.055 2951 Dihedral : 18.238 171.027 2941 Min Nonbonded Distance : 1.845 Molprobity Statistics. All-atom Clashscore : 9.76 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.31 % Favored : 96.69 % Rotamer Outliers : 3.73 % 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.90 (0.18), residues: 2054 helix: 2.24 (0.14), residues: 1251 sheet: 0.79 (0.39), residues: 181 loop : -0.65 (0.23), residues: 622 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 555 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 71 poor density : 484 time to evaluate : 2.121 Fit side-chains revert: symmetry clash outliers start: 71 outliers final: 45 residues processed: 516 average time/residue: 1.2700 time to fit residues: 737.1327 Evaluate side-chains 518 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 45 poor density : 473 time to evaluate : 2.239 Switching outliers to nearest non-outliers outliers start: 45 outliers final: 22 residues processed: 26 average time/residue: 0.8343 time to fit residues: 29.0026 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 214 random chunks: chunk 82 optimal weight: 0.3980 chunk 123 optimal weight: 0.9990 chunk 62 optimal weight: 4.9990 chunk 40 optimal weight: 1.9990 chunk 39 optimal weight: 0.9990 chunk 131 optimal weight: 10.0000 chunk 140 optimal weight: 6.9990 chunk 101 optimal weight: 2.9990 chunk 19 optimal weight: 1.9990 chunk 162 optimal weight: 7.9990 chunk 187 optimal weight: 2.9990 overall best weight: 1.2788 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 63 GLN A 106 GLN A 132 GLN ** A 142 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A1138 GLN B 132 ASN B 211 GLN B 265 ASN B1385 ASN B1393 GLN C 725 ASN D 163 GLN ** D 254 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 715 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 11 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7603 moved from start: 0.5175 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.039 18651 Z= 0.207 Angle : 0.617 14.310 25429 Z= 0.307 Chirality : 0.040 0.212 2963 Planarity : 0.004 0.055 2951 Dihedral : 18.159 170.813 2941 Min Nonbonded Distance : 1.887 Molprobity Statistics. All-atom Clashscore : 9.59 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.26 % Favored : 96.74 % Rotamer Outliers : 3.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.94 (0.18), residues: 2054 helix: 2.28 (0.14), residues: 1251 sheet: 0.74 (0.38), residues: 181 loop : -0.64 (0.24), residues: 622 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 527 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 63 poor density : 464 time to evaluate : 2.045 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 63 outliers final: 35 residues processed: 497 average time/residue: 1.2855 time to fit residues: 717.3352 Evaluate side-chains 495 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 35 poor density : 460 time to evaluate : 2.130 Switching outliers to nearest non-outliers outliers start: 35 outliers final: 25 residues processed: 11 average time/residue: 0.6277 time to fit residues: 11.4057 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 214 random chunks: chunk 197 optimal weight: 1.9990 chunk 180 optimal weight: 5.9990 chunk 192 optimal weight: 5.9990 chunk 115 optimal weight: 5.9990 chunk 83 optimal weight: 0.6980 chunk 150 optimal weight: 7.9990 chunk 58 optimal weight: 4.9990 chunk 173 optimal weight: 1.9990 chunk 181 optimal weight: 0.9990 chunk 191 optimal weight: 1.9990 chunk 126 optimal weight: 0.6980 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 63 GLN A 106 GLN A 140 ASN ** A 142 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A1022 ASN A1138 GLN A1153 ASN B 132 ASN B 211 GLN B1385 ASN B1393 GLN D 162 ASN ** D 254 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 715 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 852 GLN Total number of N/Q/H flips: 12 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7601 moved from start: 0.5243 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.047 18651 Z= 0.214 Angle : 0.636 16.583 25429 Z= 0.315 Chirality : 0.040 0.296 2963 Planarity : 0.004 0.056 2951 Dihedral : 18.116 170.964 2941 Min Nonbonded Distance : 1.890 Molprobity Statistics. All-atom Clashscore : 9.92 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.21 % Favored : 96.79 % Rotamer Outliers : 3.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.99 (0.19), residues: 2054 helix: 2.33 (0.15), residues: 1245 sheet: 0.66 (0.38), residues: 181 loop : -0.59 (0.24), residues: 628 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 530 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 57 poor density : 473 time to evaluate : 1.985 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 57 outliers final: 37 residues processed: 509 average time/residue: 1.2934 time to fit residues: 741.7510 Evaluate side-chains 503 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 37 poor density : 466 time to evaluate : 2.165 Switching outliers to nearest non-outliers outliers start: 37 outliers final: 25 residues processed: 13 average time/residue: 0.6624 time to fit residues: 13.8975 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 214 random chunks: chunk 203 optimal weight: 2.9990 chunk 124 optimal weight: 0.6980 chunk 96 optimal weight: 2.9990 chunk 141 optimal weight: 1.9990 chunk 213 optimal weight: 10.0000 chunk 196 optimal weight: 0.4980 chunk 169 optimal weight: 7.9990 chunk 17 optimal weight: 0.7980 chunk 131 optimal weight: 0.0970 chunk 104 optimal weight: 7.9990 chunk 134 optimal weight: 1.9990 overall best weight: 0.8180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 63 GLN A 106 GLN A 139 ASN ** A 142 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A1022 ASN A1138 GLN B 211 GLN B1322 ASN B1385 ASN B1392 GLN B1393 GLN ** D 254 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 852 GLN Total number of N/Q/H flips: 11 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7580 moved from start: 0.5265 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.047 18651 Z= 0.188 Angle : 0.637 16.832 25429 Z= 0.315 Chirality : 0.040 0.283 2963 Planarity : 0.004 0.055 2951 Dihedral : 18.057 171.070 2941 Min Nonbonded Distance : 1.908 Molprobity Statistics. All-atom Clashscore : 10.53 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.26 % Favored : 96.74 % Rotamer Outliers : 2.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: 2.00 (0.19), residues: 2054 helix: 2.33 (0.15), residues: 1255 sheet: 0.64 (0.39), residues: 181 loop : -0.60 (0.24), residues: 618 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4108 Ramachandran restraints generated. 2054 Oldfield, 0 Emsley, 2054 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 509 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 44 poor density : 465 time to evaluate : 2.166 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 44 outliers final: 30 residues processed: 494 average time/residue: 1.2999 time to fit residues: 720.9245 Evaluate side-chains 481 residues out of total 1907 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 30 poor density : 451 time to evaluate : 2.240 Switching outliers to nearest non-outliers outliers start: 30 outliers final: 23 residues processed: 8 average time/residue: 0.7860 time to fit residues: 10.4800 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 214 random chunks: chunk 180 optimal weight: 3.9990 chunk 52 optimal weight: 2.9990 chunk 156 optimal weight: 7.9990 chunk 25 optimal weight: 5.9990 chunk 47 optimal weight: 4.9990 chunk 170 optimal weight: 0.8980 chunk 71 optimal weight: 2.9990 chunk 174 optimal weight: 3.9990 chunk 21 optimal weight: 5.9990 chunk 31 optimal weight: 0.5980 chunk 149 optimal weight: 9.9990 overall best weight: 2.2986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 63 GLN A 106 GLN A 139 ASN ** A 142 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A1138 GLN A1153 ASN B 211 GLN B1322 ASN B1392 GLN B1393 GLN ** D 254 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** D 312 GLN D 852 GLN Total number of N/Q/H flips: 11 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3883 r_free = 0.3883 target = 0.151417 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 34)----------------| | r_work = 0.3516 r_free = 0.3516 target = 0.122824 restraints weight = 24651.756| |-----------------------------------------------------------------------------| r_work (start): 0.3507 rms_B_bonded: 1.70 r_work: 0.3489 rms_B_bonded: 1.51 restraints_weight: 0.5000 r_work: 0.3473 rms_B_bonded: 1.51 restraints_weight: 0.2500 r_work: 0.3457 rms_B_bonded: 1.59 restraints_weight: 0.1250 r_work: 0.3440 rms_B_bonded: 1.72 restraints_weight: 0.0625 r_work: 0.3422 rms_B_bonded: 1.90 restraints_weight: 0.0312 r_work: 0.3401 rms_B_bonded: 2.13 restraints_weight: 0.0156 r_work: 0.3378 rms_B_bonded: 2.42 restraints_weight: 0.0078 r_work: 0.3352 rms_B_bonded: 2.75 restraints_weight: 0.0039 r_work: 0.3323 rms_B_bonded: 3.15 restraints_weight: 0.0020 r_work: 0.3290 rms_B_bonded: 3.62 restraints_weight: 0.0010 r_work: 0.3252 rms_B_bonded: 4.18 restraints_weight: 0.0005 r_work: 0.3209 rms_B_bonded: 4.83 restraints_weight: 0.0002 r_work (final): 0.3209 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7885 moved from start: 0.5465 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.058 18651 Z= 0.305 Angle : 0.692 17.133 25429 Z= 0.342 Chirality : 0.043 0.259 2963 Planarity : 0.004 0.057 2951 Dihedral : 18.102 170.394 2941 Min Nonbonded Distance : 1.863 Molprobity Statistics. All-atom Clashscore : 10.75 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.41 % Favored : 96.59 % Rotamer Outliers : 2.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.82 (0.19), residues: 2054 helix: 2.20 (0.14), residues: 1251 sheet: 0.63 (0.39), residues: 181 loop : -0.68 (0.24), residues: 622 =============================================================================== Job complete usr+sys time: 9129.53 seconds wall clock time: 161 minutes 22.74 seconds (9682.74 seconds total)