Starting phenix.real_space_refine on Sat Aug 23 16:02:21 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/7ckm_30388/08_2025/7ckm_30388.cif Found real_map, /net/cci-nas-00/data/ceres_data/7ckm_30388/08_2025/7ckm_30388.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.37 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { model { file = "/net/cci-nas-00/data/ceres_data/7ckm_30388/08_2025/7ckm_30388.cif" } default_model = "/net/cci-nas-00/data/ceres_data/7ckm_30388/08_2025/7ckm_30388.cif" real_map_files = "/net/cci-nas-00/data/ceres_data/7ckm_30388/08_2025/7ckm_30388.map" default_real_map = "/net/cci-nas-00/data/ceres_data/7ckm_30388/08_2025/7ckm_30388.map" } resolution = 3.37 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.003 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 Mn 1 7.51 5 Zn 4 6.06 5 S 94 5.16 5 C 8358 2.51 5 N 2205 2.21 5 O 2452 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 121 residue(s): 0.02s 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: 13114 Number of models: 1 Model: "" Number of chains: 4 Chain: "A" Number of atoms: 12702 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1603, 12702 Classifications: {'peptide': 1603} Incomplete info: {'truncation_to_alanine': 29} Link IDs: {'PTRANS': 47, 'TRANS': 1555} Chain breaks: 16 Unresolved non-hydrogen bonds: 146 Unresolved non-hydrogen angles: 179 Unresolved non-hydrogen dihedrals: 130 Unresolved non-hydrogen chiralities: 1 Planarities with less than four sites: {'GLN:plan1': 1, 'ARG:plan': 7, 'TRP:plan': 1, 'ASP:plan': 2, 'GLU:plan': 7, 'HIS:plan': 2, 'PHE:plan': 3} Unresolved non-hydrogen planarities: 110 Chain: "B" Number of atoms: 407 Number of conformers: 1 Conformer: "" Number of residues, atoms: 49, 407 Classifications: {'peptide': 49} Link IDs: {'PTRANS': 2, 'TRANS': 46} Chain: "A" Number of atoms: 3 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 3 Unusual residues: {' MN': 1, ' ZN': 2} Classifications: {'undetermined': 3} Link IDs: {None: 2} Chain: "B" 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 2171 SG CYS A 284 51.305 69.871 26.726 1.00 33.30 S ATOM 2195 SG CYS A 287 54.815 70.105 24.870 1.00 35.01 S ATOM 3474 SG CYS A 470 51.367 70.972 22.964 1.00 81.64 S ATOM 11406 SG CYS A1650 41.739 16.812 43.377 1.00 82.16 S ATOM 11419 SG CYS A1652 38.652 16.222 44.857 1.00 81.31 S ATOM 11438 SG CYS A1655 40.724 20.132 43.603 1.00 64.80 S ATOM 11510 SG CYS A1664 40.880 17.040 46.700 1.00 55.65 S ATOM 12763 SG CYS B 39 23.850 80.360 54.941 1.00 42.32 S ATOM 12784 SG CYS B 42 26.136 82.388 53.549 1.00 34.63 S ATOM 12918 SG CYS B 58 23.069 83.737 54.842 1.00 57.08 S ATOM 12943 SG CYS B 61 22.984 81.129 51.433 1.00 56.47 S ATOM 12866 SG CYS B 52 16.451 73.693 57.557 1.00 73.10 S ATOM 13033 SG CYS B 72 18.354 72.080 55.250 1.00 57.48 S ATOM 13057 SG CYS B 75 18.227 70.906 58.418 1.00 46.83 S Time building chain proxies: 2.68, per 1000 atoms: 0.20 Number of scatterers: 13114 At special positions: 0 Unit cell: (122.04, 114.48, 123.12, 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 4 29.99 Mn 1 24.99 S 94 16.00 O 2452 8.00 N 2205 7.00 C 8358 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=2, symmetry=0 Simple disulfide: pdb=" SG CYS A 55 " - pdb=" SG CYS A 60 " distance=2.04 Simple disulfide: pdb=" SG CYS A1691 " - pdb=" SG CYS A1776 " 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.80 Conformation dependent library (CDL) restraints added in 518.0 milliseconds Enol-peptide restraints added in 1.2 microseconds Dynamic metal coordination Zn2+ tetrahedral coordination pdb=" ZN A2302 " pdb="ZN ZN A2302 " - pdb=" ND1 HIS A 472 " pdb="ZN ZN A2302 " - pdb=" SG CYS A 287 " pdb="ZN ZN A2302 " - pdb=" SG CYS A 284 " pdb="ZN ZN A2302 " - pdb=" SG CYS A 470 " pdb=" ZN A2303 " pdb="ZN ZN A2303 " - pdb=" SG CYS A1650 " pdb="ZN ZN A2303 " - pdb=" SG CYS A1652 " pdb="ZN ZN A2303 " - pdb=" SG CYS A1655 " pdb="ZN ZN A2303 " - pdb=" SG CYS A1664 " pdb=" ZN B 101 " pdb="ZN ZN B 101 " - pdb=" SG CYS B 39 " pdb="ZN ZN B 101 " - pdb=" SG CYS B 58 " pdb="ZN ZN B 101 " - pdb=" SG CYS B 61 " pdb="ZN ZN B 101 " - pdb=" SG CYS B 42 " pdb=" ZN B 102 " pdb="ZN ZN B 102 " - pdb=" NE2 HIS B 55 " pdb="ZN ZN B 102 " - pdb=" SG CYS B 75 " pdb="ZN ZN B 102 " - pdb=" SG CYS B 72 " pdb="ZN ZN B 102 " - pdb=" SG CYS B 52 " Number of angles added : 18 3232 Ramachandran restraints generated. 1616 Oldfield, 0 Emsley, 1616 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 3152 Finding SS restraints... Secondary structure from input PDB file: 70 helices and 10 sheets defined 52.7% alpha, 6.0% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.54 Creating SS restraints... Processing helix chain 'A' and resid 3 through 12 removed outlier: 3.751A pdb=" N GLU A 7 " --> pdb=" O GLU A 3 " (cutoff:3.500A) removed outlier: 4.065A pdb=" N LYS A 9 " --> pdb=" O VAL A 5 " (cutoff:3.500A) removed outlier: 3.711A pdb=" N GLY A 10 " --> pdb=" O GLN A 6 " (cutoff:3.500A) Processing helix chain 'A' and resid 19 through 31 removed outlier: 3.700A pdb=" N GLY A 23 " --> pdb=" O ARG A 19 " (cutoff:3.500A) removed outlier: 3.833A pdb=" N THR A 28 " --> pdb=" O HIS A 24 " (cutoff:3.500A) Processing helix chain 'A' and resid 34 through 57 removed outlier: 4.035A pdb=" N LEU A 39 " --> pdb=" O PRO A 35 " (cutoff:3.500A) removed outlier: 4.628A pdb=" N THR A 40 " --> pdb=" O SER A 36 " (cutoff:3.500A) removed outlier: 3.598A pdb=" N GLU A 41 " --> pdb=" O PRO A 37 " (cutoff:3.500A) Processing helix chain 'A' and resid 68 through 76 removed outlier: 4.412A pdb=" N ASN A 76 " --> pdb=" O ILE A 72 " (cutoff:3.500A) Processing helix chain 'A' and resid 108 through 125 removed outlier: 4.347A pdb=" N SER A 121 " --> pdb=" O LYS A 117 " (cutoff:3.500A) removed outlier: 3.893A pdb=" N LYS A 122 " --> pdb=" O ALA A 118 " (cutoff:3.500A) Processing helix chain 'A' and resid 126 through 132 removed outlier: 3.763A pdb=" N LYS A 132 " --> pdb=" O HIS A 128 " (cutoff:3.500A) Processing helix chain 'A' and resid 155 through 171 removed outlier: 4.206A pdb=" N ARG A 159 " --> pdb=" O TRP A 155 " (cutoff:3.500A) Processing helix chain 'A' and resid 182 through 190 removed outlier: 3.733A pdb=" N LEU A 187 " --> pdb=" O GLU A 183 " (cutoff:3.500A) removed outlier: 3.594A pdb=" N ILE A 188 " --> pdb=" O TYR A 184 " (cutoff:3.500A) removed outlier: 3.630A pdb=" N HIS A 189 " --> pdb=" O GLN A 185 " (cutoff:3.500A) Processing helix chain 'A' and resid 214 through 223 Processing helix chain 'A' and resid 232 through 248 removed outlier: 3.844A pdb=" N ASP A 248 " --> pdb=" O LEU A 244 " (cutoff:3.500A) Processing helix chain 'A' and resid 262 through 274 Processing helix chain 'A' and resid 284 through 289 removed outlier: 3.886A pdb=" N ASN A 288 " --> pdb=" O CYS A 284 " (cutoff:3.500A) removed outlier: 3.620A pdb=" N ASN A 289 " --> pdb=" O PRO A 285 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 284 through 289' Processing helix chain 'A' and resid 289 through 303 Processing helix chain 'A' and resid 322 through 338 Processing helix chain 'A' and resid 339 through 341 No H-bonds generated for 'chain 'A' and resid 339 through 341' Processing helix chain 'A' and resid 346 through 365 Processing helix chain 'A' and resid 389 through 391 No H-bonds generated for 'chain 'A' and resid 389 through 391' Processing helix chain 'A' and resid 396 through 405 Processing helix chain 'A' and resid 412 through 421 removed outlier: 3.520A pdb=" N CYS A 417 " --> pdb=" O TRP A 413 " (cutoff:3.500A) removed outlier: 3.617A pdb=" N LEU A 421 " --> pdb=" O CYS A 417 " (cutoff:3.500A) Processing helix chain 'A' and resid 440 through 451 removed outlier: 3.736A pdb=" N LEU A 444 " --> pdb=" O ASP A 440 " (cutoff:3.500A) Processing helix chain 'A' and resid 474 through 478 Processing helix chain 'A' and resid 481 through 501 Processing helix chain 'A' and resid 574 through 578 removed outlier: 3.561A pdb=" N ARG A 577 " --> pdb=" O ASP A 574 " (cutoff:3.500A) removed outlier: 3.503A pdb=" N PHE A 578 " --> pdb=" O PRO A 575 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 574 through 578' Processing helix chain 'A' and resid 586 through 598 Processing helix chain 'A' and resid 605 through 623 Proline residue: A 611 - end of helix Processing helix chain 'A' and resid 625 through 644 removed outlier: 3.665A pdb=" N GLN A 629 " --> pdb=" O SER A 625 " (cutoff:3.500A) removed outlier: 3.551A pdb=" N ASN A 644 " --> pdb=" O MET A 640 " (cutoff:3.500A) Processing helix chain 'A' and resid 663 through 680 removed outlier: 3.926A pdb=" N VAL A 667 " --> pdb=" O SER A 663 " (cutoff:3.500A) Processing helix chain 'A' and resid 688 through 701 Processing helix chain 'A' and resid 702 through 705 Processing helix chain 'A' and resid 712 through 732 Proline residue: A 727 - end of helix Processing helix chain 'A' and resid 745 through 760 Processing helix chain 'A' and resid 776 through 789 removed outlier: 3.824A pdb=" N ARG A 789 " --> pdb=" O SER A 785 " (cutoff:3.500A) Processing helix chain 'A' and resid 833 through 852 Processing helix chain 'A' and resid 863 through 872 removed outlier: 3.586A pdb=" N TYR A 867 " --> pdb=" O ASP A 863 " (cutoff:3.500A) Processing helix chain 'A' and resid 887 through 895 Processing helix chain 'A' and resid 898 through 921 Processing helix chain 'A' and resid 963 through 968 removed outlier: 4.429A pdb=" N VAL A 967 " --> pdb=" O SER A 963 " (cutoff:3.500A) removed outlier: 3.687A pdb=" N LYS A 968 " --> pdb=" O LEU A 964 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 963 through 968' Processing helix chain 'A' and resid 971 through 975 removed outlier: 3.541A pdb=" N LEU A 975 " --> pdb=" O PRO A 972 " (cutoff:3.500A) Processing helix chain 'A' and resid 976 through 980 removed outlier: 3.871A pdb=" N VAL A 979 " --> pdb=" O PRO A 976 " (cutoff:3.500A) Processing helix chain 'A' and resid 981 through 989 removed outlier: 4.158A pdb=" N ASP A 985 " --> pdb=" O LYS A 981 " (cutoff:3.500A) Processing helix chain 'A' and resid 1009 through 1021 removed outlier: 3.756A pdb=" N ASP A1021 " --> pdb=" O SER A1017 " (cutoff:3.500A) Processing helix chain 'A' and resid 1024 through 1035 Processing helix chain 'A' and resid 1135 through 1156 Processing helix chain 'A' and resid 1163 through 1181 Processing helix chain 'A' and resid 1197 through 1208 removed outlier: 3.823A pdb=" N PHE A1201 " --> pdb=" O SER A1197 " (cutoff:3.500A) Processing helix chain 'A' and resid 1221 through 1233 removed outlier: 3.660A pdb=" N HIS A1233 " --> pdb=" O LYS A1229 " (cutoff:3.500A) Processing helix chain 'A' and resid 1239 through 1249 Processing helix chain 'A' and resid 1264 through 1275 removed outlier: 4.679A pdb=" N GLN A1275 " --> pdb=" O HIS A1271 " (cutoff:3.500A) Processing helix chain 'A' and resid 1295 through 1319 Processing helix chain 'A' and resid 1351 through 1365 Processing helix chain 'A' and resid 1395 through 1405 removed outlier: 3.749A pdb=" N PHE A1399 " --> pdb=" O LEU A1395 " (cutoff:3.500A) Processing helix chain 'A' and resid 1411 through 1429 removed outlier: 3.914A pdb=" N ASN A1429 " --> pdb=" O GLN A1425 " (cutoff:3.500A) Processing helix chain 'A' and resid 1435 through 1450 removed outlier: 3.624A pdb=" N SER A1450 " --> pdb=" O LEU A1446 " (cutoff:3.500A) Processing helix chain 'A' and resid 1472 through 1485 removed outlier: 3.744A pdb=" N PHE A1485 " --> pdb=" O ILE A1481 " (cutoff:3.500A) Processing helix chain 'A' and resid 1489 through 1507 Processing helix chain 'A' and resid 1513 through 1520 removed outlier: 3.771A pdb=" N ASN A1517 " --> pdb=" O GLU A1513 " (cutoff:3.500A) removed outlier: 3.615A pdb=" N LEU A1518 " --> pdb=" O ASN A1514 " (cutoff:3.500A) Processing helix chain 'A' and resid 1523 through 1535 Processing helix chain 'A' and resid 1537 through 1545 removed outlier: 3.677A pdb=" N LEU A1545 " --> pdb=" O VAL A1541 " (cutoff:3.500A) Processing helix chain 'A' and resid 1579 through 1590 removed outlier: 3.839A pdb=" N LYS A1583 " --> pdb=" O PRO A1579 " (cutoff:3.500A) Processing helix chain 'A' and resid 1613 through 1625 Processing helix chain 'A' and resid 1639 through 1645 removed outlier: 3.602A pdb=" N VAL A1643 " --> pdb=" O TYR A1639 " (cutoff:3.500A) Processing helix chain 'A' and resid 1664 through 1669 removed outlier: 3.800A pdb=" N ASN A1669 " --> pdb=" O ASN A1666 " (cutoff:3.500A) Processing helix chain 'A' and resid 1678 through 1694 removed outlier: 3.827A pdb=" N GLY A1694 " --> pdb=" O ALA A1690 " (cutoff:3.500A) Processing helix chain 'A' and resid 1695 through 1699 removed outlier: 3.822A pdb=" N PHE A1698 " --> pdb=" O GLU A1695 " (cutoff:3.500A) removed outlier: 3.952A pdb=" N SER A1699 " --> pdb=" O TRP A1696 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 1695 through 1699' Processing helix chain 'A' and resid 1735 through 1744 Processing helix chain 'A' and resid 1744 through 1750 Processing helix chain 'A' and resid 1758 through 1764 removed outlier: 3.977A pdb=" N SER A1762 " --> pdb=" O ASP A1758 " (cutoff:3.500A) Processing helix chain 'A' and resid 1778 through 1783 removed outlier: 3.678A pdb=" N MET A1782 " --> pdb=" O ALA A1778 " (cutoff:3.500A) Processing helix chain 'A' and resid 1788 through 1793 removed outlier: 3.758A pdb=" N LEU A1793 " --> pdb=" O ILE A1789 " (cutoff:3.500A) Processing helix chain 'B' and resid 60 through 69 Processing sheet with id=AA1, first strand: chain 'A' and resid 90 through 94 removed outlier: 6.771A pdb=" N LEU A 98 " --> pdb=" O VAL A 138 " (cutoff:3.500A) removed outlier: 7.687A pdb=" N ILE A 140 " --> pdb=" O LEU A 98 " (cutoff:3.500A) removed outlier: 6.355A pdb=" N LEU A 100 " --> pdb=" O ILE A 140 " (cutoff:3.500A) removed outlier: 7.750A pdb=" N ASP A 142 " --> pdb=" O LEU A 100 " (cutoff:3.500A) removed outlier: 6.583A pdb=" N GLU A 102 " --> pdb=" O ASP A 142 " (cutoff:3.500A) removed outlier: 8.087A pdb=" N ARG A 144 " --> pdb=" O GLU A 102 " (cutoff:3.500A) removed outlier: 6.782A pdb=" N PHE A 104 " --> pdb=" O ARG A 144 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 258 through 260 Processing sheet with id=AA3, first strand: chain 'A' and resid 385 through 387 Processing sheet with id=AA4, first strand: chain 'A' and resid 529 through 534 removed outlier: 6.032A pdb=" N LEU A 559 " --> pdb=" O GLU A 567 " (cutoff:3.500A) removed outlier: 5.565A pdb=" N GLU A 567 " --> pdb=" O LEU A 559 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'A' and resid 1129 through 1132 removed outlier: 3.833A pdb=" N LEU A1119 " --> pdb=" O GLU A1237 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'A' and resid 1379 through 1381 Processing sheet with id=AA7, first strand: chain 'A' and resid 1548 through 1549 Processing sheet with id=AA8, first strand: chain 'A' and resid 1630 through 1631 Processing sheet with id=AA9, first strand: chain 'A' and resid 1716 through 1719 Processing sheet with id=AB1, first strand: chain 'B' and resid 50 through 51 605 hydrogen bonds defined for protein. 1752 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.76 Time building geometry restraints manager: 1.05 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.23 - 1.35: 4235 1.35 - 1.46: 2978 1.46 - 1.58: 5996 1.58 - 1.70: 0 1.70 - 1.82: 132 Bond restraints: 13341 Sorted by residual: bond pdb=" N ILE B 50 " pdb=" CA ILE B 50 " ideal model delta sigma weight residual 1.458 1.497 -0.039 1.20e-02 6.94e+03 1.04e+01 bond pdb=" N LYS B 40 " pdb=" CA LYS B 40 " ideal model delta sigma weight residual 1.460 1.497 -0.037 1.28e-02 6.10e+03 8.40e+00 bond pdb=" N LEU B 57 " pdb=" CA LEU B 57 " ideal model delta sigma weight residual 1.454 1.488 -0.034 1.21e-02 6.83e+03 7.85e+00 bond pdb=" N CYS B 58 " pdb=" CA CYS B 58 " ideal model delta sigma weight residual 1.453 1.489 -0.036 1.31e-02 5.83e+03 7.66e+00 bond pdb=" N ARG B 60 " pdb=" CA ARG B 60 " ideal model delta sigma weight residual 1.458 1.491 -0.034 1.28e-02 6.10e+03 6.86e+00 ... (remaining 13336 not shown) Histogram of bond angle deviations from ideal: 0.00 - 1.88: 17710 1.88 - 3.77: 245 3.77 - 5.65: 40 5.65 - 7.54: 2 7.54 - 9.42: 4 Bond angle restraints: 18001 Sorted by residual: angle pdb=" N ILE A1654 " pdb=" CA ILE A1654 " pdb=" C ILE A1654 " ideal model delta sigma weight residual 113.71 108.31 5.40 9.50e-01 1.11e+00 3.23e+01 angle pdb=" C LYS B 40 " pdb=" CA LYS B 40 " pdb=" CB LYS B 40 " ideal model delta sigma weight residual 109.83 115.39 -5.56 9.90e-01 1.02e+00 3.15e+01 angle pdb=" CA CYS B 39 " pdb=" C CYS B 39 " pdb=" O CYS B 39 " ideal model delta sigma weight residual 120.63 116.18 4.45 1.08e+00 8.57e-01 1.70e+01 angle pdb=" N CYS B 41 " pdb=" CA CYS B 41 " pdb=" C CYS B 41 " ideal model delta sigma weight residual 113.88 108.87 5.01 1.23e+00 6.61e-01 1.66e+01 angle pdb=" O CYS B 39 " pdb=" C CYS B 39 " pdb=" N LYS B 40 " ideal model delta sigma weight residual 122.12 126.03 -3.91 1.06e+00 8.90e-01 1.36e+01 ... (remaining 17996 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.97: 7375 17.97 - 35.95: 572 35.95 - 53.92: 127 53.92 - 71.89: 22 71.89 - 89.87: 13 Dihedral angle restraints: 8109 sinusoidal: 3274 harmonic: 4835 Sorted by residual: dihedral pdb=" CA THR A1753 " pdb=" C THR A1753 " pdb=" N PRO A1754 " pdb=" CA PRO A1754 " ideal model delta harmonic sigma weight residual 180.00 157.28 22.72 0 5.00e+00 4.00e-02 2.07e+01 dihedral pdb=" C LYS B 40 " pdb=" N LYS B 40 " pdb=" CA LYS B 40 " pdb=" CB LYS B 40 " ideal model delta harmonic sigma weight residual -122.60 -133.40 10.80 0 2.50e+00 1.60e-01 1.86e+01 dihedral pdb=" N LYS B 40 " pdb=" C LYS B 40 " pdb=" CA LYS B 40 " pdb=" CB LYS B 40 " ideal model delta harmonic sigma weight residual 122.80 133.35 -10.55 0 2.50e+00 1.60e-01 1.78e+01 ... (remaining 8106 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.069: 1897 0.069 - 0.138: 169 0.138 - 0.207: 6 0.207 - 0.275: 0 0.275 - 0.344: 1 Chirality restraints: 2073 Sorted by residual: chirality pdb=" CA LYS B 40 " pdb=" N LYS B 40 " pdb=" C LYS B 40 " pdb=" CB LYS B 40 " both_signs ideal model delta sigma weight residual False 2.51 2.17 0.34 2.00e-01 2.50e+01 2.96e+00 chirality pdb=" CA ILE B 50 " pdb=" N ILE B 50 " pdb=" C ILE B 50 " pdb=" CB ILE B 50 " both_signs ideal model delta sigma weight residual False 2.43 2.62 -0.18 2.00e-01 2.50e+01 8.48e-01 chirality pdb=" CB VAL A1465 " pdb=" CA VAL A1465 " pdb=" CG1 VAL A1465 " pdb=" CG2 VAL A1465 " both_signs ideal model delta sigma weight residual False -2.63 -2.46 -0.17 2.00e-01 2.50e+01 7.11e-01 ... (remaining 2070 not shown) Planarity restraints: 2274 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C THR A1393 " -0.029 5.00e-02 4.00e+02 4.32e-02 2.99e+00 pdb=" N PRO A1394 " 0.075 5.00e-02 4.00e+02 pdb=" CA PRO A1394 " -0.022 5.00e-02 4.00e+02 pdb=" CD PRO A1394 " -0.024 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CA CYS B 39 " -0.008 2.00e-02 2.50e+03 1.68e-02 2.81e+00 pdb=" C CYS B 39 " 0.029 2.00e-02 2.50e+03 pdb=" O CYS B 39 " -0.011 2.00e-02 2.50e+03 pdb=" N LYS B 40 " -0.010 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C ASN A 737 " 0.026 5.00e-02 4.00e+02 3.92e-02 2.45e+00 pdb=" N PRO A 738 " -0.068 5.00e-02 4.00e+02 pdb=" CA PRO A 738 " 0.020 5.00e-02 4.00e+02 pdb=" CD PRO A 738 " 0.022 5.00e-02 4.00e+02 ... (remaining 2271 not shown) Histogram of nonbonded interaction distances: 2.25 - 2.78: 2992 2.78 - 3.31: 12415 3.31 - 3.84: 19912 3.84 - 4.37: 22789 4.37 - 4.90: 40391 Nonbonded interactions: 98499 Sorted by model distance: nonbonded pdb=" O VAL A 27 " pdb=" OG SER A 31 " model vdw 2.250 3.040 nonbonded pdb=" OD2 ASP A1173 " pdb=" OH TYR A1325 " model vdw 2.277 3.040 nonbonded pdb=" O THR A1377 " pdb=" OG1 THR A1377 " model vdw 2.284 3.040 nonbonded pdb=" ND2 ASN A 788 " pdb=" OE1 GLU A1209 " model vdw 2.300 3.120 nonbonded pdb=" O ASN A 412 " pdb=" ND2 ASN A 416 " model vdw 2.303 3.120 ... (remaining 98494 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.00 ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as 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: 0.820 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.120 Check model and map are aligned: 0.050 Set scattering table: 0.020 Process input model: 12.000 Find NCS groups from input model: 0.040 Set up NCS constraints: 0.010 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.010 Load rotamer database and sin/cos tables:1.050 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 14.120 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7924 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.117 13359 Z= 0.179 Angle : 0.839 38.791 18023 Z= 0.353 Chirality : 0.041 0.344 2073 Planarity : 0.004 0.043 2274 Dihedral : 14.022 89.867 4951 Min Nonbonded Distance : 2.250 Molprobity Statistics. All-atom Clashscore : 7.21 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.60 % Favored : 91.40 % Rotamer: Outliers : 3.56 % Allowed : 10.60 % Favored : 85.84 % Cbeta Deviations : 0.06 % Peptide Plane: Cis-proline : 0.00 % 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: -2.94 (0.19), residues: 1616 helix: -0.54 (0.17), residues: 805 sheet: -1.61 (0.44), residues: 121 loop : -3.86 (0.18), residues: 690 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.005 0.000 ARG A1741 TYR 0.009 0.001 TYR A 572 PHE 0.013 0.001 PHE A 533 TRP 0.009 0.001 TRP A 155 HIS 0.004 0.001 HIS A 647 Details of bonding type rmsd covalent geometry : bond 0.00315 (13341) covalent geometry : angle 0.58106 (18001) SS BOND : bond 0.00690 ( 2) SS BOND : angle 2.25990 ( 4) hydrogen bonds : bond 0.15744 ( 605) hydrogen bonds : angle 5.45285 ( 1752) metal coordination : bond 0.06848 ( 16) metal coordination : angle 19.15173 ( 18) *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3232 Ramachandran restraints generated. 1616 Oldfield, 0 Emsley, 1616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3232 Ramachandran restraints generated. 1616 Oldfield, 0 Emsley, 1616 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 323 residues out of total 1519 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 53 poor density : 270 time to evaluate : 0.516 Fit side-chains REVERT: A 125 ASN cc_start: 0.8330 (t0) cc_final: 0.8063 (t0) REVERT: A 131 GLN cc_start: 0.7639 (tm-30) cc_final: 0.7359 (pp30) REVERT: A 183 GLU cc_start: 0.7437 (OUTLIER) cc_final: 0.7229 (mp0) REVERT: A 204 LEU cc_start: 0.7232 (OUTLIER) cc_final: 0.6648 (pp) REVERT: A 209 ARG cc_start: 0.7371 (ttt-90) cc_final: 0.7121 (ttm-80) REVERT: A 292 MET cc_start: 0.8030 (mtp) cc_final: 0.7796 (mtm) REVERT: A 606 MET cc_start: 0.8275 (ptm) cc_final: 0.8047 (ptt) REVERT: A 838 GLU cc_start: 0.7657 (tp30) cc_final: 0.6932 (tp30) REVERT: A 980 TYR cc_start: 0.6245 (OUTLIER) cc_final: 0.5567 (t80) REVERT: A 1216 ARG cc_start: 0.7644 (mmt90) cc_final: 0.7388 (mmt-90) REVERT: A 1461 GLU cc_start: 0.7530 (pt0) cc_final: 0.7068 (pt0) REVERT: A 1462 ASP cc_start: 0.8158 (OUTLIER) cc_final: 0.7821 (m-30) REVERT: B 60 ARG cc_start: 0.8038 (OUTLIER) cc_final: 0.7589 (ppt-90) outliers start: 53 outliers final: 35 residues processed: 312 average time/residue: 0.1155 time to fit residues: 50.4458 Evaluate side-chains 241 residues out of total 1519 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 40 poor density : 201 time to evaluate : 0.382 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 22 PHE Chi-restraints excluded: chain A residue 105 VAL Chi-restraints excluded: chain A residue 108 ASN Chi-restraints excluded: chain A residue 155 TRP Chi-restraints excluded: chain A residue 183 GLU Chi-restraints excluded: chain A residue 201 LEU Chi-restraints excluded: chain A residue 204 LEU Chi-restraints excluded: chain A residue 224 ILE Chi-restraints excluded: chain A residue 234 THR Chi-restraints excluded: chain A residue 381 ILE Chi-restraints excluded: chain A residue 382 VAL Chi-restraints excluded: chain A residue 421 LEU Chi-restraints excluded: chain A residue 451 LEU Chi-restraints excluded: chain A residue 564 VAL Chi-restraints excluded: chain A residue 667 VAL Chi-restraints excluded: chain A residue 735 VAL Chi-restraints excluded: chain A residue 736 VAL Chi-restraints excluded: chain A residue 772 THR Chi-restraints excluded: chain A residue 865 THR Chi-restraints excluded: chain A residue 892 LEU Chi-restraints excluded: chain A residue 974 VAL Chi-restraints excluded: chain A residue 980 TYR Chi-restraints excluded: chain A residue 1010 LEU Chi-restraints excluded: chain A residue 1162 ASN Chi-restraints excluded: chain A residue 1212 THR Chi-restraints excluded: chain A residue 1352 LEU Chi-restraints excluded: chain A residue 1377 THR Chi-restraints excluded: chain A residue 1397 THR Chi-restraints excluded: chain A residue 1426 CYS Chi-restraints excluded: chain A residue 1462 ASP Chi-restraints excluded: chain A residue 1560 LEU Chi-restraints excluded: chain A residue 1566 LEU Chi-restraints excluded: chain A residue 1578 VAL Chi-restraints excluded: chain A residue 1625 MET Chi-restraints excluded: chain A residue 1654 ILE Chi-restraints excluded: chain A residue 1731 LEU Chi-restraints excluded: chain A residue 1795 THR Chi-restraints excluded: chain B residue 57 LEU Chi-restraints excluded: chain B residue 59 LEU Chi-restraints excluded: chain B residue 60 ARG Rotamers are restrained with sigma=5.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 702, in refine_xyz self.minimization_no_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 755, in minimization_no_ncs ro = mmtbx.refinement.real_space.individual_sites.easy( 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 44, in __init__ self.weight = mmtbx.refinement.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 362, 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 100, in refine self.refined = maptbx.real_space_refinement_simple.lbfgs( File "/net/cci-filer3/home/dcliebschner/04_cryoem/phenix-2.0-5787/lib/python3.9/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-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/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-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.8727 > 50: distance: 49 - 61: 15.823 distance: 61 - 62: 12.877 distance: 62 - 63: 18.406 distance: 62 - 65: 18.867 distance: 63 - 64: 8.075 distance: 63 - 68: 12.049 distance: 65 - 66: 10.085 distance: 65 - 67: 35.128 distance: 68 - 69: 13.316 distance: 69 - 70: 29.409 distance: 69 - 72: 16.740 distance: 70 - 71: 8.941 distance: 70 - 79: 42.597 distance: 72 - 73: 28.152 distance: 73 - 74: 8.176 distance: 73 - 75: 18.604 distance: 74 - 76: 8.968 distance: 75 - 77: 15.882 distance: 76 - 78: 14.767 distance: 77 - 78: 22.091 distance: 79 - 80: 16.609 distance: 80 - 81: 38.032 distance: 80 - 83: 3.706 distance: 81 - 82: 3.579 distance: 81 - 85: 12.224 distance: 85 - 86: 16.207 distance: 86 - 87: 23.888 distance: 86 - 89: 45.886 distance: 87 - 88: 22.260 distance: 87 - 92: 33.094 distance: 89 - 90: 34.049 distance: 89 - 91: 10.201 distance: 92 - 93: 33.208 distance: 93 - 94: 25.605 distance: 93 - 96: 39.501 distance: 94 - 95: 30.165 distance: 94 - 99: 28.714 distance: 96 - 97: 38.859 distance: 99 - 100: 32.590 distance: 100 - 101: 9.628 distance: 100 - 103: 17.828 distance: 101 - 102: 39.246 distance: 101 - 108: 21.731 distance: 103 - 104: 19.418 distance: 104 - 105: 21.621 distance: 105 - 106: 4.028 distance: 106 - 107: 13.591 distance: 108 - 109: 22.711 distance: 109 - 110: 26.785 distance: 109 - 112: 33.319 distance: 110 - 111: 22.089 distance: 110 - 117: 34.199 distance: 112 - 113: 27.455 distance: 113 - 114: 17.417 distance: 114 - 115: 7.424 distance: 114 - 116: 18.974 distance: 117 - 118: 50.557 distance: 117 - 123: 42.978 distance: 118 - 119: 17.933 distance: 118 - 121: 19.654 distance: 119 - 120: 19.859 distance: 119 - 124: 19.411 distance: 121 - 122: 28.030 distance: 122 - 123: 32.899 distance: 124 - 125: 4.457 distance: 125 - 126: 5.923 distance: 125 - 128: 5.972 distance: 126 - 127: 20.484 distance: 128 - 129: 4.119 distance: 129 - 130: 10.009 distance: 130 - 131: 14.652 distance: 130 - 132: 3.233 distance: 133 - 134: 20.196 distance: 134 - 135: 21.696 distance: 134 - 137: 6.691 distance: 135 - 136: 10.744 distance: 135 - 141: 26.469 distance: 137 - 138: 6.688 distance: 138 - 139: 3.166 distance: 138 - 140: 10.393