Starting phenix.real_space_refine on Sat Aug 23 11:06:33 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/7bwa_30231/08_2025/7bwa_30231.cif Found real_map, /net/cci-nas-00/data/ceres_data/7bwa_30231/08_2025/7bwa_30231.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=4.9 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/7bwa_30231/08_2025/7bwa_30231.map" default_real_map = "/net/cci-nas-00/data/ceres_data/7bwa_30231/08_2025/7bwa_30231.map" model { file = "/net/cci-nas-00/data/ceres_data/7bwa_30231/08_2025/7bwa_30231.cif" } default_model = "/net/cci-nas-00/data/ceres_data/7bwa_30231/08_2025/7bwa_30231.cif" } resolution = 4.9 write_initial_geo_file = False refinement { macro_cycles = 10 } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs 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.004 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 94 5.16 5 C 7420 2.51 5 N 1998 2.21 5 O 2158 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 72 residue(s): 0.01s Monomer Library directory: "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/chem_data/mon_lib" Total number of atoms: 11670 Number of models: 1 Model: "" Number of chains: 4 Chain: "A" Number of atoms: 5459 Number of conformers: 1 Conformer: "" Number of residues, atoms: 674, 5459 Classifications: {'peptide': 674} Link IDs: {'PCIS': 1, 'PTRANS': 33, 'TRANS': 639} Chain breaks: 6 Chain: "C" Number of atoms: 5459 Number of conformers: 1 Conformer: "" Number of residues, atoms: 674, 5459 Classifications: {'peptide': 674} Link IDs: {'PCIS': 1, 'PTRANS': 33, 'TRANS': 639} Chain breaks: 6 Chain: "D" Number of atoms: 376 Number of conformers: 1 Conformer: "" Number of residues, atoms: 48, 376 Classifications: {'peptide': 48} Link IDs: {'PTRANS': 1, 'TRANS': 46} Chain breaks: 1 Chain: "E" Number of atoms: 376 Number of conformers: 1 Conformer: "" Number of residues, atoms: 48, 376 Classifications: {'peptide': 48} Link IDs: {'PTRANS': 1, 'TRANS': 46} Chain breaks: 1 Time building chain proxies: 2.39, per 1000 atoms: 0.20 Number of scatterers: 11670 At special positions: 0 Unit cell: (91.79, 152.07, 132.89, 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 94 16.00 O 2158 8.00 N 1998 7.00 C 7420 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=32, symmetry=0 Simple disulfide: pdb=" SG CYS A 8 " - pdb=" SG CYS A 26 " distance=2.04 Simple disulfide: pdb=" SG CYS A 192 " - pdb=" SG CYS A 201 " distance=2.03 Simple disulfide: pdb=" SG CYS A 196 " - pdb=" SG CYS A 207 " distance=2.03 Simple disulfide: pdb=" SG CYS A 208 " - pdb=" SG CYS A 216 " distance=2.02 Simple disulfide: pdb=" SG CYS A 212 " - pdb=" SG CYS A 225 " distance=2.03 Simple disulfide: pdb=" SG CYS A 228 " - pdb=" SG CYS A 237 " distance=2.03 Simple disulfide: pdb=" SG CYS A 241 " - pdb=" SG CYS A 253 " distance=2.03 Simple disulfide: pdb=" SG CYS A 259 " - pdb=" SG CYS A 284 " distance=2.04 Simple disulfide: pdb=" SG CYS A 266 " - pdb=" SG CYS A 274 " distance=2.03 Simple disulfide: pdb=" SG CYS A 288 " - pdb=" SG CYS A 301 " distance=2.03 Simple disulfide: pdb=" SG CYS A 312 " - pdb=" SG CYS A 333 " distance=2.03 Simple disulfide: pdb=" SG CYS A 435 " - pdb=" SG CYS A 468 " distance=2.02 Simple disulfide: pdb=" SG CYS A 786 " - pdb=" SG CYS A 795 " distance=2.03 Simple disulfide: pdb=" SG CYS C 8 " - pdb=" SG CYS C 26 " distance=2.04 Simple disulfide: pdb=" SG CYS C 192 " - pdb=" SG CYS C 201 " distance=2.03 Simple disulfide: pdb=" SG CYS C 196 " - pdb=" SG CYS C 207 " distance=2.03 Simple disulfide: pdb=" SG CYS C 208 " - pdb=" SG CYS C 216 " distance=2.02 Simple disulfide: pdb=" SG CYS C 212 " - pdb=" SG CYS C 225 " distance=2.03 Simple disulfide: pdb=" SG CYS C 228 " - pdb=" SG CYS C 237 " distance=2.03 Simple disulfide: pdb=" SG CYS C 241 " - pdb=" SG CYS C 253 " distance=2.03 Simple disulfide: pdb=" SG CYS C 259 " - pdb=" SG CYS C 284 " distance=2.04 Simple disulfide: pdb=" SG CYS C 266 " - pdb=" SG CYS C 274 " distance=2.03 Simple disulfide: pdb=" SG CYS C 288 " - pdb=" SG CYS C 301 " distance=2.03 Simple disulfide: pdb=" SG CYS C 312 " - pdb=" SG CYS C 333 " distance=2.03 Simple disulfide: pdb=" SG CYS C 435 " - pdb=" SG CYS C 468 " distance=2.02 Simple disulfide: pdb=" SG CYS C 786 " - pdb=" SG CYS C 795 " distance=2.03 Simple disulfide: pdb=" SG CYS D 7 " - pdb=" SG CYS D 62 " distance=2.03 Simple disulfide: pdb=" SG CYS D 19 " - pdb=" SG CYS D 75 " distance=2.03 Simple disulfide: pdb=" SG CYS D 61 " - pdb=" SG CYS D 66 " distance=2.03 Simple disulfide: pdb=" SG CYS E 7 " - pdb=" SG CYS E 62 " distance=2.03 Simple disulfide: pdb=" SG CYS E 19 " - pdb=" SG CYS E 75 " distance=2.03 Simple disulfide: pdb=" SG CYS E 61 " - pdb=" SG CYS E 66 " 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 Time building additional restraints: 0.76 Conformation dependent library (CDL) restraints added in 396.2 milliseconds Enol-peptide restraints added in 715.3 nanoseconds 2816 Ramachandran restraints generated. 1408 Oldfield, 0 Emsley, 1408 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2740 Finding SS restraints... Secondary structure from input PDB file: 34 helices and 22 sheets defined 17.9% alpha, 17.5% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.39 Creating SS restraints... Processing helix chain 'A' and resid 16 through 23 removed outlier: 4.311A pdb=" N ARG A 19 " --> pdb=" O ASN A 16 " (cutoff:3.500A) removed outlier: 4.221A pdb=" N HIS A 21 " --> pdb=" O THR A 18 " (cutoff:3.500A) removed outlier: 3.741A pdb=" N GLU A 22 " --> pdb=" O ARG A 19 " (cutoff:3.500A) Processing helix chain 'A' and resid 42 through 47 removed outlier: 4.153A pdb=" N PHE A 46 " --> pdb=" O ARG A 42 " (cutoff:3.500A) Processing helix chain 'A' and resid 139 through 143 Processing helix chain 'A' and resid 256 through 265 Processing helix chain 'A' and resid 314 through 318 Processing helix chain 'A' and resid 323 through 328 Processing helix chain 'A' and resid 329 through 331 No H-bonds generated for 'chain 'A' and resid 329 through 331' Processing helix chain 'A' and resid 350 through 358 Processing helix chain 'A' and resid 413 through 417 Processing helix chain 'A' and resid 435 through 445 Processing helix chain 'A' and resid 458 through 462 Processing helix chain 'A' and resid 638 through 642 Processing helix chain 'A' and resid 695 through 714 removed outlier: 4.182A pdb=" N SER A 699 " --> pdb=" O GLU A 695 " (cutoff:3.500A) removed outlier: 3.953A pdb=" N SER A 700 " --> pdb=" O LEU A 696 " (cutoff:3.500A) removed outlier: 3.598A pdb=" N GLU A 706 " --> pdb=" O ARG A 702 " (cutoff:3.500A) removed outlier: 3.563A pdb=" N ASP A 707 " --> pdb=" O LYS A 703 " (cutoff:3.500A) removed outlier: 3.719A pdb=" N VAL A 712 " --> pdb=" O TYR A 708 " (cutoff:3.500A) removed outlier: 3.815A pdb=" N PHE A 714 " --> pdb=" O HIS A 710 " (cutoff:3.500A) Processing helix chain 'A' and resid 809 through 813 Processing helix chain 'C' and resid 16 through 23 removed outlier: 4.311A pdb=" N ARG C 19 " --> pdb=" O ASN C 16 " (cutoff:3.500A) removed outlier: 4.221A pdb=" N HIS C 21 " --> pdb=" O THR C 18 " (cutoff:3.500A) removed outlier: 3.741A pdb=" N GLU C 22 " --> pdb=" O ARG C 19 " (cutoff:3.500A) Processing helix chain 'C' and resid 42 through 47 removed outlier: 4.154A pdb=" N PHE C 46 " --> pdb=" O ARG C 42 " (cutoff:3.500A) Processing helix chain 'C' and resid 139 through 143 Processing helix chain 'C' and resid 256 through 265 Processing helix chain 'C' and resid 314 through 318 Processing helix chain 'C' and resid 323 through 328 Processing helix chain 'C' and resid 329 through 331 No H-bonds generated for 'chain 'C' and resid 329 through 331' Processing helix chain 'C' and resid 350 through 358 Processing helix chain 'C' and resid 413 through 417 Processing helix chain 'C' and resid 435 through 445 Processing helix chain 'C' and resid 458 through 462 Processing helix chain 'C' and resid 638 through 642 Processing helix chain 'C' and resid 695 through 714 removed outlier: 4.182A pdb=" N SER C 699 " --> pdb=" O GLU C 695 " (cutoff:3.500A) removed outlier: 3.953A pdb=" N SER C 700 " --> pdb=" O LEU C 696 " (cutoff:3.500A) removed outlier: 3.598A pdb=" N GLU C 706 " --> pdb=" O ARG C 702 " (cutoff:3.500A) removed outlier: 3.564A pdb=" N ASP C 707 " --> pdb=" O LYS C 703 " (cutoff:3.500A) removed outlier: 3.718A pdb=" N VAL C 712 " --> pdb=" O TYR C 708 " (cutoff:3.500A) removed outlier: 3.815A pdb=" N PHE C 714 " --> pdb=" O HIS C 710 " (cutoff:3.500A) Processing helix chain 'C' and resid 809 through 813 Processing helix chain 'D' and resid 7 through 19 Processing helix chain 'D' and resid 57 through 62 Processing helix chain 'D' and resid 67 through 73 Processing helix chain 'E' and resid 7 through 19 Processing helix chain 'E' and resid 57 through 62 Processing helix chain 'E' and resid 67 through 73 Processing sheet with id=AA1, first strand: chain 'A' and resid 11 through 13 removed outlier: 6.538A pdb=" N MET A 11 " --> pdb=" O GLN A 34 " (cutoff:3.500A) removed outlier: 8.200A pdb=" N LEU A 36 " --> pdb=" O MET A 11 " (cutoff:3.500A) removed outlier: 5.712A pdb=" N ILE A 13 " --> pdb=" O LEU A 36 " (cutoff:3.500A) removed outlier: 6.064A pdb=" N LEU A 93 " --> pdb=" O ARG A 118 " (cutoff:3.500A) removed outlier: 7.862A pdb=" N GLU A 120 " --> pdb=" O LEU A 93 " (cutoff:3.500A) removed outlier: 6.523A pdb=" N ILE A 95 " --> pdb=" O GLU A 120 " (cutoff:3.500A) removed outlier: 6.195A pdb=" N ASN A 122 " --> pdb=" O ILE A 95 " (cutoff:3.500A) removed outlier: 6.367A pdb=" N ASN A 148 " --> pdb=" O LYS A 121 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 28 through 29 removed outlier: 3.944A pdb=" N ILE A 29 " --> pdb=" O MET A 56 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 231 through 233 Processing sheet with id=AA4, first strand: chain 'A' and resid 252 through 255 removed outlier: 7.719A pdb=" N CYS A 284 " --> pdb=" O TYR A 245 " (cutoff:3.500A) removed outlier: 6.570A pdb=" N HIS A 247 " --> pdb=" O CYS A 284 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'A' and resid 292 through 293 Processing sheet with id=AA6, first strand: chain 'A' and resid 363 through 364 removed outlier: 3.798A pdb=" N ILE A 364 " --> pdb=" O LEU A 387 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'A' and resid 368 through 371 removed outlier: 5.331A pdb=" N TYR A 398 " --> pdb=" O LYS A 425 " (cutoff:3.500A) removed outlier: 8.328A pdb=" N PHE A 427 " --> pdb=" O TYR A 398 " (cutoff:3.500A) removed outlier: 6.547A pdb=" N PHE A 400 " --> pdb=" O PHE A 427 " (cutoff:3.500A) removed outlier: 7.899A pdb=" N HIS A 429 " --> pdb=" O PHE A 400 " (cutoff:3.500A) removed outlier: 6.360A pdb=" N ALA A 402 " --> pdb=" O HIS A 429 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'A' and resid 474 through 479 removed outlier: 5.475A pdb=" N PHE A 475 " --> pdb=" O GLU A 490 " (cutoff:3.500A) removed outlier: 5.568A pdb=" N GLU A 490 " --> pdb=" O PHE A 475 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'A' and resid 532 through 533 removed outlier: 3.602A pdb=" N VAL A 532 " --> pdb=" O LEU A 505 " (cutoff:3.500A) removed outlier: 4.004A pdb=" N GLY A 502 " --> pdb=" O LEU A 569 " (cutoff:3.500A) removed outlier: 3.552A pdb=" N LEU A 569 " --> pdb=" O GLY A 502 " (cutoff:3.500A) removed outlier: 4.275A pdb=" N GLN A 561 " --> pdb=" O ALA A 510 " (cutoff:3.500A) removed outlier: 5.306A pdb=" N VAL A 566 " --> pdb=" O ASP A 584 " (cutoff:3.500A) removed outlier: 6.169A pdb=" N ASP A 584 " --> pdb=" O VAL A 566 " (cutoff:3.500A) removed outlier: 7.038A pdb=" N THR A 568 " --> pdb=" O LYS A 582 " (cutoff:3.500A) removed outlier: 7.618A pdb=" N LYS A 582 " --> pdb=" O THR A 568 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'A' and resid 599 through 604 removed outlier: 5.977A pdb=" N ASP A 600 " --> pdb=" O LYS A 616 " (cutoff:3.500A) removed outlier: 3.917A pdb=" N LYS A 616 " --> pdb=" O ASP A 600 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'A' and resid 759 through 762 Processing sheet with id=AB3, first strand: chain 'C' and resid 11 through 13 removed outlier: 6.537A pdb=" N MET C 11 " --> pdb=" O GLN C 34 " (cutoff:3.500A) removed outlier: 8.200A pdb=" N LEU C 36 " --> pdb=" O MET C 11 " (cutoff:3.500A) removed outlier: 5.713A pdb=" N ILE C 13 " --> pdb=" O LEU C 36 " (cutoff:3.500A) removed outlier: 6.064A pdb=" N LEU C 93 " --> pdb=" O ARG C 118 " (cutoff:3.500A) removed outlier: 7.863A pdb=" N GLU C 120 " --> pdb=" O LEU C 93 " (cutoff:3.500A) removed outlier: 6.523A pdb=" N ILE C 95 " --> pdb=" O GLU C 120 " (cutoff:3.500A) removed outlier: 6.194A pdb=" N ASN C 122 " --> pdb=" O ILE C 95 " (cutoff:3.500A) removed outlier: 6.368A pdb=" N ASN C 148 " --> pdb=" O LYS C 121 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'C' and resid 28 through 29 removed outlier: 3.944A pdb=" N ILE C 29 " --> pdb=" O MET C 56 " (cutoff:3.500A) Processing sheet with id=AB5, first strand: chain 'C' and resid 231 through 233 Processing sheet with id=AB6, first strand: chain 'C' and resid 252 through 255 removed outlier: 7.720A pdb=" N CYS C 284 " --> pdb=" O TYR C 245 " (cutoff:3.500A) removed outlier: 6.570A pdb=" N HIS C 247 " --> pdb=" O CYS C 284 " (cutoff:3.500A) Processing sheet with id=AB7, first strand: chain 'C' and resid 292 through 293 Processing sheet with id=AB8, first strand: chain 'C' and resid 363 through 364 removed outlier: 3.797A pdb=" N ILE C 364 " --> pdb=" O LEU C 387 " (cutoff:3.500A) Processing sheet with id=AB9, first strand: chain 'C' and resid 368 through 371 removed outlier: 5.331A pdb=" N TYR C 398 " --> pdb=" O LYS C 425 " (cutoff:3.500A) removed outlier: 8.327A pdb=" N PHE C 427 " --> pdb=" O TYR C 398 " (cutoff:3.500A) removed outlier: 6.547A pdb=" N PHE C 400 " --> pdb=" O PHE C 427 " (cutoff:3.500A) removed outlier: 7.900A pdb=" N HIS C 429 " --> pdb=" O PHE C 400 " (cutoff:3.500A) removed outlier: 6.361A pdb=" N ALA C 402 " --> pdb=" O HIS C 429 " (cutoff:3.500A) Processing sheet with id=AC1, first strand: chain 'C' and resid 474 through 479 removed outlier: 5.475A pdb=" N PHE C 475 " --> pdb=" O GLU C 490 " (cutoff:3.500A) removed outlier: 5.567A pdb=" N GLU C 490 " --> pdb=" O PHE C 475 " (cutoff:3.500A) Processing sheet with id=AC2, first strand: chain 'C' and resid 532 through 533 removed outlier: 3.602A pdb=" N VAL C 532 " --> pdb=" O LEU C 505 " (cutoff:3.500A) removed outlier: 4.004A pdb=" N GLY C 502 " --> pdb=" O LEU C 569 " (cutoff:3.500A) removed outlier: 3.553A pdb=" N LEU C 569 " --> pdb=" O GLY C 502 " (cutoff:3.500A) removed outlier: 4.276A pdb=" N GLN C 561 " --> pdb=" O ALA C 510 " (cutoff:3.500A) removed outlier: 5.306A pdb=" N VAL C 566 " --> pdb=" O ASP C 584 " (cutoff:3.500A) removed outlier: 6.169A pdb=" N ASP C 584 " --> pdb=" O VAL C 566 " (cutoff:3.500A) removed outlier: 7.038A pdb=" N THR C 568 " --> pdb=" O LYS C 582 " (cutoff:3.500A) removed outlier: 7.619A pdb=" N LYS C 582 " --> pdb=" O THR C 568 " (cutoff:3.500A) Processing sheet with id=AC3, first strand: chain 'C' and resid 599 through 604 removed outlier: 5.976A pdb=" N ASP C 600 " --> pdb=" O LYS C 616 " (cutoff:3.500A) removed outlier: 3.918A pdb=" N LYS C 616 " --> pdb=" O ASP C 600 " (cutoff:3.500A) Processing sheet with id=AC4, first strand: chain 'C' and resid 759 through 762 262 hydrogen bonds defined for protein. 636 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.85 Time building geometry restraints manager: 1.02 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: 3799 1.34 - 1.46: 2701 1.46 - 1.58: 5330 1.58 - 1.70: 0 1.70 - 1.82: 114 Bond restraints: 11944 Sorted by residual: bond pdb=" N TYR C 579 " pdb=" CA TYR C 579 " ideal model delta sigma weight residual 1.453 1.487 -0.034 1.31e-02 5.83e+03 6.90e+00 bond pdb=" CA ALA A 581 " pdb=" CB ALA A 581 " ideal model delta sigma weight residual 1.538 1.506 0.032 1.22e-02 6.72e+03 6.89e+00 bond pdb=" N TYR A 579 " pdb=" CA TYR A 579 " ideal model delta sigma weight residual 1.453 1.486 -0.033 1.31e-02 5.83e+03 6.30e+00 bond pdb=" N ARG C 576 " pdb=" CA ARG C 576 " ideal model delta sigma weight residual 1.456 1.485 -0.029 1.16e-02 7.43e+03 6.17e+00 bond pdb=" N ARG A 577 " pdb=" CA ARG A 577 " ideal model delta sigma weight residual 1.456 1.486 -0.030 1.23e-02 6.61e+03 6.09e+00 ... (remaining 11939 not shown) Histogram of bond angle deviations from ideal: 0.00 - 1.46: 15121 1.46 - 2.92: 890 2.92 - 4.38: 128 4.38 - 5.83: 27 5.83 - 7.29: 8 Bond angle restraints: 16174 Sorted by residual: angle pdb=" N GLU A 575 " pdb=" CA GLU A 575 " pdb=" C GLU A 575 " ideal model delta sigma weight residual 113.28 108.49 4.79 1.22e+00 6.72e-01 1.54e+01 angle pdb=" N GLU C 575 " pdb=" CA GLU C 575 " pdb=" C GLU C 575 " ideal model delta sigma weight residual 113.20 108.64 4.56 1.21e+00 6.83e-01 1.42e+01 angle pdb=" CA GLY C 580 " pdb=" C GLY C 580 " pdb=" O GLY C 580 " ideal model delta sigma weight residual 122.24 118.32 3.92 1.08e+00 8.57e-01 1.32e+01 angle pdb=" CA ARG C 577 " pdb=" C ARG C 577 " pdb=" O ARG C 577 " ideal model delta sigma weight residual 121.87 117.95 3.92 1.10e+00 8.26e-01 1.27e+01 angle pdb=" CA ARG C 576 " pdb=" C ARG C 576 " pdb=" O ARG C 576 " ideal model delta sigma weight residual 121.38 118.19 3.19 1.06e+00 8.90e-01 9.06e+00 ... (remaining 16169 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.95: 6337 17.95 - 35.89: 732 35.89 - 53.84: 153 53.84 - 71.79: 22 71.79 - 89.74: 8 Dihedral angle restraints: 7252 sinusoidal: 3024 harmonic: 4228 Sorted by residual: dihedral pdb=" CB CYS D 7 " pdb=" SG CYS D 7 " pdb=" SG CYS D 62 " pdb=" CB CYS D 62 " ideal model delta sinusoidal sigma weight residual -86.00 -152.59 66.59 1 1.00e+01 1.00e-02 5.79e+01 dihedral pdb=" CB CYS E 7 " pdb=" SG CYS E 7 " pdb=" SG CYS E 62 " pdb=" CB CYS E 62 " ideal model delta sinusoidal sigma weight residual -86.00 -152.58 66.58 1 1.00e+01 1.00e-02 5.79e+01 dihedral pdb=" CB CYS A 312 " pdb=" SG CYS A 312 " pdb=" SG CYS A 333 " pdb=" CB CYS A 333 " ideal model delta sinusoidal sigma weight residual 93.00 148.34 -55.34 1 1.00e+01 1.00e-02 4.14e+01 ... (remaining 7249 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.036: 1038 0.036 - 0.072: 486 0.072 - 0.109: 182 0.109 - 0.145: 48 0.145 - 0.181: 16 Chirality restraints: 1770 Sorted by residual: chirality pdb=" CA ARG A 577 " pdb=" N ARG A 577 " pdb=" C ARG A 577 " pdb=" CB ARG A 577 " both_signs ideal model delta sigma weight residual False 2.51 2.69 -0.18 2.00e-01 2.50e+01 8.18e-01 chirality pdb=" CA ARG C 577 " pdb=" N ARG C 577 " pdb=" C ARG C 577 " pdb=" CB ARG C 577 " both_signs ideal model delta sigma weight residual False 2.51 2.68 -0.17 2.00e-01 2.50e+01 6.95e-01 chirality pdb=" CA ILE C 145 " pdb=" N ILE C 145 " pdb=" C ILE C 145 " pdb=" CB ILE C 145 " both_signs ideal model delta sigma weight residual False 2.43 2.59 -0.16 2.00e-01 2.50e+01 6.30e-01 ... (remaining 1767 not shown) Planarity restraints: 2076 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA ASP A 574 " 0.010 2.00e-02 2.50e+03 2.08e-02 4.35e+00 pdb=" C ASP A 574 " -0.036 2.00e-02 2.50e+03 pdb=" O ASP A 574 " 0.014 2.00e-02 2.50e+03 pdb=" N GLU A 575 " 0.012 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA ASP C 574 " 0.008 2.00e-02 2.50e+03 1.67e-02 2.79e+00 pdb=" C ASP C 574 " -0.029 2.00e-02 2.50e+03 pdb=" O ASP C 574 " 0.011 2.00e-02 2.50e+03 pdb=" N GLU C 575 " 0.010 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C ARG C 42 " 0.022 5.00e-02 4.00e+02 3.35e-02 1.80e+00 pdb=" N PRO C 43 " -0.058 5.00e-02 4.00e+02 pdb=" CA PRO C 43 " 0.017 5.00e-02 4.00e+02 pdb=" CD PRO C 43 " 0.019 5.00e-02 4.00e+02 ... (remaining 2073 not shown) Histogram of nonbonded interaction distances: 1.97 - 2.56: 258 2.56 - 3.14: 10632 3.14 - 3.73: 20249 3.73 - 4.31: 27501 4.31 - 4.90: 40931 Nonbonded interactions: 99571 Sorted by model distance: nonbonded pdb=" O GLU C 469 " pdb=" NH2 ARG C 577 " model vdw 1.970 3.120 nonbonded pdb=" OD1 ASP C 221 " pdb=" OG1 THR C 223 " model vdw 2.122 3.040 nonbonded pdb=" OD1 ASP A 221 " pdb=" OG1 THR A 223 " model vdw 2.123 3.040 nonbonded pdb=" OH TYR C 628 " pdb=" OE2 GLU C 766 " model vdw 2.138 3.040 nonbonded pdb=" OH TYR A 628 " pdb=" OE2 GLU A 766 " model vdw 2.139 3.040 ... (remaining 99566 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.02 Found NCS groups: ncs_group { reference = chain 'A' selection = chain 'C' } ncs_group { reference = chain 'D' selection = chain 'E' } Set up NCS constraints Number of NCS constrained groups: 2 pdb_interpretation.ncs_group { reference = chain 'A' selection = chain 'C' } pdb_interpretation.ncs_group { reference = chain 'D' selection = chain 'E' } Set refine NCS operators NCS operators will be refined. 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: 0.240 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.110 Check model and map are aligned: 0.020 Set scattering table: 0.020 Process input model: 10.010 Find NCS groups from input model: 0.090 Set up NCS constraints: 0.010 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.010 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:0.970 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 11.490 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6447 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.034 11976 Z= 0.221 Angle : 0.774 7.294 16238 Z= 0.439 Chirality : 0.049 0.181 1770 Planarity : 0.004 0.034 2076 Dihedral : 15.314 89.735 4416 Min Nonbonded Distance : 1.970 Molprobity Statistics. All-atom Clashscore : 41.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 13.35 % Favored : 86.65 % Rotamer: Outliers : 0.00 % Allowed : 0.91 % Favored : 99.09 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.94 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z values with (uncertainties): Interpretation: poor |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores below are scaled independently, so they aren not related in a simple way. whole: -4.15 (0.20), residues: 1408 helix: -2.50 (0.34), residues: 166 sheet: -3.54 (0.32), residues: 178 loop : -2.89 (0.18), residues: 1064 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.004 0.001 ARG C 702 TYR 0.018 0.002 TYR C 91 PHE 0.019 0.002 PHE C 89 TRP 0.014 0.002 TRP C 632 HIS 0.007 0.002 HIS A 263 Details of bonding type rmsd covalent geometry : bond 0.00425 (11944) covalent geometry : angle 0.77212 (16174) SS BOND : bond 0.00509 ( 32) SS BOND : angle 1.09026 ( 64) hydrogen bonds : bond 0.20182 ( 262) hydrogen bonds : angle 10.08491 ( 636) *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2816 Ramachandran restraints generated. 1408 Oldfield, 0 Emsley, 1408 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2816 Ramachandran restraints generated. 1408 Oldfield, 0 Emsley, 1408 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Update NCS operators ******************** ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 117 residues out of total 618 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 117 time to evaluate : 0.141 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: symmetry clash REVERT: A 109 LEU cc_start: 0.8593 (mm) cc_final: 0.8379 (mm) REVERT: A 250 ASP cc_start: 0.7958 (m-30) cc_final: 0.7721 (t70) REVERT: A 368 LEU cc_start: 0.8239 (mm) cc_final: 0.7935 (mm) REVERT: A 442 MET cc_start: 0.8113 (mmt) cc_final: 0.7504 (mmt) REVERT: A 477 TYR cc_start: 0.8328 (m-10) cc_final: 0.7880 (m-80) REVERT: A 575 GLU cc_start: 0.5698 (pp20) cc_final: 0.5451 (pp20) REVERT: A 706 GLU cc_start: 0.8874 (mp0) cc_final: 0.8410 (mp0) outliers start: 0 outliers final: 0 residues processed: 117 average time/residue: 0.0864 time to fit residues: 12.3352 Evaluate side-chains 65 residues out of total 618 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 65 time to evaluate : 0.134 Evaluate side-chains 11 residues out of total 43 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 11 time to evaluate : 0.009 Fit side-chains REVERT: D 22 ARG cc_start: 0.7851 (mpp80) cc_final: 0.7118 (mmm160) outliers start: 0 outliers final: 0 residues processed: 11 average time/residue: 0.0495 time to fit residues: 0.5797 Evaluate side-chains 6 residues out of total 43 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 6 time to evaluate : 0.009 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 142 random chunks: chunk 98 optimal weight: 0.3980 chunk 107 optimal weight: 6.9990 chunk 10 optimal weight: 0.9980 chunk 66 optimal weight: 7.9990 chunk 130 optimal weight: 0.0040 chunk 124 optimal weight: 2.9990 chunk 103 optimal weight: 5.9990 chunk 77 optimal weight: 0.0050 chunk 122 optimal weight: 0.0870 chunk 91 optimal weight: 10.0000 chunk 55 optimal weight: 0.5980 overall best weight: 0.2184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 25 ASN ** A 199 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 209 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 282 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 406 GLN A 431 ASN ** A 610 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 3 Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Update NCS operators ******************** ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4050 r_free = 0.4050 target = 0.091584 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 48)----------------| | r_work = 0.3628 r_free = 0.3628 target = 0.073643 restraints weight = 36415.017| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 41)----------------| | r_work = 0.3675 r_free = 0.3675 target = 0.075711 restraints weight = 26351.065| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 26)----------------| | r_work = 0.3707 r_free = 0.3707 target = 0.077149 restraints weight = 20963.998| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 4; iterations = 33)----------------| | r_work = 0.3730 r_free = 0.3730 target = 0.078159 restraints weight = 17714.041| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 5; iterations = 28)----------------| | r_work = 0.3742 r_free = 0.3742 target = 0.078722 restraints weight = 15729.351| |-----------------------------------------------------------------------------| r_work (final): 0.3716 |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3446 r_free = 0.3446 target = 0.080663 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 20)----------------| | r_work = 0.3272 r_free = 0.3272 target = 0.071553 restraints weight = 2226.351| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 30)----------------| | r_work = 0.3312 r_free = 0.3312 target = 0.073409 restraints weight = 1550.158| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 27)----------------| | r_work = 0.3342 r_free = 0.3342 target = 0.074826 restraints weight = 1124.110| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 4; iterations = 28)----------------| | r_work = 0.3367 r_free = 0.3367 target = 0.075970 restraints weight = 839.162| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 5; iterations = 44)----------------| | r_work = 0.3383 r_free = 0.3383 target = 0.076751 restraints weight = 631.507| |-----------------------------------------------------------------------------| r_work (final): 0.3391 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6169 moved from start: 0.2586 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.049 11976 Z= 0.161 Angle : 0.778 8.515 16238 Z= 0.405 Chirality : 0.049 0.181 1770 Planarity : 0.005 0.060 2076 Dihedral : 6.543 21.575 1566 Min Nonbonded Distance : 2.521 Molprobity Statistics. All-atom Clashscore : 22.98 Ramachandran Plot: Outliers : 0.14 % Allowed : 11.36 % Favored : 88.49 % Rotamer: Outliers : 0.00 % Allowed : 8.47 % Favored : 91.53 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.94 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z values with (uncertainties): Interpretation: poor |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores below are scaled independently, so they aren not related in a simple way. whole: -3.64 (0.21), residues: 1408 helix: -2.24 (0.36), residues: 154 sheet: -3.19 (0.33), residues: 178 loop : -2.50 (0.19), residues: 1076 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.006 0.001 ARG C 345 TYR 0.025 0.002 TYR A 587 PHE 0.070 0.002 PHE A 701 TRP 0.011 0.002 TRP C 133 HIS 0.029 0.002 HIS E 10 Details of bonding type rmsd covalent geometry : bond 0.00346 (11944) covalent geometry : angle 0.77805 (16174) SS BOND : bond 0.00474 ( 32) SS BOND : angle 0.78715 ( 64) hydrogen bonds : bond 0.05846 ( 262) hydrogen bonds : angle 8.70098 ( 636) *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2816 Ramachandran restraints generated. 1408 Oldfield, 0 Emsley, 1408 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2816 Ramachandran restraints generated. 1408 Oldfield, 0 Emsley, 1408 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Update NCS operators ******************** ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 87 residues out of total 618 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 87 time to evaluate : 0.136 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 93 LEU cc_start: 0.6889 (tp) cc_final: 0.6684 (tp) REVERT: A 101 LEU cc_start: 0.9235 (tp) cc_final: 0.8909 (tp) REVERT: A 263 HIS cc_start: 0.8597 (t70) cc_final: 0.8110 (t70) REVERT: A 368 LEU cc_start: 0.8063 (mm) cc_final: 0.7736 (mm) REVERT: A 419 LEU cc_start: 0.8341 (tt) cc_final: 0.8064 (mt) REVERT: A 442 MET cc_start: 0.8549 (mmt) cc_final: 0.7906 (mmt) REVERT: A 504 MET cc_start: 0.9044 (mmp) cc_final: 0.8615 (mmm) REVERT: A 579 TYR cc_start: 0.8645 (m-80) cc_final: 0.8293 (m-80) REVERT: A 706 GLU cc_start: 0.8847 (mp0) cc_final: 0.8217 (pm20) outliers start: 0 outliers final: 0 residues processed: 87 average time/residue: 0.0509 time to fit residues: 5.9799 Evaluate side-chains 67 residues out of total 618 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 67 time to evaluate : 0.165 Evaluate side-chains 9 residues out of total 43 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 9 time to evaluate : 0.009 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 9 average time/residue: 0.0532 time to fit residues: 0.5085 Evaluate side-chains 6 residues out of total 43 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 6 time to evaluate : 0.016 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 142 random chunks: chunk 39 optimal weight: 4.9990 chunk 20 optimal weight: 10.0000 chunk 35 optimal weight: 7.9990 chunk 38 optimal weight: 3.9990 chunk 140 optimal weight: 0.0060 chunk 47 optimal weight: 3.9990 chunk 104 optimal weight: 0.2980 chunk 99 optimal weight: 0.1980 chunk 4 optimal weight: 0.0050 chunk 94 optimal weight: 0.0970 chunk 73 optimal weight: 0.0170 overall best weight: 0.0646 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 152 ASN A 199 HIS ** A 209 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 282 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 410 GLN A 417 HIS ** A 610 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 4 Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Update NCS operators ******************** ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4096 r_free = 0.4096 target = 0.093396 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 58)----------------| | r_work = 0.3711 r_free = 0.3711 target = 0.077075 restraints weight = 37684.869| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 49)----------------| | r_work = 0.3745 r_free = 0.3745 target = 0.078549 restraints weight = 29498.640| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 41)----------------| | r_work = 0.3767 r_free = 0.3767 target = 0.079547 restraints weight = 24548.685| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 4; iterations = 39)----------------| | r_work = 0.3784 r_free = 0.3784 target = 0.080315 restraints weight = 21537.081| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 5; iterations = 28)----------------| | r_work = 0.3794 r_free = 0.3794 target = 0.080790 restraints weight = 19504.661| |-----------------------------------------------------------------------------| r_work (final): 0.3753 |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3486 r_free = 0.3486 target = 0.082219 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 32)----------------| | r_work = 0.3310 r_free = 0.3310 target = 0.073553 restraints weight = 2109.463| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 28)----------------| | r_work = 0.3349 r_free = 0.3349 target = 0.075542 restraints weight = 1399.508| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 33)----------------| | r_work = 0.3381 r_free = 0.3381 target = 0.076971 restraints weight = 980.770| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 4; iterations = 29)----------------| | r_work = 0.3410 r_free = 0.3410 target = 0.078052 restraints weight = 709.433| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 5; iterations = 29)----------------| | r_work = 0.3429 r_free = 0.3429 target = 0.078972 restraints weight = 509.389| |-----------------------------------------------------------------------------| r_work (final): 0.3435 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6074 moved from start: 0.3713 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.046 11976 Z= 0.144 Angle : 0.736 11.446 16238 Z= 0.383 Chirality : 0.047 0.166 1770 Planarity : 0.005 0.045 2076 Dihedral : 6.121 23.570 1566 Min Nonbonded Distance : 2.461 Molprobity Statistics. All-atom Clashscore : 18.33 Ramachandran Plot: Outliers : 0.00 % Allowed : 10.09 % Favored : 89.91 % Rotamer: Outliers : 0.00 % Allowed : 5.45 % Favored : 94.55 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.94 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z values with (uncertainties): Interpretation: poor |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores below are scaled independently, so they aren not related in a simple way. whole: -3.53 (0.21), residues: 1408 helix: -2.02 (0.36), residues: 158 sheet: -2.86 (0.33), residues: 188 loop : -2.49 (0.18), residues: 1062 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.006 0.001 ARG A 702 TYR 0.018 0.002 TYR A 628 PHE 0.038 0.002 PHE A 701 TRP 0.017 0.002 TRP A 133 HIS 0.006 0.001 HIS D 10 Details of bonding type rmsd covalent geometry : bond 0.00298 (11944) covalent geometry : angle 0.73477 (16174) SS BOND : bond 0.00290 ( 32) SS BOND : angle 0.95224 ( 64) hydrogen bonds : bond 0.04705 ( 262) hydrogen bonds : angle 8.05916 ( 636) *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2816 Ramachandran restraints generated. 1408 Oldfield, 0 Emsley, 1408 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2816 Ramachandran restraints generated. 1408 Oldfield, 0 Emsley, 1408 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Update NCS operators ******************** ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 95 residues out of total 618 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 95 time to evaluate : 0.208 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 36 LEU cc_start: 0.7967 (pt) cc_final: 0.6908 (pt) REVERT: A 38 MET cc_start: 0.4238 (ttt) cc_final: 0.3343 (tpp) REVERT: A 93 LEU cc_start: 0.7176 (tp) cc_final: 0.6880 (tp) REVERT: A 109 LEU cc_start: 0.8489 (mm) cc_final: 0.8249 (mm) REVERT: A 368 LEU cc_start: 0.8602 (mm) cc_final: 0.8252 (mm) REVERT: A 442 MET cc_start: 0.8516 (mmt) cc_final: 0.7994 (mmt) REVERT: A 488 ARG cc_start: 0.7193 (mmm-85) cc_final: 0.6985 (mtt180) REVERT: A 504 MET cc_start: 0.8998 (mmp) cc_final: 0.8418 (mmm) REVERT: A 579 TYR cc_start: 0.8680 (m-80) cc_final: 0.8016 (m-80) REVERT: A 710 HIS cc_start: 0.9433 (m90) cc_final: 0.9028 (m-70) outliers start: 0 outliers final: 0 residues processed: 95 average time/residue: 0.0737 time to fit residues: 9.2252 Evaluate side-chains 64 residues out of total 618 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 64 time to evaluate : 0.218 Evaluate side-chains 8 residues out of total 43 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 8 time to evaluate : 0.015 Fit side-chains revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 8 average time/residue: 0.0774 time to fit residues: 0.6584 Evaluate side-chains 7 residues out of total 43 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 7 time to evaluate : 0.009 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/phenix/command_line/real_space_refine.py", line 9, in run_program(real_space_refine.Program) File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/iotbx/cli_parser.py", line 999, in run_program task.run() File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/phenix/programs/real_space_refine.py", line 206, in run self.rsr_result = rsr_wrappers.run_real_space_refine( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/phenix/refinement/rsr/wrappers.py", line 52, in __init__ self._refined = macro_cycle_real_space.run( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/phenix/refinement/macro_cycle_real_space.py", line 306, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/phenix/refinement/macro_cycle_real_space.py", line 704, in refine_xyz self.minimization_ncs() File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/phenix/refinement/macro_cycle_real_space.py", line 794, in minimization_ncs weight = real_space_weight.run( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/mmtbx/refinement/real_space/weight.py", line 84, in __init__ ro.refine( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 394, in refine real_space_result = refinery( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 233, in __init__ refiner.refine( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 126, in refine minimized = mmtbx.refinement.minimization_ncs_constraints.lbfgs( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/mmtbx/refinement/minimization_ncs_constraints.py", line 424, in __init__ minimizer = scitbx.lbfgs.run( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/scitbx/lbfgs/__init__.py", line 270, in run return run_c_plus_plus( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/scitbx/lbfgs/__init__.py", line 131, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/mmtbx/refinement/minimization_ncs_constraints.py", line 456, in compute_functional_and_gradients t,g = self.target_and_grads_object.target_and_gradients( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/mmtbx/refinement/minimization_ncs_constraints.py", line 242, in target_and_gradients t_restraints, g_restraints = restraints_target_and_grads( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/mmtbx/refinement/minimization_ncs_constraints.py", line 110, in restraints_target_and_grads ef = restraints_manager.energies_sites( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/cctbx/geometry_restraints/manager.py", line 1588, in energies_sites pair_proxies = self.pair_proxies( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/cctbx/geometry_restraints/manager.py", line 1491, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/site-packages/cctbx/geometry_restraints/manager.py", line 1372, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 62.2416 > 50: distance: 4 - 5: 10.403 distance: 5 - 6: 5.865 distance: 5 - 8: 4.964 distance: 6 - 7: 24.429 distance: 6 - 12: 27.933 distance: 8 - 9: 15.848 distance: 9 - 10: 17.359 distance: 9 - 11: 13.905 distance: 12 - 13: 12.074 distance: 13 - 14: 26.900 distance: 13 - 16: 22.632 distance: 14 - 15: 3.209 distance: 14 - 20: 28.599 distance: 16 - 17: 5.778 distance: 17 - 19: 8.011 distance: 20 - 21: 26.642 distance: 21 - 22: 15.571 distance: 21 - 24: 15.532 distance: 22 - 23: 11.227 distance: 22 - 28: 40.132 distance: 23 - 46: 13.771 distance: 24 - 25: 17.686 distance: 25 - 26: 25.839 distance: 25 - 27: 17.005 distance: 28 - 29: 15.081 distance: 29 - 30: 11.334 distance: 29 - 32: 8.019 distance: 30 - 31: 17.764 distance: 30 - 35: 11.680 distance: 32 - 33: 17.704 distance: 32 - 34: 8.397 distance: 35 - 36: 7.896 distance: 36 - 37: 9.012 distance: 36 - 39: 12.542 distance: 37 - 38: 4.902 distance: 37 - 46: 22.218 distance: 39 - 40: 24.477 distance: 40 - 41: 13.278 distance: 42 - 43: 3.182 distance: 46 - 47: 9.866 distance: 47 - 48: 25.036 distance: 47 - 50: 24.766 distance: 48 - 49: 14.382 distance: 48 - 54: 26.596 distance: 49 - 73: 11.617 distance: 50 - 51: 19.467 distance: 51 - 52: 17.506 distance: 51 - 53: 11.517 distance: 54 - 55: 15.361 distance: 55 - 56: 14.104 distance: 55 - 58: 8.699 distance: 56 - 57: 37.177 distance: 56 - 64: 10.504 distance: 58 - 59: 13.268 distance: 59 - 60: 19.542 distance: 59 - 61: 22.336 distance: 60 - 62: 18.323 distance: 61 - 63: 9.019 distance: 62 - 63: 20.134 distance: 64 - 65: 28.377 distance: 65 - 66: 28.530 distance: 65 - 68: 12.815 distance: 66 - 67: 44.096 distance: 66 - 73: 19.764 distance: 68 - 69: 31.203 distance: 69 - 70: 15.092 distance: 70 - 71: 23.876 distance: 70 - 72: 17.389