Starting phenix.real_space_refine on Sun Aug 24 20:39:40 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/7nbx_12262/08_2025/7nbx_12262.cif Found real_map, /net/cci-nas-00/data/ceres_data/7nbx_12262/08_2025/7nbx_12262.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=4.8 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { model { file = "/net/cci-nas-00/data/ceres_data/7nbx_12262/08_2025/7nbx_12262.cif" } default_model = "/net/cci-nas-00/data/ceres_data/7nbx_12262/08_2025/7nbx_12262.cif" real_map_files = "/net/cci-nas-00/data/ceres_data/7nbx_12262/08_2025/7nbx_12262.map" default_real_map = "/net/cci-nas-00/data/ceres_data/7nbx_12262/08_2025/7nbx_12262.map" } resolution = 4.8 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.004 sd= 0.037 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: 5 Type Number sf(0) Gaussians S 32 5.16 5 C 7251 2.51 5 N 1968 2.21 5 O 2257 1.98 5 H 11129 0.53 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 113 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: 22637 Number of models: 1 Model: "" Number of chains: 5 Chain: "A" Number of atoms: 11784 Number of conformers: 1 Conformer: "" Number of residues, atoms: 761, 11784 Classifications: {'peptide': 761} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 27, 'TRANS': 733} Chain breaks: 1 Chain: "B" Number of atoms: 5307 Number of conformers: 1 Conformer: "" Number of residues, atoms: 357, 5307 Classifications: {'peptide': 357} Link IDs: {'PTRANS': 11, 'TRANS': 345} Chain breaks: 1 Chain: "C" Number of atoms: 817 Number of conformers: 1 Conformer: "" Number of residues, atoms: 56, 817 Classifications: {'peptide': 56} Link IDs: {'PTRANS': 6, 'TRANS': 49} Chain: "D" Number of atoms: 3459 Number of conformers: 1 Conformer: "" Number of residues, atoms: 218, 3459 Classifications: {'peptide': 218} Link IDs: {'PTRANS': 12, 'TRANS': 205} Chain: "E" Number of atoms: 1270 Number of conformers: 1 Conformer: "" Number of residues, atoms: 83, 1270 Classifications: {'peptide': 83} Link IDs: {'PTRANS': 5, 'TRANS': 77} Time building chain proxies: 4.14, per 1000 atoms: 0.18 Number of scatterers: 22637 At special positions: 0 Unit cell: (117.7, 141.24, 103.79, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) S 32 16.00 O 2257 8.00 N 1968 7.00 C 7251 6.00 H 11129 1.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=1, symmetry=0 Simple disulfide: pdb=" SG CYS A 435 " - pdb=" SG CYS A 665 " 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: 1.98 Conformation dependent library (CDL) restraints added in 485.3 milliseconds Enol-peptide restraints added in 953.7 nanoseconds 2922 Ramachandran restraints generated. 1461 Oldfield, 0 Emsley, 1461 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2692 Finding SS restraints... Secondary structure from input PDB file: 25 helices and 18 sheets defined 22.0% alpha, 24.5% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.74 Creating SS restraints... Processing helix chain 'A' and resid 38 through 46 removed outlier: 3.713A pdb=" N MET A 46 " --> pdb=" O ALA A 42 " (cutoff:3.500A) Processing helix chain 'A' and resid 57 through 67 Processing helix chain 'A' and resid 106 through 118 Processing helix chain 'A' and resid 129 through 145 Processing helix chain 'A' and resid 192 through 197 Processing helix chain 'A' and resid 215 through 231 removed outlier: 3.634A pdb=" N LEU A 219 " --> pdb=" O GLN A 215 " (cutoff:3.500A) removed outlier: 3.798A pdb=" N ALA A 220 " --> pdb=" O LYS A 216 " (cutoff:3.500A) removed outlier: 3.506A pdb=" N GLU A 224 " --> pdb=" O ALA A 220 " (cutoff:3.500A) Processing helix chain 'A' and resid 280 through 288 removed outlier: 3.633A pdb=" N ILE A 284 " --> pdb=" O HIS A 280 " (cutoff:3.500A) removed outlier: 3.991A pdb=" N GLN A 286 " --> pdb=" O ALA A 282 " (cutoff:3.500A) Processing helix chain 'A' and resid 298 through 315 Processing helix chain 'A' and resid 380 through 393 removed outlier: 3.924A pdb=" N LEU A 392 " --> pdb=" O ARG A 388 " (cutoff:3.500A) removed outlier: 4.041A pdb=" N GLY A 393 " --> pdb=" O LEU A 389 " (cutoff:3.500A) Processing helix chain 'A' and resid 448 through 452 Processing helix chain 'A' and resid 542 through 553 removed outlier: 4.414A pdb=" N TRP A 546 " --> pdb=" O GLN A 542 " (cutoff:3.500A) removed outlier: 4.019A pdb=" N ARG A 547 " --> pdb=" O VAL A 543 " (cutoff:3.500A) removed outlier: 3.869A pdb=" N TYR A 548 " --> pdb=" O ALA A 544 " (cutoff:3.500A) Processing helix chain 'A' and resid 726 through 731 Processing helix chain 'D' and resid 31 through 45 removed outlier: 3.914A pdb=" N ILE D 35 " --> pdb=" O PRO D 31 " (cutoff:3.500A) Processing helix chain 'D' and resid 48 through 62 removed outlier: 3.975A pdb=" N THR D 53 " --> pdb=" O ARG D 49 " (cutoff:3.500A) Processing helix chain 'D' and resid 67 through 83 Processing helix chain 'D' and resid 84 through 99 Processing helix chain 'D' and resid 104 through 121 removed outlier: 3.674A pdb=" N VAL D 108 " --> pdb=" O ASN D 104 " (cutoff:3.500A) Processing helix chain 'D' and resid 137 through 154 Processing helix chain 'D' and resid 162 through 189 removed outlier: 3.538A pdb=" N GLU D 178 " --> pdb=" O LEU D 174 " (cutoff:3.500A) Processing helix chain 'D' and resid 192 through 205 removed outlier: 3.875A pdb=" N GLY D 200 " --> pdb=" O ASN D 196 " (cutoff:3.500A) Processing helix chain 'D' and resid 208 through 213 removed outlier: 3.773A pdb=" N ARG D 212 " --> pdb=" O THR D 208 " (cutoff:3.500A) Processing helix chain 'D' and resid 214 through 225 removed outlier: 3.717A pdb=" N MET D 218 " --> pdb=" O ALA D 214 " (cutoff:3.500A) removed outlier: 4.012A pdb=" N ARG D 223 " --> pdb=" O GLU D 219 " (cutoff:3.500A) Processing helix chain 'D' and resid 230 through 241 Processing helix chain 'E' and resid 41 through 46 removed outlier: 3.669A pdb=" N ILE E 46 " --> pdb=" O VAL E 43 " (cutoff:3.500A) Processing helix chain 'E' and resid 51 through 60 removed outlier: 3.641A pdb=" N VAL E 55 " --> pdb=" O THR E 51 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 28 through 32 Processing sheet with id=AA2, first strand: chain 'A' and resid 93 through 100 removed outlier: 6.064A pdb=" N VAL A 163 " --> pdb=" O ILE A 94 " (cutoff:3.500A) removed outlier: 4.165A pdb=" N SER A 96 " --> pdb=" O VAL A 163 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 176 through 180 removed outlier: 8.018A pdb=" N VAL A 256 " --> pdb=" O GLN A 178 " (cutoff:3.500A) removed outlier: 4.779A pdb=" N ILE A 180 " --> pdb=" O VAL A 256 " (cutoff:3.500A) removed outlier: 7.902A pdb=" N VAL A 258 " --> pdb=" O ILE A 180 " (cutoff:3.500A) removed outlier: 7.138A pdb=" N GLY A 253 " --> pdb=" O LEU A 247 " (cutoff:3.500A) removed outlier: 4.770A pdb=" N LEU A 247 " --> pdb=" O GLY A 253 " (cutoff:3.500A) removed outlier: 6.839A pdb=" N TYR A 255 " --> pdb=" O VAL A 245 " (cutoff:3.500A) removed outlier: 4.509A pdb=" N VAL A 245 " --> pdb=" O TYR A 255 " (cutoff:3.500A) removed outlier: 6.880A pdb=" N THR A 257 " --> pdb=" O THR A 243 " (cutoff:3.500A) removed outlier: 4.853A pdb=" N THR A 243 " --> pdb=" O THR A 257 " (cutoff:3.500A) removed outlier: 6.520A pdb=" N ASN A 259 " --> pdb=" O ASP A 241 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'A' and resid 271 through 272 Processing sheet with id=AA5, first strand: chain 'A' and resid 351 through 355 removed outlier: 6.776A pdb=" N LYS A 351 " --> pdb=" O VAL A 414 " (cutoff:3.500A) removed outlier: 8.339A pdb=" N TYR A 416 " --> pdb=" O LYS A 351 " (cutoff:3.500A) removed outlier: 6.604A pdb=" N ARG A 353 " --> pdb=" O TYR A 416 " (cutoff:3.500A) removed outlier: 8.084A pdb=" N VAL A 418 " --> pdb=" O ARG A 353 " (cutoff:3.500A) removed outlier: 6.419A pdb=" N GLU A 355 " --> pdb=" O VAL A 418 " (cutoff:3.500A) removed outlier: 7.122A pdb=" N GLU A 420 " --> pdb=" O GLU A 355 " (cutoff:3.500A) removed outlier: 7.487A pdb=" N ASP A 413 " --> pdb=" O THR A 402 " (cutoff:3.500A) removed outlier: 5.468A pdb=" N THR A 402 " --> pdb=" O ASP A 413 " (cutoff:3.500A) removed outlier: 6.866A pdb=" N VAL A 415 " --> pdb=" O THR A 400 " (cutoff:3.500A) removed outlier: 4.780A pdb=" N THR A 400 " --> pdb=" O VAL A 415 " (cutoff:3.500A) removed outlier: 6.757A pdb=" N LYS A 417 " --> pdb=" O VAL A 398 " (cutoff:3.500A) removed outlier: 4.652A pdb=" N VAL A 398 " --> pdb=" O LYS A 417 " (cutoff:3.500A) removed outlier: 6.425A pdb=" N LYS A 419 " --> pdb=" O GLU A 396 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'A' and resid 427 through 429 removed outlier: 3.943A pdb=" N GLY A 443 " --> pdb=" O ASN A 427 " (cutoff:3.500A) removed outlier: 5.938A pdb=" N VAL A 483 " --> pdb=" O TYR A 477 " (cutoff:3.500A) removed outlier: 5.702A pdb=" N TYR A 608 " --> pdb=" O THR A 600 " (cutoff:3.500A) removed outlier: 3.940A pdb=" N VAL A 628 " --> pdb=" O ILE A 719 " (cutoff:3.500A) removed outlier: 4.255A pdb=" N LEU A 630 " --> pdb=" O GLU A 717 " (cutoff:3.500A) removed outlier: 3.670A pdb=" N GLY A 636 " --> pdb=" O MET A 711 " (cutoff:3.500A) removed outlier: 3.778A pdb=" N THR A 735 " --> pdb=" O THR A 720 " (cutoff:3.500A) removed outlier: 6.256A pdb=" N MET A 768 " --> pdb=" O PHE A 791 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'A' and resid 671 through 672 removed outlier: 3.966A pdb=" N ALA A 672 " --> pdb=" O ASP A 703 " (cutoff:3.500A) removed outlier: 3.888A pdb=" N ASP A 703 " --> pdb=" O ALA A 672 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA7 Processing sheet with id=AA8, first strand: chain 'B' and resid 49 through 51 removed outlier: 4.530A pdb=" N LEU B 378 " --> pdb=" O ILE B 390 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'B' and resid 66 through 68 Processing sheet with id=AB1, first strand: chain 'B' and resid 80 through 81 Processing sheet with id=AB2, first strand: chain 'B' and resid 117 through 118 removed outlier: 6.688A pdb=" N VAL B 131 " --> pdb=" O GLN B 144 " (cutoff:3.500A) removed outlier: 4.499A pdb=" N GLN B 144 " --> pdb=" O VAL B 131 " (cutoff:3.500A) removed outlier: 6.852A pdb=" N ALA B 133 " --> pdb=" O ALA B 142 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'B' and resid 162 through 164 removed outlier: 6.413A pdb=" N LEU B 171 " --> pdb=" O THR B 184 " (cutoff:3.500A) removed outlier: 4.329A pdb=" N THR B 184 " --> pdb=" O LEU B 171 " (cutoff:3.500A) removed outlier: 7.087A pdb=" N ALA B 173 " --> pdb=" O LYS B 182 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'B' and resid 201 through 203 Processing sheet with id=AB5, first strand: chain 'B' and resid 252 through 253 removed outlier: 3.784A pdb=" N TRP B 279 " --> pdb=" O ALA B 269 " (cutoff:3.500A) removed outlier: 6.564A pdb=" N ASP B 271 " --> pdb=" O ILE B 277 " (cutoff:3.500A) removed outlier: 7.421A pdb=" N ILE B 277 " --> pdb=" O ASP B 271 " (cutoff:3.500A) Processing sheet with id=AB6, first strand: chain 'B' and resid 295 through 298 removed outlier: 3.537A pdb=" N LEU B 298 " --> pdb=" O MET B 306 " (cutoff:3.500A) removed outlier: 3.542A pdb=" N TRP B 317 " --> pdb=" O ALA B 307 " (cutoff:3.500A) Processing sheet with id=AB7, first strand: chain 'B' and resid 331 through 332 removed outlier: 3.927A pdb=" N VAL B 331 " --> pdb=" O VAL B 338 " (cutoff:3.500A) removed outlier: 3.868A pdb=" N VAL B 338 " --> pdb=" O VAL B 331 " (cutoff:3.500A) No H-bonds generated for sheet with id=AB7 Processing sheet with id=AB8, first strand: chain 'B' and resid 345 through 348 removed outlier: 3.757A pdb=" N GLN B 360 " --> pdb=" O LEU B 346 " (cutoff:3.500A) removed outlier: 3.830A pdb=" N ALA B 358 " --> pdb=" O TRP B 348 " (cutoff:3.500A) Processing sheet with id=AB9, first strand: chain 'E' and resid 33 through 35 456 hydrogen bonds defined for protein. 1272 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: 5.04 Time building geometry restraints manager: 2.70 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 0.83 - 1.03: 11129 1.03 - 1.22: 4 1.22 - 1.42: 5068 1.42 - 1.62: 6621 1.62 - 1.82: 62 Bond restraints: 22884 Sorted by residual: bond pdb=" N PHE A 517 " pdb=" CA PHE A 517 " ideal model delta sigma weight residual 1.456 1.430 0.026 1.04e-02 9.25e+03 6.10e+00 bond pdb=" CA ASN E 106 " pdb=" CB ASN E 106 " ideal model delta sigma weight residual 1.527 1.496 0.030 1.27e-02 6.20e+03 5.74e+00 bond pdb=" CA ASN A 709 " pdb=" CB ASN A 709 " ideal model delta sigma weight residual 1.536 1.505 0.031 1.46e-02 4.69e+03 4.62e+00 bond pdb=" C GLY D 65 " pdb=" N PRO D 66 " ideal model delta sigma weight residual 1.337 1.313 0.024 1.24e-02 6.50e+03 3.74e+00 bond pdb=" CB TRP B 143 " pdb=" CG TRP B 143 " ideal model delta sigma weight residual 1.498 1.440 0.058 3.10e-02 1.04e+03 3.45e+00 ... (remaining 22879 not shown) Histogram of bond angle deviations from ideal: 0.00 - 2.02: 39467 2.02 - 4.03: 1580 4.03 - 6.05: 108 6.05 - 8.07: 11 8.07 - 10.09: 3 Bond angle restraints: 41169 Sorted by residual: angle pdb=" N GLY B 211 " pdb=" CA GLY B 211 " pdb=" C GLY B 211 " ideal model delta sigma weight residual 111.85 104.09 7.76 1.06e+00 8.90e-01 5.37e+01 angle pdb=" C LEU B 152 " pdb=" N SER B 153 " pdb=" CA SER B 153 " ideal model delta sigma weight residual 121.54 131.63 -10.09 1.91e+00 2.74e-01 2.79e+01 angle pdb=" N VAL D 192 " pdb=" CA VAL D 192 " pdb=" CB VAL D 192 " ideal model delta sigma weight residual 111.64 107.99 3.65 8.20e-01 1.49e+00 1.98e+01 angle pdb=" N ASN A 448 " pdb=" CA ASN A 448 " pdb=" C ASN A 448 " ideal model delta sigma weight residual 110.44 105.23 5.21 1.20e+00 6.94e-01 1.89e+01 angle pdb=" N ASN E 41 " pdb=" CA ASN E 41 " pdb=" C ASN E 41 " ideal model delta sigma weight residual 114.56 109.19 5.37 1.27e+00 6.20e-01 1.79e+01 ... (remaining 41164 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.75: 9307 17.75 - 35.49: 1116 35.49 - 53.24: 279 53.24 - 70.99: 112 70.99 - 88.73: 34 Dihedral angle restraints: 10848 sinusoidal: 5851 harmonic: 4997 Sorted by residual: dihedral pdb=" CA SER A 557 " pdb=" C SER A 557 " pdb=" N THR A 558 " pdb=" CA THR A 558 " ideal model delta harmonic sigma weight residual 180.00 130.70 49.30 0 5.00e+00 4.00e-02 9.72e+01 dihedral pdb=" CA LEU B 152 " pdb=" C LEU B 152 " pdb=" N SER B 153 " pdb=" CA SER B 153 " ideal model delta harmonic sigma weight residual 180.00 -135.51 -44.49 0 5.00e+00 4.00e-02 7.92e+01 dihedral pdb=" CA PRO B 371 " pdb=" C PRO B 371 " pdb=" N VAL B 372 " pdb=" CA VAL B 372 " ideal model delta harmonic sigma weight residual 180.00 146.12 33.88 0 5.00e+00 4.00e-02 4.59e+01 ... (remaining 10845 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.048: 1219 0.048 - 0.095: 413 0.095 - 0.143: 103 0.143 - 0.191: 13 0.191 - 0.238: 2 Chirality restraints: 1750 Sorted by residual: chirality pdb=" CG LEU A 450 " pdb=" CB LEU A 450 " pdb=" CD1 LEU A 450 " pdb=" CD2 LEU A 450 " both_signs ideal model delta sigma weight residual False -2.59 -2.35 -0.24 2.00e-01 2.50e+01 1.42e+00 chirality pdb=" CB THR A 505 " pdb=" CA THR A 505 " pdb=" OG1 THR A 505 " pdb=" CG2 THR A 505 " both_signs ideal model delta sigma weight residual False 2.55 2.36 0.20 2.00e-01 2.50e+01 9.58e-01 chirality pdb=" CG LEU C 58 " pdb=" CB LEU C 58 " pdb=" CD1 LEU C 58 " pdb=" CD2 LEU C 58 " both_signs ideal model delta sigma weight residual False -2.59 -2.40 -0.19 2.00e-01 2.50e+01 9.02e-01 ... (remaining 1747 not shown) Planarity restraints: 3509 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C GLU B 370 " 0.057 5.00e-02 4.00e+02 8.52e-02 1.16e+01 pdb=" N PRO B 371 " -0.147 5.00e-02 4.00e+02 pdb=" CA PRO B 371 " 0.044 5.00e-02 4.00e+02 pdb=" CD PRO B 371 " 0.046 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C HIS A 555 " -0.053 5.00e-02 4.00e+02 8.04e-02 1.03e+01 pdb=" N PRO A 556 " 0.139 5.00e-02 4.00e+02 pdb=" CA PRO A 556 " -0.043 5.00e-02 4.00e+02 pdb=" CD PRO A 556 " -0.043 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CA ALA B 374 " -0.012 2.00e-02 2.50e+03 2.54e-02 6.46e+00 pdb=" C ALA B 374 " 0.044 2.00e-02 2.50e+03 pdb=" O ALA B 374 " -0.017 2.00e-02 2.50e+03 pdb=" N ASP B 375 " -0.015 2.00e-02 2.50e+03 ... (remaining 3506 not shown) Histogram of nonbonded interaction distances: 1.63 - 2.22: 2542 2.22 - 2.82: 50800 2.82 - 3.41: 60149 3.41 - 4.01: 78690 4.01 - 4.60: 115940 Nonbonded interactions: 308121 Sorted by model distance: nonbonded pdb=" H SER A 269 " pdb=" O VAL A 335 " model vdw 1.626 2.450 nonbonded pdb=" O MET A 539 " pdb="HD21 ASN A 563 " model vdw 1.637 2.450 nonbonded pdb=" O THR A 190 " pdb=" H ILE A 194 " model vdw 1.658 2.450 nonbonded pdb=" OD2 ASP A 410 " pdb=" H THR E 61 " model vdw 1.658 2.450 nonbonded pdb=" O SER C 73 " pdb=" HG1 THR D 161 " model vdw 1.661 2.450 ... (remaining 308116 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=0.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: 2.430 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.100 Extract box with map and model: 0.310 Check model and map are aligned: 0.070 Set scattering table: 0.070 Process input model: 20.490 Find NCS groups from input model: 0.190 Set up NCS constraints: 0.030 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: 24.750 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8249 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.008 0.062 11756 Z= 0.402 Angle : 1.043 10.086 15978 Z= 0.619 Chirality : 0.050 0.238 1750 Planarity : 0.006 0.085 2102 Dihedral : 15.897 88.734 4252 Min Nonbonded Distance : 1.950 Molprobity Statistics. All-atom Clashscore : 12.59 Ramachandran Plot: Outliers : 0.14 % Allowed : 12.25 % Favored : 87.61 % Rotamer: Outliers : 0.00 % Allowed : 4.21 % Favored : 95.79 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.28 % 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.49 (0.18), residues: 1461 helix: -2.41 (0.22), residues: 283 sheet: -2.37 (0.22), residues: 378 loop : -3.62 (0.19), residues: 800 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.007 0.001 ARG D 173 TYR 0.020 0.003 TYR A 255 PHE 0.024 0.002 PHE A 565 TRP 0.022 0.003 TRP A 376 HIS 0.011 0.003 HIS A 533 Details of bonding type rmsd covalent geometry : bond 0.00787 (11755) covalent geometry : angle 1.04290 (15976) SS BOND : bond 0.00045 ( 1) SS BOND : angle 1.50766 ( 2) hydrogen bonds : bond 0.15088 ( 456) hydrogen bonds : angle 8.56883 ( 1272) *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2922 Ramachandran restraints generated. 1461 Oldfield, 0 Emsley, 1461 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2922 Ramachandran restraints generated. 1461 Oldfield, 0 Emsley, 1461 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 153 residues out of total 1234 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 153 time to evaluate : 0.434 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: symmetry clash 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: symmetry clash 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: symmetry clash 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: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 241 ASP cc_start: 0.9127 (m-30) cc_final: 0.8801 (m-30) REVERT: A 503 ASP cc_start: 0.8686 (t0) cc_final: 0.8366 (m-30) REVERT: A 523 ASN cc_start: 0.8839 (m110) cc_final: 0.8395 (m-40) REVERT: A 652 PHE cc_start: 0.8008 (m-80) cc_final: 0.7754 (m-10) REVERT: A 671 LYS cc_start: 0.8882 (pttm) cc_final: 0.8488 (ptpp) REVERT: B 165 HIS cc_start: 0.8274 (t-90) cc_final: 0.7677 (t-90) REVERT: B 367 PHE cc_start: 0.8991 (m-80) cc_final: 0.8713 (m-80) REVERT: D 34 GLU cc_start: 0.9492 (mp0) cc_final: 0.9188 (pm20) REVERT: D 109 MET cc_start: 0.9241 (mmm) cc_final: 0.8814 (mtp) REVERT: D 148 SER cc_start: 0.9580 (m) cc_final: 0.9369 (p) REVERT: D 183 GLU cc_start: 0.9061 (mm-30) cc_final: 0.8840 (tt0) REVERT: D 219 GLU cc_start: 0.8814 (tp30) cc_final: 0.8599 (tp30) REVERT: E 36 ASN cc_start: 0.9345 (m-40) cc_final: 0.8705 (t0) REVERT: E 50 MET cc_start: 0.8702 (mtp) cc_final: 0.8082 (mtp) REVERT: E 107 LYS cc_start: 0.9239 (mmtp) cc_final: 0.8963 (mmtp) outliers start: 0 outliers final: 0 residues processed: 153 average time/residue: 0.1805 time to fit residues: 41.1151 Evaluate side-chains 113 residues out of total 1234 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 113 time to evaluate : 0.514 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 145 random chunks: chunk 98 optimal weight: 0.9980 chunk 107 optimal weight: 4.9990 chunk 10 optimal weight: 6.9990 chunk 66 optimal weight: 4.9990 chunk 130 optimal weight: 1.9990 chunk 124 optimal weight: 0.9980 chunk 103 optimal weight: 0.9990 chunk 77 optimal weight: 4.9990 chunk 122 optimal weight: 3.9990 chunk 91 optimal weight: 0.9990 chunk 55 optimal weight: 2.9990 overall best weight: 1.1986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 161 ASN A 170 GLN A 286 GLN A 323 GLN A 411 GLN A 441 GLN A 523 ASN A 540 GLN A 561 GLN A 801 GLN A 805 ASN B 40 ASN B 70 ASN B 85 ASN B 130 GLN B 144 GLN B 170 GLN B 319 GLN B 360 GLN ** B 368 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** D 69 GLN D 82 ASN E 36 ASN E 53 GLN E 54 GLN E 103 ASN Total number of N/Q/H flips: 25 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3207 r_free = 0.3207 target = 0.065020 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 53)----------------| | r_work = 0.2786 r_free = 0.2786 target = 0.048154 restraints weight = 171881.876| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 36)----------------| | r_work = 0.2861 r_free = 0.2861 target = 0.051116 restraints weight = 78669.193| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 41)----------------| | r_work = 0.2910 r_free = 0.2910 target = 0.053103 restraints weight = 45110.355| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 4; iterations = 35)----------------| | r_work = 0.2943 r_free = 0.2943 target = 0.054421 restraints weight = 29791.906| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 5; iterations = 48)----------------| | r_work = 0.2966 r_free = 0.2966 target = 0.055345 restraints weight = 22175.971| |-----------------------------------------------------------------------------| r_work (final): 0.2971 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8177 moved from start: 0.1804 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.036 11756 Z= 0.153 Angle : 0.674 6.125 15978 Z= 0.371 Chirality : 0.044 0.155 1750 Planarity : 0.005 0.071 2102 Dihedral : 7.491 44.601 1615 Min Nonbonded Distance : 2.281 Molprobity Statistics. All-atom Clashscore : 8.22 Ramachandran Plot: Outliers : 0.14 % Allowed : 10.75 % Favored : 89.12 % Rotamer: Outliers : 0.00 % Allowed : 3.65 % Favored : 96.35 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.14 % 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.56 (0.19), residues: 1461 helix: -1.00 (0.26), residues: 300 sheet: -2.04 (0.22), residues: 383 loop : -3.19 (0.20), residues: 778 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.011 0.000 ARG B 355 TYR 0.015 0.002 TYR D 177 PHE 0.016 0.001 PHE D 95 TRP 0.016 0.002 TRP A 546 HIS 0.005 0.001 HIS B 347 Details of bonding type rmsd covalent geometry : bond 0.00313 (11755) covalent geometry : angle 0.67391 (15976) SS BOND : bond 0.00034 ( 1) SS BOND : angle 1.20743 ( 2) hydrogen bonds : bond 0.04715 ( 456) hydrogen bonds : angle 7.05501 ( 1272) *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2922 Ramachandran restraints generated. 1461 Oldfield, 0 Emsley, 1461 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2922 Ramachandran restraints generated. 1461 Oldfield, 0 Emsley, 1461 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 141 residues out of total 1234 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 141 time to evaluate : 0.395 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: symmetry clash 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: symmetry clash 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: symmetry clash 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: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 369 MET cc_start: 0.9245 (mmt) cc_final: 0.9034 (mmt) REVERT: A 387 GLU cc_start: 0.9042 (mt-10) cc_final: 0.8778 (mt-10) REVERT: A 503 ASP cc_start: 0.8624 (t0) cc_final: 0.8411 (m-30) REVERT: B 165 HIS cc_start: 0.8272 (t-90) cc_final: 0.7657 (t-90) REVERT: B 369 THR cc_start: 0.9788 (p) cc_final: 0.9543 (t) REVERT: D 34 GLU cc_start: 0.9508 (mp0) cc_final: 0.9186 (pm20) REVERT: D 178 GLU cc_start: 0.9205 (mm-30) cc_final: 0.8777 (mm-30) REVERT: D 201 MET cc_start: 0.9257 (mtp) cc_final: 0.8795 (mtp) REVERT: D 219 GLU cc_start: 0.8760 (tp30) cc_final: 0.8476 (tp30) REVERT: E 50 MET cc_start: 0.8772 (mtp) cc_final: 0.8247 (mtp) outliers start: 0 outliers final: 0 residues processed: 141 average time/residue: 0.1779 time to fit residues: 38.6011 Evaluate side-chains 115 residues out of total 1234 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 115 time to evaluate : 0.429 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 145 random chunks: chunk 138 optimal weight: 2.9990 chunk 142 optimal weight: 4.9990 chunk 40 optimal weight: 9.9990 chunk 101 optimal weight: 6.9990 chunk 125 optimal weight: 7.9990 chunk 24 optimal weight: 9.9990 chunk 29 optimal weight: 4.9990 chunk 65 optimal weight: 3.9990 chunk 73 optimal weight: 0.7980 chunk 136 optimal weight: 5.9990 chunk 72 optimal weight: 4.9990 overall best weight: 3.5588 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 170 GLN A 403 GLN A 506 ASN B 324 HIS ** B 347 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** D 40 GLN D 60 ASN E 52 GLN E 83 HIS E 96 ASN Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3130 r_free = 0.3130 target = 0.061877 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 50)----------------| | r_work = 0.2694 r_free = 0.2694 target = 0.044927 restraints weight = 170378.359| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 35)----------------| | r_work = 0.2767 r_free = 0.2767 target = 0.047675 restraints weight = 79898.778| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 30)----------------| | r_work = 0.2817 r_free = 0.2817 target = 0.049552 restraints weight = 46728.486| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 4; iterations = 48)----------------| | r_work = 0.2851 r_free = 0.2851 target = 0.050856 restraints weight = 31466.584| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 5; iterations = 32)----------------| | r_work = 0.2872 r_free = 0.2872 target = 0.051692 restraints weight = 23377.262| |-----------------------------------------------------------------------------| r_work (final): 0.2874 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8257 moved from start: 0.2505 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.093 11756 Z= 0.292 Angle : 0.719 7.115 15978 Z= 0.400 Chirality : 0.043 0.175 1750 Planarity : 0.005 0.072 2102 Dihedral : 7.037 41.978 1615 Min Nonbonded Distance : 2.294 Molprobity Statistics. All-atom Clashscore : 12.33 Ramachandran Plot: Outliers : 0.14 % Allowed : 12.59 % Favored : 87.27 % Rotamer: Outliers : 0.00 % Allowed : 3.32 % Favored : 96.68 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.28 % 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.56 (0.20), residues: 1461 helix: -0.83 (0.27), residues: 299 sheet: -2.19 (0.22), residues: 430 loop : -3.19 (0.21), residues: 732 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.008 0.001 ARG D 212 TYR 0.018 0.002 TYR A 504 PHE 0.020 0.002 PHE D 95 TRP 0.017 0.002 TRP A 376 HIS 0.009 0.002 HIS B 347 Details of bonding type rmsd covalent geometry : bond 0.00559 (11755) covalent geometry : angle 0.71870 (15976) SS BOND : bond 0.00133 ( 1) SS BOND : angle 1.35928 ( 2) hydrogen bonds : bond 0.04394 ( 456) hydrogen bonds : angle 6.92218 ( 1272) *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2922 Ramachandran restraints generated. 1461 Oldfield, 0 Emsley, 1461 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2922 Ramachandran restraints generated. 1461 Oldfield, 0 Emsley, 1461 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 127 residues out of total 1234 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 127 time to evaluate : 0.580 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: symmetry clash 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: symmetry clash 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: symmetry clash 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: symmetry clash REVERT: A 369 MET cc_start: 0.9496 (mmt) cc_final: 0.9126 (mmt) REVERT: A 503 ASP cc_start: 0.8755 (t0) cc_final: 0.8380 (m-30) REVERT: A 671 LYS cc_start: 0.8721 (pttm) cc_final: 0.8223 (ptpp) REVERT: B 165 HIS cc_start: 0.8212 (t-90) cc_final: 0.7879 (t70) REVERT: C 81 ASP cc_start: 0.8547 (t0) cc_final: 0.8330 (t0) REVERT: D 34 GLU cc_start: 0.9480 (mp0) cc_final: 0.9158 (pm20) REVERT: D 178 GLU cc_start: 0.9207 (mm-30) cc_final: 0.8960 (mm-30) REVERT: D 219 GLU cc_start: 0.8748 (tp30) cc_final: 0.8148 (tp30) REVERT: E 36 ASN cc_start: 0.9374 (m110) cc_final: 0.8588 (m-40) REVERT: E 50 MET cc_start: 0.9039 (mtp) cc_final: 0.8524 (mtp) outliers start: 0 outliers final: 0 residues processed: 127 average time/residue: 0.1857 time to fit residues: 36.3670 Evaluate side-chains 110 residues out of total 1234 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 110 time to evaluate : 0.621 ------------------------------------------------------------------------------- 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 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: 71.2527 > 50: distance: 55 - 68: 34.630 distance: 68 - 69: 41.050 distance: 68 - 75: 39.352 distance: 69 - 70: 39.772 distance: 69 - 72: 57.900 distance: 69 - 76: 40.579 distance: 70 - 71: 39.123 distance: 70 - 84: 22.801 distance: 72 - 73: 68.870 distance: 72 - 74: 69.866 distance: 72 - 77: 45.640 distance: 73 - 78: 39.402 distance: 73 - 79: 40.921 distance: 73 - 80: 39.806 distance: 74 - 81: 39.827 distance: 74 - 82: 39.322 distance: 74 - 83: 40.010 distance: 84 - 85: 17.670 distance: 84 - 90: 43.897 distance: 85 - 86: 24.640 distance: 85 - 91: 4.954 distance: 86 - 87: 16.331 distance: 86 - 98: 37.042 distance: 88 - 89: 39.949 distance: 88 - 92: 40.728 distance: 88 - 93: 11.553 distance: 89 - 95: 40.418 distance: 90 - 97: 40.703 distance: 98 - 99: 65.554 distance: 98 - 112: 40.792 distance: 99 - 100: 13.472 distance: 99 - 102: 34.352 distance: 99 - 113: 53.016 distance: 100 - 101: 29.723 distance: 100 - 122: 59.625 distance: 104 - 106: 33.746 distance: 105 - 107: 28.335 distance: 105 - 108: 9.170 distance: 106 - 107: 37.563 distance: 106 - 117: 27.097 distance: 107 - 109: 36.205 distance: 108 - 110: 14.820 distance: 108 - 118: 10.807 distance: 109 - 111: 33.624 distance: 109 - 119: 23.826 distance: 110 - 111: 24.164 distance: 110 - 120: 24.328 distance: 111 - 121: 36.761 distance: 122 - 123: 13.600 distance: 122 - 136: 27.781 distance: 123 - 124: 18.690 distance: 123 - 126: 39.639 distance: 123 - 137: 14.160 distance: 124 - 125: 19.065 distance: 124 - 146: 18.706 distance: 126 - 127: 10.958 distance: 126 - 138: 5.803 distance: 127 - 128: 30.151 distance: 127 - 129: 10.421 distance: 128 - 130: 14.953 distance: 128 - 140: 31.136 distance: 129 - 131: 27.425 distance: 129 - 132: 9.892 distance: 130 - 131: 16.304 distance: 130 - 141: 35.895 distance: 131 - 133: 15.053 distance: 132 - 134: 9.317 distance: 132 - 142: 12.717 distance: 133 - 135: 15.159 distance: 133 - 143: 17.238 distance: 134 - 135: 37.112 distance: 134 - 144: 10.295 distance: 135 - 145: 7.083 distance: 146 - 147: 40.248 distance: 146 - 154: 40.244 distance: 147 - 150: 8.791 distance: 148 - 149: 56.511 distance: 150 - 151: 32.201 distance: 150 - 156: 31.877 distance: 150 - 157: 30.666 distance: 151 - 152: 5.529 distance: 151 - 153: 27.640 distance: 153 - 158: 10.536 distance: 153 - 159: 40.056 distance: 160 - 161: 40.698 distance: 161 - 162: 16.872 distance: 161 - 164: 56.551 distance: 162 - 163: 43.556 distance: 162 - 176: 23.923 distance: 164 - 166: 28.941 distance: 166 - 173: 9.558 distance: 166 - 174: 10.227 distance: 166 - 175: 30.776 distance: 176 - 177: 39.740 distance: 176 - 183: 40.113 distance: 177 - 178: 18.801 distance: 177 - 180: 39.798 distance: 177 - 184: 40.849 distance: 178 - 179: 35.239 distance: 178 - 192: 36.648 distance: 180 - 181: 40.028 distance: 180 - 182: 38.681 distance: 181 - 188: 14.092 distance: 182 - 191: 25.696 distance: 192 - 193: 41.263 distance: 192 - 196: 56.370 distance: 193 - 194: 57.654 distance: 193 - 197: 39.611 distance: 193 - 198: 31.770 distance: 194 - 195: 20.413 distance: 194 - 199: 56.432 distance: 199 - 200: 4.471 distance: 199 - 207: 40.700 distance: 200 - 201: 41.488 distance: 200 - 203: 40.342 distance: 200 - 208: 13.846 distance: 201 - 202: 39.163 distance: 201 - 211: 41.563 distance: 203 - 204: 26.221 distance: 203 - 209: 49.081 distance: 203 - 210: 56.737 distance: 204 - 205: 35.210 distance: 204 - 206: 47.709 distance: 211 - 212: 30.254 distance: 211 - 222: 29.473 distance: 213 - 235: 34.580 distance: 216 - 226: 27.905 distance: 217 - 229: 7.012 distance: 218 - 219: 13.915 distance: 218 - 230: 39.787 distance: 219 - 220: 24.524 distance: 219 - 221: 6.826 distance: 220 - 231: 3.927 distance: 220 - 232: 4.243 distance: 221 - 233: 16.296 distance: 221 - 234: 12.517