Starting phenix.real_space_refine on Thu Jul 2 20:32:19 2026 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, restraint, /net/cci-nas-00/data/ceres_data/7bzf_30252/07_2026/7bzf_30252.cif Found real_map, /net/cci-nas-00/data/ceres_data/7bzf_30252/07_2026/7bzf_30252.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.26 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { model { file = "/net/cci-nas-00/data/ceres_data/7bzf_30252/07_2026/7bzf_30252.cif" } default_model = "/net/cci-nas-00/data/ceres_data/7bzf_30252/07_2026/7bzf_30252.cif" real_map_files = "/net/cci-nas-00/data/ceres_data/7bzf_30252/07_2026/7bzf_30252.map" default_real_map = "/net/cci-nas-00/data/ceres_data/7bzf_30252/07_2026/7bzf_30252.map" restraint_files = "/net/cci-nas-00/data/ceres_data/7bzf_30252/07_2026/7bzf_30252.cif" default_restraint = "/net/cci-nas-00/data/ceres_data/7bzf_30252/07_2026/7bzf_30252.cif" } resolution = 3.26 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.000 sd= 0.059 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: 6 Type Number sf(0) Gaussians Zn 2 6.06 5 P 19 5.49 5 S 70 5.16 5 C 6245 2.51 5 N 1669 2.21 5 O 1927 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 30 residue(s): 0.01s Monomer Library directory: "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/chem_data/mon_lib" Total number of atoms: 9932 Number of models: 1 Model: "" Number of chains: 7 Chain: "A" Number of atoms: 7333 Number of conformers: 1 Conformer: "" Number of residues, atoms: 914, 7333 Classifications: {'peptide': 914} Incomplete info: {'truncation_to_alanine': 6} Link IDs: {'PCIS': 1, 'PTRANS': 29, 'TRANS': 883} Chain breaks: 1 Unresolved non-hydrogen bonds: 33 Unresolved non-hydrogen angles: 42 Unresolved non-hydrogen dihedrals: 33 Planarities with less than four sites: {'ASP:plan': 1, 'GLU:plan': 1, 'PHE:plan': 1, 'TRP:plan': 1} Unresolved non-hydrogen planarities: 22 Chain: "B" Number of atoms: 873 Number of conformers: 1 Conformer: "" Number of residues, atoms: 114, 873 Classifications: {'peptide': 114} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 5, 'TRANS': 108} Unresolved non-hydrogen bonds: 11 Unresolved non-hydrogen angles: 16 Unresolved non-hydrogen dihedrals: 13 Planarities with less than four sites: {'TRP:plan': 1} Unresolved non-hydrogen planarities: 9 Chain: "C" Number of atoms: 519 Number of conformers: 1 Conformer: "" Number of residues, atoms: 68, 519 Classifications: {'peptide': 68} Link IDs: {'TRANS': 67} Chain: "D" Number of atoms: 802 Number of conformers: 1 Conformer: "" Number of residues, atoms: 104, 802 Classifications: {'peptide': 104} Incomplete info: {'truncation_to_alanine': 2} Link IDs: {'PTRANS': 5, 'TRANS': 98} Chain breaks: 1 Unresolved non-hydrogen bonds: 7 Unresolved non-hydrogen angles: 8 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 1 Chain: "F" Number of atoms: 209 Number of conformers: 1 Conformer: "" Number of residues, atoms: 10, 209 Classifications: {'RNA': 10} Modifications used: {'rna2p_pyr': 1, 'rna3p_pur': 3, 'rna3p_pyr': 6} Link IDs: {'rna2p': 1, 'rna3p': 8} Chain: "G" Number of atoms: 194 Number of conformers: 1 Conformer: "" Number of residues, atoms: 9, 194 Classifications: {'RNA': 9} Modifications used: {'rna3p_pur': 6, 'rna3p_pyr': 3} Link IDs: {'rna3p': 8} Chain: "A" Number of atoms: 2 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 2 Unusual residues: {' ZN': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} List of CYS excluded from plausible disulfide bonds: (reason: may participate in coordination) ATOM 2450 SG CYS A 301 79.275 41.008 52.381 1.00 28.24 S ATOM 2523 SG CYS A 310 76.742 41.105 49.964 1.00 25.23 S ATOM 3925 SG CYS A 487 89.052 47.626 68.104 1.00 30.40 S ATOM 5178 SG CYS A 645 90.808 48.944 65.420 1.00 35.57 S ATOM 5184 SG CYS A 646 90.164 45.454 65.537 1.00 35.17 S Time building chain proxies: 1.89, per 1000 atoms: 0.19 Number of scatterers: 9932 At special positions: 0 Unit cell: (104.14, 103.32, 124.64, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 6 Type Number sf(0) Zn 2 29.99 S 70 16.00 P 19 15.00 O 1927 8.00 N 1669 7.00 C 6245 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.00 Amino acid : True - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 0.86 Conformation dependent library (CDL) restraints added in 221.0 milliseconds Dynamic metal coordination Zn2+ tetrahedral coordination pdb=" ZN A1001 " pdb="ZN ZN A1001 " - pdb=" SG CYS A 310 " pdb="ZN ZN A1001 " - pdb=" ND1 HIS A 295 " pdb="ZN ZN A1001 " - pdb=" SG CYS A 301 " pdb=" ZN A1002 " pdb="ZN ZN A1002 " - pdb=" ND1 HIS A 642 " pdb="ZN ZN A1002 " - pdb=" SG CYS A 487 " pdb="ZN ZN A1002 " - pdb=" SG CYS A 645 " pdb="ZN ZN A1002 " - pdb=" SG CYS A 646 " Number of angles added : 3 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2302 Finding SS restraints... Secondary structure from input PDB file: 47 helices and 11 sheets defined 43.8% alpha, 11.4% beta 9 base pairs and 14 stacking pairs defined. Time for finding SS restraints: 1.15 Creating SS restraints... Processing helix chain 'A' and resid 1 through 13 removed outlier: 3.753A pdb=" N GLY A 13 " --> pdb=" O ASN A 9 " (cutoff:3.500A) Processing helix chain 'A' and resid 61 through 63 No H-bonds generated for 'chain 'A' and resid 61 through 63' Processing helix chain 'A' and resid 76 through 88 Processing helix chain 'A' and resid 123 through 133 Processing helix chain 'A' and resid 153 through 157 Processing helix chain 'A' and resid 170 through 177 Processing helix chain 'A' and resid 178 through 198 removed outlier: 3.684A pdb=" N ALA A 195 " --> pdb=" O GLN A 191 " (cutoff:3.500A) Processing helix chain 'A' and resid 206 through 210 removed outlier: 3.508A pdb=" N ASN A 209 " --> pdb=" O THR A 206 " (cutoff:3.500A) Processing helix chain 'A' and resid 234 through 242 removed outlier: 3.627A pdb=" N TYR A 238 " --> pdb=" O VAL A 234 " (cutoff:3.500A) Processing helix chain 'A' and resid 243 through 245 No H-bonds generated for 'chain 'A' and resid 243 through 245' Processing helix chain 'A' and resid 275 through 287 Processing helix chain 'A' and resid 303 through 318 removed outlier: 4.070A pdb=" N LEU A 316 " --> pdb=" O ASN A 312 " (cutoff:3.500A) Processing helix chain 'A' and resid 367 through 377 removed outlier: 3.645A pdb=" N ASP A 377 " --> pdb=" O VAL A 373 " (cutoff:3.500A) Processing helix chain 'A' and resid 378 through 383 Processing helix chain 'A' and resid 417 through 426 removed outlier: 3.785A pdb=" N ASP A 421 " --> pdb=" O LYS A 417 " (cutoff:3.500A) removed outlier: 4.265A pdb=" N ALA A 423 " --> pdb=" O PHE A 419 " (cutoff:3.500A) Processing helix chain 'A' and resid 448 through 454 Processing helix chain 'A' and resid 455 through 459 removed outlier: 3.598A pdb=" N TYR A 458 " --> pdb=" O TYR A 455 " (cutoff:3.500A) Processing helix chain 'A' and resid 465 through 480 removed outlier: 4.346A pdb=" N PHE A 471 " --> pdb=" O ARG A 467 " (cutoff:3.500A) removed outlier: 4.827A pdb=" N GLU A 474 " --> pdb=" O LEU A 470 " (cutoff:3.500A) Processing helix chain 'A' and resid 505 through 509 Processing helix chain 'A' and resid 511 through 519 removed outlier: 3.812A pdb=" N TYR A 515 " --> pdb=" O LYS A 511 " (cutoff:3.500A) Processing helix chain 'A' and resid 520 through 532 Processing helix chain 'A' and resid 561 through 581 removed outlier: 3.640A pdb=" N THR A 565 " --> pdb=" O SER A 561 " (cutoff:3.500A) removed outlier: 4.140A pdb=" N THR A 567 " --> pdb=" O CYS A 563 " (cutoff:3.500A) removed outlier: 3.905A pdb=" N GLN A 570 " --> pdb=" O MET A 566 " (cutoff:3.500A) removed outlier: 5.671A pdb=" N LYS A 574 " --> pdb=" O GLN A 570 " (cutoff:3.500A) removed outlier: 5.786A pdb=" N LEU A 575 " --> pdb=" O PHE A 571 " (cutoff:3.500A) Processing helix chain 'A' and resid 596 through 606 Processing helix chain 'A' and resid 621 through 626 Processing helix chain 'A' and resid 627 through 640 Processing helix chain 'A' and resid 647 through 663 Processing helix chain 'A' and resid 686 through 708 removed outlier: 3.915A pdb=" N SER A 692 " --> pdb=" O ALA A 688 " (cutoff:3.500A) Processing helix chain 'A' and resid 711 through 715 removed outlier: 3.921A pdb=" N ILE A 715 " --> pdb=" O GLY A 712 " (cutoff:3.500A) Processing helix chain 'A' and resid 719 through 733 Processing helix chain 'A' and resid 738 through 752 Processing helix chain 'A' and resid 768 through 774 Processing helix chain 'A' and resid 778 through 790 removed outlier: 4.126A pdb=" N SER A 784 " --> pdb=" O LYS A 780 " (cutoff:3.500A) removed outlier: 3.855A pdb=" N VAL A 785 " --> pdb=" O ASN A 781 " (cutoff:3.500A) Processing helix chain 'A' and resid 833 through 843 Processing helix chain 'A' and resid 855 through 864 removed outlier: 3.899A pdb=" N PHE A 859 " --> pdb=" O MET A 855 " (cutoff:3.500A) removed outlier: 3.610A pdb=" N SER A 861 " --> pdb=" O GLU A 857 " (cutoff:3.500A) removed outlier: 3.832A pdb=" N LEU A 862 " --> pdb=" O ARG A 858 " (cutoff:3.500A) removed outlier: 3.561A pdb=" N ALA A 863 " --> pdb=" O PHE A 859 " (cutoff:3.500A) Processing helix chain 'A' and resid 867 through 872 Processing helix chain 'A' and resid 874 through 896 removed outlier: 3.591A pdb=" N ALA A 878 " --> pdb=" O ASN A 874 " (cutoff:3.500A) removed outlier: 3.778A pdb=" N PHE A 881 " --> pdb=" O TYR A 877 " (cutoff:3.500A) Processing helix chain 'A' and resid 922 through 926 removed outlier: 3.924A pdb=" N TYR A 925 " --> pdb=" O GLU A 922 " (cutoff:3.500A) Processing helix chain 'B' and resid 80 through 95 Processing helix chain 'B' and resid 100 through 108 removed outlier: 3.893A pdb=" N ASN B 104 " --> pdb=" O ASN B 100 " (cutoff:3.500A) Processing helix chain 'B' and resid 118 through 125 removed outlier: 4.342A pdb=" N LEU B 122 " --> pdb=" O ASN B 118 " (cutoff:3.500A) Processing helix chain 'B' and resid 134 through 142 removed outlier: 3.570A pdb=" N TYR B 138 " --> pdb=" O ASP B 134 " (cutoff:3.500A) removed outlier: 3.625A pdb=" N ASN B 140 " --> pdb=" O ASN B 136 " (cutoff:3.500A) Processing helix chain 'B' and resid 176 through 180 removed outlier: 3.689A pdb=" N LEU B 180 " --> pdb=" O SER B 177 " (cutoff:3.500A) Processing helix chain 'C' and resid 3 through 19 Processing helix chain 'C' and resid 26 through 41 Processing helix chain 'C' and resid 44 through 62 removed outlier: 4.047A pdb=" N ALA C 48 " --> pdb=" O ASP C 44 " (cutoff:3.500A) removed outlier: 3.514A pdb=" N PHE C 49 " --> pdb=" O THR C 45 " (cutoff:3.500A) removed outlier: 4.007A pdb=" N MET C 62 " --> pdb=" O VAL C 58 " (cutoff:3.500A) Processing helix chain 'D' and resid 85 through 94 Processing helix chain 'D' and resid 100 through 112 Processing sheet with id=AA1, first strand: chain 'A' and resid 19 through 20 removed outlier: 3.509A pdb=" N GLN A 57 " --> pdb=" O THR A 20 " (cutoff:3.500A) removed outlier: 7.088A pdb=" N PHE A 56 " --> pdb=" O SER A 68 " (cutoff:3.500A) removed outlier: 6.858A pdb=" N SER A 68 " --> pdb=" O PHE A 56 " (cutoff:3.500A) removed outlier: 5.386A pdb=" N GLU A 58 " --> pdb=" O ILE A 66 " (cutoff:3.500A) removed outlier: 3.506A pdb=" N ILE A 66 " --> pdb=" O GLU A 58 " (cutoff:3.500A) removed outlier: 3.517A pdb=" N ASN A 64 " --> pdb=" O ASP A 60 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 35 through 38 removed outlier: 3.505A pdb=" N ILE A 37 " --> pdb=" O GLY A 44 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 102 through 104 removed outlier: 4.282A pdb=" N PHE A 102 " --> pdb=" O HIS A 113 " (cutoff:3.500A) removed outlier: 4.255A pdb=" N HIS A 113 " --> pdb=" O PHE A 102 " (cutoff:3.500A) removed outlier: 3.914A pdb=" N PHE A 104 " --> pdb=" O VAL A 111 " (cutoff:3.500A) removed outlier: 3.634A pdb=" N VAL A 111 " --> pdb=" O PHE A 104 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA3 Processing sheet with id=AA4, first strand: chain 'A' and resid 223 through 224 removed outlier: 6.527A pdb=" N VAL A 233 " --> pdb=" O ILE A 201 " (cutoff:3.500A) removed outlier: 6.831A pdb=" N GLY A 203 " --> pdb=" O VAL A 233 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'A' and resid 353 through 354 removed outlier: 3.570A pdb=" N VAL A 342 " --> pdb=" O ARG A 331 " (cutoff:3.500A) removed outlier: 6.367A pdb=" N LEU A 329 " --> pdb=" O THR A 344 " (cutoff:3.500A) removed outlier: 5.340A pdb=" N TYR A 346 " --> pdb=" O GLY A 327 " (cutoff:3.500A) removed outlier: 7.163A pdb=" N GLY A 327 " --> pdb=" O TYR A 346 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'A' and resid 556 through 559 removed outlier: 3.554A pdb=" N MET A 542 " --> pdb=" O MET A 668 " (cutoff:3.500A) removed outlier: 4.130A pdb=" N ALA A 399 " --> pdb=" O LEU A 673 " (cutoff:3.500A) removed outlier: 6.292A pdb=" N VAL A 675 " --> pdb=" O SER A 397 " (cutoff:3.500A) removed outlier: 6.651A pdb=" N SER A 397 " --> pdb=" O VAL A 675 " (cutoff:3.500A) removed outlier: 4.815A pdb=" N VAL A 398 " --> pdb=" O LEU A 388 " (cutoff:3.500A) removed outlier: 3.951A pdb=" N LEU A 388 " --> pdb=" O VAL A 398 " (cutoff:3.500A) removed outlier: 6.384A pdb=" N LEU A 387 " --> pdb=" O MET B 129 " (cutoff:3.500A) removed outlier: 7.791A pdb=" N VAL B 131 " --> pdb=" O LEU A 387 " (cutoff:3.500A) removed outlier: 6.683A pdb=" N LEU A 389 " --> pdb=" O VAL B 131 " (cutoff:3.500A) removed outlier: 4.041A pdb=" N THR B 187 " --> pdb=" O GLN B 158 " (cutoff:3.500A) removed outlier: 4.796A pdb=" N GLN B 158 " --> pdb=" O THR B 187 " (cutoff:3.500A) removed outlier: 6.758A pdb=" N LEU B 189 " --> pdb=" O ILE B 156 " (cutoff:3.500A) removed outlier: 5.954A pdb=" N ILE B 156 " --> pdb=" O LEU B 189 " (cutoff:3.500A) removed outlier: 6.375A pdb=" N ALA B 191 " --> pdb=" O TRP B 154 " (cutoff:3.500A) removed outlier: 5.606A pdb=" N TRP B 154 " --> pdb=" O ALA B 191 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'A' and resid 612 through 616 removed outlier: 3.525A pdb=" N PHE A 766 " --> pdb=" O HIS A 613 " (cutoff:3.500A) removed outlier: 3.521A pdb=" N SER A 754 " --> pdb=" O CYS A 765 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'A' and resid 612 through 616 removed outlier: 3.525A pdb=" N PHE A 766 " --> pdb=" O HIS A 613 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'A' and resid 816 through 821 Processing sheet with id=AB1, first strand: chain 'D' and resid 115 through 116 removed outlier: 3.597A pdb=" N THR D 187 " --> pdb=" O GLN D 158 " (cutoff:3.500A) removed outlier: 3.884A pdb=" N GLN D 158 " --> pdb=" O THR D 187 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'D' and resid 142 through 143 379 hydrogen bonds defined for protein. 1065 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 23 hydrogen bonds 46 hydrogen bond angles 0 basepair planarities 9 basepair parallelities 14 stacking parallelities Total time for adding SS restraints: 1.64 Time building geometry restraints manager: 1.08 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: 3113 1.34 - 1.46: 1710 1.46 - 1.57: 5220 1.57 - 1.69: 36 1.69 - 1.81: 105 Bond restraints: 10184 Sorted by residual: bond pdb=" C TRP B 182 " pdb=" N PRO B 183 " ideal model delta sigma weight residual 1.335 1.350 -0.015 8.70e-03 1.32e+04 3.10e+00 bond pdb=" CA SER A 709 " pdb=" C SER A 709 " ideal model delta sigma weight residual 1.523 1.499 0.024 1.80e-02 3.09e+03 1.80e+00 bond pdb=" N ILE A 847 " pdb=" CA ILE A 847 " ideal model delta sigma weight residual 1.474 1.456 0.018 1.57e-02 4.06e+03 1.33e+00 bond pdb=" CA VAL A 844 " pdb=" CB VAL A 844 " ideal model delta sigma weight residual 1.550 1.539 0.012 1.03e-02 9.43e+03 1.28e+00 bond pdb=" CA VAL D 167 " pdb=" CB VAL D 167 " ideal model delta sigma weight residual 1.530 1.542 -0.012 1.05e-02 9.07e+03 1.24e+00 ... (remaining 10179 not shown) Histogram of bond angle deviations from ideal: 0.00 - 2.10: 13723 2.10 - 4.19: 172 4.19 - 6.29: 19 6.29 - 8.38: 1 8.38 - 10.48: 1 Bond angle restraints: 13916 Sorted by residual: angle pdb=" C LEU A 119 " pdb=" N THR A 120 " pdb=" CA THR A 120 " ideal model delta sigma weight residual 121.54 127.94 -6.40 1.91e+00 2.74e-01 1.12e+01 angle pdb=" C TYR A 787 " pdb=" N TYR A 788 " pdb=" CA TYR A 788 " ideal model delta sigma weight residual 122.26 117.24 5.02 1.73e+00 3.34e-01 8.43e+00 angle pdb=" C TRP A 916 " pdb=" N GLU A 917 " pdb=" CA GLU A 917 " ideal model delta sigma weight residual 121.95 132.43 -10.48 3.74e+00 7.15e-02 7.84e+00 angle pdb=" C ASP A 846 " pdb=" N ILE A 847 " pdb=" CA ILE A 847 " ideal model delta sigma weight residual 121.65 119.11 2.54 9.40e-01 1.13e+00 7.33e+00 angle pdb=" N ILE A 847 " pdb=" CA ILE A 847 " pdb=" C ILE A 847 " ideal model delta sigma weight residual 106.21 109.10 -2.89 1.07e+00 8.73e-01 7.30e+00 ... (remaining 13911 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 28.37: 5821 28.37 - 56.73: 193 56.73 - 85.10: 12 85.10 - 113.46: 3 113.46 - 141.83: 1 Dihedral angle restraints: 6030 sinusoidal: 2471 harmonic: 3559 Sorted by residual: dihedral pdb=" CA ILE A 106 " pdb=" C ILE A 106 " pdb=" N ASP A 107 " pdb=" CA ASP A 107 " ideal model delta harmonic sigma weight residual 180.00 151.29 28.71 0 5.00e+00 4.00e-02 3.30e+01 dihedral pdb=" CA ASP A 107 " pdb=" C ASP A 107 " pdb=" N GLY A 108 " pdb=" CA GLY A 108 " ideal model delta harmonic sigma weight residual -180.00 -152.19 -27.81 0 5.00e+00 4.00e-02 3.09e+01 dihedral pdb=" CA ILE A 757 " pdb=" C ILE A 757 " pdb=" N LEU A 758 " pdb=" CA LEU A 758 " ideal model delta harmonic sigma weight residual 180.00 154.74 25.26 0 5.00e+00 4.00e-02 2.55e+01 ... (remaining 6027 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.056: 1311 0.056 - 0.111: 263 0.111 - 0.167: 17 0.167 - 0.222: 1 0.222 - 0.278: 1 Chirality restraints: 1593 Sorted by residual: chirality pdb=" CB THR A 120 " pdb=" CA THR A 120 " pdb=" OG1 THR A 120 " pdb=" CG2 THR A 120 " both_signs ideal model delta sigma weight residual False 2.55 2.27 0.28 2.00e-01 2.50e+01 1.93e+00 chirality pdb=" C3' G F -3 " pdb=" C4' G F -3 " pdb=" O3' G F -3 " pdb=" C2' G F -3 " both_signs ideal model delta sigma weight residual False -2.48 -2.29 -0.19 2.00e-01 2.50e+01 8.72e-01 chirality pdb=" CA TRP D 182 " pdb=" N TRP D 182 " pdb=" C TRP D 182 " pdb=" CB TRP D 182 " both_signs ideal model delta sigma weight residual False 2.51 2.35 0.16 2.00e-01 2.50e+01 6.61e-01 ... (remaining 1590 not shown) Planarity restraints: 1712 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C SER B 177 " -0.026 5.00e-02 4.00e+02 3.85e-02 2.38e+00 pdb=" N PRO B 178 " 0.067 5.00e-02 4.00e+02 pdb=" CA PRO B 178 " -0.019 5.00e-02 4.00e+02 pdb=" CD PRO B 178 " -0.021 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C TRP D 182 " 0.026 5.00e-02 4.00e+02 3.84e-02 2.36e+00 pdb=" N PRO D 183 " -0.066 5.00e-02 4.00e+02 pdb=" CA PRO D 183 " 0.019 5.00e-02 4.00e+02 pdb=" CD PRO D 183 " 0.022 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C1' G F -1 " -0.019 2.00e-02 2.50e+03 7.91e-03 1.88e+00 pdb=" N9 G F -1 " 0.018 2.00e-02 2.50e+03 pdb=" C8 G F -1 " 0.003 2.00e-02 2.50e+03 pdb=" N7 G F -1 " -0.000 2.00e-02 2.50e+03 pdb=" C5 G F -1 " -0.000 2.00e-02 2.50e+03 pdb=" C6 G F -1 " -0.002 2.00e-02 2.50e+03 pdb=" O6 G F -1 " -0.005 2.00e-02 2.50e+03 pdb=" N1 G F -1 " -0.002 2.00e-02 2.50e+03 pdb=" C2 G F -1 " 0.001 2.00e-02 2.50e+03 pdb=" N2 G F -1 " 0.001 2.00e-02 2.50e+03 pdb=" N3 G F -1 " 0.003 2.00e-02 2.50e+03 pdb=" C4 G F -1 " 0.003 2.00e-02 2.50e+03 ... (remaining 1709 not shown) Histogram of nonbonded interaction distances: 2.20 - 2.74: 964 2.74 - 3.28: 10318 3.28 - 3.82: 16222 3.82 - 4.36: 20084 4.36 - 4.90: 33800 Nonbonded interactions: 81388 Sorted by model distance: nonbonded pdb=" OG SER A 239 " pdb=" OD1 ASP A 465 " model vdw 2.199 3.040 nonbonded pdb=" OH TYR A 294 " pdb=" OD1 ASN A 314 " model vdw 2.233 3.040 nonbonded pdb=" OD2 ASP A 390 " pdb=" OH TYR A 674 " model vdw 2.237 3.040 nonbonded pdb=" OG SER A 433 " pdb=" O VAL A 435 " model vdw 2.254 3.040 nonbonded pdb=" OH TYR A 420 " pdb=" OG SER C 4 " model vdw 2.259 3.040 ... (remaining 81383 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 'B' and (resid 84 through 122 or (resid 127 and (name N or name CA or nam \ e C or name O or name CB )) or resid 128 through 188 or (resid 189 and (name N o \ r name CA or name C or name O or name CB )) or resid 190 through 191)) selection = (chain 'D' and (resid 84 through 181 or (resid 182 and (name N or name CA or nam \ e C or name O or name CB )) or resid 183 through 191)) } 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: 2.510 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.200 Check model and map are aligned: 0.040 Set scattering table: 0.020 Process input model: 10.990 Find NCS groups from input model: 0.060 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:1.240 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 15.080 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7716 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.012 0.692 10191 Z= 0.404 Angle : 0.726 28.251 13919 Z= 0.360 Chirality : 0.043 0.278 1593 Planarity : 0.003 0.039 1712 Dihedral : 13.936 141.826 3728 Min Nonbonded Distance : 2.199 Molprobity Statistics. All-atom Clashscore : 11.62 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.50 % Favored : 91.50 % Rotamer: Outliers : 0.00 % Allowed : 4.75 % Favored : 95.25 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.50 % 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 are not related in a simple way. whole: -4.42 (0.19), residues: 1188 helix: -1.96 (0.19), residues: 500 sheet: -3.57 (0.36), residues: 133 loop : -3.74 (0.20), residues: 555 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.002 0.000 ARG A 18 TYR 0.013 0.001 TYR A 921 PHE 0.015 0.001 PHE A 753 TRP 0.021 0.001 TRP D 182 HIS 0.005 0.001 HIS A 113 Details of bonding type rmsd/Z covalent geometry : bond 0.00344 / 0.16 (10184) covalent geometry : angle 0.61030 / 0.35 (13916) hydrogen bonds : bond 0.12949 / 8.89 ( 399) hydrogen bonds : angle 6.20992 / 4.36 ( 1111) metal coordination : bond 0.43123 / 25.49 ( 7) metal coordination : angle 26.81661 / 14.49 ( 3) *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 167 residues out of total 1062 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 167 time to evaluate : 0.380 Fit side-chains revert: symmetry clash REVERT: A 452 ASP cc_start: 0.7776 (t70) cc_final: 0.7492 (t70) REVERT: D 106 ILE cc_start: 0.5964 (mm) cc_final: 0.5570 (mt) outliers start: 0 outliers final: 0 residues processed: 167 average time/residue: 0.0999 time to fit residues: 22.5420 Evaluate side-chains 120 residues out of total 1062 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 120 time to evaluate : 0.258 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 108 optimal weight: 0.9980 chunk 49 optimal weight: 2.9990 chunk 97 optimal weight: 0.0000 chunk 113 optimal weight: 5.9990 chunk 53 optimal weight: 2.9990 chunk 5 optimal weight: 10.0000 chunk 33 optimal weight: 0.0770 chunk 65 optimal weight: 0.9980 chunk 62 optimal weight: 6.9990 chunk 51 optimal weight: 6.9990 chunk 117 optimal weight: 5.9990 overall best weight: 1.0144 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... A 75 HIS A 81 GLN A 113 HIS A 133 HIS A 312 ASN A 356 ASN A 439 HIS A 572 HIS A 650 HIS A 789 GLN B 108 ASN B 168 GLN C 36 HIS C 63 GLN D 104 ASN Total number of N/Q/H flips: 15 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4030 r_free = 0.4030 target = 0.159813 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 43)----------------| | r_work = 0.3739 r_free = 0.3739 target = 0.135717 restraints weight = 15874.287| |-----------------------------------------------------------------------------| r_work (start): 0.3737 rms_B_bonded: 2.45 r_work: 0.3642 rms_B_bonded: 3.00 restraints_weight: 0.5000 r_work (final): 0.3642 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7714 moved from start: 0.1252 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.030 10191 Z= 0.140 Angle : 0.595 7.759 13919 Z= 0.311 Chirality : 0.043 0.236 1593 Planarity : 0.004 0.059 1712 Dihedral : 8.089 145.719 1530 Min Nonbonded Distance : 2.395 Molprobity Statistics. All-atom Clashscore : 10.18 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.58 % Favored : 92.42 % Rotamer: Outliers : 1.14 % Allowed : 11.49 % Favored : 87.37 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.50 % 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 are not related in a simple way. whole: -3.14 (0.22), residues: 1188 helix: -0.77 (0.22), residues: 491 sheet: -2.84 (0.40), residues: 127 loop : -3.16 (0.22), residues: 570 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.006 0.000 ARG A 331 TYR 0.015 0.001 TYR A 887 PHE 0.016 0.001 PHE A 741 TRP 0.012 0.001 TRP D 182 HIS 0.004 0.001 HIS A 295 Details of bonding type rmsd/Z covalent geometry : bond 0.00319 / 0.14 (10184) covalent geometry : angle 0.58792 / 0.31 (13916) hydrogen bonds : bond 0.03886 / 2.63 ( 399) hydrogen bonds : angle 5.07387 / 3.59 ( 1111) metal coordination : bond 0.01134 / 0.62 ( 7) metal coordination : angle 6.28750 / 3.44 ( 3) *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 151 residues out of total 1062 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 12 poor density : 139 time to evaluate : 0.226 Fit side-chains revert: symmetry clash REVERT: A 197 ARG cc_start: 0.5771 (tpp-160) cc_final: 0.5510 (tpp-160) REVERT: A 452 ASP cc_start: 0.8476 (t70) cc_final: 0.8107 (t70) REVERT: A 500 LYS cc_start: 0.7580 (mmtt) cc_final: 0.7379 (mttm) REVERT: A 583 ARG cc_start: 0.8676 (mpt90) cc_final: 0.8395 (mmt180) REVERT: A 601 MET cc_start: 0.7606 (mtp) cc_final: 0.7391 (mtm) REVERT: A 818 MET cc_start: 0.8672 (ttp) cc_final: 0.8451 (ttt) REVERT: A 821 LYS cc_start: 0.8782 (tttm) cc_final: 0.8518 (tptp) REVERT: A 855 MET cc_start: 0.4960 (tpp) cc_final: 0.4571 (tpp) REVERT: A 919 GLU cc_start: 0.8029 (mp0) cc_final: 0.7746 (mp0) REVERT: D 94 MET cc_start: 0.6409 (mtt) cc_final: 0.6039 (mtp) REVERT: D 106 ILE cc_start: 0.5592 (mm) cc_final: 0.5241 (mt) outliers start: 12 outliers final: 9 residues processed: 144 average time/residue: 0.0858 time to fit residues: 17.2719 Evaluate side-chains 130 residues out of total 1062 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 9 poor density : 121 time to evaluate : 0.273 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 120 THR Chi-restraints excluded: chain A residue 123 THR Chi-restraints excluded: chain A residue 223 ILE Chi-restraints excluded: chain A residue 608 ASP Chi-restraints excluded: chain A residue 662 VAL Chi-restraints excluded: chain A residue 672 SER Chi-restraints excluded: chain A residue 805 LEU Chi-restraints excluded: chain B residue 164 SER Chi-restraints excluded: chain D residue 153 LEU Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 66 optimal weight: 5.9990 chunk 23 optimal weight: 1.9990 chunk 74 optimal weight: 3.9990 chunk 82 optimal weight: 0.2980 chunk 71 optimal weight: 5.9990 chunk 56 optimal weight: 0.5980 chunk 28 optimal weight: 0.8980 chunk 52 optimal weight: 0.6980 chunk 60 optimal weight: 0.3980 chunk 81 optimal weight: 1.9990 chunk 114 optimal weight: 3.9990 overall best weight: 0.5780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... ** A 57 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 75 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 752 HIS C 31 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4059 r_free = 0.4059 target = 0.162555 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 44)----------------| | r_work = 0.3770 r_free = 0.3770 target = 0.138273 restraints weight = 15854.881| |-----------------------------------------------------------------------------| r_work (start): 0.3773 rms_B_bonded: 2.47 r_work: 0.3678 rms_B_bonded: 3.03 restraints_weight: 0.5000 r_work (final): 0.3678 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7651 moved from start: 0.1612 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.048 10191 Z= 0.112 Angle : 0.536 5.930 13919 Z= 0.283 Chirality : 0.041 0.235 1593 Planarity : 0.003 0.053 1712 Dihedral : 7.885 146.659 1530 Min Nonbonded Distance : 2.524 Molprobity Statistics. All-atom Clashscore : 9.15 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.57 % Favored : 93.43 % Rotamer: Outliers : 1.71 % Allowed : 14.25 % Favored : 84.05 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.50 % 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 are not related in a simple way. whole: -2.43 (0.23), residues: 1188 helix: -0.16 (0.23), residues: 491 sheet: -2.21 (0.44), residues: 124 loop : -2.86 (0.23), residues: 573 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.002 0.000 ARG A 132 TYR 0.015 0.001 TYR A 887 PHE 0.015 0.001 PHE A 741 TRP 0.011 0.001 TRP D 182 HIS 0.004 0.001 HIS A 309 Details of bonding type rmsd/Z covalent geometry : bond 0.00251 / 0.11 (10184) covalent geometry : angle 0.53280 / 0.28 (13916) hydrogen bonds : bond 0.03305 / 2.21 ( 399) hydrogen bonds : angle 4.73673 / 3.36 ( 1111) metal coordination : bond 0.00649 / 0.46 ( 7) metal coordination : angle 4.07965 / 2.26 ( 3) *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 151 residues out of total 1062 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 18 poor density : 133 time to evaluate : 0.362 Fit side-chains REVERT: A 197 ARG cc_start: 0.5643 (tpp-160) cc_final: 0.5366 (tpp-160) REVERT: A 452 ASP cc_start: 0.8383 (t70) cc_final: 0.7993 (t70) REVERT: A 583 ARG cc_start: 0.8641 (mpt90) cc_final: 0.8323 (mmt180) REVERT: A 821 LYS cc_start: 0.8735 (tttm) cc_final: 0.8511 (tptp) REVERT: A 919 GLU cc_start: 0.8026 (mp0) cc_final: 0.7645 (mp0) REVERT: D 94 MET cc_start: 0.6258 (mtt) cc_final: 0.5927 (mtp) REVERT: D 106 ILE cc_start: 0.5502 (mm) cc_final: 0.5148 (mt) outliers start: 18 outliers final: 11 residues processed: 145 average time/residue: 0.0775 time to fit residues: 16.2301 Evaluate side-chains 132 residues out of total 1062 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 11 poor density : 121 time to evaluate : 0.256 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 120 THR Chi-restraints excluded: chain A residue 123 THR Chi-restraints excluded: chain A residue 608 ASP Chi-restraints excluded: chain A residue 662 VAL Chi-restraints excluded: chain A residue 681 SER Chi-restraints excluded: chain A residue 799 CYS Chi-restraints excluded: chain A residue 805 LEU Chi-restraints excluded: chain A residue 820 VAL Chi-restraints excluded: chain A residue 928 HIS Chi-restraints excluded: chain B residue 164 SER Chi-restraints excluded: chain D residue 114 CYS Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 119 random chunks: chunk 1 optimal weight: 0.9990 chunk 17 optimal weight: 0.9980 chunk 105 optimal weight: 2.9990 chunk 100 optimal weight: 6.9990 chunk 86 optimal weight: 2.9990 chunk 107 optimal weight: 0.8980 chunk 75 optimal weight: 1.9990 chunk 8 optimal weight: 0.6980 chunk 48 optimal weight: 0.7980 chunk 34 optimal weight: 4.9990 chunk 101 optimal weight: 3.9990 overall best weight: 0.8782 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... ** A 57 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 75 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** A 168 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4061 r_free = 0.4061 target = 0.162791 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 61)----------------| | r_work = 0.3771 r_free = 0.3771 target = 0.138330 restraints weight = 15903.210| |-----------------------------------------------------------------------------| r_work (start): 0.3770 rms_B_bonded: 2.48 r_work: 0.3673 rms_B_bonded: 3.05 restraints_weight: 0.5000 r_work (final): 0.3673 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7677 moved from start: 0.1905 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.025 10191 Z= 0.120 Angle : 0.531 5.656 13919 Z= 0.279 Chirality : 0.041 0.238 1593 Planarity : 0.003 0.047 1712 Dihedral : 7.732 145.545 1530 Min Nonbonded Distance : 2.505 Molprobity Statistics. All-atom Clashscore : 8.79 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.82 % Favored : 93.18 % Rotamer: Outliers : 1.99 % Allowed : 14.72 % Favored : 83.29 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 2.50 % 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 are not related in a simple way. whole: -1.98 (0.24), residues: 1188 helix: 0.22 (0.24), residues: 488 sheet: -1.94 (0.46), residues: 124 loop : -2.62 (0.23), residues: 576 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.002 0.000 ARG A 365 TYR 0.014 0.001 TYR A 887 PHE 0.015 0.001 PHE A 741 TRP 0.009 0.001 TRP D 182 HIS 0.004 0.001 HIS A 309 Details of bonding type rmsd/Z covalent geometry : bond 0.00270 / 0.12 (10184) covalent geometry : angle 0.52838 / 0.28 (13916) hydrogen bonds : bond 0.03242 / 2.16 ( 399) hydrogen bonds : angle 4.59909 / 3.27 ( 1111) metal coordination : bond 0.00602 / 0.41 ( 7) metal coordination : angle 3.39390 / 1.92 ( 3) *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2376 Ramachandran restraints generated. 1188 Oldfield, 0 Emsley, 1188 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 148 residues out of total 1062 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 21 poor density : 127 time to evaluate : 0.343 Fit side-chains revert: symmetry clash REVERT: A 102 PHE cc_start: 0.6422 (p90) cc_final: 0.5976 (p90) REVERT: A 160 LYS cc_start: 0.7705 (tppt) cc_final: 0.7132 (ptmt) REVERT: A 452 ASP cc_start: 0.8364 (t70) cc_final: 0.7976 (t70) REVERT: A 552 ASN cc_start: 0.7671 (p0) cc_final: 0.7092 (t0) REVERT: A 583 ARG cc_start: 0.8634 (mpt90) cc_final: 0.8318 (mmt180) REVERT: A 608 ASP cc_start: 0.7659 (OUTLIER) cc_final: 0.7412 (m-30) REVERT: A 880 VAL cc_start: 0.8175 (t) cc_final: 0.7892 (t) REVERT: A 919 GLU cc_start: 0.8043 (mp0) cc_final: 0.7635 (mp0) REVERT: B 80 ARG cc_start: 0.8095 (OUTLIER) cc_final: 0.7675 (mtm-85) REVERT: D 94 MET cc_start: 0.6242 (mtt) cc_final: 0.5962 (mtp) REVERT: D 106 ILE cc_start: 0.5555 (mm) cc_final: 0.5179 (mt) outliers start: 21 outliers final: 14 residues processed: 140 average time/residue: 0.0825 time to fit residues: 16.4640 Evaluate side-chains 140 residues out of total 1062 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 16 poor density : 124 time to evaluate : 0.356 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 120 THR Chi-restraints excluded: chain A residue 123 THR Chi-restraints excluded: chain A residue 223 ILE Chi-restraints excluded: chain A residue 350 GLU Chi-restraints excluded: chain A residue 608 ASP Chi-restraints excluded: chain A residue 644 THR Chi-restraints excluded: chain A residue 662 VAL Chi-restraints excluded: chain A residue 681 SER Chi-restraints excluded: chain A residue 799 CYS Chi-restraints excluded: chain A residue 805 LEU Chi-restraints excluded: chain A residue 820 VAL Chi-restraints excluded: chain A residue 928 HIS Chi-restraints excluded: chain B residue 80 ARG Chi-restraints excluded: chain B residue 164 SER Chi-restraints excluded: chain D residue 101 ASP Chi-restraints excluded: chain D residue 114 CYS Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/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-dev-6108/lib/python3.11/site-packages/iotbx/cli_parser.py", line 987, in run_program task.run() File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/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-dev-6108/lib/python3.11/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-dev-6108/lib/python3.11/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-dev-6108/lib/python3.11/site-packages/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/phenix/refinement/macro_cycle_real_space.py", line 702, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/phenix/refinement/macro_cycle_real_space.py", line 755, in minimization_no_ncs ro = mmtbx.refinement.real_space.individual_sites.easy( ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 44, in __init__ self.weight = mmtbx.refinement.real_space.weight.run( ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/mmtbx/refinement/real_space/weight.py", line 84, in __init__ ro.refine( File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 362, in refine real_space_result = refinery( ^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 233, in __init__ refiner.refine( File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 100, in refine self.refined = maptbx.real_space_refinement_simple.lbfgs( ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/cctbx/maptbx/real_space_refinement_simple.py", line 168, in __init__ O.minimizer = scitbx.lbfgs.run( ^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/scitbx/lbfgs/__init__.py", line 270, in run return run_c_plus_plus( ^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/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-dev-6108/lib/python3.11/site-packages/cctbx/maptbx/real_space_refinement_simple.py", line 243, in compute_functional_and_gradients gr_e = O.geometry_restraints_manager.energies_sites( ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/cctbx/geometry_restraints/manager.py", line 1602, in energies_sites pair_proxies = self.pair_proxies( ^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/cctbx/geometry_restraints/manager.py", line 1505, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/cctbx/geometry_restraints/manager.py", line 1386, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 67.6342 > 50: distance: 81 - 87: 23.655 distance: 87 - 88: 4.517 distance: 88 - 89: 25.938 distance: 88 - 91: 27.903 distance: 89 - 90: 9.018 distance: 89 - 99: 43.644 distance: 91 - 92: 12.730 distance: 92 - 93: 13.123 distance: 92 - 94: 20.735 distance: 93 - 95: 10.201 distance: 94 - 96: 11.862 distance: 95 - 97: 7.851 distance: 96 - 97: 10.230 distance: 99 - 100: 48.191 distance: 100 - 101: 18.814 distance: 100 - 103: 22.529 distance: 101 - 102: 26.352 distance: 101 - 104: 38.585 distance: 104 - 105: 14.729 distance: 105 - 106: 35.442 distance: 105 - 108: 26.399 distance: 106 - 107: 33.145 distance: 106 - 113: 46.528 distance: 108 - 109: 26.994 distance: 109 - 110: 21.115 distance: 110 - 111: 39.085 distance: 110 - 112: 31.232 distance: 113 - 114: 16.040 distance: 113 - 119: 30.528 distance: 114 - 115: 48.870 distance: 114 - 117: 11.259 distance: 115 - 120: 67.617 distance: 117 - 118: 56.464 distance: 118 - 119: 44.332 distance: 120 - 121: 12.922 distance: 121 - 122: 14.091 distance: 121 - 124: 33.132 distance: 122 - 123: 7.245 distance: 122 - 129: 17.675 distance: 124 - 125: 7.871 distance: 125 - 126: 18.042 distance: 126 - 127: 9.273 distance: 126 - 128: 6.090 distance: 129 - 130: 19.024 distance: 130 - 131: 27.515 distance: 130 - 133: 21.024 distance: 131 - 132: 35.552 distance: 131 - 140: 13.991 distance: 133 - 134: 4.919 distance: 134 - 135: 6.924 distance: 134 - 136: 7.357 distance: 135 - 137: 5.438 distance: 136 - 138: 7.257 distance: 137 - 139: 3.334 distance: 138 - 139: 5.263 distance: 141 - 142: 7.811 distance: 141 - 144: 5.357 distance: 142 - 143: 11.010 distance: 145 - 146: 3.198 distance: 145 - 147: 6.941 distance: 148 - 150: 5.604 distance: 149 - 150: 5.969 distance: 152 - 153: 7.107 distance: 153 - 154: 15.518 distance: 153 - 156: 14.273 distance: 154 - 155: 13.330 distance: 154 - 161: 8.368 distance: 156 - 157: 7.883 distance: 157 - 158: 3.596 distance: 158 - 159: 4.305 distance: 158 - 160: 5.710 distance: 161 - 162: 11.270 distance: 162 - 163: 16.944 distance: 162 - 165: 20.253 distance: 163 - 164: 6.779 distance: 163 - 166: 11.990 distance: 164 - 185: 10.087 distance: 166 - 167: 32.087 distance: 167 - 168: 8.345 distance: 167 - 170: 18.448 distance: 168 - 169: 24.320 distance: 168 - 174: 3.855 distance: 170 - 171: 23.617 distance: 171 - 172: 6.385 distance: 172 - 173: 9.718