Starting phenix.real_space_refine (version: 1.21rc1) on Tue Jul 4 04:26:21 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8erl_28554/07_2023/8erl_28554.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8erl_28554/07_2023/8erl_28554.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.9 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8erl_28554/07_2023/8erl_28554.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8erl_28554/07_2023/8erl_28554.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8erl_28554/07_2023/8erl_28554.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8erl_28554/07_2023/8erl_28554.pdb" } resolution = 3.9 write_initial_geo_file = False refinement { macro_cycles = 10 } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= -0.000 sd= 0.008 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 36 5.16 5 C 4208 2.51 5 N 1130 2.21 5 O 1226 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped None Time to flip residues: 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5013/modules/chem_data/mon_lib" Total number of atoms: 6600 Number of models: 1 Model: "" Number of chains: 2 Chain: "A" Number of atoms: 3300 Number of conformers: 1 Conformer: "" Number of residues, atoms: 415, 3300 Classifications: {'peptide': 415} Link IDs: {'PTRANS': 20, 'TRANS': 394} Chain breaks: 1 Chain: "B" Number of atoms: 3300 Number of conformers: 1 Conformer: "" Number of residues, atoms: 415, 3300 Classifications: {'peptide': 415} Link IDs: {'PTRANS': 20, 'TRANS': 394} Chain breaks: 1 Time building chain proxies: 3.68, per 1000 atoms: 0.56 Number of scatterers: 6600 At special positions: 0 Unit cell: (76.56, 157.52, 58.96, 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 36 16.00 O 1226 8.00 N 1130 7.00 C 4208 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=8, symmetry=0 Simple disulfide: pdb=" SG CYS A 57 " - pdb=" SG CYS A 70 " distance=2.03 Simple disulfide: pdb=" SG CYS A 294 " - pdb=" SG CYS A 305 " distance=2.03 Simple disulfide: pdb=" SG CYS A 308 " - pdb=" SG CYS A 313 " distance=2.03 Simple disulfide: pdb=" SG CYS A 448 " - pdb=" SG CYS A 468 " distance=2.03 Simple disulfide: pdb=" SG CYS B 57 " - pdb=" SG CYS B 70 " distance=2.03 Simple disulfide: pdb=" SG CYS B 294 " - pdb=" SG CYS B 305 " distance=2.03 Simple disulfide: pdb=" SG CYS B 308 " - pdb=" SG CYS B 313 " distance=2.03 Simple disulfide: pdb=" SG CYS B 448 " - pdb=" SG CYS B 468 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 2.37 Conformation dependent library (CDL) restraints added in 1.1 seconds 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1548 Finding SS restraints... Secondary structure from input PDB file: 20 helices and 6 sheets defined 23.7% alpha, 17.8% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.38 Creating SS restraints... Processing helix chain 'A' and resid 94 through 106 Processing helix chain 'A' and resid 116 through 120 removed outlier: 3.826A pdb=" N ALA A 120 " --> pdb=" O LEU A 117 " (cutoff:3.500A) Processing helix chain 'A' and resid 123 through 149 removed outlier: 4.393A pdb=" N LYS A 132 " --> pdb=" O ALA A 128 " (cutoff:3.500A) removed outlier: 4.373A pdb=" N LEU A 133 " --> pdb=" O GLY A 129 " (cutoff:3.500A) removed outlier: 3.929A pdb=" N MET A 145 " --> pdb=" O PHE A 141 " (cutoff:3.500A) removed outlier: 3.874A pdb=" N ALA A 146 " --> pdb=" O MET A 142 " (cutoff:3.500A) removed outlier: 4.190A pdb=" N ASP A 147 " --> pdb=" O ASN A 143 " (cutoff:3.500A) Processing helix chain 'A' and resid 163 through 175 Processing helix chain 'A' and resid 196 through 200 Processing helix chain 'A' and resid 202 through 206 removed outlier: 3.762A pdb=" N ALA A 206 " --> pdb=" O PRO A 203 " (cutoff:3.500A) Processing helix chain 'A' and resid 239 through 241 No H-bonds generated for 'chain 'A' and resid 239 through 241' Processing helix chain 'A' and resid 270 through 284 Processing helix chain 'A' and resid 296 through 303 Processing helix chain 'A' and resid 451 through 455 removed outlier: 3.964A pdb=" N SER A 454 " --> pdb=" O GLU A 451 " (cutoff:3.500A) removed outlier: 3.667A pdb=" N TYR A 455 " --> pdb=" O LYS A 452 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 451 through 455' Processing helix chain 'B' and resid 94 through 106 Processing helix chain 'B' and resid 116 through 120 removed outlier: 3.826A pdb=" N ALA B 120 " --> pdb=" O LEU B 117 " (cutoff:3.500A) Processing helix chain 'B' and resid 123 through 149 removed outlier: 4.393A pdb=" N LYS B 132 " --> pdb=" O ALA B 128 " (cutoff:3.500A) removed outlier: 4.372A pdb=" N LEU B 133 " --> pdb=" O GLY B 129 " (cutoff:3.500A) removed outlier: 3.929A pdb=" N MET B 145 " --> pdb=" O PHE B 141 " (cutoff:3.500A) removed outlier: 3.874A pdb=" N ALA B 146 " --> pdb=" O MET B 142 " (cutoff:3.500A) removed outlier: 4.190A pdb=" N ASP B 147 " --> pdb=" O ASN B 143 " (cutoff:3.500A) Processing helix chain 'B' and resid 163 through 175 Processing helix chain 'B' and resid 196 through 200 Processing helix chain 'B' and resid 202 through 206 removed outlier: 3.761A pdb=" N ALA B 206 " --> pdb=" O PRO B 203 " (cutoff:3.500A) Processing helix chain 'B' and resid 239 through 241 No H-bonds generated for 'chain 'B' and resid 239 through 241' Processing helix chain 'B' and resid 271 through 284 removed outlier: 3.696A pdb=" N PHE B 278 " --> pdb=" O SER B 274 " (cutoff:3.500A) Processing helix chain 'B' and resid 296 through 303 Processing helix chain 'B' and resid 451 through 455 removed outlier: 3.964A pdb=" N SER B 454 " --> pdb=" O GLU B 451 " (cutoff:3.500A) removed outlier: 3.666A pdb=" N TYR B 455 " --> pdb=" O LYS B 452 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 451 through 455' Processing sheet with id=AA1, first strand: chain 'A' and resid 44 through 47 removed outlier: 4.704A pdb=" N VAL A 111 " --> pdb=" O ARG A 47 " (cutoff:3.500A) removed outlier: 6.591A pdb=" N LYS A 77 " --> pdb=" O ASN A 110 " (cutoff:3.500A) removed outlier: 7.750A pdb=" N ILE A 112 " --> pdb=" O LYS A 77 " (cutoff:3.500A) removed outlier: 7.165A pdb=" N PHE A 79 " --> pdb=" O ILE A 112 " (cutoff:3.500A) removed outlier: 9.158A pdb=" N VAL A 114 " --> pdb=" O PHE A 79 " (cutoff:3.500A) removed outlier: 8.072A pdb=" N VAL A 81 " --> pdb=" O VAL A 114 " (cutoff:3.500A) removed outlier: 4.260A pdb=" N TYR A 161 " --> pdb=" O ILE A 82 " (cutoff:3.500A) removed outlier: 8.342A pdb=" N THR A 183 " --> pdb=" O VAL A 156 " (cutoff:3.500A) removed outlier: 6.945A pdb=" N LEU A 158 " --> pdb=" O THR A 183 " (cutoff:3.500A) removed outlier: 7.483A pdb=" N LEU A 185 " --> pdb=" O LEU A 158 " (cutoff:3.500A) removed outlier: 6.930A pdb=" N GLY A 160 " --> pdb=" O LEU A 185 " (cutoff:3.500A) removed outlier: 5.039A pdb=" N SER A 162 " --> pdb=" O PRO A 187 " (cutoff:3.500A) removed outlier: 3.619A pdb=" N LEU A 212 " --> pdb=" O GLY A 184 " (cutoff:3.500A) removed outlier: 3.689A pdb=" N ASP A 186 " --> pdb=" O LEU A 212 " (cutoff:3.500A) removed outlier: 6.892A pdb=" N VAL A 209 " --> pdb=" O VAL A 233 " (cutoff:3.500A) removed outlier: 8.663A pdb=" N ILE A 235 " --> pdb=" O VAL A 209 " (cutoff:3.500A) removed outlier: 7.396A pdb=" N VAL A 211 " --> pdb=" O ILE A 235 " (cutoff:3.500A) removed outlier: 5.610A pdb=" N ASP A 234 " --> pdb=" O MET A 331 " (cutoff:3.500A) removed outlier: 7.269A pdb=" N LEU A 333 " --> pdb=" O ASP A 234 " (cutoff:3.500A) removed outlier: 5.997A pdb=" N TYR A 236 " --> pdb=" O LEU A 333 " (cutoff:3.500A) removed outlier: 4.100A pdb=" N LYS A 330 " --> pdb=" O CYS A 294 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 395 through 398 removed outlier: 4.239A pdb=" N VAL A 436 " --> pdb=" O VAL A 445 " (cutoff:3.500A) removed outlier: 4.451A pdb=" N VAL A 445 " --> pdb=" O VAL A 436 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 375 through 381 removed outlier: 5.731A pdb=" N GLU A 375 " --> pdb=" O GLY A 371 " (cutoff:3.500A) removed outlier: 7.332A pdb=" N GLY A 371 " --> pdb=" O GLU A 375 " (cutoff:3.500A) removed outlier: 5.578A pdb=" N GLU A 377 " --> pdb=" O LEU A 369 " (cutoff:3.500A) removed outlier: 5.270A pdb=" N LEU A 369 " --> pdb=" O GLU A 377 " (cutoff:3.500A) removed outlier: 3.613A pdb=" N ILE A 367 " --> pdb=" O ILE A 379 " (cutoff:3.500A) removed outlier: 4.041A pdb=" N PHE A 365 " --> pdb=" O PHE A 381 " (cutoff:3.500A) removed outlier: 3.856A pdb=" N TYR A 370 " --> pdb=" O MET A 407 " (cutoff:3.500A) removed outlier: 4.152A pdb=" N MET A 407 " --> pdb=" O TYR A 370 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'B' and resid 44 through 47 removed outlier: 4.705A pdb=" N VAL B 111 " --> pdb=" O ARG B 47 " (cutoff:3.500A) removed outlier: 6.591A pdb=" N LYS B 77 " --> pdb=" O ASN B 110 " (cutoff:3.500A) removed outlier: 7.751A pdb=" N ILE B 112 " --> pdb=" O LYS B 77 " (cutoff:3.500A) removed outlier: 7.165A pdb=" N PHE B 79 " --> pdb=" O ILE B 112 " (cutoff:3.500A) removed outlier: 9.156A pdb=" N VAL B 114 " --> pdb=" O PHE B 79 " (cutoff:3.500A) removed outlier: 8.073A pdb=" N VAL B 81 " --> pdb=" O VAL B 114 " (cutoff:3.500A) removed outlier: 4.260A pdb=" N TYR B 161 " --> pdb=" O ILE B 82 " (cutoff:3.500A) removed outlier: 8.342A pdb=" N THR B 183 " --> pdb=" O VAL B 156 " (cutoff:3.500A) removed outlier: 6.945A pdb=" N LEU B 158 " --> pdb=" O THR B 183 " (cutoff:3.500A) removed outlier: 7.483A pdb=" N LEU B 185 " --> pdb=" O LEU B 158 " (cutoff:3.500A) removed outlier: 6.931A pdb=" N GLY B 160 " --> pdb=" O LEU B 185 " (cutoff:3.500A) removed outlier: 5.039A pdb=" N SER B 162 " --> pdb=" O PRO B 187 " (cutoff:3.500A) removed outlier: 3.619A pdb=" N LEU B 212 " --> pdb=" O GLY B 184 " (cutoff:3.500A) removed outlier: 3.688A pdb=" N ASP B 186 " --> pdb=" O LEU B 212 " (cutoff:3.500A) removed outlier: 6.893A pdb=" N VAL B 209 " --> pdb=" O VAL B 233 " (cutoff:3.500A) removed outlier: 8.663A pdb=" N ILE B 235 " --> pdb=" O VAL B 209 " (cutoff:3.500A) removed outlier: 7.396A pdb=" N VAL B 211 " --> pdb=" O ILE B 235 " (cutoff:3.500A) removed outlier: 5.611A pdb=" N ASP B 234 " --> pdb=" O MET B 331 " (cutoff:3.500A) removed outlier: 7.269A pdb=" N LEU B 333 " --> pdb=" O ASP B 234 " (cutoff:3.500A) removed outlier: 5.997A pdb=" N TYR B 236 " --> pdb=" O LEU B 333 " (cutoff:3.500A) removed outlier: 4.101A pdb=" N LYS B 330 " --> pdb=" O CYS B 294 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'B' and resid 395 through 398 removed outlier: 4.238A pdb=" N VAL B 436 " --> pdb=" O VAL B 445 " (cutoff:3.500A) removed outlier: 4.452A pdb=" N VAL B 445 " --> pdb=" O VAL B 436 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'B' and resid 375 through 381 removed outlier: 5.730A pdb=" N GLU B 375 " --> pdb=" O GLY B 371 " (cutoff:3.500A) removed outlier: 7.331A pdb=" N GLY B 371 " --> pdb=" O GLU B 375 " (cutoff:3.500A) removed outlier: 5.578A pdb=" N GLU B 377 " --> pdb=" O LEU B 369 " (cutoff:3.500A) removed outlier: 5.271A pdb=" N LEU B 369 " --> pdb=" O GLU B 377 " (cutoff:3.500A) removed outlier: 3.613A pdb=" N ILE B 367 " --> pdb=" O ILE B 379 " (cutoff:3.500A) removed outlier: 4.041A pdb=" N PHE B 365 " --> pdb=" O PHE B 381 " (cutoff:3.500A) removed outlier: 3.857A pdb=" N TYR B 370 " --> pdb=" O MET B 407 " (cutoff:3.500A) removed outlier: 4.152A pdb=" N MET B 407 " --> pdb=" O TYR B 370 " (cutoff:3.500A) 172 hydrogen bonds defined for protein. 468 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: 1.68 Time building geometry restraints manager: 2.60 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.34: 2071 1.34 - 1.46: 1655 1.46 - 1.59: 2998 1.59 - 1.71: 0 1.71 - 1.83: 52 Bond restraints: 6776 Sorted by residual: bond pdb=" C VAL A 95 " pdb=" N PRO A 96 " ideal model delta sigma weight residual 1.336 1.361 -0.025 1.20e-02 6.94e+03 4.38e+00 bond pdb=" C VAL B 95 " pdb=" N PRO B 96 " ideal model delta sigma weight residual 1.336 1.361 -0.025 1.20e-02 6.94e+03 4.22e+00 bond pdb=" CB TRP A 412 " pdb=" CG TRP A 412 " ideal model delta sigma weight residual 1.498 1.551 -0.053 3.10e-02 1.04e+03 2.89e+00 bond pdb=" CB TRP B 412 " pdb=" CG TRP B 412 " ideal model delta sigma weight residual 1.498 1.550 -0.052 3.10e-02 1.04e+03 2.79e+00 bond pdb=" CB PHE B 79 " pdb=" CG PHE B 79 " ideal model delta sigma weight residual 1.502 1.539 -0.037 2.30e-02 1.89e+03 2.59e+00 ... (remaining 6771 not shown) Histogram of bond angle deviations from ideal: 95.02 - 103.13: 52 103.13 - 111.25: 2472 111.25 - 119.37: 3080 119.37 - 127.48: 3460 127.48 - 135.60: 112 Bond angle restraints: 9176 Sorted by residual: angle pdb=" CA LEU B 102 " pdb=" CB LEU B 102 " pdb=" CG LEU B 102 " ideal model delta sigma weight residual 116.30 134.78 -18.48 3.50e+00 8.16e-02 2.79e+01 angle pdb=" CA LEU A 102 " pdb=" CB LEU A 102 " pdb=" CG LEU A 102 " ideal model delta sigma weight residual 116.30 134.76 -18.46 3.50e+00 8.16e-02 2.78e+01 angle pdb=" CG1 ILE A 320 " pdb=" CB ILE A 320 " pdb=" CG2 ILE A 320 " ideal model delta sigma weight residual 110.70 95.02 15.68 3.00e+00 1.11e-01 2.73e+01 angle pdb=" CG1 ILE B 320 " pdb=" CB ILE B 320 " pdb=" CG2 ILE B 320 " ideal model delta sigma weight residual 110.70 95.04 15.66 3.00e+00 1.11e-01 2.72e+01 angle pdb=" CA LEU B 277 " pdb=" CB LEU B 277 " pdb=" CG LEU B 277 " ideal model delta sigma weight residual 116.30 133.94 -17.64 3.50e+00 8.16e-02 2.54e+01 ... (remaining 9171 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.79: 3541 17.79 - 35.57: 392 35.57 - 53.35: 84 53.35 - 71.14: 12 71.14 - 88.92: 9 Dihedral angle restraints: 4038 sinusoidal: 1626 harmonic: 2412 Sorted by residual: dihedral pdb=" CB CYS B 448 " pdb=" SG CYS B 448 " pdb=" SG CYS B 468 " pdb=" CB CYS B 468 " ideal model delta sinusoidal sigma weight residual -86.00 -172.57 86.57 1 1.00e+01 1.00e-02 9.03e+01 dihedral pdb=" CB CYS A 448 " pdb=" SG CYS A 448 " pdb=" SG CYS A 468 " pdb=" CB CYS A 468 " ideal model delta sinusoidal sigma weight residual -86.00 -172.53 86.53 1 1.00e+01 1.00e-02 9.02e+01 dihedral pdb=" CB CYS B 294 " pdb=" SG CYS B 294 " pdb=" SG CYS B 305 " pdb=" CB CYS B 305 " ideal model delta sinusoidal sigma weight residual 93.00 41.28 51.72 1 1.00e+01 1.00e-02 3.65e+01 ... (remaining 4035 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.062: 715 0.062 - 0.124: 196 0.124 - 0.186: 57 0.186 - 0.248: 8 0.248 - 0.311: 4 Chirality restraints: 980 Sorted by residual: chirality pdb=" CG LEU B 333 " pdb=" CB LEU B 333 " pdb=" CD1 LEU B 333 " pdb=" CD2 LEU B 333 " both_signs ideal model delta sigma weight residual False -2.59 -2.90 0.31 2.00e-01 2.50e+01 2.41e+00 chirality pdb=" CG LEU A 333 " pdb=" CB LEU A 333 " pdb=" CD1 LEU A 333 " pdb=" CD2 LEU A 333 " both_signs ideal model delta sigma weight residual False -2.59 -2.90 0.31 2.00e-01 2.50e+01 2.41e+00 chirality pdb=" CB ILE B 320 " pdb=" CA ILE B 320 " pdb=" CG1 ILE B 320 " pdb=" CG2 ILE B 320 " both_signs ideal model delta sigma weight residual False 2.64 2.93 -0.28 2.00e-01 2.50e+01 2.02e+00 ... (remaining 977 not shown) Planarity restraints: 1174 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA ARG B 273 " 0.014 2.00e-02 2.50e+03 2.87e-02 8.23e+00 pdb=" C ARG B 273 " -0.050 2.00e-02 2.50e+03 pdb=" O ARG B 273 " 0.019 2.00e-02 2.50e+03 pdb=" N SER B 274 " 0.017 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE B 79 " 0.016 2.00e-02 2.50e+03 2.02e-02 7.17e+00 pdb=" CG PHE B 79 " -0.047 2.00e-02 2.50e+03 pdb=" CD1 PHE B 79 " 0.013 2.00e-02 2.50e+03 pdb=" CD2 PHE B 79 " 0.015 2.00e-02 2.50e+03 pdb=" CE1 PHE B 79 " 0.002 2.00e-02 2.50e+03 pdb=" CE2 PHE B 79 " -0.001 2.00e-02 2.50e+03 pdb=" CZ PHE B 79 " 0.001 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE A 79 " -0.016 2.00e-02 2.50e+03 2.02e-02 7.12e+00 pdb=" CG PHE A 79 " 0.047 2.00e-02 2.50e+03 pdb=" CD1 PHE A 79 " -0.013 2.00e-02 2.50e+03 pdb=" CD2 PHE A 79 " -0.015 2.00e-02 2.50e+03 pdb=" CE1 PHE A 79 " -0.002 2.00e-02 2.50e+03 pdb=" CE2 PHE A 79 " 0.001 2.00e-02 2.50e+03 pdb=" CZ PHE A 79 " -0.001 2.00e-02 2.50e+03 ... (remaining 1171 not shown) Histogram of nonbonded interaction distances: 2.30 - 2.82: 1885 2.82 - 3.34: 5895 3.34 - 3.86: 11272 3.86 - 4.38: 12283 4.38 - 4.90: 20881 Nonbonded interactions: 52216 Sorted by model distance: nonbonded pdb=" NE2 GLN B 457 " pdb=" O GLY B 459 " model vdw 2.296 2.520 nonbonded pdb=" NE2 GLN A 457 " pdb=" O GLY A 459 " model vdw 2.296 2.520 nonbonded pdb=" O ALA B 128 " pdb=" OG1 THR B 131 " model vdw 2.299 2.440 nonbonded pdb=" O ALA A 128 " pdb=" OG1 THR A 131 " model vdw 2.299 2.440 nonbonded pdb=" O ASN B 73 " pdb=" OG SER B 76 " model vdw 2.337 2.440 ... (remaining 52211 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = chain 'A' selection = chain 'B' } 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 group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 1.330 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: 7.130 Check model and map are aligned: 0.100 Set scattering table: 0.070 Process input model: 19.830 Find NCS groups from input model: 0.310 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.000 Load rotamer database and sin/cos tables:3.050 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 31.860 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7750 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.005 0.053 6776 Z= 0.358 Angle : 1.277 18.484 9176 Z= 0.623 Chirality : 0.065 0.311 980 Planarity : 0.006 0.042 1174 Dihedral : 15.420 88.923 2466 Min Nonbonded Distance : 2.296 Molprobity Statistics. All-atom Clashscore : 4.84 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.49 % Favored : 90.51 % Rotamer Outliers : 0.27 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.82 (0.24), residues: 822 helix: -2.63 (0.29), residues: 156 sheet: -1.66 (0.35), residues: 206 loop : -3.01 (0.23), residues: 460 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 92 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 90 time to evaluate : 0.853 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 2 outliers final: 1 residues processed: 92 average time/residue: 0.1820 time to fit residues: 23.0619 Evaluate side-chains 86 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 85 time to evaluate : 0.797 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.3362 time to fit residues: 1.4940 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 82 random chunks: chunk 69 optimal weight: 1.9990 chunk 62 optimal weight: 0.9980 chunk 34 optimal weight: 0.0970 chunk 21 optimal weight: 5.9990 chunk 41 optimal weight: 10.0000 chunk 33 optimal weight: 7.9990 chunk 64 optimal weight: 0.5980 chunk 24 optimal weight: 2.9990 chunk 39 optimal weight: 0.6980 chunk 47 optimal weight: 2.9990 chunk 74 optimal weight: 0.9980 overall best weight: 0.6778 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 191 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 191 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 271 HIS Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7801 moved from start: 0.1371 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.023 6776 Z= 0.152 Angle : 0.672 12.225 9176 Z= 0.328 Chirality : 0.044 0.194 980 Planarity : 0.004 0.067 1174 Dihedral : 4.905 20.234 898 Min Nonbonded Distance : 2.146 Molprobity Statistics. All-atom Clashscore : 10.44 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.69 % Favored : 93.31 % Rotamer Outliers : 2.47 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.14 (0.26), residues: 822 helix: -1.70 (0.33), residues: 168 sheet: -1.52 (0.36), residues: 206 loop : -2.56 (0.25), residues: 448 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 94 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 76 time to evaluate : 0.791 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 18 outliers final: 8 residues processed: 90 average time/residue: 0.1508 time to fit residues: 19.5677 Evaluate side-chains 78 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 70 time to evaluate : 0.859 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 0 residues processed: 8 average time/residue: 0.0688 time to fit residues: 2.0322 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 82 random chunks: chunk 41 optimal weight: 10.0000 chunk 23 optimal weight: 4.9990 chunk 61 optimal weight: 2.9990 chunk 50 optimal weight: 0.9990 chunk 20 optimal weight: 4.9990 chunk 74 optimal weight: 0.1980 chunk 80 optimal weight: 0.8980 chunk 66 optimal weight: 6.9990 chunk 73 optimal weight: 0.9990 chunk 25 optimal weight: 0.4980 chunk 59 optimal weight: 0.7980 overall best weight: 0.6782 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** B 191 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7814 moved from start: 0.1693 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.022 6776 Z= 0.141 Angle : 0.634 12.146 9176 Z= 0.307 Chirality : 0.044 0.155 980 Planarity : 0.003 0.044 1174 Dihedral : 4.632 16.263 898 Min Nonbonded Distance : 2.106 Molprobity Statistics. All-atom Clashscore : 11.67 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.81 % Favored : 93.19 % Rotamer Outliers : 1.92 % 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.77 (0.26), residues: 822 helix: -1.08 (0.36), residues: 164 sheet: -1.21 (0.37), residues: 186 loop : -2.48 (0.25), residues: 472 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 84 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 70 time to evaluate : 0.752 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 14 outliers final: 10 residues processed: 81 average time/residue: 0.1357 time to fit residues: 16.2840 Evaluate side-chains 80 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 10 poor density : 70 time to evaluate : 0.747 Switching outliers to nearest non-outliers outliers start: 10 outliers final: 0 residues processed: 10 average time/residue: 0.0669 time to fit residues: 2.1615 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 82 random chunks: chunk 73 optimal weight: 2.9990 chunk 56 optimal weight: 9.9990 chunk 38 optimal weight: 1.9990 chunk 8 optimal weight: 10.0000 chunk 35 optimal weight: 6.9990 chunk 50 optimal weight: 0.8980 chunk 74 optimal weight: 0.8980 chunk 79 optimal weight: 1.9990 chunk 39 optimal weight: 9.9990 chunk 70 optimal weight: 1.9990 chunk 21 optimal weight: 0.9980 overall best weight: 1.3584 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 58 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7880 moved from start: 0.2504 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.026 6776 Z= 0.191 Angle : 0.645 12.236 9176 Z= 0.316 Chirality : 0.044 0.145 980 Planarity : 0.003 0.031 1174 Dihedral : 4.667 17.585 898 Min Nonbonded Distance : 2.070 Molprobity Statistics. All-atom Clashscore : 17.04 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.91 % Favored : 92.09 % Rotamer Outliers : 2.47 % 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.67 (0.27), residues: 822 helix: -0.63 (0.38), residues: 164 sheet: -1.48 (0.35), residues: 208 loop : -2.42 (0.25), residues: 450 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 85 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 67 time to evaluate : 0.805 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 18 outliers final: 13 residues processed: 76 average time/residue: 0.1400 time to fit residues: 15.9336 Evaluate side-chains 75 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 62 time to evaluate : 0.748 Switching outliers to nearest non-outliers outliers start: 13 outliers final: 0 residues processed: 13 average time/residue: 0.0750 time to fit residues: 2.7101 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 82 random chunks: chunk 66 optimal weight: 5.9990 chunk 44 optimal weight: 5.9990 chunk 1 optimal weight: 6.9990 chunk 59 optimal weight: 9.9990 chunk 32 optimal weight: 0.0000 chunk 67 optimal weight: 7.9990 chunk 54 optimal weight: 5.9990 chunk 0 optimal weight: 9.9990 chunk 40 optimal weight: 3.9990 chunk 71 optimal weight: 4.9990 chunk 20 optimal weight: 4.9990 overall best weight: 3.9992 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 271 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 469 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 469 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8002 moved from start: 0.4275 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.007 0.045 6776 Z= 0.449 Angle : 0.849 9.834 9176 Z= 0.445 Chirality : 0.048 0.203 980 Planarity : 0.005 0.036 1174 Dihedral : 5.982 20.388 898 Min Nonbonded Distance : 2.016 Molprobity Statistics. All-atom Clashscore : 37.38 Ramachandran Plot: Outliers : 0.00 % Allowed : 11.80 % Favored : 88.20 % Rotamer Outliers : 4.67 % 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.81 (0.27), residues: 822 helix: -0.35 (0.43), residues: 152 sheet: -1.68 (0.38), residues: 188 loop : -2.56 (0.25), residues: 482 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 104 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 34 poor density : 70 time to evaluate : 0.812 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 34 outliers final: 19 residues processed: 99 average time/residue: 0.1447 time to fit residues: 20.6282 Evaluate side-chains 85 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 19 poor density : 66 time to evaluate : 0.786 Switching outliers to nearest non-outliers outliers start: 19 outliers final: 0 residues processed: 19 average time/residue: 0.0709 time to fit residues: 3.3676 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 82 random chunks: chunk 26 optimal weight: 0.0050 chunk 71 optimal weight: 6.9990 chunk 15 optimal weight: 0.8980 chunk 46 optimal weight: 5.9990 chunk 19 optimal weight: 2.9990 chunk 79 optimal weight: 2.9990 chunk 65 optimal weight: 1.9990 chunk 36 optimal weight: 7.9990 chunk 6 optimal weight: 0.6980 chunk 41 optimal weight: 10.0000 chunk 76 optimal weight: 0.9990 overall best weight: 0.9198 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 271 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 469 HIS ** B 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 469 HIS Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7918 moved from start: 0.3873 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.028 6776 Z= 0.179 Angle : 0.711 13.901 9176 Z= 0.341 Chirality : 0.045 0.130 980 Planarity : 0.003 0.032 1174 Dihedral : 5.122 20.320 898 Min Nonbonded Distance : 2.086 Molprobity Statistics. All-atom Clashscore : 20.03 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.00 % Favored : 91.00 % Rotamer Outliers : 2.06 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.81 (0.27), residues: 822 helix: -0.33 (0.42), residues: 152 sheet: -2.19 (0.39), residues: 170 loop : -2.33 (0.25), residues: 500 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 90 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 75 time to evaluate : 0.845 Fit side-chains revert: symmetry clash outliers start: 15 outliers final: 8 residues processed: 85 average time/residue: 0.1474 time to fit residues: 18.2629 Evaluate side-chains 76 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 68 time to evaluate : 0.802 Switching outliers to nearest non-outliers outliers start: 8 outliers final: 0 residues processed: 8 average time/residue: 0.0736 time to fit residues: 2.0735 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 82 random chunks: chunk 8 optimal weight: 10.0000 chunk 45 optimal weight: 4.9990 chunk 57 optimal weight: 4.9990 chunk 44 optimal weight: 2.9990 chunk 66 optimal weight: 4.9990 chunk 79 optimal weight: 6.9990 chunk 49 optimal weight: 5.9990 chunk 48 optimal weight: 4.9990 chunk 36 optimal weight: 7.9990 chunk 31 optimal weight: 0.9980 chunk 47 optimal weight: 1.9990 overall best weight: 3.1988 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 271 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7992 moved from start: 0.4582 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.006 0.036 6776 Z= 0.369 Angle : 0.793 10.794 9176 Z= 0.402 Chirality : 0.048 0.153 980 Planarity : 0.004 0.038 1174 Dihedral : 5.679 19.799 898 Min Nonbonded Distance : 2.044 Molprobity Statistics. All-atom Clashscore : 34.39 Ramachandran Plot: Outliers : 0.00 % Allowed : 12.53 % Favored : 87.47 % Rotamer Outliers : 2.75 % 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.88 (0.28), residues: 822 helix: -0.38 (0.43), residues: 152 sheet: -1.90 (0.42), residues: 154 loop : -2.53 (0.25), residues: 516 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 85 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 20 poor density : 65 time to evaluate : 0.801 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 20 outliers final: 15 residues processed: 82 average time/residue: 0.1520 time to fit residues: 18.2308 Evaluate side-chains 76 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 61 time to evaluate : 0.829 Switching outliers to nearest non-outliers outliers start: 15 outliers final: 0 residues processed: 15 average time/residue: 0.0759 time to fit residues: 3.0343 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 82 random chunks: chunk 23 optimal weight: 1.9990 chunk 15 optimal weight: 4.9990 chunk 50 optimal weight: 0.6980 chunk 53 optimal weight: 0.9980 chunk 39 optimal weight: 0.0870 chunk 7 optimal weight: 4.9990 chunk 62 optimal weight: 4.9990 chunk 71 optimal weight: 8.9990 chunk 75 optimal weight: 3.9990 chunk 69 optimal weight: 0.7980 chunk 73 optimal weight: 0.4980 overall best weight: 0.6158 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 271 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7907 moved from start: 0.4126 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.025 6776 Z= 0.168 Angle : 0.731 14.962 9176 Z= 0.344 Chirality : 0.045 0.140 980 Planarity : 0.003 0.032 1174 Dihedral : 4.990 19.900 898 Min Nonbonded Distance : 2.104 Molprobity Statistics. All-atom Clashscore : 20.03 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.88 % Favored : 91.12 % Rotamer Outliers : 0.82 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.69 (0.28), residues: 822 helix: -0.08 (0.43), residues: 152 sheet: -2.14 (0.36), residues: 192 loop : -2.27 (0.26), residues: 478 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 83 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 77 time to evaluate : 0.827 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 6 outliers final: 4 residues processed: 78 average time/residue: 0.1447 time to fit residues: 16.7405 Evaluate side-chains 72 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 68 time to evaluate : 0.817 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.0653 time to fit residues: 1.5532 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 82 random chunks: chunk 75 optimal weight: 0.7980 chunk 44 optimal weight: 0.4980 chunk 32 optimal weight: 0.6980 chunk 57 optimal weight: 3.9990 chunk 22 optimal weight: 0.7980 chunk 66 optimal weight: 5.9990 chunk 69 optimal weight: 3.9990 chunk 73 optimal weight: 0.9990 chunk 48 optimal weight: 0.8980 chunk 77 optimal weight: 4.9990 chunk 47 optimal weight: 0.8980 overall best weight: 0.7380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 271 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7908 moved from start: 0.4086 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.029 6776 Z= 0.171 Angle : 0.742 15.559 9176 Z= 0.345 Chirality : 0.045 0.133 980 Planarity : 0.003 0.030 1174 Dihedral : 4.809 19.289 898 Min Nonbonded Distance : 2.102 Molprobity Statistics. All-atom Clashscore : 20.49 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.12 % Favored : 90.88 % Rotamer Outliers : 0.55 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.65 (0.28), residues: 822 helix: -0.02 (0.42), residues: 152 sheet: -2.09 (0.37), residues: 178 loop : -2.28 (0.26), residues: 492 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 74 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 70 time to evaluate : 0.816 Fit side-chains revert: symmetry clash outliers start: 4 outliers final: 3 residues processed: 70 average time/residue: 0.1464 time to fit residues: 15.2703 Evaluate side-chains 71 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 68 time to evaluate : 0.838 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0992 time to fit residues: 1.5520 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 82 random chunks: chunk 36 optimal weight: 9.9990 chunk 54 optimal weight: 10.0000 chunk 81 optimal weight: 3.9990 chunk 75 optimal weight: 2.9990 chunk 65 optimal weight: 3.9990 chunk 6 optimal weight: 3.9990 chunk 50 optimal weight: 0.8980 chunk 39 optimal weight: 2.9990 chunk 51 optimal weight: 2.9990 chunk 69 optimal weight: 0.8980 chunk 19 optimal weight: 0.7980 overall best weight: 1.7184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 271 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7935 moved from start: 0.4396 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.033 6776 Z= 0.228 Angle : 0.742 13.991 9176 Z= 0.357 Chirality : 0.044 0.134 980 Planarity : 0.003 0.032 1174 Dihedral : 4.943 18.513 898 Min Nonbonded Distance : 2.094 Molprobity Statistics. All-atom Clashscore : 23.79 Ramachandran Plot: Outliers : 0.00 % Allowed : 10.58 % Favored : 89.42 % Rotamer Outliers : 0.55 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.67 (0.28), residues: 822 helix: -0.16 (0.43), residues: 152 sheet: -1.99 (0.39), residues: 158 loop : -2.31 (0.26), residues: 512 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1644 Ramachandran restraints generated. 822 Oldfield, 0 Emsley, 822 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 70 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 66 time to evaluate : 0.812 Fit side-chains revert: symmetry clash outliers start: 4 outliers final: 4 residues processed: 67 average time/residue: 0.1429 time to fit residues: 14.5114 Evaluate side-chains 69 residues out of total 728 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 65 time to evaluate : 0.817 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.0766 time to fit residues: 1.6012 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 82 random chunks: chunk 59 optimal weight: 9.9990 chunk 9 optimal weight: 6.9990 chunk 18 optimal weight: 6.9990 chunk 65 optimal weight: 0.0570 chunk 27 optimal weight: 4.9990 chunk 66 optimal weight: 0.1980 chunk 8 optimal weight: 1.9990 chunk 12 optimal weight: 0.9980 chunk 57 optimal weight: 0.8980 chunk 3 optimal weight: 3.9990 chunk 47 optimal weight: 2.9990 overall best weight: 0.8300 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 363 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3791 r_free = 0.3791 target = 0.067593 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 29)----------------| | r_work = 0.3355 r_free = 0.3355 target = 0.050331 restraints weight = 31196.658| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 32)----------------| | r_work = 0.3412 r_free = 0.3412 target = 0.052436 restraints weight = 19989.870| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 24)----------------| | r_work = 0.3452 r_free = 0.3452 target = 0.053922 restraints weight = 14776.772| |-----------------------------------------------------------------------------| r_work (final): 0.3421 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7974 moved from start: 0.4303 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.028 6776 Z= 0.172 Angle : 0.742 15.514 9176 Z= 0.348 Chirality : 0.045 0.133 980 Planarity : 0.003 0.033 1174 Dihedral : 4.767 18.821 898 Min Nonbonded Distance : 2.101 Molprobity Statistics. All-atom Clashscore : 21.34 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.49 % Favored : 90.51 % Rotamer Outliers : 0.14 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.63 (0.28), residues: 822 helix: -0.12 (0.43), residues: 152 sheet: -2.15 (0.38), residues: 178 loop : -2.20 (0.26), residues: 492 =============================================================================== Job complete usr+sys time: 1478.98 seconds wall clock time: 27 minutes 57.02 seconds (1677.02 seconds total)