Starting phenix.real_space_refine on Fri Dec 8 22:29:00 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7whn_32505/12_2023/7whn_32505.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7whn_32505/12_2023/7whn_32505.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.3 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7whn_32505/12_2023/7whn_32505.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7whn_32505/12_2023/7whn_32505.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7whn_32505/12_2023/7whn_32505.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7whn_32505/12_2023/7whn_32505.pdb" } resolution = 3.3 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= -0.000 sd= 0.005 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 4 Type Number sf(0) Gaussians S 39 5.16 5 C 5816 2.51 5 N 1629 2.21 5 O 1800 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A ARG 45": "NH1" <-> "NH2" Residue "A ARG 47": "NH1" <-> "NH2" Residue "A ARG 57": "NH1" <-> "NH2" Residue "A ARG 62": "NH1" <-> "NH2" Residue "A GLU 71": "OE1" <-> "OE2" Residue "A GLU 91": "OE1" <-> "OE2" Residue "A ARG 129": "NH1" <-> "NH2" Residue "A GLU 145": "OE1" <-> "OE2" Residue "A ARG 166": "NH1" <-> "NH2" Residue "A GLU 184": "OE1" <-> "OE2" Residue "A ARG 198": "NH1" <-> "NH2" Residue "A GLU 210": "OE1" <-> "OE2" Residue "A ARG 219": "NH1" <-> "NH2" Residue "A ARG 247": "NH1" <-> "NH2" Residue "A ARG 277": "NH1" <-> "NH2" Residue "A ARG 280": "NH1" <-> "NH2" Residue "A ARG 312": "NH1" <-> "NH2" Residue "A GLU 339": "OE1" <-> "OE2" Residue "A ARG 353": "NH1" <-> "NH2" Residue "A ARG 382": "NH1" <-> "NH2" Residue "A ARG 386": "NH1" <-> "NH2" Residue "A ARG 437": "NH1" <-> "NH2" Residue "A ARG 467": "NH1" <-> "NH2" Residue "A ARG 471": "NH1" <-> "NH2" Residue "A ARG 518": "NH1" <-> "NH2" Residue "A ARG 521": "NH1" <-> "NH2" Residue "A ARG 524": "NH1" <-> "NH2" Residue "A ARG 551": "NH1" <-> "NH2" Residue "A ARG 563": "NH1" <-> "NH2" Residue "A ARG 575": "NH1" <-> "NH2" Residue "A ARG 645": "NH1" <-> "NH2" Residue "A GLU 667": "OE1" <-> "OE2" Residue "A GLU 718": "OE1" <-> "OE2" Residue "A ARG 729": "NH1" <-> "NH2" Residue "A GLU 768": "OE1" <-> "OE2" Residue "A GLU 800": "OE1" <-> "OE2" Residue "A TYR 816": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 878": "OE1" <-> "OE2" Residue "A ARG 890": "NH1" <-> "NH2" Residue "A GLU 896": "OE1" <-> "OE2" Residue "A ARG 925": "NH1" <-> "NH2" Residue "A ARG 960": "NH1" <-> "NH2" Residue "A GLU 971": "OE1" <-> "OE2" Residue "A ARG 1080": "NH1" <-> "NH2" Residue "A ARG 1104": "NH1" <-> "NH2" Residue "A GLU 1119": "OE1" <-> "OE2" Residue "A ARG 1147": "NH1" <-> "NH2" Time to flip residues: 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5172/modules/chem_data/mon_lib" Total number of atoms: 9284 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 9284 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1170, 9284 Classifications: {'peptide': 1170} Link IDs: {'PCIS': 4, 'PTRANS': 38, 'TRANS': 1127} Time building chain proxies: 5.49, per 1000 atoms: 0.59 Number of scatterers: 9284 At special positions: 0 Unit cell: (92.394, 98.766, 210.276, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 4 Type Number sf(0) S 39 16.00 O 1800 8.00 N 1629 7.00 C 5816 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 3.23 Conformation dependent library (CDL) restraints added in 1.9 seconds 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2208 Finding SS restraints... Secondary structure from input PDB file: 31 helices and 9 sheets defined 35.3% alpha, 21.9% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.20 Creating SS restraints... Processing helix chain 'A' and resid 28 through 53 Processing helix chain 'A' and resid 58 through 74 Processing helix chain 'A' and resid 86 through 97 Processing helix chain 'A' and resid 241 through 249 Processing helix chain 'A' and resid 306 through 338 removed outlier: 3.910A pdb=" N HIS A 317 " --> pdb=" O ALA A 313 " (cutoff:3.500A) removed outlier: 3.596A pdb=" N SER A 338 " --> pdb=" O SER A 334 " (cutoff:3.500A) Processing helix chain 'A' and resid 349 through 372 removed outlier: 3.806A pdb=" N ARG A 353 " --> pdb=" O ASP A 349 " (cutoff:3.500A) Processing helix chain 'A' and resid 373 through 377 Processing helix chain 'A' and resid 652 through 658 Processing helix chain 'A' and resid 682 through 688 removed outlier: 3.688A pdb=" N ASN A 688 " --> pdb=" O TYR A 684 " (cutoff:3.500A) Processing helix chain 'A' and resid 745 through 754 Processing helix chain 'A' and resid 756 through 772 removed outlier: 4.306A pdb=" N GLY A 770 " --> pdb=" O SER A 766 " (cutoff:3.500A) removed outlier: 3.941A pdb=" N ARG A 771 " --> pdb=" O GLN A 767 " (cutoff:3.500A) Processing helix chain 'A' and resid 774 through 777 removed outlier: 3.627A pdb=" N LEU A 777 " --> pdb=" O ASP A 774 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 774 through 777' Processing helix chain 'A' and resid 778 through 796 removed outlier: 3.611A pdb=" N SER A 796 " --> pdb=" O MET A 792 " (cutoff:3.500A) Processing helix chain 'A' and resid 800 through 815 removed outlier: 3.690A pdb=" N ILE A 804 " --> pdb=" O GLU A 800 " (cutoff:3.500A) Processing helix chain 'A' and resid 820 through 835 Processing helix chain 'A' and resid 835 through 850 removed outlier: 3.578A pdb=" N THR A 839 " --> pdb=" O ALA A 835 " (cutoff:3.500A) Processing helix chain 'A' and resid 852 through 861 removed outlier: 4.090A pdb=" N GLY A 861 " --> pdb=" O THR A 857 " (cutoff:3.500A) Processing helix chain 'A' and resid 867 through 885 Processing helix chain 'A' and resid 887 through 891 Processing helix chain 'A' and resid 897 through 901 Processing helix chain 'A' and resid 904 through 909 Processing helix chain 'A' and resid 909 through 923 removed outlier: 3.606A pdb=" N PHE A 922 " --> pdb=" O MET A 918 " (cutoff:3.500A) Processing helix chain 'A' and resid 932 through 938 removed outlier: 3.873A pdb=" N ARG A 938 " --> pdb=" O ASP A 934 " (cutoff:3.500A) Processing helix chain 'A' and resid 979 through 994 Processing helix chain 'A' and resid 1036 through 1040 Processing helix chain 'A' and resid 1074 through 1087 Processing helix chain 'A' and resid 1094 through 1107 removed outlier: 3.528A pdb=" N VAL A1106 " --> pdb=" O ILE A1102 " (cutoff:3.500A) removed outlier: 3.709A pdb=" N ASP A1107 " --> pdb=" O SER A1103 " (cutoff:3.500A) Processing helix chain 'A' and resid 1123 through 1136 removed outlier: 3.593A pdb=" N LEU A1133 " --> pdb=" O MET A1129 " (cutoff:3.500A) Processing helix chain 'A' and resid 1141 through 1146 Processing helix chain 'A' and resid 1148 through 1161 removed outlier: 4.119A pdb=" N PHE A1152 " --> pdb=" O ASN A1148 " (cutoff:3.500A) removed outlier: 4.032A pdb=" N VAL A1153 " --> pdb=" O CYS A1149 " (cutoff:3.500A) Processing helix chain 'A' and resid 1170 through 1194 removed outlier: 4.322A pdb=" N ALA A1174 " --> pdb=" O ASP A1170 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 83 through 85 removed outlier: 3.949A pdb=" N ASN A 105 " --> pdb=" O GLU A 295 " (cutoff:3.500A) removed outlier: 3.720A pdb=" N VAL A 106 " --> pdb=" O LEU A 192 " (cutoff:3.500A) removed outlier: 6.498A pdb=" N GLU A 231 " --> pdb=" O LEU A 263 " (cutoff:3.500A) removed outlier: 4.521A pdb=" N LEU A 263 " --> pdb=" O GLU A 231 " (cutoff:3.500A) removed outlier: 6.203A pdb=" N LEU A 233 " --> pdb=" O LEU A 261 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 171 through 181 removed outlier: 6.455A pdb=" N ASP A 155 " --> pdb=" O TYR A 149 " (cutoff:3.500A) removed outlier: 4.115A pdb=" N TYR A 149 " --> pdb=" O ASP A 155 " (cutoff:3.500A) removed outlier: 6.673A pdb=" N GLU A 157 " --> pdb=" O LEU A 147 " (cutoff:3.500A) removed outlier: 4.174A pdb=" N GLN A 124 " --> pdb=" O ILE A 144 " (cutoff:3.500A) removed outlier: 6.695A pdb=" N ARG A 125 " --> pdb=" O ILE A 284 " (cutoff:3.500A) removed outlier: 4.513A pdb=" N ILE A 284 " --> pdb=" O ARG A 125 " (cutoff:3.500A) removed outlier: 6.717A pdb=" N PHE A 127 " --> pdb=" O GLU A 282 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 213 through 214 removed outlier: 4.133A pdb=" N ASN A1053 " --> pdb=" O ILE A1062 " (cutoff:3.500A) removed outlier: 7.073A pdb=" N ILE A1062 " --> pdb=" O ASN A1053 " (cutoff:3.500A) removed outlier: 6.763A pdb=" N LEU A 300 " --> pdb=" O ILE A1063 " (cutoff:3.500A) removed outlier: 8.134A pdb=" N ARG A1065 " --> pdb=" O LEU A 300 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'A' and resid 213 through 214 removed outlier: 6.731A pdb=" N SER A1016 " --> pdb=" O GLU A1012 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'A' and resid 379 through 384 removed outlier: 4.160A pdb=" N LEU A 379 " --> pdb=" O TYR A 740 " (cutoff:3.500A) removed outlier: 3.572A pdb=" N ALA A 736 " --> pdb=" O ALA A 383 " (cutoff:3.500A) removed outlier: 3.694A pdb=" N PHE A 394 " --> pdb=" O LEU A 726 " (cutoff:3.500A) removed outlier: 3.958A pdb=" N PHE A 728 " --> pdb=" O GLY A 392 " (cutoff:3.500A) removed outlier: 4.206A pdb=" N GLY A 392 " --> pdb=" O PHE A 728 " (cutoff:3.500A) removed outlier: 3.651A pdb=" N ILE A 730 " --> pdb=" O THR A 390 " (cutoff:3.500A) removed outlier: 5.244A pdb=" N GLY A 732 " --> pdb=" O ASP A 388 " (cutoff:3.500A) removed outlier: 4.840A pdb=" N ASP A 388 " --> pdb=" O GLY A 732 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'A' and resid 379 through 384 removed outlier: 4.160A pdb=" N LEU A 379 " --> pdb=" O TYR A 740 " (cutoff:3.500A) removed outlier: 3.572A pdb=" N ALA A 736 " --> pdb=" O ALA A 383 " (cutoff:3.500A) removed outlier: 3.762A pdb=" N ARG A 729 " --> pdb=" O GLU A 567 " (cutoff:3.500A) removed outlier: 4.583A pdb=" N GLU A 567 " --> pdb=" O ARG A 729 " (cutoff:3.500A) removed outlier: 5.069A pdb=" N HIS A 731 " --> pdb=" O SER A 565 " (cutoff:3.500A) removed outlier: 4.781A pdb=" N SER A 565 " --> pdb=" O HIS A 731 " (cutoff:3.500A) removed outlier: 3.952A pdb=" N ASN A 733 " --> pdb=" O ARG A 563 " (cutoff:3.500A) removed outlier: 4.110A pdb=" N ARG A 737 " --> pdb=" O ASP A 559 " (cutoff:3.500A) removed outlier: 5.747A pdb=" N ASP A 559 " --> pdb=" O ARG A 737 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'A' and resid 407 through 408 removed outlier: 3.613A pdb=" N TYR A 407 " --> pdb=" O PHE A 440 " (cutoff:3.500A) removed outlier: 3.775A pdb=" N ILE A 439 " --> pdb=" O PHE A 676 " (cutoff:3.500A) removed outlier: 3.623A pdb=" N LYS A 441 " --> pdb=" O VAL A 678 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'A' and resid 451 through 453 removed outlier: 3.819A pdb=" N PHE A 468 " --> pdb=" O HIS A 452 " (cutoff:3.500A) removed outlier: 3.638A pdb=" N ARG A 467 " --> pdb=" O GLU A 545 " (cutoff:3.500A) removed outlier: 4.920A pdb=" N GLU A 545 " --> pdb=" O ARG A 467 " (cutoff:3.500A) removed outlier: 7.227A pdb=" N SER A 469 " --> pdb=" O VAL A 543 " (cutoff:3.500A) removed outlier: 5.797A pdb=" N VAL A 543 " --> pdb=" O SER A 469 " (cutoff:3.500A) removed outlier: 3.689A pdb=" N GLU A 545 " --> pdb=" O ILE A 602 " (cutoff:3.500A) removed outlier: 6.246A pdb=" N ILE A 602 " --> pdb=" O GLU A 545 " (cutoff:3.500A) removed outlier: 6.332A pdb=" N LEU A 547 " --> pdb=" O CYS A 600 " (cutoff:3.500A) removed outlier: 7.748A pdb=" N CYS A 600 " --> pdb=" O LEU A 547 " (cutoff:3.500A) removed outlier: 6.009A pdb=" N GLY A 549 " --> pdb=" O LEU A 598 " (cutoff:3.500A) removed outlier: 7.241A pdb=" N LEU A 598 " --> pdb=" O GLY A 549 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'A' and resid 864 through 865 removed outlier: 6.633A pdb=" N ASN A 864 " --> pdb=" O GLU A1119 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA9 404 hydrogen bonds defined for protein. 1164 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 0 stacking parallelities Total time for adding SS restraints: 3.31 Time building geometry restraints manager: 3.56 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.21 - 1.33: 2732 1.33 - 1.45: 1440 1.45 - 1.57: 5209 1.57 - 1.69: 0 1.69 - 1.81: 74 Bond restraints: 9455 Sorted by residual: bond pdb=" CB ASN A1030 " pdb=" CG ASN A1030 " ideal model delta sigma weight residual 1.516 1.472 0.044 2.50e-02 1.60e+03 3.11e+00 bond pdb=" C GLU A 679 " pdb=" N LEU A 680 " ideal model delta sigma weight residual 1.325 1.345 -0.020 1.25e-02 6.40e+03 2.67e+00 bond pdb=" CG GLN A1179 " pdb=" CD GLN A1179 " ideal model delta sigma weight residual 1.516 1.477 0.039 2.50e-02 1.60e+03 2.42e+00 bond pdb=" CB ASP A 182 " pdb=" CG ASP A 182 " ideal model delta sigma weight residual 1.516 1.477 0.039 2.50e-02 1.60e+03 2.38e+00 bond pdb=" CB ASP A 756 " pdb=" CG ASP A 756 " ideal model delta sigma weight residual 1.516 1.478 0.038 2.50e-02 1.60e+03 2.34e+00 ... (remaining 9450 not shown) Histogram of bond angle deviations from ideal: 97.44 - 104.78: 138 104.78 - 112.13: 4725 112.13 - 119.47: 3115 119.47 - 126.81: 4764 126.81 - 134.16: 86 Bond angle restraints: 12828 Sorted by residual: angle pdb=" N GLY A 815 " pdb=" CA GLY A 815 " pdb=" C GLY A 815 " ideal model delta sigma weight residual 114.92 120.28 -5.36 1.30e+00 5.92e-01 1.70e+01 angle pdb=" CA LEU A 777 " pdb=" C LEU A 777 " pdb=" N ASN A 778 " ideal model delta sigma weight residual 118.27 112.71 5.56 1.59e+00 3.96e-01 1.22e+01 angle pdb=" C LEU A 240 " pdb=" N HIS A 241 " pdb=" CA HIS A 241 " ideal model delta sigma weight residual 122.21 117.99 4.22 1.48e+00 4.57e-01 8.11e+00 angle pdb=" N VAL A 532 " pdb=" CA VAL A 532 " pdb=" C VAL A 532 " ideal model delta sigma weight residual 107.71 111.90 -4.19 1.58e+00 4.01e-01 7.02e+00 angle pdb=" N ASP A 629 " pdb=" CA ASP A 629 " pdb=" C ASP A 629 " ideal model delta sigma weight residual 111.36 114.08 -2.72 1.09e+00 8.42e-01 6.24e+00 ... (remaining 12823 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 16.96: 5139 16.96 - 33.91: 489 33.91 - 50.86: 77 50.86 - 67.82: 18 67.82 - 84.77: 7 Dihedral angle restraints: 5730 sinusoidal: 2307 harmonic: 3423 Sorted by residual: dihedral pdb=" CA ILE A 730 " pdb=" C ILE A 730 " pdb=" N HIS A 731 " pdb=" CA HIS A 731 " ideal model delta harmonic sigma weight residual 180.00 156.63 23.37 0 5.00e+00 4.00e-02 2.18e+01 dihedral pdb=" CA GLN A 755 " pdb=" C GLN A 755 " pdb=" N ASP A 756 " pdb=" CA ASP A 756 " ideal model delta harmonic sigma weight residual -180.00 -159.28 -20.72 0 5.00e+00 4.00e-02 1.72e+01 dihedral pdb=" CA HIS A 136 " pdb=" C HIS A 136 " pdb=" N PHE A 137 " pdb=" CA PHE A 137 " ideal model delta harmonic sigma weight residual -180.00 -160.34 -19.66 0 5.00e+00 4.00e-02 1.55e+01 ... (remaining 5727 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.040: 831 0.040 - 0.081: 448 0.081 - 0.121: 151 0.121 - 0.162: 32 0.162 - 0.202: 2 Chirality restraints: 1464 Sorted by residual: chirality pdb=" CB THR A1007 " pdb=" CA THR A1007 " pdb=" OG1 THR A1007 " pdb=" CG2 THR A1007 " both_signs ideal model delta sigma weight residual False 2.55 2.35 0.20 2.00e-01 2.50e+01 1.02e+00 chirality pdb=" CG LEU A 240 " pdb=" CB LEU A 240 " pdb=" CD1 LEU A 240 " pdb=" CD2 LEU A 240 " both_signs ideal model delta sigma weight residual False -2.59 -2.43 -0.16 2.00e-01 2.50e+01 6.76e-01 chirality pdb=" CB ILE A1063 " pdb=" CA ILE A1063 " pdb=" CG1 ILE A1063 " pdb=" CG2 ILE A1063 " both_signs ideal model delta sigma weight residual False 2.64 2.49 0.16 2.00e-01 2.50e+01 6.16e-01 ... (remaining 1461 not shown) Planarity restraints: 1674 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C LEU A 952 " -0.053 5.00e-02 4.00e+02 8.08e-02 1.04e+01 pdb=" N PRO A 953 " 0.140 5.00e-02 4.00e+02 pdb=" CA PRO A 953 " -0.042 5.00e-02 4.00e+02 pdb=" CD PRO A 953 " -0.044 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CB TYR A 234 " -0.011 2.00e-02 2.50e+03 1.26e-02 3.19e+00 pdb=" CG TYR A 234 " 0.031 2.00e-02 2.50e+03 pdb=" CD1 TYR A 234 " -0.012 2.00e-02 2.50e+03 pdb=" CD2 TYR A 234 " -0.004 2.00e-02 2.50e+03 pdb=" CE1 TYR A 234 " 0.003 2.00e-02 2.50e+03 pdb=" CE2 TYR A 234 " -0.006 2.00e-02 2.50e+03 pdb=" CZ TYR A 234 " -0.000 2.00e-02 2.50e+03 pdb=" OH TYR A 234 " 0.000 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C CYS A 600 " 0.029 5.00e-02 4.00e+02 4.42e-02 3.13e+00 pdb=" N PRO A 601 " -0.077 5.00e-02 4.00e+02 pdb=" CA PRO A 601 " 0.023 5.00e-02 4.00e+02 pdb=" CD PRO A 601 " 0.025 5.00e-02 4.00e+02 ... (remaining 1671 not shown) Histogram of nonbonded interaction distances: 2.15 - 2.70: 280 2.70 - 3.25: 8660 3.25 - 3.80: 13899 3.80 - 4.35: 17952 4.35 - 4.90: 30950 Nonbonded interactions: 71741 Sorted by model distance: nonbonded pdb=" OH TYR A 176 " pdb=" OG SER A 237 " model vdw 2.147 2.440 nonbonded pdb=" O SER A 614 " pdb=" NH2 ARG A 698 " model vdw 2.207 2.520 nonbonded pdb=" O ILE A 546 " pdb=" OG1 THR A 599 " model vdw 2.247 2.440 nonbonded pdb=" NE2 GLN A 94 " pdb=" O LEU A 276 " model vdw 2.249 2.520 nonbonded pdb=" O ARG A 721 " pdb=" N THR A 725 " model vdw 2.270 2.520 ... (remaining 71736 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=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: 3.120 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: 8.000 Check model and map are aligned: 0.130 Set scattering table: 0.090 Process input model: 28.280 Find NCS groups from input model: 0.270 Set up NCS constraints: 0.030 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.010 Load rotamer database and sin/cos tables:2.910 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 42.850 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8054 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.057 9455 Z= 0.577 Angle : 0.800 7.926 12828 Z= 0.470 Chirality : 0.053 0.202 1464 Planarity : 0.005 0.081 1674 Dihedral : 13.811 84.774 3522 Min Nonbonded Distance : 2.147 Molprobity Statistics. All-atom Clashscore : 7.64 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.54 % Favored : 95.46 % Rotamer: Outliers : 0.10 % Allowed : 0.10 % Favored : 99.81 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % 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.29 (0.22), residues: 1168 helix: 0.45 (0.26), residues: 373 sheet: -0.68 (0.32), residues: 241 loop : -1.87 (0.22), residues: 554 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.017 0.002 TRP A1099 HIS 0.007 0.002 HIS A 317 PHE 0.027 0.003 PHE A1140 TYR 0.031 0.003 TYR A 234 ARG 0.014 0.001 ARG A 219 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 105 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 104 time to evaluate : 1.094 Fit side-chains outliers start: 1 outliers final: 0 residues processed: 104 average time/residue: 0.3284 time to fit residues: 46.0772 Evaluate side-chains 95 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 95 time to evaluate : 0.989 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.5702 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 117 random chunks: chunk 98 optimal weight: 4.9990 chunk 88 optimal weight: 0.9980 chunk 49 optimal weight: 0.0010 chunk 30 optimal weight: 0.9990 chunk 59 optimal weight: 20.0000 chunk 47 optimal weight: 2.9990 chunk 91 optimal weight: 0.9980 chunk 35 optimal weight: 10.0000 chunk 55 optimal weight: 2.9990 chunk 68 optimal weight: 0.9990 chunk 106 optimal weight: 0.8980 overall best weight: 0.7788 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 335 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 452 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 767 GLN A 907 HIS A 924 ASN A 927 ASN ** A 993 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1030 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A1167 ASN A1179 GLN Total number of N/Q/H flips: 6 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7951 moved from start: 0.1161 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.032 9455 Z= 0.148 Angle : 0.494 5.833 12828 Z= 0.269 Chirality : 0.043 0.146 1464 Planarity : 0.004 0.063 1674 Dihedral : 4.783 17.489 1298 Min Nonbonded Distance : 2.189 Molprobity Statistics. All-atom Clashscore : 6.01 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.40 % Favored : 97.60 % Rotamer: Outliers : 0.29 % Allowed : 4.91 % Favored : 94.80 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.35 (0.24), residues: 1168 helix: 1.64 (0.27), residues: 383 sheet: -0.52 (0.33), residues: 231 loop : -1.49 (0.23), residues: 554 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A1099 HIS 0.003 0.001 HIS A 317 PHE 0.014 0.001 PHE A 128 TYR 0.016 0.001 TYR A1050 ARG 0.004 0.000 ARG A 381 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 116 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 113 time to evaluate : 1.120 Fit side-chains revert: symmetry clash outliers start: 3 outliers final: 3 residues processed: 113 average time/residue: 0.2652 time to fit residues: 40.6151 Evaluate side-chains 104 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 101 time to evaluate : 1.190 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0915 time to fit residues: 2.0573 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 117 random chunks: chunk 59 optimal weight: 20.0000 chunk 32 optimal weight: 9.9990 chunk 88 optimal weight: 0.6980 chunk 72 optimal weight: 0.9980 chunk 29 optimal weight: 7.9990 chunk 106 optimal weight: 6.9990 chunk 114 optimal weight: 6.9990 chunk 94 optimal weight: 4.9990 chunk 105 optimal weight: 1.9990 chunk 36 optimal weight: 6.9990 chunk 85 optimal weight: 0.9990 overall best weight: 1.9386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 329 GLN ** A 452 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 767 GLN A 864 ASN ** A 993 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1030 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1053 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A1179 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7985 moved from start: 0.1263 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.036 9455 Z= 0.224 Angle : 0.500 5.798 12828 Z= 0.269 Chirality : 0.044 0.141 1464 Planarity : 0.004 0.058 1674 Dihedral : 4.627 16.853 1298 Min Nonbonded Distance : 2.216 Molprobity Statistics. All-atom Clashscore : 5.15 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.91 % Favored : 97.09 % Rotamer: Outliers : 0.96 % Allowed : 7.03 % Favored : 92.01 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.15 (0.24), residues: 1168 helix: 1.80 (0.27), residues: 388 sheet: -0.20 (0.32), residues: 237 loop : -1.51 (0.24), residues: 543 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.001 TRP A1099 HIS 0.005 0.001 HIS A 317 PHE 0.016 0.001 PHE A1003 TYR 0.015 0.001 TYR A 234 ARG 0.003 0.000 ARG A 219 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 116 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 106 time to evaluate : 1.152 Fit side-chains revert: symmetry clash outliers start: 10 outliers final: 5 residues processed: 111 average time/residue: 0.2662 time to fit residues: 39.9722 Evaluate side-chains 108 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 103 time to evaluate : 0.967 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.0965 time to fit residues: 2.3181 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 117 random chunks: chunk 105 optimal weight: 0.7980 chunk 80 optimal weight: 6.9990 chunk 55 optimal weight: 7.9990 chunk 11 optimal weight: 2.9990 chunk 50 optimal weight: 2.9990 chunk 71 optimal weight: 5.9990 chunk 106 optimal weight: 3.9990 chunk 113 optimal weight: 0.7980 chunk 101 optimal weight: 0.7980 chunk 30 optimal weight: 0.8980 chunk 94 optimal weight: 8.9990 overall best weight: 1.2582 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 452 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 754 HIS A 767 GLN ** A 993 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A1053 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7961 moved from start: 0.1528 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.034 9455 Z= 0.165 Angle : 0.458 5.676 12828 Z= 0.247 Chirality : 0.042 0.138 1464 Planarity : 0.004 0.054 1674 Dihedral : 4.370 15.192 1298 Min Nonbonded Distance : 2.251 Molprobity Statistics. All-atom Clashscore : 5.58 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.74 % Favored : 97.26 % Rotamer: Outliers : 1.06 % Allowed : 8.57 % Favored : 90.38 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.13 (0.24), residues: 1168 helix: 2.07 (0.27), residues: 388 sheet: 0.04 (0.33), residues: 222 loop : -1.39 (0.24), residues: 558 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A1099 HIS 0.004 0.001 HIS A 317 PHE 0.015 0.001 PHE A1003 TYR 0.014 0.001 TYR A1050 ARG 0.004 0.000 ARG A 219 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 119 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 108 time to evaluate : 1.228 Fit side-chains revert: symmetry clash outliers start: 11 outliers final: 2 residues processed: 116 average time/residue: 0.2947 time to fit residues: 48.4288 Evaluate side-chains 104 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 102 time to evaluate : 1.144 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0853 time to fit residues: 1.7656 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 117 random chunks: chunk 64 optimal weight: 0.8980 chunk 1 optimal weight: 0.9990 chunk 84 optimal weight: 0.8980 chunk 46 optimal weight: 7.9990 chunk 96 optimal weight: 8.9990 chunk 78 optimal weight: 0.9990 chunk 0 optimal weight: 0.8980 chunk 57 optimal weight: 20.0000 chunk 101 optimal weight: 1.9990 chunk 28 optimal weight: 0.8980 chunk 38 optimal weight: 8.9990 overall best weight: 0.9182 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 329 GLN ** A 452 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 767 GLN A 993 ASN A1053 ASN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7937 moved from start: 0.1815 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.030 9455 Z= 0.138 Angle : 0.432 5.585 12828 Z= 0.233 Chirality : 0.042 0.134 1464 Planarity : 0.003 0.052 1674 Dihedral : 4.105 13.920 1298 Min Nonbonded Distance : 2.282 Molprobity Statistics. All-atom Clashscore : 5.31 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.48 % Favored : 97.52 % Rotamer: Outliers : 0.96 % Allowed : 9.82 % Favored : 89.22 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.45 (0.24), residues: 1168 helix: 2.35 (0.27), residues: 385 sheet: 0.15 (0.33), residues: 231 loop : -1.18 (0.24), residues: 552 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP A1099 HIS 0.003 0.000 HIS A 317 PHE 0.014 0.001 PHE A1003 TYR 0.013 0.001 TYR A1050 ARG 0.005 0.000 ARG A 219 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 116 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 106 time to evaluate : 0.994 Fit side-chains revert: symmetry clash outliers start: 10 outliers final: 1 residues processed: 114 average time/residue: 0.2473 time to fit residues: 38.8077 Evaluate side-chains 104 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 103 time to evaluate : 1.033 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0838 time to fit residues: 1.6355 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 117 random chunks: chunk 101 optimal weight: 4.9990 chunk 22 optimal weight: 2.9990 chunk 66 optimal weight: 6.9990 chunk 27 optimal weight: 0.3980 chunk 113 optimal weight: 0.6980 chunk 93 optimal weight: 2.9990 chunk 52 optimal weight: 5.9990 chunk 9 optimal weight: 0.3980 chunk 37 optimal weight: 30.0000 chunk 59 optimal weight: 6.9990 chunk 109 optimal weight: 5.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 452 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 767 GLN A 924 ASN A 927 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7957 moved from start: 0.1837 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 9455 Z= 0.181 Angle : 0.462 7.539 12828 Z= 0.245 Chirality : 0.042 0.135 1464 Planarity : 0.004 0.051 1674 Dihedral : 4.146 14.554 1298 Min Nonbonded Distance : 2.247 Molprobity Statistics. All-atom Clashscore : 5.42 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.57 % Favored : 97.43 % Rotamer: Outliers : 1.25 % Allowed : 10.30 % Favored : 88.45 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.55 (0.24), residues: 1168 helix: 2.47 (0.27), residues: 379 sheet: 0.16 (0.33), residues: 229 loop : -1.08 (0.24), residues: 560 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A1099 HIS 0.004 0.001 HIS A 317 PHE 0.015 0.001 PHE A1003 TYR 0.014 0.001 TYR A1050 ARG 0.002 0.000 ARG A 341 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 118 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 105 time to evaluate : 1.265 Fit side-chains revert: symmetry clash outliers start: 13 outliers final: 6 residues processed: 111 average time/residue: 0.2788 time to fit residues: 41.9161 Evaluate side-chains 109 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 103 time to evaluate : 1.034 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 0 residues processed: 6 average time/residue: 0.1057 time to fit residues: 2.6733 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 117 random chunks: chunk 12 optimal weight: 6.9990 chunk 64 optimal weight: 2.9990 chunk 82 optimal weight: 0.9990 chunk 95 optimal weight: 0.9990 chunk 63 optimal weight: 10.0000 chunk 112 optimal weight: 2.9990 chunk 70 optimal weight: 10.0000 chunk 68 optimal weight: 20.0000 chunk 52 optimal weight: 0.7980 chunk 69 optimal weight: 6.9990 chunk 45 optimal weight: 20.0000 overall best weight: 1.7588 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 452 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 767 GLN A1053 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7965 moved from start: 0.1855 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.035 9455 Z= 0.201 Angle : 0.465 7.096 12828 Z= 0.246 Chirality : 0.043 0.136 1464 Planarity : 0.004 0.050 1674 Dihedral : 4.177 15.177 1298 Min Nonbonded Distance : 2.260 Molprobity Statistics. All-atom Clashscore : 5.36 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.91 % Favored : 97.09 % Rotamer: Outliers : 0.48 % Allowed : 10.97 % Favored : 88.55 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.56 (0.24), residues: 1168 helix: 2.43 (0.27), residues: 379 sheet: 0.16 (0.33), residues: 229 loop : -1.03 (0.24), residues: 560 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A1099 HIS 0.005 0.001 HIS A 317 PHE 0.016 0.001 PHE A1003 TYR 0.015 0.001 TYR A1050 ARG 0.002 0.000 ARG A 341 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 109 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 104 time to evaluate : 1.062 Fit side-chains revert: symmetry clash outliers start: 5 outliers final: 0 residues processed: 107 average time/residue: 0.2799 time to fit residues: 40.5578 Evaluate side-chains 104 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 104 time to evaluate : 1.091 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.5317 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 117 random chunks: chunk 67 optimal weight: 0.0770 chunk 33 optimal weight: 3.9990 chunk 22 optimal weight: 0.0020 chunk 21 optimal weight: 0.9990 chunk 71 optimal weight: 5.9990 chunk 76 optimal weight: 2.9990 chunk 55 optimal weight: 7.9990 chunk 10 optimal weight: 0.9990 chunk 88 optimal weight: 0.5980 chunk 102 optimal weight: 0.8980 chunk 108 optimal weight: 0.9980 overall best weight: 0.5146 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 138 GLN ** A 452 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 767 GLN A1053 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7917 moved from start: 0.2135 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.026 9455 Z= 0.110 Angle : 0.428 6.617 12828 Z= 0.227 Chirality : 0.041 0.133 1464 Planarity : 0.003 0.049 1674 Dihedral : 3.935 13.797 1298 Min Nonbonded Distance : 2.267 Molprobity Statistics. All-atom Clashscore : 5.42 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.23 % Favored : 97.77 % Rotamer: Outliers : 0.67 % Allowed : 10.78 % Favored : 88.55 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.78 (0.25), residues: 1168 helix: 2.67 (0.27), residues: 379 sheet: -0.11 (0.34), residues: 236 loop : -0.74 (0.25), residues: 553 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP A1099 HIS 0.003 0.000 HIS A 317 PHE 0.013 0.001 PHE A 128 TYR 0.014 0.001 TYR A1050 ARG 0.005 0.000 ARG A1181 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 112 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 105 time to evaluate : 1.017 Fit side-chains revert: symmetry clash outliers start: 7 outliers final: 3 residues processed: 108 average time/residue: 0.2636 time to fit residues: 38.5170 Evaluate side-chains 104 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 101 time to evaluate : 1.029 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0863 time to fit residues: 1.9404 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 117 random chunks: chunk 98 optimal weight: 0.9990 chunk 105 optimal weight: 2.9990 chunk 108 optimal weight: 9.9990 chunk 63 optimal weight: 10.0000 chunk 45 optimal weight: 20.0000 chunk 82 optimal weight: 0.8980 chunk 32 optimal weight: 9.9990 chunk 94 optimal weight: 8.9990 chunk 99 optimal weight: 2.9990 chunk 104 optimal weight: 0.8980 chunk 69 optimal weight: 30.0000 overall best weight: 1.7586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 452 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 767 GLN ** A1030 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7960 moved from start: 0.2015 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.034 9455 Z= 0.204 Angle : 0.470 6.871 12828 Z= 0.249 Chirality : 0.043 0.262 1464 Planarity : 0.004 0.049 1674 Dihedral : 4.091 15.294 1298 Min Nonbonded Distance : 2.241 Molprobity Statistics. All-atom Clashscore : 5.31 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.91 % Favored : 97.09 % Rotamer: Outliers : 0.10 % Allowed : 11.55 % Favored : 88.35 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.68 (0.25), residues: 1168 helix: 2.51 (0.27), residues: 379 sheet: -0.00 (0.33), residues: 241 loop : -0.83 (0.25), residues: 548 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A1099 HIS 0.005 0.001 HIS A 317 PHE 0.016 0.001 PHE A1003 TYR 0.013 0.001 TYR A1050 ARG 0.004 0.000 ARG A1181 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 107 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 106 time to evaluate : 1.066 Fit side-chains revert: symmetry clash outliers start: 1 outliers final: 0 residues processed: 107 average time/residue: 0.2778 time to fit residues: 40.3357 Evaluate side-chains 105 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 105 time to evaluate : 1.065 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.5872 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 117 random chunks: chunk 111 optimal weight: 3.9990 chunk 67 optimal weight: 0.0050 chunk 52 optimal weight: 0.9990 chunk 77 optimal weight: 0.9990 chunk 116 optimal weight: 7.9990 chunk 107 optimal weight: 0.9990 chunk 92 optimal weight: 0.9990 chunk 9 optimal weight: 7.9990 chunk 71 optimal weight: 5.9990 chunk 56 optimal weight: 7.9990 chunk 73 optimal weight: 0.5980 overall best weight: 0.7200 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 452 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 535 GLN A 767 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7916 moved from start: 0.2284 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.027 9455 Z= 0.122 Angle : 0.428 6.279 12828 Z= 0.228 Chirality : 0.042 0.273 1464 Planarity : 0.003 0.048 1674 Dihedral : 3.885 14.358 1298 Min Nonbonded Distance : 2.270 Molprobity Statistics. All-atom Clashscore : 5.69 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.48 % Favored : 97.52 % Rotamer: Outliers : 0.29 % Allowed : 11.74 % Favored : 87.97 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.85 (0.25), residues: 1168 helix: 2.68 (0.27), residues: 379 sheet: -0.07 (0.34), residues: 236 loop : -0.65 (0.25), residues: 553 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.009 0.001 TRP A1099 HIS 0.003 0.000 HIS A 317 PHE 0.013 0.001 PHE A 128 TYR 0.012 0.001 TYR A1050 ARG 0.002 0.000 ARG A 341 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2336 Ramachandran restraints generated. 1168 Oldfield, 0 Emsley, 1168 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 106 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 103 time to evaluate : 0.989 Fit side-chains revert: symmetry clash outliers start: 3 outliers final: 0 residues processed: 104 average time/residue: 0.2713 time to fit residues: 38.2049 Evaluate side-chains 103 residues out of total 1039 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 103 time to evaluate : 1.017 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.4304 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 117 random chunks: chunk 98 optimal weight: 1.9990 chunk 28 optimal weight: 0.9980 chunk 85 optimal weight: 6.9990 chunk 13 optimal weight: 1.9990 chunk 25 optimal weight: 5.9990 chunk 92 optimal weight: 0.0980 chunk 38 optimal weight: 0.0670 chunk 95 optimal weight: 1.9990 chunk 11 optimal weight: 3.9990 chunk 17 optimal weight: 0.9990 chunk 81 optimal weight: 3.9990 overall best weight: 0.8322 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 438 ASN ** A 452 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3952 r_free = 0.3952 target = 0.188697 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 37)----------------| | r_work = 0.3341 r_free = 0.3341 target = 0.130033 restraints weight = 10346.241| |-----------------------------------------------------------------------------| r_work (start): 0.3275 rms_B_bonded: 1.97 r_work: 0.3020 rms_B_bonded: 1.95 restraints_weight: 0.5000 r_work: 0.2911 rms_B_bonded: 3.30 restraints_weight: 0.2500 r_work (final): 0.2911 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7907 moved from start: 0.2344 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.026 9455 Z= 0.129 Angle : 0.427 6.513 12828 Z= 0.226 Chirality : 0.042 0.261 1464 Planarity : 0.003 0.048 1674 Dihedral : 3.822 14.001 1298 Min Nonbonded Distance : 2.281 Molprobity Statistics. All-atom Clashscore : 5.47 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.57 % Favored : 97.43 % Rotamer: Outliers : 0.19 % Allowed : 11.84 % Favored : 87.97 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 9.52 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.94 (0.25), residues: 1168 helix: 2.71 (0.27), residues: 379 sheet: -0.02 (0.34), residues: 236 loop : -0.56 (0.25), residues: 553 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP A1099 HIS 0.004 0.000 HIS A 317 PHE 0.013 0.001 PHE A1003 TYR 0.012 0.001 TYR A1050 ARG 0.002 0.000 ARG A 341 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 2248.45 seconds wall clock time: 41 minutes 26.21 seconds (2486.21 seconds total)