Starting phenix.real_space_refine on Tue Feb 11 20:26:18 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/7kmb_22922/02_2025/7kmb_22922.cif Found real_map, /net/cci-nas-00/data/ceres_data/7kmb_22922/02_2025/7kmb_22922.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.39 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/7kmb_22922/02_2025/7kmb_22922.map" default_real_map = "/net/cci-nas-00/data/ceres_data/7kmb_22922/02_2025/7kmb_22922.map" model { file = "/net/cci-nas-00/data/ceres_data/7kmb_22922/02_2025/7kmb_22922.cif" } default_model = "/net/cci-nas-00/data/ceres_data/7kmb_22922/02_2025/7kmb_22922.cif" } resolution = 3.39 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.001 sd= 0.031 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 37 5.16 5 C 4159 2.51 5 N 1067 2.21 5 O 1246 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 51 residue(s): 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/chem_data/mon_lib" Total number of atoms: 6509 Number of models: 1 Model: "" Number of chains: 5 Chain: "F" Number of atoms: 4862 Number of conformers: 1 Conformer: "" Number of residues, atoms: 596, 4862 Classifications: {'peptide': 596} Link IDs: {'PTRANS': 27, 'TRANS': 568} Chain: "G" Number of atoms: 1521 Number of conformers: 1 Conformer: "" Number of residues, atoms: 192, 1521 Classifications: {'peptide': 192} Link IDs: {'PTRANS': 10, 'TRANS': 181} Chain: "A" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Chain: "B" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Chain: "F" Number of atoms: 70 Number of conformers: 1 Conformer: "" Number of residues, atoms: 5, 70 Unusual residues: {'NAG': 5} Classifications: {'undetermined': 5} Link IDs: {None: 4} Unresolved non-hydrogen bonds: 5 Unresolved non-hydrogen angles: 10 Unresolved non-hydrogen dihedrals: 15 Unresolved non-hydrogen chiralities: 5 Time building chain proxies: 6.36, per 1000 atoms: 0.98 Number of scatterers: 6509 At special positions: 0 Unit cell: (99.452, 75.118, 111.09, 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 37 16.00 O 1246 8.00 N 1067 7.00 C 4159 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=7, symmetry=0 Simple disulfide: pdb=" SG CYS F 133 " - pdb=" SG CYS F 141 " distance=2.03 Simple disulfide: pdb=" SG CYS F 344 " - pdb=" SG CYS F 361 " distance=2.03 Simple disulfide: pdb=" SG CYS F 530 " - pdb=" SG CYS F 542 " distance=2.03 Simple disulfide: pdb=" SG CYS G 336 " - pdb=" SG CYS G 361 " distance=2.03 Simple disulfide: pdb=" SG CYS G 379 " - pdb=" SG CYS G 432 " distance=2.03 Simple disulfide: pdb=" SG CYS G 391 " - pdb=" SG CYS G 525 " distance=2.03 Simple disulfide: pdb=" SG CYS G 480 " - pdb=" SG CYS G 488 " distance=2.03 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 Links applied BETA1-4 " NAG A 1 " - " NAG A 2 " " NAG B 1 " - " NAG B 2 " NAG-ASN " NAG A 1 " - " ASN F 546 " " NAG B 1 " - " ASN G 343 " " NAG F 701 " - " ASN F 322 " " NAG F 702 " - " ASN F 432 " " NAG F 703 " - " ASN F 103 " " NAG F 704 " - " ASN F 90 " " NAG F 705 " - " ASN F 53 " Time building additional restraints: 1.75 Conformation dependent library (CDL) restraints added in 1.1 seconds 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1482 Finding SS restraints... Secondary structure from input PDB file: 36 helices and 4 sheets defined 52.8% alpha, 4.8% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.62 Creating SS restraints... Processing helix chain 'F' and resid 20 through 51 removed outlier: 3.788A pdb=" N HIS F 34 " --> pdb=" O ASP F 30 " (cutoff:3.500A) removed outlier: 3.634A pdb=" N ASP F 38 " --> pdb=" O HIS F 34 " (cutoff:3.500A) Processing helix chain 'F' and resid 55 through 82 removed outlier: 3.796A pdb=" N VAL F 59 " --> pdb=" O THR F 55 " (cutoff:3.500A) removed outlier: 3.965A pdb=" N GLN F 60 " --> pdb=" O GLU F 56 " (cutoff:3.500A) removed outlier: 3.524A pdb=" N ALA F 65 " --> pdb=" O ASN F 61 " (cutoff:3.500A) removed outlier: 3.583A pdb=" N THR F 78 " --> pdb=" O LYS F 74 " (cutoff:3.500A) removed outlier: 4.275A pdb=" N MET F 82 " --> pdb=" O THR F 78 " (cutoff:3.500A) Processing helix chain 'F' and resid 91 through 101 removed outlier: 3.714A pdb=" N GLN F 96 " --> pdb=" O THR F 92 " (cutoff:3.500A) Processing helix chain 'F' and resid 109 through 128 removed outlier: 4.551A pdb=" N ARG F 115 " --> pdb=" O ASP F 111 " (cutoff:3.500A) removed outlier: 4.016A pdb=" N LEU F 116 " --> pdb=" O LYS F 112 " (cutoff:3.500A) removed outlier: 3.975A pdb=" N ASN F 117 " --> pdb=" O SER F 113 " (cutoff:3.500A) removed outlier: 3.502A pdb=" N THR F 118 " --> pdb=" O LYS F 114 " (cutoff:3.500A) removed outlier: 3.588A pdb=" N ILE F 119 " --> pdb=" O ARG F 115 " (cutoff:3.500A) Processing helix chain 'F' and resid 143 through 154 removed outlier: 4.223A pdb=" N ASN F 149 " --> pdb=" O GLU F 145 " (cutoff:3.500A) removed outlier: 4.741A pdb=" N GLU F 150 " --> pdb=" O PRO F 146 " (cutoff:3.500A) removed outlier: 3.980A pdb=" N ILE F 151 " --> pdb=" O GLY F 147 " (cutoff:3.500A) removed outlier: 3.772A pdb=" N MET F 152 " --> pdb=" O LEU F 148 " (cutoff:3.500A) Processing helix chain 'F' and resid 157 through 170 removed outlier: 3.751A pdb=" N TRP F 168 " --> pdb=" O ALA F 164 " (cutoff:3.500A) removed outlier: 3.507A pdb=" N SER F 170 " --> pdb=" O GLU F 166 " (cutoff:3.500A) Processing helix chain 'F' and resid 172 through 175 removed outlier: 3.608A pdb=" N GLN F 175 " --> pdb=" O VAL F 172 " (cutoff:3.500A) No H-bonds generated for 'chain 'F' and resid 172 through 175' Processing helix chain 'F' and resid 176 through 194 removed outlier: 3.607A pdb=" N GLU F 181 " --> pdb=" O ARG F 177 " (cutoff:3.500A) removed outlier: 3.655A pdb=" N GLU F 182 " --> pdb=" O PRO F 178 " (cutoff:3.500A) removed outlier: 3.753A pdb=" N LEU F 186 " --> pdb=" O GLU F 182 " (cutoff:3.500A) removed outlier: 4.038A pdb=" N ASN F 188 " --> pdb=" O VAL F 184 " (cutoff:3.500A) removed outlier: 3.992A pdb=" N GLU F 189 " --> pdb=" O VAL F 185 " (cutoff:3.500A) removed outlier: 4.041A pdb=" N ALA F 191 " --> pdb=" O LYS F 187 " (cutoff:3.500A) removed outlier: 4.332A pdb=" N ARG F 192 " --> pdb=" O ASN F 188 " (cutoff:3.500A) removed outlier: 3.527A pdb=" N ALA F 193 " --> pdb=" O GLU F 189 " (cutoff:3.500A) Processing helix chain 'F' and resid 198 through 204 removed outlier: 3.569A pdb=" N ARG F 204 " --> pdb=" O GLY F 200 " (cutoff:3.500A) Processing helix chain 'F' and resid 220 through 249 removed outlier: 3.707A pdb=" N GLU F 224 " --> pdb=" O GLY F 220 " (cutoff:3.500A) removed outlier: 4.002A pdb=" N GLU F 231 " --> pdb=" O GLU F 227 " (cutoff:3.500A) removed outlier: 4.892A pdb=" N GLU F 232 " --> pdb=" O HIS F 228 " (cutoff:3.500A) removed outlier: 3.628A pdb=" N ILE F 233 " --> pdb=" O THR F 229 " (cutoff:3.500A) Proline residue: F 235 - end of helix removed outlier: 3.778A pdb=" N HIS F 239 " --> pdb=" O PRO F 235 " (cutoff:3.500A) removed outlier: 3.685A pdb=" N ALA F 246 " --> pdb=" O ALA F 242 " (cutoff:3.500A) removed outlier: 3.730A pdb=" N MET F 249 " --> pdb=" O ARG F 245 " (cutoff:3.500A) Processing helix chain 'F' and resid 265 through 267 No H-bonds generated for 'chain 'F' and resid 265 through 267' Processing helix chain 'F' and resid 275 through 279 Processing helix chain 'F' and resid 293 through 300 removed outlier: 3.967A pdb=" N VAL F 298 " --> pdb=" O THR F 294 " (cutoff:3.500A) removed outlier: 3.729A pdb=" N ASP F 299 " --> pdb=" O ASP F 295 " (cutoff:3.500A) removed outlier: 3.521A pdb=" N GLN F 300 " --> pdb=" O ALA F 296 " (cutoff:3.500A) Processing helix chain 'F' and resid 303 through 319 removed outlier: 3.543A pdb=" N PHE F 308 " --> pdb=" O ALA F 304 " (cutoff:3.500A) removed outlier: 3.635A pdb=" N ALA F 311 " --> pdb=" O ILE F 307 " (cutoff:3.500A) Processing helix chain 'F' and resid 324 through 329 removed outlier: 4.140A pdb=" N TRP F 328 " --> pdb=" O THR F 324 " (cutoff:3.500A) removed outlier: 4.409A pdb=" N GLU F 329 " --> pdb=" O GLN F 325 " (cutoff:3.500A) No H-bonds generated for 'chain 'F' and resid 324 through 329' Processing helix chain 'F' and resid 365 through 383 removed outlier: 3.925A pdb=" N HIS F 378 " --> pdb=" O HIS F 374 " (cutoff:3.500A) removed outlier: 3.798A pdb=" N ILE F 379 " --> pdb=" O GLU F 375 " (cutoff:3.500A) Processing helix chain 'F' and resid 399 through 413 removed outlier: 4.040A pdb=" N GLY F 405 " --> pdb=" O HIS F 401 " (cutoff:3.500A) removed outlier: 3.559A pdb=" N GLU F 406 " --> pdb=" O GLU F 402 " (cutoff:3.500A) removed outlier: 4.405A pdb=" N MET F 408 " --> pdb=" O VAL F 404 " (cutoff:3.500A) removed outlier: 4.694A pdb=" N SER F 409 " --> pdb=" O GLY F 405 " (cutoff:3.500A) Processing helix chain 'F' and resid 414 through 422 Processing helix chain 'F' and resid 432 through 447 removed outlier: 4.082A pdb=" N PHE F 438 " --> pdb=" O THR F 434 " (cutoff:3.500A) removed outlier: 3.690A pdb=" N THR F 445 " --> pdb=" O LYS F 441 " (cutoff:3.500A) Processing helix chain 'F' and resid 449 through 466 Processing helix chain 'F' and resid 469 through 471 No H-bonds generated for 'chain 'F' and resid 469 through 471' Processing helix chain 'F' and resid 472 through 484 removed outlier: 3.703A pdb=" N LYS F 476 " --> pdb=" O GLN F 472 " (cutoff:3.500A) removed outlier: 3.576A pdb=" N GLU F 479 " --> pdb=" O LYS F 475 " (cutoff:3.500A) Processing helix chain 'F' and resid 498 through 502 Processing helix chain 'F' and resid 503 through 508 Processing helix chain 'F' and resid 512 through 523 removed outlier: 4.074A pdb=" N TYR F 516 " --> pdb=" O PHE F 512 " (cutoff:3.500A) removed outlier: 3.767A pdb=" N THR F 519 " --> pdb=" O TYR F 515 " (cutoff:3.500A) removed outlier: 3.744A pdb=" N LEU F 520 " --> pdb=" O TYR F 516 " (cutoff:3.500A) Processing helix chain 'F' and resid 523 through 533 Processing helix chain 'F' and resid 547 through 559 removed outlier: 3.966A pdb=" N MET F 557 " --> pdb=" O LYS F 553 " (cutoff:3.500A) Processing helix chain 'F' and resid 565 through 572 removed outlier: 3.632A pdb=" N GLU F 571 " --> pdb=" O THR F 567 " (cutoff:3.500A) Processing helix chain 'F' and resid 581 through 586 removed outlier: 4.132A pdb=" N LEU F 585 " --> pdb=" O VAL F 581 " (cutoff:3.500A) removed outlier: 3.733A pdb=" N ASN F 586 " --> pdb=" O ARG F 582 " (cutoff:3.500A) No H-bonds generated for 'chain 'F' and resid 581 through 586' Processing helix chain 'F' and resid 588 through 599 removed outlier: 3.636A pdb=" N TRP F 594 " --> pdb=" O PRO F 590 " (cutoff:3.500A) Processing helix chain 'G' and resid 338 through 343 Processing helix chain 'G' and resid 365 through 370 removed outlier: 3.675A pdb=" N TYR G 369 " --> pdb=" O TYR G 365 " (cutoff:3.500A) removed outlier: 3.776A pdb=" N ASN G 370 " --> pdb=" O SER G 366 " (cutoff:3.500A) No H-bonds generated for 'chain 'G' and resid 365 through 370' Processing helix chain 'G' and resid 385 through 390 removed outlier: 3.978A pdb=" N ASP G 389 " --> pdb=" O THR G 385 " (cutoff:3.500A) Processing helix chain 'G' and resid 405 through 410 removed outlier: 4.508A pdb=" N ARG G 408 " --> pdb=" O ASP G 405 " (cutoff:3.500A) Processing helix chain 'G' and resid 417 through 422 Processing helix chain 'G' and resid 502 through 505 Processing sheet with id=AA1, first strand: chain 'F' and resid 262 through 263 removed outlier: 5.664A pdb=" N LEU F 262 " --> pdb=" O VAL F 488 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA1 Processing sheet with id=AA2, first strand: chain 'F' and resid 347 through 352 removed outlier: 6.023A pdb=" N ASP F 355 " --> pdb=" O LEU F 351 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'G' and resid 354 through 356 removed outlier: 3.604A pdb=" N ASN G 354 " --> pdb=" O SER G 399 " (cutoff:3.500A) removed outlier: 3.508A pdb=" N VAL G 510 " --> pdb=" O PHE G 400 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'G' and resid 452 through 454 236 hydrogen bonds defined for protein. 675 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.83 Time building geometry restraints manager: 2.03 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.35: 2068 1.35 - 1.47: 1771 1.47 - 1.59: 2792 1.59 - 1.72: 0 1.72 - 1.84: 58 Bond restraints: 6689 Sorted by residual: bond pdb=" C1 NAG A 1 " pdb=" O5 NAG A 1 " ideal model delta sigma weight residual 1.406 1.447 -0.041 2.00e-02 2.50e+03 4.17e+00 bond pdb=" C1 NAG F 705 " pdb=" O5 NAG F 705 " ideal model delta sigma weight residual 1.406 1.446 -0.040 2.00e-02 2.50e+03 4.07e+00 bond pdb=" C1 NAG F 703 " pdb=" O5 NAG F 703 " ideal model delta sigma weight residual 1.406 1.445 -0.039 2.00e-02 2.50e+03 3.81e+00 bond pdb=" C1 NAG F 702 " pdb=" O5 NAG F 702 " ideal model delta sigma weight residual 1.406 1.445 -0.039 2.00e-02 2.50e+03 3.74e+00 bond pdb=" C1 NAG B 2 " pdb=" O5 NAG B 2 " ideal model delta sigma weight residual 1.406 1.444 -0.038 2.00e-02 2.50e+03 3.57e+00 ... (remaining 6684 not shown) Histogram of bond angle deviations from ideal: 0.00 - 2.29: 8914 2.29 - 4.57: 143 4.57 - 6.86: 19 6.86 - 9.14: 13 9.14 - 11.43: 2 Bond angle restraints: 9091 Sorted by residual: angle pdb=" CB MET F 152 " pdb=" CG MET F 152 " pdb=" SD MET F 152 " ideal model delta sigma weight residual 112.70 124.13 -11.43 3.00e+00 1.11e-01 1.45e+01 angle pdb=" CA MET F 557 " pdb=" CB MET F 557 " pdb=" CG MET F 557 " ideal model delta sigma weight residual 114.10 121.21 -7.11 2.00e+00 2.50e-01 1.27e+01 angle pdb=" CB MET F 557 " pdb=" CG MET F 557 " pdb=" SD MET F 557 " ideal model delta sigma weight residual 112.70 123.33 -10.63 3.00e+00 1.11e-01 1.26e+01 angle pdb=" N ASP F 499 " pdb=" CA ASP F 499 " pdb=" C ASP F 499 " ideal model delta sigma weight residual 109.81 117.51 -7.70 2.21e+00 2.05e-01 1.21e+01 angle pdb=" CB GLU F 238 " pdb=" CG GLU F 238 " pdb=" CD GLU F 238 " ideal model delta sigma weight residual 112.60 118.18 -5.58 1.70e+00 3.46e-01 1.08e+01 ... (remaining 9086 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.56: 3596 17.56 - 35.12: 373 35.12 - 52.68: 70 52.68 - 70.24: 17 70.24 - 87.80: 5 Dihedral angle restraints: 4061 sinusoidal: 1750 harmonic: 2311 Sorted by residual: dihedral pdb=" CB CYS G 336 " pdb=" SG CYS G 336 " pdb=" SG CYS G 361 " pdb=" CB CYS G 361 " ideal model delta sinusoidal sigma weight residual 93.00 26.69 66.31 1 1.00e+01 1.00e-02 5.74e+01 dihedral pdb=" CB CYS F 530 " pdb=" SG CYS F 530 " pdb=" SG CYS F 542 " pdb=" CB CYS F 542 " ideal model delta sinusoidal sigma weight residual -86.00 -122.53 36.53 1 1.00e+01 1.00e-02 1.89e+01 dihedral pdb=" SG CYS G 391 " pdb=" CB CYS G 525 " pdb=" SG CYS G 525 " pdb=" CA CYS G 525 " ideal model delta sinusoidal sigma weight residual 79.00 16.22 62.78 1 2.00e+01 2.50e-03 1.30e+01 ... (remaining 4058 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.036: 594 0.036 - 0.072: 267 0.072 - 0.108: 84 0.108 - 0.144: 25 0.144 - 0.180: 7 Chirality restraints: 977 Sorted by residual: chirality pdb=" CB VAL F 364 " pdb=" CA VAL F 364 " pdb=" CG1 VAL F 364 " pdb=" CG2 VAL F 364 " both_signs ideal model delta sigma weight residual False -2.63 -2.81 0.18 2.00e-01 2.50e+01 8.07e-01 chirality pdb=" CA ASP F 499 " pdb=" N ASP F 499 " pdb=" C ASP F 499 " pdb=" CB ASP F 499 " both_signs ideal model delta sigma weight residual False 2.51 2.33 0.18 2.00e-01 2.50e+01 7.82e-01 chirality pdb=" CA ILE F 468 " pdb=" N ILE F 468 " pdb=" C ILE F 468 " pdb=" CB ILE F 468 " both_signs ideal model delta sigma weight residual False 2.43 2.61 -0.18 2.00e-01 2.50e+01 7.81e-01 ... (remaining 974 not shown) Planarity restraints: 1174 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA VAL G 395 " -0.012 2.00e-02 2.50e+03 2.42e-02 5.85e+00 pdb=" C VAL G 395 " 0.042 2.00e-02 2.50e+03 pdb=" O VAL G 395 " -0.016 2.00e-02 2.50e+03 pdb=" N TYR G 396 " -0.014 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA VAL F 491 " -0.008 2.00e-02 2.50e+03 1.63e-02 2.67e+00 pdb=" C VAL F 491 " 0.028 2.00e-02 2.50e+03 pdb=" O VAL F 491 " -0.010 2.00e-02 2.50e+03 pdb=" N PRO F 492 " -0.010 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C VAL F 491 " -0.026 5.00e-02 4.00e+02 3.95e-02 2.49e+00 pdb=" N PRO F 492 " 0.068 5.00e-02 4.00e+02 pdb=" CA PRO F 492 " -0.020 5.00e-02 4.00e+02 pdb=" CD PRO F 492 " -0.022 5.00e-02 4.00e+02 ... (remaining 1171 not shown) Histogram of nonbonded interaction distances: 2.30 - 2.82: 2005 2.82 - 3.34: 5801 3.34 - 3.86: 10775 3.86 - 4.38: 12487 4.38 - 4.90: 21006 Nonbonded interactions: 52074 Sorted by model distance: nonbonded pdb=" O ALA F 501 " pdb=" OG SER F 507 " model vdw 2.298 3.040 nonbonded pdb=" O ILE F 544 " pdb=" OG SER F 547 " model vdw 2.300 3.040 nonbonded pdb=" O THR F 118 " pdb=" OG1 THR F 122 " model vdw 2.303 3.040 nonbonded pdb=" OG SER G 438 " pdb=" O PRO G 507 " model vdw 2.305 3.040 nonbonded pdb=" O GLU F 406 " pdb=" OG SER F 409 " model vdw 2.342 3.040 ... (remaining 52069 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 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 individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 8.320 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.000 Extract box with map and model: 0.000 Check model and map are aligned: 0.000 Set scattering table: 0.000 Process input model: 20.510 Find NCS groups from input model: 0.070 Set up NCS constraints: 0.020 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:2.980 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 31.900 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7863 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.004 0.041 6689 Z= 0.236 Angle : 0.741 11.430 9091 Z= 0.367 Chirality : 0.047 0.180 977 Planarity : 0.005 0.039 1167 Dihedral : 14.633 87.805 2558 Min Nonbonded Distance : 2.298 Molprobity Statistics. All-atom Clashscore : 3.86 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.85 % Favored : 95.15 % Rotamer: Outliers : 0.00 % Allowed : 0.29 % Favored : 99.71 % 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.57 (0.24), residues: 784 helix: -3.19 (0.17), residues: 348 sheet: -0.40 (0.79), residues: 41 loop : -1.70 (0.28), residues: 395 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP F 478 HIS 0.006 0.001 HIS F 540 PHE 0.009 0.001 PHE G 486 TYR 0.018 0.001 TYR G 396 ARG 0.006 0.001 ARG G 466 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 86 residues out of total 691 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 86 time to evaluate : 0.701 Fit side-chains revert: symmetry clash revert: symmetry clash REVERT: F 121 ASN cc_start: 0.8004 (m110) cc_final: 0.7788 (m110) outliers start: 0 outliers final: 0 residues processed: 86 average time/residue: 0.1697 time to fit residues: 20.3267 Evaluate side-chains 67 residues out of total 691 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 67 time to evaluate : 0.755 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 65 optimal weight: 2.9990 chunk 59 optimal weight: 0.9980 chunk 32 optimal weight: 2.9990 chunk 20 optimal weight: 0.5980 chunk 39 optimal weight: 5.9990 chunk 31 optimal weight: 5.9990 chunk 61 optimal weight: 1.9990 chunk 23 optimal weight: 1.9990 chunk 37 optimal weight: 0.0020 chunk 45 optimal weight: 1.9990 chunk 70 optimal weight: 0.0870 overall best weight: 0.7368 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: F 188 ASN F 373 HIS Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3631 r_free = 0.3631 target = 0.094614 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 36)----------------| | r_work = 0.3426 r_free = 0.3426 target = 0.083071 restraints weight = 17399.119| |-----------------------------------------------------------------------------| r_work (start): 0.3421 rms_B_bonded: 2.94 r_work: 0.3314 rms_B_bonded: 3.47 restraints_weight: 0.5000 r_work (final): 0.3314 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8134 moved from start: 0.1496 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.040 6689 Z= 0.201 Angle : 0.600 8.223 9091 Z= 0.308 Chirality : 0.043 0.148 977 Planarity : 0.004 0.031 1167 Dihedral : 4.914 45.720 1045 Min Nonbonded Distance : 2.438 Molprobity Statistics. All-atom Clashscore : 6.38 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.34 % Favored : 95.66 % Rotamer: Outliers : 0.72 % Allowed : 10.27 % Favored : 89.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.11 (0.27), residues: 784 helix: -1.21 (0.25), residues: 342 sheet: -0.33 (0.70), residues: 42 loop : -1.66 (0.28), residues: 400 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.001 TRP F 477 HIS 0.012 0.001 HIS F 540 PHE 0.011 0.001 PHE F 369 TYR 0.022 0.001 TYR F 183 ARG 0.003 0.000 ARG G 466 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 77 residues out of total 691 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 72 time to evaluate : 0.702 Fit side-chains revert: symmetry clash revert: symmetry clash REVERT: F 441 LYS cc_start: 0.8694 (tmtt) cc_final: 0.8452 (tmtt) outliers start: 5 outliers final: 2 residues processed: 74 average time/residue: 0.1677 time to fit residues: 17.3033 Evaluate side-chains 69 residues out of total 691 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 2 poor density : 67 time to evaluate : 0.680 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain F residue 82 MET Chi-restraints excluded: chain F residue 248 LEU Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 54 optimal weight: 0.9980 chunk 73 optimal weight: 0.5980 chunk 60 optimal weight: 0.6980 chunk 55 optimal weight: 0.0050 chunk 49 optimal weight: 2.9990 chunk 77 optimal weight: 3.9990 chunk 47 optimal weight: 4.9990 chunk 31 optimal weight: 0.0070 chunk 10 optimal weight: 1.9990 chunk 28 optimal weight: 0.9990 chunk 61 optimal weight: 0.9990 overall best weight: 0.4612 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3634 r_free = 0.3634 target = 0.094787 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 36)----------------| | r_work = 0.3430 r_free = 0.3430 target = 0.083271 restraints weight = 17520.698| |-----------------------------------------------------------------------------| r_work (start): 0.3426 rms_B_bonded: 2.96 r_work: 0.3319 rms_B_bonded: 3.46 restraints_weight: 0.5000 r_work (final): 0.3319 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8125 moved from start: 0.1828 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.034 6689 Z= 0.152 Angle : 0.528 6.563 9091 Z= 0.267 Chirality : 0.041 0.146 977 Planarity : 0.003 0.029 1167 Dihedral : 4.350 36.895 1045 Min Nonbonded Distance : 2.458 Molprobity Statistics. All-atom Clashscore : 6.77 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.70 % Favored : 96.30 % Rotamer: Outliers : 0.72 % Allowed : 13.17 % Favored : 86.11 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.39 (0.28), residues: 784 helix: -0.34 (0.27), residues: 344 sheet: -0.55 (0.69), residues: 42 loop : -1.47 (0.29), residues: 398 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP F 69 HIS 0.009 0.001 HIS F 540 PHE 0.010 0.001 PHE F 523 TYR 0.017 0.001 TYR F 183 ARG 0.004 0.000 ARG G 466 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 75 residues out of total 691 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 70 time to evaluate : 0.758 Fit side-chains revert: symmetry clash REVERT: F 332 MET cc_start: 0.6986 (ppp) cc_final: 0.6583 (ptt) REVERT: F 441 LYS cc_start: 0.8705 (tmtt) cc_final: 0.8467 (tmtt) outliers start: 5 outliers final: 3 residues processed: 73 average time/residue: 0.1729 time to fit residues: 17.5364 Evaluate side-chains 65 residues out of total 691 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 62 time to evaluate : 0.925 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain F residue 82 MET Chi-restraints excluded: chain F residue 276 THR Chi-restraints excluded: chain F residue 450 LEU Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 36 optimal weight: 3.9990 chunk 61 optimal weight: 0.6980 chunk 19 optimal weight: 0.6980 chunk 41 optimal weight: 2.9990 chunk 33 optimal weight: 0.6980 chunk 17 optimal weight: 3.9990 chunk 35 optimal weight: 1.9990 chunk 18 optimal weight: 0.4980 chunk 26 optimal weight: 0.2980 chunk 21 optimal weight: 0.9980 chunk 52 optimal weight: 6.9990 overall best weight: 0.5780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3625 r_free = 0.3625 target = 0.094195 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 29)----------------| | r_work = 0.3421 r_free = 0.3421 target = 0.082693 restraints weight = 17642.582| |-----------------------------------------------------------------------------| r_work (start): 0.3418 rms_B_bonded: 2.94 r_work: 0.3314 rms_B_bonded: 3.42 restraints_weight: 0.5000 r_work (final): 0.3314 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8140 moved from start: 0.2062 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 6689 Z= 0.169 Angle : 0.529 6.467 9091 Z= 0.269 Chirality : 0.042 0.147 977 Planarity : 0.003 0.028 1167 Dihedral : 4.183 28.004 1045 Min Nonbonded Distance : 2.478 Molprobity Statistics. All-atom Clashscore : 6.22 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.46 % Favored : 95.54 % Rotamer: Outliers : 0.58 % Allowed : 14.76 % Favored : 84.66 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.11 (0.29), residues: 784 helix: 0.00 (0.28), residues: 344 sheet: -0.46 (0.70), residues: 41 loop : -1.43 (0.29), residues: 399 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.001 TRP F 610 HIS 0.009 0.001 HIS F 540 PHE 0.011 0.001 PHE F 504 TYR 0.017 0.001 TYR F 183 ARG 0.004 0.000 ARG G 466 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 69 residues out of total 691 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 4 poor density : 65 time to evaluate : 0.711 Fit side-chains revert: symmetry clash revert: symmetry clash REVERT: F 332 MET cc_start: 0.6947 (ppp) cc_final: 0.6651 (ptt) REVERT: F 441 LYS cc_start: 0.8698 (tmtt) cc_final: 0.8461 (tmtt) outliers start: 4 outliers final: 3 residues processed: 66 average time/residue: 0.1636 time to fit residues: 15.4361 Evaluate side-chains 64 residues out of total 691 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 61 time to evaluate : 0.732 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain F residue 82 MET Chi-restraints excluded: chain F residue 276 THR Chi-restraints excluded: chain F residue 408 MET Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 78 random chunks: chunk 27 optimal weight: 0.8980 chunk 54 optimal weight: 0.5980 chunk 43 optimal weight: 2.9990 chunk 16 optimal weight: 0.5980 chunk 42 optimal weight: 0.9990 chunk 31 optimal weight: 0.8980 chunk 2 optimal weight: 0.7980 chunk 26 optimal weight: 0.9980 chunk 77 optimal weight: 3.9990 chunk 20 optimal weight: 0.6980 chunk 17 optimal weight: 0.6980 overall best weight: 0.6780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: F 175 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3624 r_free = 0.3624 target = 0.093662 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 44)----------------| | r_work = 0.3425 r_free = 0.3425 target = 0.082250 restraints weight = 17552.494| |-----------------------------------------------------------------------------| r_work (start): 0.3413 rms_B_bonded: 2.89 r_work: 0.3310 rms_B_bonded: 3.38 restraints_weight: 0.5000 r_work (final): 0.3310 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8161 moved from start: 0.2262 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 6689 Z= 0.179 Angle : 0.534 6.435 9091 Z= 0.271 Chirality : 0.042 0.145 977 Planarity : 0.003 0.029 1167 Dihedral : 4.096 21.867 1045 Min Nonbonded Distance : 2.471 Molprobity Statistics. All-atom Clashscore : 6.93 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.08 % Favored : 95.92 % Rotamer: Outliers : 1.30 % Allowed : 15.77 % Favored : 82.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: -0.96 (0.29), residues: 784 helix: 0.16 (0.28), residues: 352 sheet: -0.44 (0.72), residues: 41 loop : -1.40 (0.30), residues: 391 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.009 0.001 TRP F 610 HIS 0.009 0.001 HIS F 540 PHE 0.011 0.001 PHE F 504 TYR 0.016 0.001 TYR F 183 ARG 0.003 0.000 ARG F 460 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1568 Ramachandran restraints generated. 784 Oldfield, 0 Emsley, 784 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 75 residues out of total 691 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 9 poor density : 66 time to evaluate : 0.689 Fit side-chains revert: symmetry clash revert: symmetry clash REVERT: F 60 GLN cc_start: 0.7998 (tp40) cc_final: 0.7645 (tp40) REVERT: F 116 LEU cc_start: 0.7972 (tt) cc_final: 0.7686 (mt) REVERT: F 332 MET cc_start: 0.7153 (ppp) cc_final: 0.6857 (ptt) REVERT: G 430 THR cc_start: 0.8297 (p) cc_final: 0.8087 (p) outliers start: 9 outliers final: 6 residues processed: 72 average time/residue: 0.1840 time to fit residues: 18.4948 Evaluate side-chains 68 residues out of total 691 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 6 poor density : 62 time to evaluate : 0.760 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain F residue 82 MET Chi-restraints excluded: chain F residue 276 THR Chi-restraints excluded: chain F residue 408 MET Chi-restraints excluded: chain F residue 570 LEU Chi-restraints excluded: chain G residue 402 ILE Chi-restraints excluded: chain G residue 510 VAL Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 8, in run_program(real_space_refine.Program) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/iotbx/cli_parser.py", line 946, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/phenix/phenix/programs/real_space_refine.py", line 210, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 306, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 702, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 767, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 51, in __init__ gradients_method = gradients_method) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/mmtbx/refinement/real_space/weight.py", line 87, in __init__ rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 368, in refine rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 112, in refine states_collector = self.states_accumulator) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 174, in __init__ exception_handling_params=lbfgs_exception_handling_params) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 131, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 245, in compute_functional_and_gradients compute_gradients=True) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1525, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1427, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5592/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1308, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 55.3022 > 50: distance: 77 - 81: 7.403 distance: 81 - 82: 12.208 distance: 81 - 87: 5.099 distance: 82 - 83: 15.626 distance: 82 - 85: 17.040 distance: 83 - 84: 15.418 distance: 83 - 88: 21.689 distance: 85 - 86: 12.214 distance: 86 - 87: 16.395 distance: 88 - 89: 18.134 distance: 89 - 90: 4.507 distance: 89 - 92: 7.701 distance: 90 - 91: 6.749 distance: 90 - 96: 5.235 distance: 92 - 93: 15.060 distance: 92 - 94: 14.283 distance: 93 - 95: 9.191 distance: 96 - 97: 13.334 distance: 97 - 98: 21.323 distance: 98 - 99: 24.970 distance: 98 - 100: 14.089 distance: 100 - 101: 11.369 distance: 101 - 102: 6.934 distance: 101 - 104: 16.933 distance: 102 - 103: 7.238 distance: 102 - 106: 41.064 distance: 104 - 105: 28.664 distance: 106 - 107: 23.276 distance: 107 - 108: 13.967 distance: 107 - 110: 28.335 distance: 108 - 109: 6.667 distance: 108 - 114: 9.707 distance: 110 - 111: 20.271 distance: 111 - 112: 20.299 distance: 111 - 113: 29.237 distance: 114 - 115: 5.536 distance: 114 - 120: 14.040 distance: 115 - 116: 10.926 distance: 115 - 118: 12.051 distance: 116 - 117: 31.651 distance: 116 - 121: 24.426 distance: 118 - 119: 11.636 distance: 119 - 120: 17.613 distance: 121 - 122: 8.102 distance: 122 - 123: 7.002 distance: 122 - 125: 9.282 distance: 123 - 124: 4.100 distance: 123 - 126: 6.686 distance: 126 - 127: 14.287 distance: 127 - 128: 4.706 distance: 127 - 130: 12.603 distance: 128 - 129: 17.415 distance: 128 - 136: 5.163 distance: 130 - 131: 4.090 distance: 131 - 132: 9.994 distance: 131 - 133: 7.980 distance: 132 - 134: 9.981 distance: 133 - 135: 11.030 distance: 134 - 135: 6.061 distance: 136 - 137: 10.376 distance: 137 - 140: 8.841 distance: 138 - 139: 21.595 distance: 138 - 144: 15.615 distance: 140 - 141: 11.630 distance: 141 - 142: 26.283 distance: 141 - 143: 19.061 distance: 144 - 145: 8.585 distance: 145 - 146: 9.830 distance: 145 - 148: 10.189 distance: 146 - 147: 25.306 distance: 146 - 152: 6.896 distance: 148 - 149: 16.395 distance: 149 - 150: 12.338 distance: 149 - 151: 8.133