Starting phenix.real_space_refine on Mon Jun 9 15:51:04 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/9kbd_62222/06_2025/9kbd_62222.cif Found real_map, /net/cci-nas-00/data/ceres_data/9kbd_62222/06_2025/9kbd_62222.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.7 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/9kbd_62222/06_2025/9kbd_62222.map" default_real_map = "/net/cci-nas-00/data/ceres_data/9kbd_62222/06_2025/9kbd_62222.map" model { file = "/net/cci-nas-00/data/ceres_data/9kbd_62222/06_2025/9kbd_62222.cif" } default_model = "/net/cci-nas-00/data/ceres_data/9kbd_62222/06_2025/9kbd_62222.cif" } resolution = 3.7 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.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 56 5.16 5 C 6824 2.51 5 N 1813 2.21 5 O 2035 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 5 residue(s): 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/chem_data/mon_lib" Total number of atoms: 10728 Number of models: 1 Model: "" Number of chains: 4 Chain: "B" Number of atoms: 1281 Number of conformers: 1 Conformer: "" Number of residues, atoms: 160, 1281 Classifications: {'peptide': 160} Link IDs: {'PTRANS': 9, 'TRANS': 150} Chain: "C" Number of atoms: 3277 Number of conformers: 1 Conformer: "" Number of residues, atoms: 405, 3277 Classifications: {'peptide': 405} Link IDs: {'PCIS': 2, 'PTRANS': 19, 'TRANS': 383} Chain: "R" Number of atoms: 144 Number of conformers: 1 Conformer: "" Number of residues, atoms: 16, 144 Classifications: {'peptide': 16} Link IDs: {'TRANS': 15} Chain: "A" Number of atoms: 6026 Number of conformers: 1 Conformer: "" Number of residues, atoms: 737, 6026 Classifications: {'peptide': 737} Link IDs: {'PTRANS': 13, 'TRANS': 723} Chain breaks: 1 Time building chain proxies: 6.97, per 1000 atoms: 0.65 Number of scatterers: 10728 At special positions: 0 Unit cell: (82, 127.1, 179.58, 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 56 16.00 O 2035 8.00 N 1813 7.00 C 6824 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.00 Amino acid : True - 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.31 Conformation dependent library (CDL) restraints added in 1.4 seconds 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2528 Finding SS restraints... Secondary structure from input PDB file: 57 helices and 14 sheets defined 58.0% alpha, 9.3% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.10 Creating SS restraints... Processing helix chain 'B' and resid 18 through 22 Processing helix chain 'B' and resid 24 through 34 removed outlier: 3.600A pdb=" N LEU B 34 " --> pdb=" O MET B 30 " (cutoff:3.500A) Processing helix chain 'B' and resid 51 through 65 Processing helix chain 'B' and resid 81 through 85 removed outlier: 3.607A pdb=" N ASP B 84 " --> pdb=" O ARG B 81 " (cutoff:3.500A) removed outlier: 4.039A pdb=" N ILE B 85 " --> pdb=" O THR B 82 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 81 through 85' Processing helix chain 'B' and resid 86 through 93 Processing helix chain 'B' and resid 96 through 110 Processing helix chain 'B' and resid 112 through 126 removed outlier: 3.536A pdb=" N MET B 126 " --> pdb=" O THR B 122 " (cutoff:3.500A) Processing helix chain 'B' and resid 131 through 140 Processing helix chain 'B' and resid 146 through 158 removed outlier: 4.306A pdb=" N GLU B 156 " --> pdb=" O GLN B 152 " (cutoff:3.500A) removed outlier: 3.574A pdb=" N GLN B 158 " --> pdb=" O ARG B 154 " (cutoff:3.500A) Processing helix chain 'B' and resid 159 through 161 No H-bonds generated for 'chain 'B' and resid 159 through 161' Processing helix chain 'C' and resid 17 through 27 removed outlier: 3.603A pdb=" N PHE C 27 " --> pdb=" O LEU C 23 " (cutoff:3.500A) Processing helix chain 'C' and resid 29 through 36 Processing helix chain 'C' and resid 37 through 39 No H-bonds generated for 'chain 'C' and resid 37 through 39' Processing helix chain 'C' and resid 40 through 48 Processing helix chain 'C' and resid 50 through 54 Processing helix chain 'C' and resid 64 through 72 Processing helix chain 'C' and resid 74 through 87 removed outlier: 3.579A pdb=" N ASP C 85 " --> pdb=" O ASP C 81 " (cutoff:3.500A) Processing helix chain 'C' and resid 89 through 111 removed outlier: 4.541A pdb=" N ALA C 94 " --> pdb=" O ILE C 90 " (cutoff:3.500A) removed outlier: 5.796A pdb=" N ALA C 95 " --> pdb=" O ASP C 91 " (cutoff:3.500A) removed outlier: 4.122A pdb=" N GLU C 108 " --> pdb=" O LYS C 104 " (cutoff:3.500A) Proline residue: C 109 - end of helix Processing helix chain 'C' and resid 111 through 118 removed outlier: 3.670A pdb=" N VAL C 115 " --> pdb=" O CYS C 111 " (cutoff:3.500A) Processing helix chain 'C' and resid 123 through 135 removed outlier: 3.568A pdb=" N LEU C 127 " --> pdb=" O ARG C 123 " (cutoff:3.500A) Processing helix chain 'C' and resid 139 through 147 removed outlier: 3.617A pdb=" N SER C 145 " --> pdb=" O ASP C 141 " (cutoff:3.500A) removed outlier: 3.670A pdb=" N TYR C 146 " --> pdb=" O TYR C 142 " (cutoff:3.500A) Processing helix chain 'C' and resid 176 through 183 Processing helix chain 'C' and resid 244 through 259 removed outlier: 3.544A pdb=" N GLY C 259 " --> pdb=" O ASN C 255 " (cutoff:3.500A) Processing helix chain 'C' and resid 294 through 297 Processing helix chain 'C' and resid 409 through 414 removed outlier: 3.705A pdb=" N LEU C 413 " --> pdb=" O SER C 409 " (cutoff:3.500A) removed outlier: 3.517A pdb=" N GLU C 414 " --> pdb=" O ILE C 410 " (cutoff:3.500A) No H-bonds generated for 'chain 'C' and resid 409 through 414' Processing helix chain 'A' and resid 20 through 33 Processing helix chain 'A' and resid 38 through 54 removed outlier: 3.529A pdb=" N TYR A 52 " --> pdb=" O HIS A 48 " (cutoff:3.500A) Processing helix chain 'A' and resid 85 through 105 removed outlier: 3.640A pdb=" N LYS A 105 " --> pdb=" O THR A 101 " (cutoff:3.500A) Processing helix chain 'A' and resid 112 through 137 removed outlier: 3.629A pdb=" N GLY A 135 " --> pdb=" O LYS A 131 " (cutoff:3.500A) Processing helix chain 'A' and resid 139 through 144 removed outlier: 3.667A pdb=" N TRP A 144 " --> pdb=" O LEU A 140 " (cutoff:3.500A) Processing helix chain 'A' and resid 144 through 151 removed outlier: 3.708A pdb=" N GLU A 151 " --> pdb=" O ARG A 147 " (cutoff:3.500A) Processing helix chain 'A' and resid 158 through 171 Processing helix chain 'A' and resid 175 through 192 removed outlier: 3.676A pdb=" N VAL A 179 " --> pdb=" O LEU A 175 " (cutoff:3.500A) Processing helix chain 'A' and resid 197 through 211 removed outlier: 3.947A pdb=" N ILE A 201 " --> pdb=" O ASN A 197 " (cutoff:3.500A) Processing helix chain 'A' and resid 227 through 232 Processing helix chain 'A' and resid 232 through 252 removed outlier: 3.932A pdb=" N GLU A 241 " --> pdb=" O LEU A 237 " (cutoff:3.500A) removed outlier: 4.072A pdb=" N ARG A 242 " --> pdb=" O ALA A 238 " (cutoff:3.500A) Processing helix chain 'A' and resid 256 through 274 removed outlier: 3.775A pdb=" N LYS A 262 " --> pdb=" O THR A 258 " (cutoff:3.500A) Processing helix chain 'A' and resid 282 through 296 removed outlier: 4.720A pdb=" N GLU A 286 " --> pdb=" O SER A 282 " (cutoff:3.500A) Processing helix chain 'A' and resid 299 through 314 removed outlier: 3.575A pdb=" N PHE A 303 " --> pdb=" O HIS A 299 " (cutoff:3.500A) removed outlier: 3.891A pdb=" N THR A 305 " --> pdb=" O GLU A 301 " (cutoff:3.500A) removed outlier: 3.776A pdb=" N GLN A 308 " --> pdb=" O HIS A 304 " (cutoff:3.500A) removed outlier: 4.527A pdb=" N ASN A 309 " --> pdb=" O THR A 305 " (cutoff:3.500A) Processing helix chain 'A' and resid 315 through 327 Processing helix chain 'A' and resid 333 through 356 removed outlier: 3.619A pdb=" N LYS A 338 " --> pdb=" O GLY A 334 " (cutoff:3.500A) Processing helix chain 'A' and resid 357 through 362 removed outlier: 4.287A pdb=" N ASN A 361 " --> pdb=" O GLU A 357 " (cutoff:3.500A) Processing helix chain 'A' and resid 362 through 383 removed outlier: 3.688A pdb=" N LEU A 380 " --> pdb=" O LYS A 376 " (cutoff:3.500A) removed outlier: 4.270A pdb=" N VAL A 381 " --> pdb=" O TYR A 377 " (cutoff:3.500A) Processing helix chain 'A' and resid 388 through 404 removed outlier: 4.045A pdb=" N ASN A 404 " --> pdb=" O GLY A 400 " (cutoff:3.500A) Processing helix chain 'A' and resid 406 through 412 Processing helix chain 'A' and resid 416 through 431 removed outlier: 3.546A pdb=" N LEU A 430 " --> pdb=" O CYS A 426 " (cutoff:3.500A) Processing helix chain 'A' and resid 438 through 454 Processing helix chain 'A' and resid 458 through 475 removed outlier: 3.873A pdb=" N PHE A 462 " --> pdb=" O ASP A 458 " (cutoff:3.500A) Processing helix chain 'A' and resid 481 through 496 Processing helix chain 'A' and resid 501 through 528 Processing helix chain 'A' and resid 556 through 574 Processing helix chain 'A' and resid 605 through 617 Processing helix chain 'A' and resid 622 through 630 Processing helix chain 'A' and resid 632 through 646 Processing helix chain 'A' and resid 688 through 723 removed outlier: 4.076A pdb=" N GLN A 692 " --> pdb=" O MET A 688 " (cutoff:3.500A) removed outlier: 4.518A pdb=" N LYS A 693 " --> pdb=" O LYS A 689 " (cutoff:3.500A) removed outlier: 3.816A pdb=" N GLN A 694 " --> pdb=" O THR A 690 " (cutoff:3.500A) Processing helix chain 'A' and resid 727 through 739 Processing helix chain 'A' and resid 745 through 758 Processing sheet with id=AA1, first strand: chain 'B' and resid 14 through 17 removed outlier: 8.954A pdb=" N VAL B 45 " --> pdb=" O SER B 3 " (cutoff:3.500A) removed outlier: 6.965A pdb=" N LYS B 5 " --> pdb=" O VAL B 45 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'C' and resid 162 through 163 Processing sheet with id=AA3, first strand: chain 'C' and resid 173 through 174 Processing sheet with id=AA4, first strand: chain 'C' and resid 193 through 195 Processing sheet with id=AA5, first strand: chain 'C' and resid 221 through 222 removed outlier: 7.128A pdb=" N VAL C 221 " --> pdb=" O ILE C 241 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'C' and resid 263 through 264 Processing sheet with id=AA7, first strand: chain 'C' and resid 284 through 285 removed outlier: 3.518A pdb=" N GLU C 313 " --> pdb=" O THR C 285 " (cutoff:3.500A) removed outlier: 3.521A pdb=" N TYR C 363 " --> pdb=" O ILE C 310 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'C' and resid 291 through 292 Processing sheet with id=AA9, first strand: chain 'C' and resid 340 through 341 removed outlier: 3.766A pdb=" N GLU C 377 " --> pdb=" O ARG C 332 " (cutoff:3.500A) removed outlier: 6.911A pdb=" N TYR C 379 " --> pdb=" O ARG C 329 " (cutoff:3.500A) removed outlier: 5.097A pdb=" N ARG C 329 " --> pdb=" O TYR C 379 " (cutoff:3.500A) removed outlier: 6.743A pdb=" N THR C 381 " --> pdb=" O ASP C 327 " (cutoff:3.500A) removed outlier: 3.540A pdb=" N CYS C 324 " --> pdb=" O ILE C 355 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'C' and resid 340 through 341 removed outlier: 3.766A pdb=" N GLU C 377 " --> pdb=" O ARG C 332 " (cutoff:3.500A) removed outlier: 3.985A pdb=" N GLY C 378 " --> pdb=" O ILE C 396 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'A' and resid 534 through 541 removed outlier: 3.819A pdb=" N ALA R 31 " --> pdb=" O GLN A 538 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'A' and resid 534 through 541 removed outlier: 3.819A pdb=" N ALA R 31 " --> pdb=" O GLN A 538 " (cutoff:3.500A) removed outlier: 3.577A pdb=" N GLU R 23 " --> pdb=" O VAL A 591 " (cutoff:3.500A) removed outlier: 7.046A pdb=" N GLU A 589 " --> pdb=" O LYS R 25 " (cutoff:3.500A) removed outlier: 4.934A pdb=" N TRP R 27 " --> pdb=" O LYS A 587 " (cutoff:3.500A) removed outlier: 7.203A pdb=" N LYS A 587 " --> pdb=" O TRP R 27 " (cutoff:3.500A) removed outlier: 3.734A pdb=" N LEU A 601 " --> pdb=" O LEU A 590 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'A' and resid 620 through 621 Processing sheet with id=AB5, first strand: chain 'A' and resid 725 through 726 removed outlier: 4.208A pdb=" N LEU A 725 " --> pdb=" O TYR A 772 " (cutoff:3.500A) removed outlier: 4.276A pdb=" N TYR A 772 " --> pdb=" O LEU A 725 " (cutoff:3.500A) No H-bonds generated for sheet with id=AB5 559 hydrogen bonds defined for protein. 1632 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.43 Time building geometry restraints manager: 2.81 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: 3501 1.34 - 1.46: 1890 1.46 - 1.58: 5472 1.58 - 1.69: 0 1.69 - 1.81: 83 Bond restraints: 10946 Sorted by residual: bond pdb=" N TRP B 159 " pdb=" CA TRP B 159 " ideal model delta sigma weight residual 1.457 1.485 -0.027 1.35e-02 5.49e+03 4.13e+00 bond pdb=" CB VAL A 183 " pdb=" CG1 VAL A 183 " ideal model delta sigma weight residual 1.521 1.466 0.055 3.30e-02 9.18e+02 2.78e+00 bond pdb=" N GLU B 161 " pdb=" CA GLU B 161 " ideal model delta sigma weight residual 1.458 1.487 -0.029 1.90e-02 2.77e+03 2.27e+00 bond pdb=" CB LYS A 647 " pdb=" CG LYS A 647 " ideal model delta sigma weight residual 1.520 1.557 -0.037 3.00e-02 1.11e+03 1.53e+00 bond pdb=" N CYS B 160 " pdb=" CA CYS B 160 " ideal model delta sigma weight residual 1.460 1.486 -0.026 2.20e-02 2.07e+03 1.45e+00 ... (remaining 10941 not shown) Histogram of bond angle deviations from ideal: 0.00 - 1.81: 14502 1.81 - 3.62: 236 3.62 - 5.43: 33 5.43 - 7.24: 13 7.24 - 9.05: 3 Bond angle restraints: 14787 Sorted by residual: angle pdb=" CA LYS A 647 " pdb=" CB LYS A 647 " pdb=" CG LYS A 647 " ideal model delta sigma weight residual 114.10 121.92 -7.82 2.00e+00 2.50e-01 1.53e+01 angle pdb=" CB ARG A 274 " pdb=" CG ARG A 274 " pdb=" CD ARG A 274 " ideal model delta sigma weight residual 111.30 119.44 -8.14 2.30e+00 1.89e-01 1.25e+01 angle pdb=" C VAL C 348 " pdb=" N VAL C 349 " pdb=" CA VAL C 349 " ideal model delta sigma weight residual 121.97 128.27 -6.30 1.80e+00 3.09e-01 1.22e+01 angle pdb=" C ALA A 412 " pdb=" N GLN A 413 " pdb=" CA GLN A 413 " ideal model delta sigma weight residual 121.54 127.69 -6.15 1.91e+00 2.74e-01 1.04e+01 angle pdb=" C ARG A 722 " pdb=" N LYS A 723 " pdb=" CA LYS A 723 " ideal model delta sigma weight residual 121.54 127.41 -5.87 1.91e+00 2.74e-01 9.44e+00 ... (remaining 14782 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.84: 5873 17.84 - 35.67: 650 35.67 - 53.51: 98 53.51 - 71.34: 26 71.34 - 89.18: 16 Dihedral angle restraints: 6663 sinusoidal: 2762 harmonic: 3901 Sorted by residual: dihedral pdb=" CA THR A 249 " pdb=" C THR A 249 " pdb=" N GLU A 250 " pdb=" CA GLU A 250 " ideal model delta harmonic sigma weight residual 180.00 153.33 26.67 0 5.00e+00 4.00e-02 2.85e+01 dihedral pdb=" CA ARG A 274 " pdb=" C ARG A 274 " pdb=" N VAL A 275 " pdb=" CA VAL A 275 " ideal model delta harmonic sigma weight residual 180.00 157.96 22.04 0 5.00e+00 4.00e-02 1.94e+01 dihedral pdb=" CA GLU A 271 " pdb=" C GLU A 271 " pdb=" N GLN A 272 " pdb=" CA GLN A 272 " ideal model delta harmonic sigma weight residual 180.00 160.37 19.63 0 5.00e+00 4.00e-02 1.54e+01 ... (remaining 6660 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.056: 1373 0.056 - 0.111: 228 0.111 - 0.167: 24 0.167 - 0.223: 3 0.223 - 0.278: 1 Chirality restraints: 1629 Sorted by residual: chirality pdb=" CB VAL C 349 " pdb=" CA VAL C 349 " pdb=" CG1 VAL C 349 " pdb=" CG2 VAL C 349 " both_signs ideal model delta sigma weight residual False -2.63 -2.35 -0.28 2.00e-01 2.50e+01 1.94e+00 chirality pdb=" CB VAL C 286 " pdb=" CA VAL C 286 " pdb=" CG1 VAL C 286 " pdb=" CG2 VAL C 286 " both_signs ideal model delta sigma weight residual False -2.63 -2.41 -0.22 2.00e-01 2.50e+01 1.16e+00 chirality pdb=" CA TYR A 244 " pdb=" N TYR A 244 " pdb=" C TYR A 244 " pdb=" CB TYR A 244 " both_signs ideal model delta sigma weight residual False 2.51 2.30 0.21 2.00e-01 2.50e+01 1.07e+00 ... (remaining 1626 not shown) Planarity restraints: 1898 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB TRP A 581 " -0.012 2.00e-02 2.50e+03 1.24e-02 3.86e+00 pdb=" CG TRP A 581 " 0.032 2.00e-02 2.50e+03 pdb=" CD1 TRP A 581 " -0.018 2.00e-02 2.50e+03 pdb=" CD2 TRP A 581 " -0.000 2.00e-02 2.50e+03 pdb=" NE1 TRP A 581 " 0.002 2.00e-02 2.50e+03 pdb=" CE2 TRP A 581 " 0.003 2.00e-02 2.50e+03 pdb=" CE3 TRP A 581 " -0.005 2.00e-02 2.50e+03 pdb=" CZ2 TRP A 581 " -0.001 2.00e-02 2.50e+03 pdb=" CZ3 TRP A 581 " 0.001 2.00e-02 2.50e+03 pdb=" CH2 TRP A 581 " -0.003 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C LYS A 664 " 0.029 5.00e-02 4.00e+02 4.44e-02 3.16e+00 pdb=" N PRO A 665 " -0.077 5.00e-02 4.00e+02 pdb=" CA PRO A 665 " 0.023 5.00e-02 4.00e+02 pdb=" CD PRO A 665 " 0.025 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C VAL A 746 " -0.028 5.00e-02 4.00e+02 4.15e-02 2.75e+00 pdb=" N PRO A 747 " 0.072 5.00e-02 4.00e+02 pdb=" CA PRO A 747 " -0.021 5.00e-02 4.00e+02 pdb=" CD PRO A 747 " -0.024 5.00e-02 4.00e+02 ... (remaining 1895 not shown) Histogram of nonbonded interaction distances: 2.24 - 2.77: 1982 2.77 - 3.30: 10536 3.30 - 3.83: 18492 3.83 - 4.37: 20851 4.37 - 4.90: 35536 Nonbonded interactions: 87397 Sorted by model distance: nonbonded pdb=" OG SER A 514 " pdb=" O VAL A 539 " model vdw 2.236 3.040 nonbonded pdb=" OH TYR B 109 " pdb=" OG1 THR A 47 " model vdw 2.241 3.040 nonbonded pdb=" O LYS A 459 " pdb=" OH TYR A 500 " model vdw 2.272 3.040 nonbonded pdb=" OH TYR A 260 " pdb=" O LEU A 295 " model vdw 2.289 3.040 nonbonded pdb=" OG1 THR A 552 " pdb=" O GLN A 630 " model vdw 2.289 3.040 ... (remaining 87392 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.00 ========== 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 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: 4.910 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: 0.450 Check model and map are aligned: 0.070 Set scattering table: 0.100 Process input model: 26.720 Find NCS groups from input model: 0.100 Set up NCS constraints: 0.040 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: 35.320 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7464 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.003 0.055 10946 Z= 0.151 Angle : 0.601 9.054 14787 Z= 0.318 Chirality : 0.042 0.278 1629 Planarity : 0.004 0.044 1898 Dihedral : 15.304 89.177 4135 Min Nonbonded Distance : 2.236 Molprobity Statistics. All-atom Clashscore : 7.25 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.05 % Favored : 95.95 % Rotamer: Outliers : 0.17 % Allowed : 0.17 % Favored : 99.66 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % 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.05 (0.22), residues: 1308 helix: 0.03 (0.18), residues: 714 sheet: -2.96 (0.47), residues: 76 loop : -1.24 (0.26), residues: 518 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.032 0.002 TRP A 581 HIS 0.007 0.001 HIS A 343 PHE 0.019 0.001 PHE A 566 TYR 0.026 0.002 TYR A 244 ARG 0.010 0.000 ARG A 274 Details of bonding type rmsd hydrogen bonds : bond 0.17836 ( 547) hydrogen bonds : angle 6.85546 ( 1632) covalent geometry : bond 0.00333 (10946) covalent geometry : angle 0.60062 (14787) *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 257 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 2 poor density : 255 time to evaluate : 1.164 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: B 130 LYS cc_start: 0.8150 (mmtp) cc_final: 0.7888 (mmtm) REVERT: C 73 GLN cc_start: 0.7756 (tt0) cc_final: 0.7491 (mt0) REVERT: C 80 ILE cc_start: 0.8371 (mt) cc_final: 0.8142 (mt) REVERT: C 137 LYS cc_start: 0.8141 (mmtt) cc_final: 0.7875 (ptpp) REVERT: C 325 GLN cc_start: 0.6942 (tm-30) cc_final: 0.6544 (tm-30) REVERT: A 346 ASN cc_start: 0.8209 (t0) cc_final: 0.7981 (t0) REVERT: A 450 MET cc_start: 0.7819 (tpp) cc_final: 0.7371 (mmp) REVERT: A 468 LYS cc_start: 0.8521 (mtpp) cc_final: 0.8190 (ttmm) outliers start: 2 outliers final: 0 residues processed: 256 average time/residue: 0.2767 time to fit residues: 97.3250 Evaluate side-chains 147 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 147 time to evaluate : 1.191 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 130 random chunks: chunk 109 optimal weight: 3.9990 chunk 98 optimal weight: 0.6980 chunk 54 optimal weight: 7.9990 chunk 33 optimal weight: 10.0000 chunk 66 optimal weight: 0.5980 chunk 52 optimal weight: 9.9990 chunk 101 optimal weight: 2.9990 chunk 39 optimal weight: 1.9990 chunk 61 optimal weight: 2.9990 chunk 75 optimal weight: 3.9990 chunk 118 optimal weight: 20.0000 overall best weight: 1.8586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 7 GLN B 90 GLN C 56 HIS C 72 ASN C 168 HIS A 361 ASN A 509 GLN ** A 584 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 641 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 7 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4360 r_free = 0.4360 target = 0.202610 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 66)----------------| | r_work = 0.4022 r_free = 0.4022 target = 0.170669 restraints weight = 14727.128| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 24)----------------| | r_work = 0.4018 r_free = 0.4018 target = 0.170954 restraints weight = 35310.405| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 24)----------------| | r_work = 0.4026 r_free = 0.4026 target = 0.172128 restraints weight = 29794.871| |-----------------------------------------------------------------------------| r_work (final): 0.4021 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7146 moved from start: 0.2357 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.007 0.065 10946 Z= 0.289 Angle : 0.692 8.345 14787 Z= 0.367 Chirality : 0.045 0.163 1629 Planarity : 0.005 0.054 1898 Dihedral : 4.333 26.188 1439 Min Nonbonded Distance : 2.453 Molprobity Statistics. All-atom Clashscore : 11.65 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.36 % Favored : 95.64 % Rotamer: Outliers : 1.68 % Allowed : 11.70 % Favored : 86.62 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % 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.64 (0.23), residues: 1308 helix: 0.34 (0.19), residues: 727 sheet: -2.77 (0.48), residues: 76 loop : -1.04 (0.28), residues: 505 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.020 0.002 TRP B 159 HIS 0.007 0.001 HIS A 343 PHE 0.026 0.002 PHE A 402 TYR 0.022 0.002 TYR C 363 ARG 0.006 0.001 ARG A 147 Details of bonding type rmsd hydrogen bonds : bond 0.04482 ( 547) hydrogen bonds : angle 5.17230 ( 1632) covalent geometry : bond 0.00679 (10946) covalent geometry : angle 0.69209 (14787) *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 168 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 20 poor density : 148 time to evaluate : 1.159 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: B 136 ARG cc_start: 0.8817 (mtm-85) cc_final: 0.8217 (mtm110) REVERT: A 157 TYR cc_start: 0.8427 (m-80) cc_final: 0.7881 (m-80) REVERT: A 468 LYS cc_start: 0.7939 (mtpp) cc_final: 0.7598 (ttmm) REVERT: A 508 PHE cc_start: 0.7736 (OUTLIER) cc_final: 0.7504 (m-80) outliers start: 20 outliers final: 12 residues processed: 160 average time/residue: 0.2221 time to fit residues: 50.5802 Evaluate side-chains 139 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 13 poor density : 126 time to evaluate : 1.215 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 3 SER Chi-restraints excluded: chain B residue 26 THR Chi-restraints excluded: chain B residue 95 VAL Chi-restraints excluded: chain B residue 131 THR Chi-restraints excluded: chain C residue 12 THR Chi-restraints excluded: chain C residue 307 THR Chi-restraints excluded: chain C residue 405 THR Chi-restraints excluded: chain A residue 115 VAL Chi-restraints excluded: chain A residue 341 GLU Chi-restraints excluded: chain A residue 392 VAL Chi-restraints excluded: chain A residue 451 VAL Chi-restraints excluded: chain A residue 508 PHE Chi-restraints excluded: chain A residue 511 ILE Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 130 random chunks: chunk 122 optimal weight: 7.9990 chunk 45 optimal weight: 5.9990 chunk 67 optimal weight: 0.5980 chunk 37 optimal weight: 0.9980 chunk 48 optimal weight: 10.0000 chunk 79 optimal weight: 0.5980 chunk 72 optimal weight: 0.0980 chunk 68 optimal weight: 0.1980 chunk 74 optimal weight: 0.2980 chunk 81 optimal weight: 0.7980 chunk 11 optimal weight: 0.5980 overall best weight: 0.3580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 152 GLN A 35 GLN A 361 ASN A 374 HIS ** A 584 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 641 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 694 GLN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4437 r_free = 0.4437 target = 0.210166 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 38)----------------| | r_work = 0.4148 r_free = 0.4148 target = 0.181509 restraints weight = 14814.375| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 26)----------------| | r_work = 0.4154 r_free = 0.4154 target = 0.183361 restraints weight = 31525.431| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 24)----------------| | r_work = 0.4154 r_free = 0.4154 target = 0.183612 restraints weight = 22375.747| |-----------------------------------------------------------------------------| r_work (final): 0.4145 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6952 moved from start: 0.2694 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.035 10946 Z= 0.122 Angle : 0.530 6.656 14787 Z= 0.281 Chirality : 0.039 0.159 1629 Planarity : 0.004 0.058 1898 Dihedral : 4.161 25.359 1439 Min Nonbonded Distance : 2.522 Molprobity Statistics. All-atom Clashscore : 10.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.20 % Favored : 95.80 % Rotamer: Outliers : 1.18 % Allowed : 14.14 % Favored : 84.68 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % 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.05 (0.23), residues: 1308 helix: 0.89 (0.19), residues: 727 sheet: -2.32 (0.49), residues: 85 loop : -0.89 (0.29), residues: 496 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.014 0.001 TRP A 21 HIS 0.008 0.001 HIS A 143 PHE 0.023 0.001 PHE C 184 TYR 0.023 0.001 TYR A 260 ARG 0.011 0.001 ARG C 170 Details of bonding type rmsd hydrogen bonds : bond 0.03670 ( 547) hydrogen bonds : angle 4.61626 ( 1632) covalent geometry : bond 0.00260 (10946) covalent geometry : angle 0.53018 (14787) *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 165 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 14 poor density : 151 time to evaluate : 3.173 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: B 159 TRP cc_start: 0.7855 (p-90) cc_final: 0.7381 (p-90) REVERT: A 346 ASN cc_start: 0.8327 (t0) cc_final: 0.8090 (t0) REVERT: A 468 LYS cc_start: 0.8089 (mtpp) cc_final: 0.7740 (ttmm) REVERT: A 574 HIS cc_start: 0.6740 (p90) cc_final: 0.6535 (p-80) outliers start: 14 outliers final: 6 residues processed: 159 average time/residue: 0.2747 time to fit residues: 64.1714 Evaluate side-chains 141 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 6 poor density : 135 time to evaluate : 1.681 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain C residue 156 VAL Chi-restraints excluded: chain C residue 177 VAL Chi-restraints excluded: chain A residue 37 MET Chi-restraints excluded: chain A residue 392 VAL Chi-restraints excluded: chain A residue 451 VAL Chi-restraints excluded: chain A residue 511 ILE Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 130 random chunks: chunk 38 optimal weight: 0.0000 chunk 47 optimal weight: 9.9990 chunk 35 optimal weight: 5.9990 chunk 102 optimal weight: 0.9990 chunk 65 optimal weight: 0.9990 chunk 49 optimal weight: 0.5980 chunk 103 optimal weight: 1.9990 chunk 17 optimal weight: 0.9990 chunk 52 optimal weight: 9.9990 chunk 112 optimal weight: 8.9990 chunk 50 optimal weight: 9.9990 overall best weight: 0.7190 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 152 GLN C 400 HIS A 30 GLN A 361 ASN ** A 584 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 641 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4418 r_free = 0.4418 target = 0.208149 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 58)----------------| | r_work = 0.4110 r_free = 0.4110 target = 0.178190 restraints weight = 14589.074| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 24)----------------| | r_work = 0.4103 r_free = 0.4103 target = 0.178622 restraints weight = 31229.043| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 26)----------------| | r_work = 0.4112 r_free = 0.4112 target = 0.179821 restraints weight = 27574.612| |-----------------------------------------------------------------------------| r_work (final): 0.4103 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7012 moved from start: 0.2973 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.036 10946 Z= 0.143 Angle : 0.539 7.193 14787 Z= 0.282 Chirality : 0.040 0.146 1629 Planarity : 0.004 0.056 1898 Dihedral : 4.062 25.329 1439 Min Nonbonded Distance : 2.509 Molprobity Statistics. All-atom Clashscore : 10.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.05 % Favored : 95.95 % Rotamer: Outliers : 1.68 % Allowed : 16.92 % Favored : 81.40 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % 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.08 (0.23), residues: 1308 helix: 1.01 (0.19), residues: 731 sheet: -2.17 (0.50), residues: 85 loop : -0.91 (0.29), residues: 492 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.001 TRP A 21 HIS 0.007 0.001 HIS A 143 PHE 0.018 0.001 PHE C 184 TYR 0.015 0.001 TYR A 278 ARG 0.003 0.000 ARG C 187 Details of bonding type rmsd hydrogen bonds : bond 0.03534 ( 547) hydrogen bonds : angle 4.48462 ( 1632) covalent geometry : bond 0.00329 (10946) covalent geometry : angle 0.53911 (14787) *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 159 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 20 poor density : 139 time to evaluate : 1.200 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: B 159 TRP cc_start: 0.7752 (p-90) cc_final: 0.7406 (p-90) REVERT: A 261 MET cc_start: 0.7434 (mtm) cc_final: 0.6643 (mtm) REVERT: A 346 ASN cc_start: 0.8401 (t0) cc_final: 0.8150 (t0) REVERT: A 468 LYS cc_start: 0.7887 (mtpp) cc_final: 0.7594 (ttmm) outliers start: 20 outliers final: 11 residues processed: 151 average time/residue: 0.2244 time to fit residues: 49.8154 Evaluate side-chains 141 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 11 poor density : 130 time to evaluate : 1.359 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 95 VAL Chi-restraints excluded: chain B residue 141 ILE Chi-restraints excluded: chain C residue 156 VAL Chi-restraints excluded: chain C residue 177 VAL Chi-restraints excluded: chain C residue 204 LEU Chi-restraints excluded: chain C residue 368 THR Chi-restraints excluded: chain A residue 37 MET Chi-restraints excluded: chain A residue 115 VAL Chi-restraints excluded: chain A residue 157 TYR Chi-restraints excluded: chain A residue 392 VAL Chi-restraints excluded: chain A residue 451 VAL Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 130 random chunks: chunk 95 optimal weight: 0.8980 chunk 42 optimal weight: 0.0050 chunk 68 optimal weight: 2.9990 chunk 52 optimal weight: 0.9990 chunk 113 optimal weight: 1.9990 chunk 50 optimal weight: 7.9990 chunk 27 optimal weight: 10.0000 chunk 79 optimal weight: 1.9990 chunk 99 optimal weight: 0.5980 chunk 60 optimal weight: 0.3980 chunk 91 optimal weight: 0.9990 overall best weight: 0.5796 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: C 152 GLN ** C 267 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** R 28 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 361 ASN A 518 ASN ** A 641 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4301 r_free = 0.4301 target = 0.198107 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 42)----------------| | r_work = 0.4005 r_free = 0.4005 target = 0.169261 restraints weight = 15086.735| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 46)----------------| | r_work = 0.3973 r_free = 0.3973 target = 0.168363 restraints weight = 33398.774| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 30)----------------| | r_work = 0.3949 r_free = 0.3949 target = 0.166096 restraints weight = 20924.753| |-----------------------------------------------------------------------------| r_work (final): 0.3936 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7106 moved from start: 0.3195 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.064 10946 Z= 0.133 Angle : 0.546 14.115 14787 Z= 0.281 Chirality : 0.039 0.144 1629 Planarity : 0.004 0.053 1898 Dihedral : 4.077 31.586 1439 Min Nonbonded Distance : 2.466 Molprobity Statistics. All-atom Clashscore : 9.92 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.20 % Favored : 95.80 % Rotamer: Outliers : 1.85 % Allowed : 18.10 % Favored : 80.05 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.08 % 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: 1308 helix: 1.07 (0.19), residues: 734 sheet: -2.17 (0.48), residues: 99 loop : -0.84 (0.29), residues: 475 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.027 0.002 TRP C 248 HIS 0.008 0.001 HIS A 143 PHE 0.034 0.002 PHE C 222 TYR 0.016 0.001 TYR A 260 ARG 0.003 0.000 ARG C 187 Details of bonding type rmsd hydrogen bonds : bond 0.03433 ( 547) hydrogen bonds : angle 4.38976 ( 1632) covalent geometry : bond 0.00304 (10946) covalent geometry : angle 0.54573 (14787) *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 165 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 22 poor density : 143 time to evaluate : 1.206 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 346 ASN cc_start: 0.8367 (t0) cc_final: 0.8116 (t0) REVERT: A 468 LYS cc_start: 0.7909 (mtpp) cc_final: 0.7678 (ttmm) outliers start: 22 outliers final: 13 residues processed: 158 average time/residue: 0.2331 time to fit residues: 54.6747 Evaluate side-chains 143 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 13 poor density : 130 time to evaluate : 1.634 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 95 VAL Chi-restraints excluded: chain B residue 141 ILE Chi-restraints excluded: chain C residue 156 VAL Chi-restraints excluded: chain C residue 177 VAL Chi-restraints excluded: chain C residue 204 LEU Chi-restraints excluded: chain C residue 247 ASP Chi-restraints excluded: chain C residue 307 THR Chi-restraints excluded: chain C residue 368 THR Chi-restraints excluded: chain C residue 375 TYR Chi-restraints excluded: chain A residue 37 MET Chi-restraints excluded: chain A residue 207 SER Chi-restraints excluded: chain A residue 392 VAL Chi-restraints excluded: chain A residue 451 VAL Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 130 random chunks: chunk 39 optimal weight: 3.9990 chunk 96 optimal weight: 0.0040 chunk 21 optimal weight: 1.9990 chunk 52 optimal weight: 0.7980 chunk 115 optimal weight: 20.0000 chunk 111 optimal weight: 0.7980 chunk 16 optimal weight: 0.2980 chunk 90 optimal weight: 1.9990 chunk 31 optimal weight: 6.9990 chunk 122 optimal weight: 10.0000 chunk 128 optimal weight: 4.9990 overall best weight: 0.7794 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 152 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 361 ASN ** A 641 GLN ** 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.4286 r_free = 0.4286 target = 0.196981 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 39)----------------| | r_work = 0.3984 r_free = 0.3984 target = 0.167603 restraints weight = 14841.249| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 37)----------------| | r_work = 0.3967 r_free = 0.3967 target = 0.167788 restraints weight = 33349.776| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 30)----------------| | r_work = 0.3951 r_free = 0.3951 target = 0.166380 restraints weight = 20436.450| |-----------------------------------------------------------------------------| r_work (final): 0.3947 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7083 moved from start: 0.3358 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.036 10946 Z= 0.148 Angle : 0.562 8.723 14787 Z= 0.291 Chirality : 0.040 0.156 1629 Planarity : 0.003 0.049 1898 Dihedral : 4.042 26.111 1439 Min Nonbonded Distance : 2.387 Molprobity Statistics. All-atom Clashscore : 11.23 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.43 % Favored : 95.57 % Rotamer: Outliers : 1.52 % Allowed : 20.03 % Favored : 78.45 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % 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.27 (0.24), residues: 1308 helix: 1.14 (0.19), residues: 729 sheet: -2.17 (0.56), residues: 77 loop : -0.79 (0.28), residues: 502 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.015 0.001 TRP C 248 HIS 0.008 0.001 HIS A 143 PHE 0.032 0.002 PHE C 399 TYR 0.014 0.001 TYR C 89 ARG 0.003 0.000 ARG C 88 Details of bonding type rmsd hydrogen bonds : bond 0.03413 ( 547) hydrogen bonds : angle 4.37098 ( 1632) covalent geometry : bond 0.00345 (10946) covalent geometry : angle 0.56183 (14787) *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 155 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 18 poor density : 137 time to evaluate : 1.180 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 244 TYR cc_start: 0.8540 (m-10) cc_final: 0.8010 (m-10) REVERT: A 346 ASN cc_start: 0.8470 (t0) cc_final: 0.8230 (t0) REVERT: A 468 LYS cc_start: 0.7965 (mtpp) cc_final: 0.7677 (ttmm) outliers start: 18 outliers final: 14 residues processed: 150 average time/residue: 0.2068 time to fit residues: 45.7379 Evaluate side-chains 137 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 14 poor density : 123 time to evaluate : 1.199 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 95 VAL Chi-restraints excluded: chain B residue 141 ILE Chi-restraints excluded: chain C residue 156 VAL Chi-restraints excluded: chain C residue 177 VAL Chi-restraints excluded: chain C residue 204 LEU Chi-restraints excluded: chain C residue 307 THR Chi-restraints excluded: chain C residue 368 THR Chi-restraints excluded: chain C residue 375 TYR Chi-restraints excluded: chain C residue 405 THR Chi-restraints excluded: chain A residue 37 MET Chi-restraints excluded: chain A residue 115 VAL Chi-restraints excluded: chain A residue 207 SER Chi-restraints excluded: chain A residue 392 VAL Chi-restraints excluded: chain A residue 451 VAL Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 130 random chunks: chunk 55 optimal weight: 9.9990 chunk 51 optimal weight: 8.9990 chunk 60 optimal weight: 0.5980 chunk 123 optimal weight: 0.9990 chunk 49 optimal weight: 20.0000 chunk 75 optimal weight: 3.9990 chunk 35 optimal weight: 6.9990 chunk 6 optimal weight: 0.8980 chunk 121 optimal weight: 9.9990 chunk 36 optimal weight: 2.9990 chunk 88 optimal weight: 1.9990 overall best weight: 1.4986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 152 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 201 HIS ** A 641 GLN ** 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.4241 r_free = 0.4241 target = 0.192507 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 41)----------------| | r_work = 0.3935 r_free = 0.3935 target = 0.163359 restraints weight = 15376.882| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 40)----------------| | r_work = 0.3911 r_free = 0.3911 target = 0.163201 restraints weight = 35599.388| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 28)----------------| | r_work = 0.3872 r_free = 0.3872 target = 0.159767 restraints weight = 20924.867| |-----------------------------------------------------------------------------| r_work (final): 0.3868 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7235 moved from start: 0.3575 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.056 10946 Z= 0.231 Angle : 0.634 8.690 14787 Z= 0.327 Chirality : 0.043 0.144 1629 Planarity : 0.004 0.053 1898 Dihedral : 4.149 25.920 1439 Min Nonbonded Distance : 2.419 Molprobity Statistics. All-atom Clashscore : 13.29 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.51 % Favored : 95.49 % Rotamer: Outliers : 2.27 % Allowed : 20.88 % Favored : 76.85 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % 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.07 (0.23), residues: 1308 helix: 0.95 (0.19), residues: 728 sheet: -2.21 (0.56), residues: 77 loop : -0.84 (0.28), residues: 503 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.013 0.002 TRP C 248 HIS 0.009 0.002 HIS A 143 PHE 0.026 0.002 PHE A 402 TYR 0.018 0.002 TYR A 278 ARG 0.004 0.000 ARG C 55 Details of bonding type rmsd hydrogen bonds : bond 0.03720 ( 547) hydrogen bonds : angle 4.52586 ( 1632) covalent geometry : bond 0.00552 (10946) covalent geometry : angle 0.63402 (14787) *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 153 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 27 poor density : 126 time to evaluate : 1.174 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: C 34 ILE cc_start: 0.8768 (OUTLIER) cc_final: 0.8545 (mt) REVERT: A 157 TYR cc_start: 0.8829 (OUTLIER) cc_final: 0.7954 (m-80) REVERT: A 468 LYS cc_start: 0.8030 (mtpp) cc_final: 0.7655 (ttmm) outliers start: 27 outliers final: 19 residues processed: 146 average time/residue: 0.2021 time to fit residues: 44.5208 Evaluate side-chains 139 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 21 poor density : 118 time to evaluate : 1.201 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 94 LYS Chi-restraints excluded: chain B residue 95 VAL Chi-restraints excluded: chain B residue 131 THR Chi-restraints excluded: chain B residue 141 ILE Chi-restraints excluded: chain C residue 33 LEU Chi-restraints excluded: chain C residue 34 ILE Chi-restraints excluded: chain C residue 156 VAL Chi-restraints excluded: chain C residue 177 VAL Chi-restraints excluded: chain C residue 204 LEU Chi-restraints excluded: chain C residue 307 THR Chi-restraints excluded: chain C residue 368 THR Chi-restraints excluded: chain C residue 375 TYR Chi-restraints excluded: chain C residue 405 THR Chi-restraints excluded: chain A residue 37 MET Chi-restraints excluded: chain A residue 115 VAL Chi-restraints excluded: chain A residue 157 TYR Chi-restraints excluded: chain A residue 207 SER Chi-restraints excluded: chain A residue 392 VAL Chi-restraints excluded: chain A residue 402 PHE Chi-restraints excluded: chain A residue 451 VAL Chi-restraints excluded: chain A residue 508 PHE Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 130 random chunks: chunk 76 optimal weight: 0.6980 chunk 99 optimal weight: 0.8980 chunk 69 optimal weight: 0.6980 chunk 86 optimal weight: 1.9990 chunk 4 optimal weight: 1.9990 chunk 52 optimal weight: 9.9990 chunk 79 optimal weight: 0.5980 chunk 34 optimal weight: 0.9980 chunk 3 optimal weight: 0.6980 chunk 45 optimal weight: 0.7980 chunk 127 optimal weight: 0.8980 overall best weight: 0.6980 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 152 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 361 ASN ** A 641 GLN ** 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.4281 r_free = 0.4281 target = 0.196540 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 42)----------------| | r_work = 0.3977 r_free = 0.3977 target = 0.166987 restraints weight = 14784.357| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 26)----------------| | r_work = 0.3982 r_free = 0.3982 target = 0.168539 restraints weight = 34556.488| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 27)----------------| | r_work = 0.3974 r_free = 0.3974 target = 0.167927 restraints weight = 20422.642| |-----------------------------------------------------------------------------| r_work (final): 0.3974 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7029 moved from start: 0.3666 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.035 10946 Z= 0.143 Angle : 0.576 9.999 14787 Z= 0.297 Chirality : 0.040 0.142 1629 Planarity : 0.004 0.056 1898 Dihedral : 4.075 25.096 1439 Min Nonbonded Distance : 2.425 Molprobity Statistics. All-atom Clashscore : 12.59 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.20 % Favored : 95.80 % Rotamer: Outliers : 1.85 % Allowed : 20.96 % Favored : 77.19 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % 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.18 (0.23), residues: 1308 helix: 1.03 (0.19), residues: 734 sheet: -1.91 (0.56), residues: 75 loop : -0.85 (0.28), residues: 499 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.017 0.001 TRP C 248 HIS 0.010 0.001 HIS A 143 PHE 0.021 0.002 PHE C 399 TYR 0.015 0.001 TYR A 278 ARG 0.003 0.000 ARG C 332 Details of bonding type rmsd hydrogen bonds : bond 0.03434 ( 547) hydrogen bonds : angle 4.36878 ( 1632) covalent geometry : bond 0.00333 (10946) covalent geometry : angle 0.57642 (14787) *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 150 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 22 poor density : 128 time to evaluate : 1.245 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 157 TYR cc_start: 0.8778 (OUTLIER) cc_final: 0.7961 (m-80) REVERT: A 244 TYR cc_start: 0.8581 (m-10) cc_final: 0.7989 (m-10) REVERT: A 468 LYS cc_start: 0.7884 (mtpp) cc_final: 0.7587 (ttmm) outliers start: 22 outliers final: 16 residues processed: 142 average time/residue: 0.2306 time to fit residues: 49.3088 Evaluate side-chains 137 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 17 poor density : 120 time to evaluate : 1.508 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 94 LYS Chi-restraints excluded: chain B residue 95 VAL Chi-restraints excluded: chain B residue 131 THR Chi-restraints excluded: chain B residue 141 ILE Chi-restraints excluded: chain C residue 33 LEU Chi-restraints excluded: chain C residue 156 VAL Chi-restraints excluded: chain C residue 177 VAL Chi-restraints excluded: chain C residue 204 LEU Chi-restraints excluded: chain C residue 307 THR Chi-restraints excluded: chain C residue 368 THR Chi-restraints excluded: chain C residue 375 TYR Chi-restraints excluded: chain C residue 405 THR Chi-restraints excluded: chain A residue 157 TYR Chi-restraints excluded: chain A residue 207 SER Chi-restraints excluded: chain A residue 392 VAL Chi-restraints excluded: chain A residue 451 VAL Chi-restraints excluded: chain A residue 508 PHE Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 130 random chunks: chunk 1 optimal weight: 5.9990 chunk 129 optimal weight: 0.0770 chunk 15 optimal weight: 2.9990 chunk 53 optimal weight: 0.5980 chunk 111 optimal weight: 3.9990 chunk 107 optimal weight: 1.9990 chunk 0 optimal weight: 3.9990 chunk 56 optimal weight: 0.9990 chunk 60 optimal weight: 0.9990 chunk 70 optimal weight: 0.5980 chunk 6 optimal weight: 0.1980 overall best weight: 0.4940 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 152 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 641 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** A 694 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4304 r_free = 0.4304 target = 0.198557 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 56)----------------| | r_work = 0.3970 r_free = 0.3970 target = 0.167149 restraints weight = 14672.600| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 29)----------------| | r_work = 0.3952 r_free = 0.3952 target = 0.166389 restraints weight = 29638.067| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 26)----------------| | r_work = 0.3962 r_free = 0.3962 target = 0.167373 restraints weight = 29312.650| |-----------------------------------------------------------------------------| r_work (final): 0.3959 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7056 moved from start: 0.3774 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.036 10946 Z= 0.125 Angle : 0.565 10.951 14787 Z= 0.290 Chirality : 0.040 0.144 1629 Planarity : 0.003 0.054 1898 Dihedral : 4.047 24.971 1439 Min Nonbonded Distance : 2.415 Molprobity Statistics. All-atom Clashscore : 11.84 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.05 % Favored : 95.95 % Rotamer: Outliers : 1.68 % Allowed : 21.30 % Favored : 77.02 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % 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.26 (0.24), residues: 1308 helix: 1.15 (0.19), residues: 732 sheet: -2.29 (0.56), residues: 77 loop : -0.81 (0.28), residues: 499 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.016 0.002 TRP C 248 HIS 0.009 0.001 HIS A 143 PHE 0.019 0.001 PHE C 399 TYR 0.015 0.001 TYR A 52 ARG 0.003 0.000 ARG A 506 Details of bonding type rmsd hydrogen bonds : bond 0.03337 ( 547) hydrogen bonds : angle 4.32111 ( 1632) covalent geometry : bond 0.00286 (10946) covalent geometry : angle 0.56548 (14787) *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 145 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 20 poor density : 125 time to evaluate : 1.245 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: R 27 TRP cc_start: 0.7705 (t60) cc_final: 0.6549 (t60) REVERT: A 157 TYR cc_start: 0.8741 (OUTLIER) cc_final: 0.7944 (m-80) REVERT: A 244 TYR cc_start: 0.8438 (m-10) cc_final: 0.7912 (m-10) REVERT: A 468 LYS cc_start: 0.7785 (mtpp) cc_final: 0.7524 (ttmm) outliers start: 20 outliers final: 18 residues processed: 139 average time/residue: 0.2619 time to fit residues: 55.2244 Evaluate side-chains 135 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 19 poor density : 116 time to evaluate : 1.266 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 94 LYS Chi-restraints excluded: chain B residue 95 VAL Chi-restraints excluded: chain B residue 131 THR Chi-restraints excluded: chain B residue 141 ILE Chi-restraints excluded: chain C residue 156 VAL Chi-restraints excluded: chain C residue 177 VAL Chi-restraints excluded: chain C residue 204 LEU Chi-restraints excluded: chain C residue 307 THR Chi-restraints excluded: chain C residue 368 THR Chi-restraints excluded: chain C residue 375 TYR Chi-restraints excluded: chain C residue 405 THR Chi-restraints excluded: chain A residue 115 VAL Chi-restraints excluded: chain A residue 157 TYR Chi-restraints excluded: chain A residue 207 SER Chi-restraints excluded: chain A residue 275 VAL Chi-restraints excluded: chain A residue 392 VAL Chi-restraints excluded: chain A residue 451 VAL Chi-restraints excluded: chain A residue 507 MET Chi-restraints excluded: chain A residue 508 PHE Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 130 random chunks: chunk 86 optimal weight: 1.9990 chunk 112 optimal weight: 30.0000 chunk 95 optimal weight: 0.4980 chunk 87 optimal weight: 0.5980 chunk 33 optimal weight: 5.9990 chunk 70 optimal weight: 0.3980 chunk 106 optimal weight: 6.9990 chunk 114 optimal weight: 8.9990 chunk 91 optimal weight: 0.9990 chunk 30 optimal weight: 0.0070 chunk 40 optimal weight: 0.8980 overall best weight: 0.4798 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 152 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 267 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 641 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.4313 r_free = 0.4313 target = 0.199532 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 67)----------------| | r_work = 0.3969 r_free = 0.3969 target = 0.167416 restraints weight = 14675.807| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 26)----------------| | r_work = 0.3932 r_free = 0.3932 target = 0.164803 restraints weight = 28652.843| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 26)----------------| | r_work = 0.3943 r_free = 0.3943 target = 0.165944 restraints weight = 30078.147| |-----------------------------------------------------------------------------| r_work (final): 0.3929 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7110 moved from start: 0.3916 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.041 10946 Z= 0.126 Angle : 0.573 11.599 14787 Z= 0.294 Chirality : 0.040 0.143 1629 Planarity : 0.004 0.053 1898 Dihedral : 3.996 25.082 1439 Min Nonbonded Distance : 2.414 Molprobity Statistics. All-atom Clashscore : 11.56 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.28 % Favored : 95.72 % Rotamer: Outliers : 1.60 % Allowed : 21.38 % Favored : 77.02 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % 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.27 (0.24), residues: 1308 helix: 1.22 (0.19), residues: 729 sheet: -2.36 (0.48), residues: 95 loop : -0.80 (0.29), residues: 484 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.059 0.002 TRP B 159 HIS 0.009 0.001 HIS A 143 PHE 0.020 0.002 PHE C 399 TYR 0.017 0.001 TYR C 361 ARG 0.002 0.000 ARG C 55 Details of bonding type rmsd hydrogen bonds : bond 0.03281 ( 547) hydrogen bonds : angle 4.23502 ( 1632) covalent geometry : bond 0.00290 (10946) covalent geometry : angle 0.57293 (14787) ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2616 Ramachandran restraints generated. 1308 Oldfield, 0 Emsley, 1308 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 138 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 19 poor density : 119 time to evaluate : 1.190 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: R 27 TRP cc_start: 0.7815 (t60) cc_final: 0.6638 (t60) REVERT: A 157 TYR cc_start: 0.8722 (OUTLIER) cc_final: 0.7842 (m-80) REVERT: A 244 TYR cc_start: 0.8428 (m-10) cc_final: 0.7898 (m-10) REVERT: A 468 LYS cc_start: 0.7787 (mtpp) cc_final: 0.7557 (ttmm) outliers start: 19 outliers final: 17 residues processed: 132 average time/residue: 0.2016 time to fit residues: 39.9545 Evaluate side-chains 131 residues out of total 1188 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 18 poor density : 113 time to evaluate : 1.182 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain B residue 94 LYS Chi-restraints excluded: chain B residue 95 VAL Chi-restraints excluded: chain B residue 131 THR Chi-restraints excluded: chain B residue 141 ILE Chi-restraints excluded: chain C residue 156 VAL Chi-restraints excluded: chain C residue 177 VAL Chi-restraints excluded: chain C residue 307 THR Chi-restraints excluded: chain C residue 368 THR Chi-restraints excluded: chain C residue 375 TYR Chi-restraints excluded: chain C residue 405 THR Chi-restraints excluded: chain A residue 115 VAL Chi-restraints excluded: chain A residue 157 TYR Chi-restraints excluded: chain A residue 207 SER Chi-restraints excluded: chain A residue 275 VAL Chi-restraints excluded: chain A residue 392 VAL Chi-restraints excluded: chain A residue 451 VAL Chi-restraints excluded: chain A residue 507 MET Chi-restraints excluded: chain A residue 508 PHE Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 130 random chunks: chunk 23 optimal weight: 0.7980 chunk 64 optimal weight: 0.9990 chunk 119 optimal weight: 7.9990 chunk 39 optimal weight: 0.9980 chunk 25 optimal weight: 2.9990 chunk 94 optimal weight: 0.0770 chunk 10 optimal weight: 0.0980 chunk 126 optimal weight: 9.9990 chunk 128 optimal weight: 0.5980 chunk 45 optimal weight: 0.9990 chunk 101 optimal weight: 6.9990 overall best weight: 0.5138 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 152 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 641 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.4315 r_free = 0.4315 target = 0.199536 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 46)----------------| | r_work = 0.3993 r_free = 0.3993 target = 0.168809 restraints weight = 14745.069| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 31)----------------| | r_work = 0.3970 r_free = 0.3970 target = 0.167166 restraints weight = 29896.561| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 28)----------------| | r_work = 0.3984 r_free = 0.3984 target = 0.168571 restraints weight = 24764.697| |-----------------------------------------------------------------------------| r_work (final): 0.3974 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7025 moved from start: 0.3987 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.036 10946 Z= 0.126 Angle : 0.570 11.023 14787 Z= 0.293 Chirality : 0.040 0.140 1629 Planarity : 0.004 0.051 1898 Dihedral : 3.968 25.120 1439 Min Nonbonded Distance : 2.395 Molprobity Statistics. All-atom Clashscore : 11.28 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.05 % Favored : 95.95 % Rotamer: Outliers : 1.43 % Allowed : 21.63 % Favored : 76.94 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 4.65 % 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.34 (0.24), residues: 1308 helix: 1.24 (0.19), residues: 731 sheet: -2.51 (0.54), residues: 77 loop : -0.78 (0.28), residues: 500 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.022 0.002 TRP B 159 HIS 0.009 0.001 HIS A 143 PHE 0.033 0.002 PHE C 399 TYR 0.016 0.001 TYR C 361 ARG 0.003 0.000 ARG C 55 Details of bonding type rmsd hydrogen bonds : bond 0.03214 ( 547) hydrogen bonds : angle 4.24408 ( 1632) covalent geometry : bond 0.00289 (10946) covalent geometry : angle 0.56970 (14787) Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 3665.30 seconds wall clock time: 65 minutes 56.01 seconds (3956.01 seconds total)