Starting phenix.real_space_refine on Tue Nov 18 12:39:12 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/8wo4_37678/11_2025/8wo4_37678.cif Found real_map, /net/cci-nas-00/data/ceres_data/8wo4_37678/11_2025/8wo4_37678.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.8 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/8wo4_37678/11_2025/8wo4_37678.map" default_real_map = "/net/cci-nas-00/data/ceres_data/8wo4_37678/11_2025/8wo4_37678.map" model { file = "/net/cci-nas-00/data/ceres_data/8wo4_37678/11_2025/8wo4_37678.cif" } default_model = "/net/cci-nas-00/data/ceres_data/8wo4_37678/11_2025/8wo4_37678.cif" } resolution = 2.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.001 sd= 0.050 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 4 Type Number sf(0) Gaussians S 74 5.16 5 C 8888 2.51 5 N 2456 2.21 5 O 2770 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 96 residue(s): 0.01s Monomer Library directory: "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-2.0-5881/lib/python3.9/site-packages/chem_data/mon_lib" Total number of atoms: 14188 Number of models: 1 Model: "" Number of chains: 10 Chain: "A" Number of atoms: 1652 Number of conformers: 1 Conformer: "" Number of residues, atoms: 215, 1652 Classifications: {'peptide': 215} Link IDs: {'PCIS': 3, 'PTRANS': 8, 'TRANS': 203} Chain: "B" Number of atoms: 1689 Number of conformers: 1 Conformer: "" Number of residues, atoms: 228, 1689 Classifications: {'peptide': 228} Link IDs: {'PCIS': 2, 'PTRANS': 9, 'TRANS': 216} Chain: "C" Number of atoms: 1652 Number of conformers: 1 Conformer: "" Number of residues, atoms: 215, 1652 Classifications: {'peptide': 215} Link IDs: {'PCIS': 3, 'PTRANS': 8, 'TRANS': 203} Chain breaks: 1 Chain: "D" Number of atoms: 1689 Number of conformers: 1 Conformer: "" Number of residues, atoms: 228, 1689 Classifications: {'peptide': 228} Link IDs: {'PCIS': 2, 'PTRANS': 9, 'TRANS': 216} Chain: "E" Number of atoms: 2817 Number of conformers: 1 Conformer: "" Number of residues, atoms: 352, 2817 Classifications: {'peptide': 352} Link IDs: {'PCIS': 2, 'PTRANS': 17, 'TRANS': 332} Chain: "F" Number of atoms: 2817 Number of conformers: 1 Conformer: "" Number of residues, atoms: 352, 2817 Classifications: {'peptide': 352} Link IDs: {'PCIS': 2, 'PTRANS': 17, 'TRANS': 332} Chain: "G" Number of atoms: 922 Number of conformers: 1 Conformer: "" Number of residues, atoms: 120, 922 Classifications: {'peptide': 120} Link IDs: {'PTRANS': 3, 'TRANS': 116} Chain: "H" Number of atoms: 922 Number of conformers: 1 Conformer: "" Number of residues, atoms: 120, 922 Classifications: {'peptide': 120} Link IDs: {'PTRANS': 3, 'TRANS': 116} Chain: "E" Number of atoms: 14 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 14 Unusual residues: {'NAG': 1} Classifications: {'undetermined': 1} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen dihedrals: 3 Unresolved non-hydrogen chiralities: 1 Chain: "F" Number of atoms: 14 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 14 Unusual residues: {'NAG': 1} Classifications: {'undetermined': 1} Unresolved non-hydrogen bonds: 1 Unresolved non-hydrogen angles: 2 Unresolved non-hydrogen dihedrals: 3 Unresolved non-hydrogen chiralities: 1 Time building chain proxies: 3.58, per 1000 atoms: 0.25 Number of scatterers: 14188 At special positions: 0 Unit cell: (89.08, 150.28, 132.6, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 4 Type Number sf(0) S 74 16.00 O 2770 8.00 N 2456 7.00 C 8888 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=22, symmetry=0 Simple disulfide: pdb=" SG CYS A 24 " - pdb=" SG CYS A 89 " distance=2.03 Simple disulfide: pdb=" SG CYS A 135 " - pdb=" SG CYS A 195 " distance=2.03 Simple disulfide: pdb=" SG CYS A 215 " - pdb=" SG CYS B 223 " distance=2.03 Simple disulfide: pdb=" SG CYS B 23 " - pdb=" SG CYS B 96 " distance=2.03 Simple disulfide: pdb=" SG CYS B 147 " - pdb=" SG CYS B 203 " distance=2.03 Simple disulfide: pdb=" SG CYS C 24 " - pdb=" SG CYS C 89 " distance=2.04 Simple disulfide: pdb=" SG CYS C 135 " - pdb=" SG CYS C 195 " distance=2.03 Simple disulfide: pdb=" SG CYS C 215 " - pdb=" SG CYS D 223 " distance=2.03 Simple disulfide: pdb=" SG CYS D 23 " - pdb=" SG CYS D 96 " distance=2.03 Simple disulfide: pdb=" SG CYS D 147 " - pdb=" SG CYS D 203 " distance=2.03 Simple disulfide: pdb=" SG CYS E 4 " - pdb=" SG CYS F 4 " distance=2.05 Simple disulfide: pdb=" SG CYS E 55 " - pdb=" SG CYS E 143 " distance=2.03 Simple disulfide: pdb=" SG CYS E 179 " - pdb=" SG CYS E 223 " distance=2.03 Simple disulfide: pdb=" SG CYS E 280 " - pdb=" SG CYS E 329 " distance=2.03 Simple disulfide: pdb=" SG CYS E 291 " - pdb=" SG CYS E 312 " distance=2.03 Simple disulfide: pdb=" SG CYS E 313 " - pdb=" SG CYS E 316 " distance=2.04 Simple disulfide: pdb=" SG CYS F 4 " - pdb=" SG CYS F 15 " distance=2.04 Simple disulfide: pdb=" SG CYS F 55 " - pdb=" SG CYS F 143 " distance=2.03 Simple disulfide: pdb=" SG CYS F 179 " - pdb=" SG CYS F 223 " distance=2.03 Simple disulfide: pdb=" SG CYS F 280 " - pdb=" SG CYS F 329 " distance=2.03 Simple disulfide: pdb=" SG CYS F 291 " - pdb=" SG CYS F 312 " distance=2.03 Simple disulfide: pdb=" SG CYS F 313 " - pdb=" SG CYS F 316 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.00 Amino acid : True - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Links applied NAG-ASN " NAG E 401 " - " ASN E 207 " " NAG F 401 " - " ASN F 207 " Time building additional restraints: 1.06 Conformation dependent library (CDL) restraints added in 558.2 milliseconds 3620 Ramachandran restraints generated. 1810 Oldfield, 0 Emsley, 1810 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 3352 Finding SS restraints... Secondary structure from input PDB file: 39 helices and 35 sheets defined 12.8% alpha, 36.9% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.65 Creating SS restraints... Processing helix chain 'A' and resid 80 through 84 removed outlier: 4.190A pdb=" N PHE A 84 " --> pdb=" O PRO A 81 " (cutoff:3.500A) Processing helix chain 'A' and resid 122 through 128 removed outlier: 3.714A pdb=" N SER A 128 " --> pdb=" O GLU A 124 " (cutoff:3.500A) Processing helix chain 'A' and resid 184 through 189 Processing helix chain 'B' and resid 62 through 65 Processing helix chain 'B' and resid 87 through 91 Processing helix chain 'B' and resid 134 through 136 No H-bonds generated for 'chain 'B' and resid 134 through 136' Processing helix chain 'B' and resid 163 through 165 No H-bonds generated for 'chain 'B' and resid 163 through 165' Processing helix chain 'B' and resid 194 through 196 No H-bonds generated for 'chain 'B' and resid 194 through 196' Processing helix chain 'B' and resid 224 through 228 removed outlier: 3.882A pdb=" N HIS B 227 " --> pdb=" O ASP B 224 " (cutoff:3.500A) Processing helix chain 'C' and resid 80 through 84 removed outlier: 4.033A pdb=" N PHE C 84 " --> pdb=" O PRO C 81 " (cutoff:3.500A) Processing helix chain 'C' and resid 122 through 128 removed outlier: 3.623A pdb=" N SER C 128 " --> pdb=" O GLU C 124 " (cutoff:3.500A) Processing helix chain 'C' and resid 184 through 190 removed outlier: 3.587A pdb=" N LYS C 189 " --> pdb=" O ALA C 185 " (cutoff:3.500A) removed outlier: 3.916A pdb=" N HIS C 190 " --> pdb=" O ASP C 186 " (cutoff:3.500A) Processing helix chain 'D' and resid 62 through 65 Processing helix chain 'D' and resid 87 through 91 removed outlier: 3.569A pdb=" N ASP D 90 " --> pdb=" O ARG D 87 " (cutoff:3.500A) Processing helix chain 'D' and resid 134 through 136 No H-bonds generated for 'chain 'D' and resid 134 through 136' Processing helix chain 'D' and resid 163 through 165 No H-bonds generated for 'chain 'D' and resid 163 through 165' Processing helix chain 'D' and resid 194 through 196 No H-bonds generated for 'chain 'D' and resid 194 through 196' Processing helix chain 'D' and resid 224 through 228 removed outlier: 4.029A pdb=" N HIS D 227 " --> pdb=" O ASP D 224 " (cutoff:3.500A) Processing helix chain 'E' and resid 24 through 28 Processing helix chain 'E' and resid 38 through 53 removed outlier: 3.775A pdb=" N ALA E 44 " --> pdb=" O ARG E 40 " (cutoff:3.500A) Processing helix chain 'E' and resid 61 through 80 removed outlier: 5.995A pdb=" N GLY E 73 " --> pdb=" O LYS E 69 " (cutoff:3.500A) removed outlier: 5.872A pdb=" N GLU E 74 " --> pdb=" O SER E 70 " (cutoff:3.500A) Processing helix chain 'E' and resid 114 through 118 Processing helix chain 'E' and resid 139 through 143 Processing helix chain 'E' and resid 144 through 146 No H-bonds generated for 'chain 'E' and resid 144 through 146' Processing helix chain 'E' and resid 180 through 183 removed outlier: 4.056A pdb=" N VAL E 183 " --> pdb=" O ASP E 180 " (cutoff:3.500A) No H-bonds generated for 'chain 'E' and resid 180 through 183' Processing helix chain 'E' and resid 244 through 248 Processing helix chain 'E' and resid 341 through 345 Processing helix chain 'F' and resid 38 through 52 removed outlier: 3.635A pdb=" N GLU F 52 " --> pdb=" O GLN F 48 " (cutoff:3.500A) Processing helix chain 'F' and resid 61 through 82 removed outlier: 3.575A pdb=" N ASN F 65 " --> pdb=" O SER F 61 " (cutoff:3.500A) removed outlier: 5.052A pdb=" N GLY F 73 " --> pdb=" O LYS F 69 " (cutoff:3.500A) removed outlier: 5.861A pdb=" N GLU F 74 " --> pdb=" O SER F 70 " (cutoff:3.500A) Processing helix chain 'F' and resid 144 through 146 No H-bonds generated for 'chain 'F' and resid 144 through 146' Processing helix chain 'F' and resid 180 through 183 removed outlier: 3.972A pdb=" N VAL F 183 " --> pdb=" O ASP F 180 " (cutoff:3.500A) No H-bonds generated for 'chain 'F' and resid 180 through 183' Processing helix chain 'F' and resid 237 through 241 removed outlier: 3.517A pdb=" N LEU F 241 " --> pdb=" O GLU F 238 " (cutoff:3.500A) Processing helix chain 'F' and resid 244 through 248 removed outlier: 3.504A pdb=" N ALA F 248 " --> pdb=" O LYS F 245 " (cutoff:3.500A) Processing helix chain 'F' and resid 252 through 256 Processing helix chain 'F' and resid 341 through 345 removed outlier: 3.549A pdb=" N ASN F 344 " --> pdb=" O PRO F 341 " (cutoff:3.500A) Processing helix chain 'G' and resid 28 through 32 removed outlier: 3.575A pdb=" N ARG G 31 " --> pdb=" O THR G 28 " (cutoff:3.500A) Processing helix chain 'G' and resid 83 through 87 Processing helix chain 'H' and resid 28 through 32 Processing helix chain 'H' and resid 83 through 87 Processing sheet with id=AA1, first strand: chain 'A' and resid 5 through 8 removed outlier: 3.694A pdb=" N VAL A 20 " --> pdb=" O ILE A 76 " (cutoff:3.500A) removed outlier: 4.149A pdb=" N ASP A 71 " --> pdb=" O SER A 68 " (cutoff:3.500A) removed outlier: 3.570A pdb=" N SER A 68 " --> pdb=" O ASP A 71 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 11 through 14 removed outlier: 6.181A pdb=" N LEU A 12 " --> pdb=" O GLU A 106 " (cutoff:3.500A) removed outlier: 6.627A pdb=" N LEU A 34 " --> pdb=" O TYR A 50 " (cutoff:3.500A) removed outlier: 4.462A pdb=" N TYR A 50 " --> pdb=" O LEU A 34 " (cutoff:3.500A) removed outlier: 6.598A pdb=" N TRP A 36 " --> pdb=" O LEU A 48 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 115 through 119 removed outlier: 5.715A pdb=" N TYR A 174 " --> pdb=" O ASN A 139 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'A' and resid 154 through 155 Processing sheet with id=AA5, first strand: chain 'B' and resid 4 through 8 removed outlier: 3.572A pdb=" N GLN B 4 " --> pdb=" O SER B 26 " (cutoff:3.500A) removed outlier: 3.587A pdb=" N SER B 26 " --> pdb=" O GLN B 4 " (cutoff:3.500A) removed outlier: 3.585A pdb=" N THR B 78 " --> pdb=" O ASP B 73 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'B' and resid 12 through 13 removed outlier: 3.522A pdb=" N ALA B 92 " --> pdb=" O VAL B 116 " (cutoff:3.500A) removed outlier: 6.755A pdb=" N MET B 35 " --> pdb=" O VAL B 51 " (cutoff:3.500A) removed outlier: 4.388A pdb=" N VAL B 51 " --> pdb=" O MET B 35 " (cutoff:3.500A) removed outlier: 6.737A pdb=" N TRP B 37 " --> pdb=" O VAL B 49 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 12 through 13 removed outlier: 3.522A pdb=" N ALA B 92 " --> pdb=" O VAL B 116 " (cutoff:3.500A) removed outlier: 3.935A pdb=" N ILE B 109 " --> pdb=" O ARG B 98 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'B' and resid 127 through 131 removed outlier: 6.312A pdb=" N TYR B 183 " --> pdb=" O ASP B 151 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'B' and resid 138 through 139 removed outlier: 6.312A pdb=" N TYR B 183 " --> pdb=" O ASP B 151 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'B' and resid 158 through 161 removed outlier: 4.265A pdb=" N TYR B 201 " --> pdb=" O VAL B 218 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'C' and resid 5 through 8 Processing sheet with id=AB3, first strand: chain 'C' and resid 11 through 13 removed outlier: 6.545A pdb=" N LEU C 34 " --> pdb=" O TYR C 50 " (cutoff:3.500A) removed outlier: 4.439A pdb=" N TYR C 50 " --> pdb=" O LEU C 34 " (cutoff:3.500A) removed outlier: 6.646A pdb=" N TRP C 36 " --> pdb=" O LEU C 48 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'C' and resid 115 through 119 removed outlier: 5.935A pdb=" N TYR C 174 " --> pdb=" O ASN C 139 " (cutoff:3.500A) Processing sheet with id=AB5, first strand: chain 'C' and resid 154 through 155 Processing sheet with id=AB6, first strand: chain 'D' and resid 5 through 8 Processing sheet with id=AB7, first strand: chain 'D' and resid 11 through 13 removed outlier: 6.222A pdb=" N GLY D 11 " --> pdb=" O THR D 117 " (cutoff:3.500A) removed outlier: 3.622A pdb=" N SER D 36 " --> pdb=" O ALA D 97 " (cutoff:3.500A) removed outlier: 6.788A pdb=" N MET D 35 " --> pdb=" O VAL D 51 " (cutoff:3.500A) removed outlier: 4.391A pdb=" N VAL D 51 " --> pdb=" O MET D 35 " (cutoff:3.500A) removed outlier: 6.551A pdb=" N TRP D 37 " --> pdb=" O VAL D 49 " (cutoff:3.500A) Processing sheet with id=AB8, first strand: chain 'D' and resid 11 through 13 removed outlier: 6.222A pdb=" N GLY D 11 " --> pdb=" O THR D 117 " (cutoff:3.500A) Processing sheet with id=AB9, first strand: chain 'D' and resid 127 through 131 removed outlier: 6.232A pdb=" N TYR D 183 " --> pdb=" O ASP D 151 " (cutoff:3.500A) Processing sheet with id=AC1, first strand: chain 'D' and resid 138 through 139 removed outlier: 6.232A pdb=" N TYR D 183 " --> pdb=" O ASP D 151 " (cutoff:3.500A) Processing sheet with id=AC2, first strand: chain 'D' and resid 158 through 161 removed outlier: 4.531A pdb=" N TYR D 201 " --> pdb=" O VAL D 218 " (cutoff:3.500A) Processing sheet with id=AC3, first strand: chain 'E' and resid 3 through 6 removed outlier: 3.973A pdb=" N VAL F 2 " --> pdb=" O VAL E 6 " (cutoff:3.500A) removed outlier: 3.617A pdb=" N ARG F 14 " --> pdb=" O SER F 5 " (cutoff:3.500A) Processing sheet with id=AC4, first strand: chain 'E' and resid 87 through 90 removed outlier: 6.467A pdb=" N THR E 87 " --> pdb=" O PHE E 133 " (cutoff:3.500A) removed outlier: 6.976A pdb=" N VAL E 135 " --> pdb=" O THR E 87 " (cutoff:3.500A) removed outlier: 6.105A pdb=" N VAL E 89 " --> pdb=" O VAL E 135 " (cutoff:3.500A) removed outlier: 7.419A pdb=" N GLY E 56 " --> pdb=" O ALA E 149 " (cutoff:3.500A) Processing sheet with id=AC5, first strand: chain 'E' and resid 153 through 157 removed outlier: 6.982A pdb=" N TRP E 168 " --> pdb=" O GLU E 156 " (cutoff:3.500A) Processing sheet with id=AC6, first strand: chain 'E' and resid 329 through 331 removed outlier: 6.559A pdb=" N GLU E 203 " --> pdb=" O LYS E 213 " (cutoff:3.500A) removed outlier: 4.263A pdb=" N ALA E 215 " --> pdb=" O TRP E 201 " (cutoff:3.500A) removed outlier: 6.773A pdb=" N TRP E 201 " --> pdb=" O ALA E 215 " (cutoff:3.500A) removed outlier: 3.522A pdb=" N GLU F 192 " --> pdb=" O LYS F 189 " (cutoff:3.500A) removed outlier: 5.556A pdb=" N GLU F 205 " --> pdb=" O LEU F 212 " (cutoff:3.500A) removed outlier: 5.247A pdb=" N LEU F 212 " --> pdb=" O GLU F 205 " (cutoff:3.500A) removed outlier: 9.411A pdb=" N LYS F 213 " --> pdb=" O PHE F 277 " (cutoff:3.500A) removed outlier: 4.857A pdb=" N PHE F 277 " --> pdb=" O LYS F 213 " (cutoff:3.500A) Processing sheet with id=AC7, first strand: chain 'E' and resid 284 through 285 removed outlier: 6.275A pdb=" N LYS E 284 " --> pdb=" O TRP E 311 " (cutoff:3.500A) Processing sheet with id=AC8, first strand: chain 'F' and resid 56 through 57 removed outlier: 7.176A pdb=" N GLY F 56 " --> pdb=" O ALA F 149 " (cutoff:3.500A) No H-bonds generated for sheet with id=AC8 Processing sheet with id=AC9, first strand: chain 'F' and resid 87 through 91 removed outlier: 6.500A pdb=" N THR F 87 " --> pdb=" O PHE F 133 " (cutoff:3.500A) removed outlier: 7.041A pdb=" N VAL F 135 " --> pdb=" O THR F 87 " (cutoff:3.500A) removed outlier: 6.203A pdb=" N VAL F 89 " --> pdb=" O VAL F 135 " (cutoff:3.500A) No H-bonds generated for sheet with id=AC9 Processing sheet with id=AD1, first strand: chain 'F' and resid 284 through 287 Processing sheet with id=AD2, first strand: chain 'G' and resid 3 through 7 removed outlier: 3.522A pdb=" N CYS G 22 " --> pdb=" O VAL G 78 " (cutoff:3.500A) removed outlier: 3.922A pdb=" N VAL G 78 " --> pdb=" O CYS G 22 " (cutoff:3.500A) removed outlier: 3.641A pdb=" N LEU G 18 " --> pdb=" O MET G 82 " (cutoff:3.500A) Processing sheet with id=AD3, first strand: chain 'G' and resid 3 through 7 removed outlier: 3.522A pdb=" N CYS G 22 " --> pdb=" O VAL G 78 " (cutoff:3.500A) removed outlier: 3.922A pdb=" N VAL G 78 " --> pdb=" O CYS G 22 " (cutoff:3.500A) removed outlier: 3.515A pdb=" N THR G 77 " --> pdb=" O ASP G 72 " (cutoff:3.500A) Processing sheet with id=AD4, first strand: chain 'G' and resid 57 through 59 removed outlier: 5.547A pdb=" N PHE G 47 " --> pdb=" O ARG G 38 " (cutoff:3.500A) removed outlier: 5.926A pdb=" N ARG G 38 " --> pdb=" O PHE G 47 " (cutoff:3.500A) removed outlier: 4.001A pdb=" N ALA G 33 " --> pdb=" O ASP G 95 " (cutoff:3.500A) removed outlier: 6.795A pdb=" N CYS G 92 " --> pdb=" O TRP G 103 " (cutoff:3.500A) removed outlier: 4.586A pdb=" N TRP G 103 " --> pdb=" O CYS G 92 " (cutoff:3.500A) removed outlier: 7.500A pdb=" N ILE G 94 " --> pdb=" O ASP G 101 " (cutoff:3.500A) removed outlier: 3.725A pdb=" N ASP G 101 " --> pdb=" O ILE G 94 " (cutoff:3.500A) Processing sheet with id=AD5, first strand: chain 'G' and resid 57 through 59 removed outlier: 5.547A pdb=" N PHE G 47 " --> pdb=" O ARG G 38 " (cutoff:3.500A) removed outlier: 5.926A pdb=" N ARG G 38 " --> pdb=" O PHE G 47 " (cutoff:3.500A) removed outlier: 4.001A pdb=" N ALA G 33 " --> pdb=" O ASP G 95 " (cutoff:3.500A) removed outlier: 3.583A pdb=" N THR G 107 " --> pdb=" O TYR G 90 " (cutoff:3.500A) Processing sheet with id=AD6, first strand: chain 'H' and resid 3 through 7 removed outlier: 3.669A pdb=" N GLN H 5 " --> pdb=" O ALA H 23 " (cutoff:3.500A) Processing sheet with id=AD7, first strand: chain 'H' and resid 57 through 59 removed outlier: 3.534A pdb=" N PHE H 58 " --> pdb=" O VAL H 50 " (cutoff:3.500A) removed outlier: 5.464A pdb=" N PHE H 47 " --> pdb=" O ARG H 38 " (cutoff:3.500A) removed outlier: 5.052A pdb=" N ARG H 38 " --> pdb=" O PHE H 47 " (cutoff:3.500A) removed outlier: 3.547A pdb=" N ALA H 49 " --> pdb=" O TRP H 36 " (cutoff:3.500A) removed outlier: 4.182A pdb=" N ALA H 33 " --> pdb=" O ASP H 95 " (cutoff:3.500A) removed outlier: 4.440A pdb=" N TYR H 102 " --> pdb=" O ILE H 94 " (cutoff:3.500A) Processing sheet with id=AD8, first strand: chain 'H' and resid 57 through 59 removed outlier: 3.534A pdb=" N PHE H 58 " --> pdb=" O VAL H 50 " (cutoff:3.500A) removed outlier: 5.464A pdb=" N PHE H 47 " --> pdb=" O ARG H 38 " (cutoff:3.500A) removed outlier: 5.052A pdb=" N ARG H 38 " --> pdb=" O PHE H 47 " (cutoff:3.500A) removed outlier: 3.547A pdb=" N ALA H 49 " --> pdb=" O TRP H 36 " (cutoff:3.500A) removed outlier: 4.182A pdb=" N ALA H 33 " --> pdb=" O ASP H 95 " (cutoff:3.500A) removed outlier: 3.573A pdb=" N THR H 107 " --> pdb=" O TYR H 90 " (cutoff:3.500A) 525 hydrogen bonds defined for protein. 1323 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: 2.57 Time building geometry restraints manager: 1.47 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.57: 14430 1.57 - 1.93: 100 1.93 - 2.28: 0 2.28 - 2.63: 0 2.63 - 2.98: 1 Bond restraints: 14531 Sorted by residual: bond pdb=" C ARG A 109 " pdb=" N THR A 110 " ideal model delta sigma weight residual 1.329 2.983 -1.655 1.44e-02 4.82e+03 1.32e+04 bond pdb=" N THR A 110 " pdb=" CA THR A 110 " ideal model delta sigma weight residual 1.453 1.488 -0.035 1.22e-02 6.72e+03 8.17e+00 bond pdb=" N THR C 110 " pdb=" CA THR C 110 " ideal model delta sigma weight residual 1.458 1.492 -0.034 1.90e-02 2.77e+03 3.21e+00 bond pdb=" N ASP E 1 " pdb=" CA ASP E 1 " ideal model delta sigma weight residual 1.458 1.492 -0.034 1.90e-02 2.77e+03 3.19e+00 bond pdb=" N SER B 1 " pdb=" CA SER B 1 " ideal model delta sigma weight residual 1.458 1.491 -0.033 1.90e-02 2.77e+03 2.96e+00 ... (remaining 14526 not shown) Histogram of bond angle deviations from ideal: 0.00 - 6.40: 19708 6.40 - 12.81: 3 12.81 - 19.21: 0 19.21 - 25.62: 1 25.62 - 32.02: 1 Bond angle restraints: 19713 Sorted by residual: angle pdb=" O ARG A 109 " pdb=" C ARG A 109 " pdb=" N THR A 110 " ideal model delta sigma weight residual 123.33 101.03 22.30 1.13e+00 7.83e-01 3.89e+02 angle pdb=" C ARG A 109 " pdb=" N THR A 110 " pdb=" CA THR A 110 " ideal model delta sigma weight residual 122.03 90.01 32.02 2.01e+00 2.48e-01 2.54e+02 angle pdb=" C ASP F 136 " pdb=" CA ASP F 136 " pdb=" CB ASP F 136 " ideal model delta sigma weight residual 117.23 111.52 5.71 1.36e+00 5.41e-01 1.76e+01 angle pdb=" C1 NAG F 401 " pdb=" O5 NAG F 401 " pdb=" C5 NAG F 401 " ideal model delta sigma weight residual 113.21 124.37 -11.16 3.00e+00 1.11e-01 1.38e+01 angle pdb=" CA CYS E 4 " pdb=" CB CYS E 4 " pdb=" SG CYS E 4 " ideal model delta sigma weight residual 114.40 120.84 -6.44 2.30e+00 1.89e-01 7.84e+00 ... (remaining 19708 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.81: 8073 17.81 - 35.62: 465 35.62 - 53.44: 120 53.44 - 71.25: 35 71.25 - 89.06: 8 Dihedral angle restraints: 8701 sinusoidal: 3454 harmonic: 5247 Sorted by residual: dihedral pdb=" CB CYS E 4 " pdb=" SG CYS E 4 " pdb=" SG CYS F 4 " pdb=" CB CYS F 4 " ideal model delta sinusoidal sigma weight residual -86.00 -15.55 -70.45 1 1.00e+01 1.00e-02 6.39e+01 dihedral pdb=" CA LYS E 33 " pdb=" C LYS E 33 " pdb=" N TYR E 34 " pdb=" CA TYR E 34 " ideal model delta harmonic sigma weight residual -180.00 -154.02 -25.98 0 5.00e+00 4.00e-02 2.70e+01 dihedral pdb=" CA TRP E 28 " pdb=" C TRP E 28 " pdb=" N ARG E 29 " pdb=" CA ARG E 29 " ideal model delta harmonic sigma weight residual -180.00 -154.52 -25.48 0 5.00e+00 4.00e-02 2.60e+01 ... (remaining 8698 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.040: 1596 0.040 - 0.080: 356 0.080 - 0.121: 160 0.121 - 0.161: 23 0.161 - 0.201: 1 Chirality restraints: 2136 Sorted by residual: chirality pdb=" C1 NAG F 401 " pdb=" ND2 ASN F 207 " pdb=" C2 NAG F 401 " pdb=" O5 NAG F 401 " both_signs ideal model delta sigma weight residual False -2.40 -2.20 -0.20 2.00e-01 2.50e+01 1.01e+00 chirality pdb=" CA VAL H 2 " pdb=" N VAL H 2 " pdb=" C VAL H 2 " pdb=" CB VAL H 2 " both_signs ideal model delta sigma weight residual False 2.44 2.58 -0.14 2.00e-01 2.50e+01 4.81e-01 chirality pdb=" CA ILE C 22 " pdb=" N ILE C 22 " pdb=" C ILE C 22 " pdb=" CB ILE C 22 " both_signs ideal model delta sigma weight residual False 2.43 2.57 -0.14 2.00e-01 2.50e+01 4.75e-01 ... (remaining 2133 not shown) Planarity restraints: 2537 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA ARG A 109 " -0.133 2.00e-02 2.50e+03 1.99e-01 3.95e+02 pdb=" C ARG A 109 " 0.337 2.00e-02 2.50e+03 pdb=" O ARG A 109 " -0.157 2.00e-02 2.50e+03 pdb=" N THR A 110 " -0.046 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C GLU D 155 " 0.032 5.00e-02 4.00e+02 4.80e-02 3.69e+00 pdb=" N PRO D 156 " -0.083 5.00e-02 4.00e+02 pdb=" CA PRO D 156 " 0.027 5.00e-02 4.00e+02 pdb=" CD PRO D 156 " 0.024 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C GLU B 155 " 0.031 5.00e-02 4.00e+02 4.66e-02 3.47e+00 pdb=" N PRO B 156 " -0.081 5.00e-02 4.00e+02 pdb=" CA PRO B 156 " 0.026 5.00e-02 4.00e+02 pdb=" CD PRO B 156 " 0.023 5.00e-02 4.00e+02 ... (remaining 2534 not shown) Histogram of nonbonded interaction distances: 1.96 - 2.54: 96 2.54 - 3.13: 10829 3.13 - 3.72: 20068 3.72 - 4.31: 28770 4.31 - 4.90: 49852 Nonbonded interactions: 109615 Sorted by model distance: nonbonded pdb=" O GLY F 159 " pdb=" OG SER F 166 " model vdw 1.955 3.040 nonbonded pdb=" OD1 ASP E 234 " pdb=" OG1 THR F 233 " model vdw 1.974 3.040 nonbonded pdb=" NE2 GLN A 101 " pdb=" O GLY A 102 " model vdw 2.003 3.120 nonbonded pdb=" OD1 ASP B 224 " pdb=" OG1 THR B 226 " model vdw 2.228 3.040 nonbonded pdb=" OD1 ASP C 18 " pdb=" N ARG C 19 " model vdw 2.228 3.120 ... (remaining 109610 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.01 Found NCS groups: ncs_group { reference = chain 'A' selection = chain 'C' } ncs_group { reference = chain 'B' selection = chain 'D' } ncs_group { reference = chain 'E' selection = chain 'F' } ncs_group { reference = chain 'G' selection = chain 'H' } 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: 19.700 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.010 Extract box with map and model: 0.350 Check model and map are aligned: 0.050 Set scattering table: 0.040 Process input model: 14.390 Find NCS groups from input model: 0.210 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.000 Load rotamer database and sin/cos tables:9.360 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 44.140 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6216 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.014 1.655 14555 Z= 0.961 Angle : 0.604 32.023 19763 Z= 0.342 Chirality : 0.042 0.201 2136 Planarity : 0.005 0.199 2535 Dihedral : 12.778 89.062 5283 Min Nonbonded Distance : 1.955 Molprobity Statistics. All-atom Clashscore : 5.67 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.81 % Favored : 96.19 % Rotamer: Outliers : 0.96 % Allowed : 8.39 % Favored : 90.65 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 15.91 % 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: 0.01 (0.20), residues: 1810 helix: 0.49 (0.51), residues: 105 sheet: 0.58 (0.22), residues: 589 loop : -0.32 (0.19), residues: 1116 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.004 0.000 ARG A 212 TYR 0.022 0.001 TYR H 91 PHE 0.009 0.001 PHE G 99 TRP 0.015 0.001 TRP D 99 HIS 0.003 0.001 HIS E 195 Details of bonding type rmsd covalent geometry : bond 0.01397 (14531) covalent geometry : angle 0.59732 (19713) SS BOND : bond 0.00623 ( 22) SS BOND : angle 1.47973 ( 44) hydrogen bonds : bond 0.24879 ( 469) hydrogen bonds : angle 9.62659 ( 1323) link_NAG-ASN : bond 0.00311 ( 2) link_NAG-ASN : angle 3.55065 ( 6) *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3620 Ramachandran restraints generated. 1810 Oldfield, 0 Emsley, 1810 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3620 Ramachandran restraints generated. 1810 Oldfield, 0 Emsley, 1810 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 148 residues out of total 1562 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 15 poor density : 133 time to evaluate : 0.600 Fit side-chains REVERT: D 65 LYS cc_start: 0.6649 (ttpt) cc_final: 0.6342 (tptt) REVERT: F 138 ASP cc_start: 0.7366 (t0) cc_final: 0.6812 (t0) outliers start: 15 outliers final: 2 residues processed: 143 average time/residue: 0.6455 time to fit residues: 100.9762 Evaluate side-chains 103 residues out of total 1562 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 2 poor density : 101 time to evaluate : 0.550 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 105 VAL Chi-restraints excluded: chain F residue 4 CYS Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 180 random chunks: chunk 98 optimal weight: 0.0370 chunk 107 optimal weight: 4.9990 chunk 10 optimal weight: 0.9980 chunk 66 optimal weight: 0.0370 chunk 130 optimal weight: 4.9990 chunk 124 optimal weight: 6.9990 chunk 103 optimal weight: 3.9990 chunk 77 optimal weight: 0.0020 chunk 122 optimal weight: 0.0670 chunk 91 optimal weight: 5.9990 chunk 149 optimal weight: 1.9990 overall best weight: 0.2282 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 35 ASN A 56 GLN A 156 GLN A 159 ASN C 91 GLN C 211 ASN ** F 35 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 6 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4784 r_free = 0.4784 target = 0.243972 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 32)----------------| | r_work = 0.4451 r_free = 0.4451 target = 0.200010 restraints weight = 17387.159| |-----------------------------------------------------------------------------| r_work (start): 0.4365 rms_B_bonded: 3.38 r_work: 0.4047 rms_B_bonded: 3.45 restraints_weight: 0.5000 r_work (final): 0.4047 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6659 moved from start: 0.1277 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.036 14555 Z= 0.114 Angle : 0.592 9.207 19763 Z= 0.307 Chirality : 0.042 0.158 2136 Planarity : 0.004 0.051 2535 Dihedral : 4.850 57.499 2033 Min Nonbonded Distance : 2.410 Molprobity Statistics. All-atom Clashscore : 3.69 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.59 % Favored : 97.41 % Rotamer: Outliers : 1.92 % Allowed : 9.48 % Favored : 88.60 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 15.91 % 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: 0.33 (0.20), residues: 1812 helix: 0.68 (0.52), residues: 110 sheet: 0.85 (0.22), residues: 606 loop : -0.15 (0.20), residues: 1096 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.005 0.000 ARG A 109 TYR 0.022 0.001 TYR H 91 PHE 0.023 0.001 PHE D 68 TRP 0.016 0.001 TRP H 36 HIS 0.005 0.001 HIS E 195 Details of bonding type rmsd covalent geometry : bond 0.00245 (14531) covalent geometry : angle 0.58389 (19713) SS BOND : bond 0.00542 ( 22) SS BOND : angle 1.45893 ( 44) hydrogen bonds : bond 0.04032 ( 469) hydrogen bonds : angle 6.52251 ( 1323) link_NAG-ASN : bond 0.00323 ( 2) link_NAG-ASN : angle 4.51863 ( 6) *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3624 Ramachandran restraints generated. 1812 Oldfield, 0 Emsley, 1812 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3624 Ramachandran restraints generated. 1812 Oldfield, 0 Emsley, 1812 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 155 residues out of total 1562 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 30 poor density : 125 time to evaluate : 0.572 Fit side-chains REVERT: C 76 ILE cc_start: 0.6087 (OUTLIER) cc_final: 0.5372 (tp) REVERT: C 148 GLN cc_start: 0.6058 (OUTLIER) cc_final: 0.5855 (mm-40) REVERT: D 65 LYS cc_start: 0.6596 (ttpt) cc_final: 0.6079 (tptt) REVERT: F 314 ARG cc_start: 0.6600 (ttm-80) cc_final: 0.6042 (ttm-80) outliers start: 30 outliers final: 13 residues processed: 148 average time/residue: 0.6222 time to fit residues: 101.3061 Evaluate side-chains 125 residues out of total 1562 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 15 poor density : 110 time to evaluate : 0.537 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 21 THR Chi-restraints excluded: chain B residue 13 ILE Chi-restraints excluded: chain C residue 76 ILE Chi-restraints excluded: chain C residue 148 GLN Chi-restraints excluded: chain C residue 197 VAL Chi-restraints excluded: chain F residue 1 ASP Chi-restraints excluded: chain F residue 4 CYS Chi-restraints excluded: chain F residue 67 MET Chi-restraints excluded: chain F residue 71 VAL Chi-restraints excluded: chain F residue 313 CYS Chi-restraints excluded: chain F residue 316 CYS Chi-restraints excluded: chain F residue 329 CYS Chi-restraints excluded: chain G residue 29 ILE Chi-restraints excluded: chain H residue 36 TRP Chi-restraints excluded: chain H residue 89 VAL Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 180 random chunks: chunk 7 optimal weight: 4.9990 chunk 165 optimal weight: 20.0000 chunk 24 optimal weight: 3.9990 chunk 160 optimal weight: 1.9990 chunk 94 optimal weight: 3.9990 chunk 11 optimal weight: 0.9980 chunk 81 optimal weight: 7.9990 chunk 103 optimal weight: 2.9990 chunk 26 optimal weight: 0.0870 chunk 149 optimal weight: 0.8980 chunk 42 optimal weight: 10.0000 overall best weight: 1.3962 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4747 r_free = 0.4747 target = 0.238744 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 58)----------------| | r_work = 0.4372 r_free = 0.4372 target = 0.191794 restraints weight = 17355.034| |-----------------------------------------------------------------------------| r_work (start): 0.4271 rms_B_bonded: 4.03 r_work: 0.3971 rms_B_bonded: 3.17 restraints_weight: 0.5000 r_work (final): 0.3971 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6776 moved from start: 0.1717 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.070 14555 Z= 0.153 Angle : 0.636 10.353 19763 Z= 0.329 Chirality : 0.044 0.181 2136 Planarity : 0.005 0.049 2535 Dihedral : 5.087 58.302 2032 Min Nonbonded Distance : 2.379 Molprobity Statistics. All-atom Clashscore : 5.24 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.97 % Favored : 96.03 % Rotamer: Outliers : 3.33 % Allowed : 10.50 % Favored : 86.17 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 15.91 % 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: 0.04 (0.20), residues: 1812 helix: 0.04 (0.49), residues: 116 sheet: 0.71 (0.22), residues: 596 loop : -0.32 (0.19), residues: 1100 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.003 0.001 ARG B 67 TYR 0.019 0.002 TYR H 91 PHE 0.024 0.002 PHE D 68 TRP 0.031 0.002 TRP D 99 HIS 0.006 0.001 HIS E 195 Details of bonding type rmsd covalent geometry : bond 0.00345 (14531) covalent geometry : angle 0.62618 (19713) SS BOND : bond 0.00930 ( 22) SS BOND : angle 1.85626 ( 44) hydrogen bonds : bond 0.04313 ( 469) hydrogen bonds : angle 6.01489 ( 1323) link_NAG-ASN : bond 0.00351 ( 2) link_NAG-ASN : angle 4.37658 ( 6) *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3624 Ramachandran restraints generated. 1812 Oldfield, 0 Emsley, 1812 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3624 Ramachandran restraints generated. 1812 Oldfield, 0 Emsley, 1812 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 170 residues out of total 1562 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 52 poor density : 118 time to evaluate : 0.566 Fit side-chains REVERT: C 76 ILE cc_start: 0.6132 (OUTLIER) cc_final: 0.5301 (tp) REVERT: C 78 SER cc_start: 0.6706 (p) cc_final: 0.6322 (t) REVERT: D 58 THR cc_start: 0.7997 (p) cc_final: 0.7692 (p) REVERT: D 65 LYS cc_start: 0.6893 (ttpt) cc_final: 0.6282 (tptt) REVERT: E 11 LYS cc_start: 0.6117 (tppt) cc_final: 0.4818 (mppt) REVERT: E 85 GLN cc_start: 0.7972 (OUTLIER) cc_final: 0.7664 (mm-40) REVERT: F 85 GLN cc_start: 0.6738 (tm-30) cc_final: 0.6020 (tp40) REVERT: F 133 PHE cc_start: 0.7586 (t80) cc_final: 0.7310 (t80) REVERT: F 138 ASP cc_start: 0.7552 (t0) cc_final: 0.6992 (t0) outliers start: 52 outliers final: 25 residues processed: 154 average time/residue: 0.6136 time to fit residues: 103.9902 Evaluate side-chains 140 residues out of total 1562 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 27 poor density : 113 time to evaluate : 0.360 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 13 SER Chi-restraints excluded: chain A residue 21 THR Chi-restraints excluded: chain B residue 13 ILE Chi-restraints excluded: chain B residue 30 VAL Chi-restraints excluded: chain B residue 47 GLU Chi-restraints excluded: chain B residue 69 THR Chi-restraints excluded: chain C residue 21 THR Chi-restraints excluded: chain C residue 76 ILE Chi-restraints excluded: chain C residue 80 GLN Chi-restraints excluded: chain C residue 197 VAL Chi-restraints excluded: chain D residue 19 LEU Chi-restraints excluded: chain D residue 159 VAL Chi-restraints excluded: chain E residue 85 GLN Chi-restraints excluded: chain E residue 138 ASP Chi-restraints excluded: chain E residue 171 VAL Chi-restraints excluded: chain E residue 293 THR Chi-restraints excluded: chain F residue 1 ASP Chi-restraints excluded: chain F residue 57 ILE Chi-restraints excluded: chain F residue 71 VAL Chi-restraints excluded: chain F residue 87 THR Chi-restraints excluded: chain F residue 206 LYS Chi-restraints excluded: chain F residue 313 CYS Chi-restraints excluded: chain F residue 329 CYS Chi-restraints excluded: chain G residue 29 ILE Chi-restraints excluded: chain G residue 48 VAL Chi-restraints excluded: chain H residue 36 TRP Chi-restraints excluded: chain H residue 89 VAL Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 180 random chunks: chunk 89 optimal weight: 0.8980 chunk 67 optimal weight: 3.9990 chunk 143 optimal weight: 2.9990 chunk 75 optimal weight: 0.6980 chunk 20 optimal weight: 7.9990 chunk 176 optimal weight: 10.0000 chunk 163 optimal weight: 20.0000 chunk 3 optimal weight: 3.9990 chunk 151 optimal weight: 8.9990 chunk 133 optimal weight: 6.9990 chunk 81 optimal weight: 10.0000 overall best weight: 2.5186 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 138 ASN B 204 ASN C 167 GLN D 204 ASN ** F 164 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4699 r_free = 0.4699 target = 0.234089 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 61)----------------| | r_work = 0.4107 r_free = 0.4107 target = 0.171136 restraints weight = 17551.097| |-----------------------------------------------------------------------------| r_work (start): 0.4082 rms_B_bonded: 2.21 r_work: 0.3983 rms_B_bonded: 2.61 restraints_weight: 0.5000 r_work: 0.3868 rms_B_bonded: 4.24 restraints_weight: 0.2500 r_work (final): 0.3868 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6970 moved from start: 0.2335 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.099 14555 Z= 0.223 Angle : 0.733 10.227 19763 Z= 0.382 Chirality : 0.047 0.185 2136 Planarity : 0.006 0.056 2535 Dihedral : 5.570 58.318 2031 Min Nonbonded Distance : 2.341 Molprobity Statistics. All-atom Clashscore : 6.96 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.36 % Favored : 95.64 % Rotamer: Outliers : 3.71 % Allowed : 12.87 % Favored : 83.42 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 15.91 % 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: -0.37 (0.20), residues: 1812 helix: -0.35 (0.49), residues: 105 sheet: 0.53 (0.21), residues: 623 loop : -0.69 (0.19), residues: 1084 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.005 0.001 ARG A 212 TYR 0.026 0.002 TYR F 32 PHE 0.045 0.003 PHE D 68 TRP 0.051 0.003 TRP D 99 HIS 0.008 0.002 HIS E 195 Details of bonding type rmsd covalent geometry : bond 0.00519 (14531) covalent geometry : angle 0.72060 (19713) SS BOND : bond 0.01223 ( 22) SS BOND : angle 2.22805 ( 44) hydrogen bonds : bond 0.04875 ( 469) hydrogen bonds : angle 6.11274 ( 1323) link_NAG-ASN : bond 0.00525 ( 2) link_NAG-ASN : angle 5.00156 ( 6) *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3624 Ramachandran restraints generated. 1812 Oldfield, 0 Emsley, 1812 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3624 Ramachandran restraints generated. 1812 Oldfield, 0 Emsley, 1812 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 182 residues out of total 1562 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 58 poor density : 124 time to evaluate : 0.554 Fit side-chains REVERT: A 111 VAL cc_start: 0.5651 (m) cc_final: 0.5334 (t) REVERT: A 124 GLU cc_start: 0.4903 (tm-30) cc_final: 0.4472 (tm-30) REVERT: A 127 LYS cc_start: 0.6327 (OUTLIER) cc_final: 0.5804 (tttm) REVERT: C 76 ILE cc_start: 0.6384 (OUTLIER) cc_final: 0.5571 (tp) REVERT: C 78 SER cc_start: 0.6753 (p) cc_final: 0.6465 (t) REVERT: C 101 GLN cc_start: 0.8057 (OUTLIER) cc_final: 0.7575 (mt0) REVERT: E 11 LYS cc_start: 0.6180 (tppt) cc_final: 0.5184 (mppt) REVERT: E 85 GLN cc_start: 0.8061 (OUTLIER) cc_final: 0.7815 (mm-40) REVERT: F 47 LYS cc_start: 0.6310 (OUTLIER) cc_final: 0.5834 (ttmm) outliers start: 58 outliers final: 30 residues processed: 167 average time/residue: 0.6060 time to fit residues: 110.9226 Evaluate side-chains 143 residues out of total 1562 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 35 poor density : 108 time to evaluate : 0.597 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 12 LEU Chi-restraints excluded: chain A residue 13 SER Chi-restraints excluded: chain A residue 21 THR Chi-restraints excluded: chain A residue 127 LYS Chi-restraints excluded: chain A residue 173 THR Chi-restraints excluded: chain B residue 13 ILE Chi-restraints excluded: chain B residue 30 VAL Chi-restraints excluded: chain B residue 47 GLU Chi-restraints excluded: chain B residue 69 THR Chi-restraints excluded: chain B residue 95 TYR Chi-restraints excluded: chain C residue 21 THR Chi-restraints excluded: chain C residue 76 ILE Chi-restraints excluded: chain C residue 80 GLN Chi-restraints excluded: chain C residue 101 GLN Chi-restraints excluded: chain C residue 197 VAL Chi-restraints excluded: chain D residue 19 LEU Chi-restraints excluded: chain D residue 159 VAL Chi-restraints excluded: chain E residue 85 GLN Chi-restraints excluded: chain E residue 87 THR Chi-restraints excluded: chain E residue 138 ASP Chi-restraints excluded: chain E residue 171 VAL Chi-restraints excluded: chain E residue 205 GLU Chi-restraints excluded: chain E residue 217 LEU Chi-restraints excluded: chain E residue 264 VAL Chi-restraints excluded: chain E residue 293 THR Chi-restraints excluded: chain F residue 47 LYS Chi-restraints excluded: chain F residue 57 ILE Chi-restraints excluded: chain F residue 71 VAL Chi-restraints excluded: chain F residue 87 THR Chi-restraints excluded: chain F residue 313 CYS Chi-restraints excluded: chain F residue 329 CYS Chi-restraints excluded: chain G residue 29 ILE Chi-restraints excluded: chain G residue 48 VAL Chi-restraints excluded: chain H residue 36 TRP Chi-restraints excluded: chain H residue 89 VAL Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 180 random chunks: chunk 104 optimal weight: 2.9990 chunk 134 optimal weight: 5.9990 chunk 93 optimal weight: 0.6980 chunk 76 optimal weight: 5.9990 chunk 12 optimal weight: 9.9990 chunk 87 optimal weight: 0.8980 chunk 100 optimal weight: 2.9990 chunk 5 optimal weight: 0.9980 chunk 97 optimal weight: 1.9990 chunk 172 optimal weight: 3.9990 chunk 44 optimal weight: 4.9990 overall best weight: 1.5184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 138 ASN ** A 200 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** E 254 HIS ** F 164 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.4715 r_free = 0.4715 target = 0.235566 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 41)----------------| | r_work = 0.4200 r_free = 0.4200 target = 0.178833 restraints weight = 17564.216| |-----------------------------------------------------------------------------| r_work (start): 0.4141 rms_B_bonded: 3.03 r_work: 0.3966 rms_B_bonded: 2.76 restraints_weight: 0.5000 r_work: 0.3850 rms_B_bonded: 4.47 restraints_weight: 0.2500 r_work (final): 0.3850 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.6940 moved from start: 0.2570 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.069 14555 Z= 0.156 Angle : 0.632 8.897 19763 Z= 0.329 Chirality : 0.043 0.166 2136 Planarity : 0.005 0.056 2535 Dihedral : 5.251 58.037 2031 Min Nonbonded Distance : 2.394 Molprobity Statistics. All-atom Clashscore : 6.35 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.42 % Favored : 95.58 % Rotamer: Outliers : 3.71 % Allowed : 14.21 % Favored : 82.07 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 15.91 % 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: -0.31 (0.20), residues: 1812 helix: 0.15 (0.52), residues: 103 sheet: 0.53 (0.21), residues: 630 loop : -0.69 (0.19), residues: 1079 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.006 0.000 ARG F 314 TYR 0.020 0.002 TYR F 32 PHE 0.035 0.002 PHE D 68 TRP 0.039 0.002 TRP D 99 HIS 0.007 0.001 HIS E 195 Details of bonding type rmsd covalent geometry : bond 0.00358 (14531) covalent geometry : angle 0.62254 (19713) SS BOND : bond 0.00861 ( 22) SS BOND : angle 1.81354 ( 44) hydrogen bonds : bond 0.04053 ( 469) hydrogen bonds : angle 5.79648 ( 1323) link_NAG-ASN : bond 0.00443 ( 2) link_NAG-ASN : angle 4.33153 ( 6) *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3624 Ramachandran restraints generated. 1812 Oldfield, 0 Emsley, 1812 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3624 Ramachandran restraints generated. 1812 Oldfield, 0 Emsley, 1812 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 177 residues out of total 1562 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 58 poor density : 119 time to evaluate : 0.533 Fit side-chains REVERT: A 19 ARG cc_start: 0.6563 (mmm-85) cc_final: 0.6316 (tpt-90) REVERT: A 111 VAL cc_start: 0.5675 (m) cc_final: 0.5352 (t) REVERT: C 76 ILE cc_start: 0.6378 (OUTLIER) cc_final: 0.5575 (tp) REVERT: C 78 SER cc_start: 0.6732 (p) cc_final: 0.6386 (t) REVERT: C 101 GLN cc_start: 0.8003 (OUTLIER) cc_final: 0.7523 (mt0) REVERT: E 11 LYS cc_start: 0.6222 (tppt) cc_final: 0.5229 (mppt) REVERT: F 21 ILE cc_start: 0.7889 (OUTLIER) cc_final: 0.6964 (mp) REVERT: F 314 ARG cc_start: 0.6735 (ttm-80) cc_final: 0.6405 (ttm-80) REVERT: H 82 MET cc_start: 0.0702 (OUTLIER) cc_final: 0.0410 (mtt) outliers start: 58 outliers final: 34 residues processed: 160 average time/residue: 0.5852 time to fit residues: 103.1544 Evaluate side-chains 150 residues out of total 1562 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 38 poor density : 112 time to evaluate : 0.538 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 12 LEU Chi-restraints excluded: chain A residue 13 SER Chi-restraints excluded: chain A residue 21 THR Chi-restraints excluded: chain A residue 79 LEU Chi-restraints excluded: chain A residue 173 THR Chi-restraints excluded: chain B residue 13 ILE Chi-restraints excluded: chain B residue 30 VAL Chi-restraints excluded: chain B residue 47 GLU Chi-restraints excluded: chain B residue 69 THR Chi-restraints excluded: chain C residue 20 VAL Chi-restraints excluded: chain C residue 21 THR Chi-restraints excluded: chain C residue 76 ILE Chi-restraints excluded: chain C residue 80 GLN Chi-restraints excluded: chain C residue 101 GLN Chi-restraints excluded: chain C residue 197 VAL Chi-restraints excluded: chain D residue 19 LEU Chi-restraints excluded: chain D residue 47 GLU Chi-restraints excluded: chain D residue 159 VAL Chi-restraints excluded: chain E residue 6 VAL Chi-restraints excluded: chain E residue 87 THR Chi-restraints excluded: chain E residue 171 VAL Chi-restraints excluded: chain E residue 264 VAL Chi-restraints excluded: chain E residue 293 THR Chi-restraints excluded: chain F residue 1 ASP Chi-restraints excluded: chain F residue 21 ILE Chi-restraints excluded: chain F residue 57 ILE Chi-restraints excluded: chain F residue 58 SER Chi-restraints excluded: chain F residue 71 VAL Chi-restraints excluded: chain F residue 87 THR Chi-restraints excluded: chain F residue 140 LEU Chi-restraints excluded: chain F residue 313 CYS Chi-restraints excluded: chain F residue 329 CYS Chi-restraints excluded: chain G residue 29 ILE Chi-restraints excluded: chain G residue 44 GLU Chi-restraints excluded: chain G residue 48 VAL Chi-restraints excluded: chain H residue 36 TRP Chi-restraints excluded: chain H residue 82 MET Chi-restraints excluded: chain H residue 89 VAL Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-2.0-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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-5881/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: 68.9924 > 50: distance: 15 - 48: 29.763 distance: 18 - 45: 13.715 distance: 33 - 38: 4.911 distance: 38 - 39: 5.708 distance: 39 - 40: 9.297 distance: 39 - 42: 6.651 distance: 40 - 41: 18.368 distance: 40 - 45: 20.018 distance: 42 - 43: 12.611 distance: 42 - 44: 16.720 distance: 45 - 46: 18.976 distance: 46 - 47: 15.361 distance: 46 - 49: 11.702 distance: 47 - 48: 13.115 distance: 47 - 52: 42.282 distance: 49 - 50: 24.489 distance: 49 - 51: 15.683 distance: 52 - 53: 13.019 distance: 52 - 58: 13.126 distance: 53 - 54: 45.365 distance: 53 - 56: 27.672 distance: 54 - 59: 42.546 distance: 56 - 57: 24.551 distance: 57 - 58: 39.410 distance: 59 - 60: 12.423 distance: 60 - 61: 38.056 distance: 60 - 63: 19.765 distance: 61 - 62: 22.324 distance: 61 - 65: 8.350 distance: 63 - 64: 22.506 distance: 65 - 66: 15.072 distance: 66 - 67: 13.770 distance: 66 - 69: 15.220 distance: 67 - 68: 34.157 distance: 67 - 71: 42.364 distance: 69 - 70: 33.797 distance: 71 - 72: 16.866 distance: 72 - 73: 19.032 distance: 72 - 75: 18.509 distance: 73 - 74: 27.922 distance: 73 - 77: 12.744 distance: 75 - 76: 4.833 distance: 77 - 78: 20.126 distance: 78 - 79: 9.899 distance: 78 - 81: 43.625 distance: 79 - 80: 14.304 distance: 79 - 85: 16.969 distance: 81 - 82: 17.202 distance: 82 - 83: 22.374 distance: 82 - 84: 10.095 distance: 85 - 86: 45.746 distance: 86 - 87: 10.667 distance: 87 - 88: 31.778 distance: 87 - 89: 28.923 distance: 89 - 90: 16.656 distance: 90 - 91: 11.449 distance: 90 - 93: 23.861 distance: 91 - 92: 15.453 distance: 91 - 96: 17.494 distance: 93 - 94: 16.449 distance: 93 - 95: 32.741 distance: 96 - 97: 22.801 distance: 97 - 98: 16.886 distance: 97 - 100: 37.695 distance: 98 - 99: 10.899 distance: 98 - 105: 30.806 distance: 100 - 101: 29.326 distance: 101 - 102: 4.152 distance: 102 - 103: 10.169 distance: 102 - 104: 21.209