Starting phenix.real_space_refine on Sat Jul 4 11:04:34 2026 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, restraint, /net/cci-nas-00/data/ceres_data/6wge_21663/07_2026/6wge_21663.cif Found real_map, /net/cci-nas-00/data/ceres_data/6wge_21663/07_2026/6wge_21663.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.9 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { model { file = "/net/cci-nas-00/data/ceres_data/6wge_21663/07_2026/6wge_21663.cif" } default_model = "/net/cci-nas-00/data/ceres_data/6wge_21663/07_2026/6wge_21663.cif" restraint_files = "/net/cci-nas-00/data/ceres_data/6wge_21663/07_2026/6wge_21663.cif" default_restraint = "/net/cci-nas-00/data/ceres_data/6wge_21663/07_2026/6wge_21663.cif" real_map_files = "/net/cci-nas-00/data/ceres_data/6wge_21663/07_2026/6wge_21663.map" default_real_map = "/net/cci-nas-00/data/ceres_data/6wge_21663/07_2026/6wge_21663.map" } resolution = 3.9 write_initial_geo_file = False refinement { macro_cycles = 10 } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= 0.000 sd= 0.006 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 6 Type Number sf(0) Gaussians P 92 5.49 5 Mg 2 5.21 5 S 106 5.16 5 C 12574 2.51 5 N 3470 2.21 5 O 3950 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 133 residue(s): 0.02s Monomer Library directory: "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/chem_data/mon_lib" Total number of atoms: 20194 Number of models: 1 Model: "" Number of chains: 8 Chain: "A" Number of atoms: 3099 Number of conformers: 1 Conformer: "" Number of residues, atoms: 390, 3099 Classifications: {'peptide': 390} Link IDs: {'PTRANS': 14, 'TRANS': 375} Chain breaks: 1 Chain: "B" Number of atoms: 4008 Number of conformers: 1 Conformer: "" Number of residues, atoms: 490, 4008 Classifications: {'peptide': 490} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 11, 'TRANS': 478} Chain breaks: 2 Chain: "C" Number of atoms: 1395 Number of conformers: 1 Conformer: "" Number of residues, atoms: 174, 1395 Classifications: {'peptide': 174} Link IDs: {'PTRANS': 6, 'TRANS': 167} Chain breaks: 2 Chain: "E" Number of atoms: 9865 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1226, 9865 Classifications: {'peptide': 1226} Link IDs: {'PTRANS': 35, 'TRANS': 1190} Chain breaks: 11 Chain: "F" Number of atoms: 903 Number of conformers: 1 Conformer: "" Number of residues, atoms: 43, 903 Classifications: {'DNA': 43} Link IDs: {'rna3p': 42} Chain: "G" Number of atoms: 860 Number of conformers: 1 Conformer: "" Number of residues, atoms: 43, 860 Classifications: {'DNA': 43} Link IDs: {'rna3p': 42} Chain: "A" Number of atoms: 32 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 32 Unusual residues: {' MG': 1, 'ANP': 1} Classifications: {'undetermined': 2} Link IDs: {None: 1} Chain: "B" Number of atoms: 32 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 32 Unusual residues: {' MG': 1, 'ANP': 1} Classifications: {'undetermined': 2} Link IDs: {None: 1} Time building chain proxies: 3.78, per 1000 atoms: 0.19 Number of scatterers: 20194 At special positions: 0 Unit cell: (147.4, 142.04, 172.86, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 6 Type Number sf(0) S 106 16.00 P 92 15.00 Mg 2 11.99 O 3950 8.00 N 3470 7.00 C 12574 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.00 Amino acid : True - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 1.79 Conformation dependent library (CDL) restraints added in 757.2 milliseconds 4480 Ramachandran restraints generated. 2240 Oldfield, 0 Emsley, 2240 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 4376 Finding SS restraints... Secondary structure from input PDB file: 97 helices and 10 sheets defined 68.1% alpha, 6.0% beta 35 base pairs and 84 stacking pairs defined. Time for finding SS restraints: 2.06 Creating SS restraints... Processing helix chain 'A' and resid 38 through 49 removed outlier: 3.639A pdb=" N MET A 42 " --> pdb=" O LYS A 38 " (cutoff:3.500A) removed outlier: 3.658A pdb=" N LEU A 49 " --> pdb=" O ILE A 45 " (cutoff:3.500A) Processing helix chain 'A' and resid 52 through 57 removed outlier: 3.625A pdb=" N ARG A 57 " --> pdb=" O SER A 54 " (cutoff:3.500A) Processing helix chain 'A' and resid 113 through 125 Processing helix chain 'A' and resid 138 through 145 removed outlier: 3.716A pdb=" N SER A 142 " --> pdb=" O GLY A 138 " (cutoff:3.500A) removed outlier: 4.061A pdb=" N ALA A 144 " --> pdb=" O VAL A 140 " (cutoff:3.500A) removed outlier: 3.824A pdb=" N MET A 145 " --> pdb=" O GLU A 141 " (cutoff:3.500A) Processing helix chain 'A' and resid 147 through 159 removed outlier: 3.755A pdb=" N SER A 159 " --> pdb=" O PHE A 155 " (cutoff:3.500A) Processing helix chain 'A' and resid 164 through 199 removed outlier: 3.593A pdb=" N TYR A 168 " --> pdb=" O LEU A 164 " (cutoff:3.500A) removed outlier: 3.595A pdb=" N ALA A 199 " --> pdb=" O GLU A 195 " (cutoff:3.500A) Processing helix chain 'A' and resid 1035 through 1089 removed outlier: 3.914A pdb=" N GLN A1057 " --> pdb=" O LYS A1053 " (cutoff:3.500A) removed outlier: 4.326A pdb=" N ALA A1058 " --> pdb=" O LYS A1054 " (cutoff:3.500A) Processing helix chain 'A' and resid 1105 through 1109 Processing helix chain 'A' and resid 1129 through 1148 removed outlier: 3.679A pdb=" N VAL A1135 " --> pdb=" O GLY A1131 " (cutoff:3.500A) removed outlier: 3.682A pdb=" N LYS A1148 " --> pdb=" O ILE A1144 " (cutoff:3.500A) Processing helix chain 'A' and resid 1163 through 1179 removed outlier: 4.337A pdb=" N THR A1179 " --> pdb=" O LYS A1175 " (cutoff:3.500A) Processing helix chain 'A' and resid 1190 through 1195 removed outlier: 3.700A pdb=" N TYR A1194 " --> pdb=" O LYS A1190 " (cutoff:3.500A) Processing helix chain 'B' and resid 38 through 50 removed outlier: 3.898A pdb=" N PHE B 42 " --> pdb=" O LYS B 38 " (cutoff:3.500A) removed outlier: 3.536A pdb=" N LEU B 49 " --> pdb=" O ILE B 45 " (cutoff:3.500A) Processing helix chain 'B' and resid 51 through 54 Processing helix chain 'B' and resid 57 through 63 Processing helix chain 'B' and resid 119 through 129 Processing helix chain 'B' and resid 144 through 148 Processing helix chain 'B' and resid 151 through 164 removed outlier: 3.614A pdb=" N ARG B 155 " --> pdb=" O PRO B 151 " (cutoff:3.500A) Processing helix chain 'B' and resid 166 through 236 removed outlier: 4.342A pdb=" N SER B 177 " --> pdb=" O GLU B 173 " (cutoff:3.500A) removed outlier: 3.981A pdb=" N LEU B 178 " --> pdb=" O GLU B 174 " (cutoff:3.500A) Processing helix chain 'B' and resid 935 through 947 removed outlier: 4.172A pdb=" N GLY B 947 " --> pdb=" O ILE B 943 " (cutoff:3.500A) Processing helix chain 'B' and resid 963 through 979 Processing helix chain 'B' and resid 985 through 1048 Processing helix chain 'B' and resid 1116 through 1135 Processing helix chain 'B' and resid 1144 through 1148 removed outlier: 3.769A pdb=" N GLN B1147 " --> pdb=" O GLN B1144 " (cutoff:3.500A) removed outlier: 4.535A pdb=" N ALA B1148 " --> pdb=" O ILE B1145 " (cutoff:3.500A) No H-bonds generated for 'chain 'B' and resid 1144 through 1148' Processing helix chain 'B' and resid 1150 through 1165 Processing helix chain 'B' and resid 1176 through 1182 removed outlier: 3.799A pdb=" N LEU B1179 " --> pdb=" O ARG B1176 " (cutoff:3.500A) removed outlier: 4.260A pdb=" N SER B1182 " --> pdb=" O LEU B1179 " (cutoff:3.500A) Processing helix chain 'B' and resid 1202 through 1211 removed outlier: 3.793A pdb=" N GLU B1211 " --> pdb=" O LYS B1207 " (cutoff:3.500A) Processing helix chain 'C' and resid 14 through 23 Processing helix chain 'C' and resid 28 through 35 removed outlier: 3.906A pdb=" N PHE C 33 " --> pdb=" O LYS C 29 " (cutoff:3.500A) removed outlier: 4.440A pdb=" N GLU C 34 " --> pdb=" O ALA C 30 " (cutoff:3.500A) Processing helix chain 'C' and resid 37 through 46 removed outlier: 4.198A pdb=" N SER C 43 " --> pdb=" O SER C 39 " (cutoff:3.500A) removed outlier: 3.576A pdb=" N SER C 46 " --> pdb=" O GLU C 42 " (cutoff:3.500A) Processing helix chain 'C' and resid 52 through 86 removed outlier: 3.566A pdb=" N LYS C 86 " --> pdb=" O PHE C 82 " (cutoff:3.500A) Processing helix chain 'C' and resid 559 through 572 Processing helix chain 'C' and resid 581 through 585 Processing helix chain 'C' and resid 589 through 606 removed outlier: 3.772A pdb=" N GLN C 606 " --> pdb=" O VAL C 602 " (cutoff:3.500A) Processing helix chain 'E' and resid 1198 through 1216 removed outlier: 3.537A pdb=" N MET E1216 " --> pdb=" O ASN E1212 " (cutoff:3.500A) Processing helix chain 'E' and resid 1234 through 1247 removed outlier: 3.551A pdb=" N LEU E1238 " --> pdb=" O GLY E1234 " (cutoff:3.500A) removed outlier: 3.937A pdb=" N LYS E1247 " --> pdb=" O SER E1243 " (cutoff:3.500A) Processing helix chain 'E' and resid 1248 through 1253 removed outlier: 3.892A pdb=" N GLY E1252 " --> pdb=" O ILE E1248 " (cutoff:3.500A) removed outlier: 3.936A pdb=" N ILE E1253 " --> pdb=" O LYS E1249 " (cutoff:3.500A) No H-bonds generated for 'chain 'E' and resid 1248 through 1253' Processing helix chain 'E' and resid 1257 through 1273 Processing helix chain 'E' and resid 1294 through 1319 removed outlier: 3.761A pdb=" N MET E1318 " --> pdb=" O THR E1314 " (cutoff:3.500A) Processing helix chain 'E' and resid 1329 through 1353 removed outlier: 3.962A pdb=" N ILE E1333 " --> pdb=" O ILE E1329 " (cutoff:3.500A) removed outlier: 5.621A pdb=" N TYR E1349 " --> pdb=" O GLN E1345 " (cutoff:3.500A) Proline residue: E1350 - end of helix Processing helix chain 'E' and resid 1381 through 1404 removed outlier: 3.546A pdb=" N GLN E1404 " --> pdb=" O LEU E1400 " (cutoff:3.500A) Processing helix chain 'E' and resid 1407 through 1420 removed outlier: 3.846A pdb=" N GLN E1413 " --> pdb=" O THR E1409 " (cutoff:3.500A) Processing helix chain 'E' and resid 1421 through 1424 removed outlier: 3.712A pdb=" N VAL E1424 " --> pdb=" O PRO E1421 " (cutoff:3.500A) No H-bonds generated for 'chain 'E' and resid 1421 through 1424' Processing helix chain 'E' and resid 1427 through 1445 removed outlier: 3.597A pdb=" N TYR E1445 " --> pdb=" O VAL E1441 " (cutoff:3.500A) Processing helix chain 'E' and resid 1448 through 1458 Processing helix chain 'E' and resid 1489 through 1501 removed outlier: 3.516A pdb=" N VAL E1501 " --> pdb=" O LEU E1497 " (cutoff:3.500A) Processing helix chain 'E' and resid 1525 through 1551 removed outlier: 4.505A pdb=" N SER E1551 " --> pdb=" O LYS E1547 " (cutoff:3.500A) Processing helix chain 'E' and resid 1558 through 1572 Processing helix chain 'E' and resid 1580 through 1597 Processing helix chain 'E' and resid 1603 through 1628 Processing helix chain 'E' and resid 1647 through 1666 removed outlier: 3.590A pdb=" N GLU E1664 " --> pdb=" O ASP E1660 " (cutoff:3.500A) Processing helix chain 'E' and resid 1668 through 1690 removed outlier: 3.976A pdb=" N SER E1672 " --> pdb=" O SER E1668 " (cutoff:3.500A) Processing helix chain 'E' and resid 1709 through 1728 Processing helix chain 'E' and resid 1749 through 1764 Processing helix chain 'E' and resid 1766 through 1768 No H-bonds generated for 'chain 'E' and resid 1766 through 1768' Processing helix chain 'E' and resid 1769 through 1781 removed outlier: 3.545A pdb=" N LEU E1781 " --> pdb=" O ILE E1777 " (cutoff:3.500A) Processing helix chain 'E' and resid 1785 through 1803 removed outlier: 3.535A pdb=" N ARG E1789 " --> pdb=" O ALA E1785 " (cutoff:3.500A) Processing helix chain 'E' and resid 1809 through 1821 Processing helix chain 'E' and resid 1824 through 1842 Processing helix chain 'E' and resid 1849 through 1855 Processing helix chain 'E' and resid 1856 through 1858 No H-bonds generated for 'chain 'E' and resid 1856 through 1858' Processing helix chain 'E' and resid 1862 through 1879 Processing helix chain 'E' and resid 1884 through 1895 Processing helix chain 'E' and resid 1896 through 1898 No H-bonds generated for 'chain 'E' and resid 1896 through 1898' Processing helix chain 'E' and resid 1900 through 1914 Processing helix chain 'E' and resid 1922 through 1941 removed outlier: 3.965A pdb=" N MET E1926 " --> pdb=" O ASP E1922 " (cutoff:3.500A) Processing helix chain 'E' and resid 1945 through 1958 Processing helix chain 'E' and resid 1965 through 1986 removed outlier: 3.835A pdb=" N ALA E1970 " --> pdb=" O PRO E1966 " (cutoff:3.500A) removed outlier: 3.955A pdb=" N GLU E1986 " --> pdb=" O ILE E1982 " (cutoff:3.500A) Processing helix chain 'E' and resid 1999 through 2016 removed outlier: 3.518A pdb=" N ARG E2016 " --> pdb=" O PHE E2012 " (cutoff:3.500A) Processing helix chain 'E' and resid 2017 through 2021 Processing helix chain 'E' and resid 2023 through 2033 Proline residue: E2029 - end of helix removed outlier: 3.824A pdb=" N THR E2032 " --> pdb=" O GLN E2028 " (cutoff:3.500A) removed outlier: 4.116A pdb=" N THR E2033 " --> pdb=" O PRO E2029 " (cutoff:3.500A) Processing helix chain 'E' and resid 2037 through 2055 Processing helix chain 'E' and resid 2056 through 2058 No H-bonds generated for 'chain 'E' and resid 2056 through 2058' Processing helix chain 'E' and resid 2062 through 2078 Processing helix chain 'E' and resid 2081 through 2099 removed outlier: 4.043A pdb=" N CYS E2091 " --> pdb=" O HIS E2087 " (cutoff:3.500A) Processing helix chain 'E' and resid 2102 through 2126 removed outlier: 4.465A pdb=" N VAL E2106 " --> pdb=" O ASN E2102 " (cutoff:3.500A) removed outlier: 4.181A pdb=" N GLU E2126 " --> pdb=" O SER E2122 " (cutoff:3.500A) Processing helix chain 'E' and resid 2132 through 2154 Proline residue: E2138 - end of helix Processing helix chain 'E' and resid 2158 through 2163 removed outlier: 4.294A pdb=" N LYS E2162 " --> pdb=" O LEU E2158 " (cutoff:3.500A) removed outlier: 3.881A pdb=" N GLY E2163 " --> pdb=" O GLU E2159 " (cutoff:3.500A) No H-bonds generated for 'chain 'E' and resid 2158 through 2163' Processing helix chain 'E' and resid 2168 through 2182 removed outlier: 4.148A pdb=" N LYS E2182 " --> pdb=" O MET E2178 " (cutoff:3.500A) Processing helix chain 'E' and resid 2186 through 2203 removed outlier: 3.801A pdb=" N GLN E2202 " --> pdb=" O PHE E2198 " (cutoff:3.500A) Processing helix chain 'E' and resid 2204 through 2208 Processing helix chain 'E' and resid 2209 through 2220 Processing helix chain 'E' and resid 2227 through 2253 removed outlier: 3.602A pdb=" N ILE E2231 " --> pdb=" O VAL E2227 " (cutoff:3.500A) Processing helix chain 'E' and resid 2258 through 2262 Processing helix chain 'E' and resid 2273 through 2289 removed outlier: 4.300A pdb=" N MET E2279 " --> pdb=" O SER E2275 " (cutoff:3.500A) removed outlier: 3.651A pdb=" N GLN E2280 " --> pdb=" O SER E2276 " (cutoff:3.500A) removed outlier: 3.691A pdb=" N LEU E2281 " --> pdb=" O SER E2277 " (cutoff:3.500A) removed outlier: 3.703A pdb=" N TYR E2282 " --> pdb=" O ILE E2278 " (cutoff:3.500A) removed outlier: 5.780A pdb=" N LYS E2284 " --> pdb=" O GLN E2280 " (cutoff:3.500A) removed outlier: 6.194A pdb=" N GLN E2285 " --> pdb=" O LEU E2281 " (cutoff:3.500A) Processing helix chain 'E' and resid 2295 through 2310 removed outlier: 3.750A pdb=" N ASN E2310 " --> pdb=" O ALA E2306 " (cutoff:3.500A) Processing helix chain 'E' and resid 2320 through 2326 Processing helix chain 'E' and resid 2331 through 2349 removed outlier: 3.712A pdb=" N ARG E2335 " --> pdb=" O GLU E2331 " (cutoff:3.500A) removed outlier: 3.774A pdb=" N LYS E2348 " --> pdb=" O GLU E2344 " (cutoff:3.500A) removed outlier: 3.732A pdb=" N TYR E2349 " --> pdb=" O ILE E2345 " (cutoff:3.500A) Processing helix chain 'E' and resid 2351 through 2356 removed outlier: 3.723A pdb=" N MET E2355 " --> pdb=" O GLY E2351 " (cutoff:3.500A) removed outlier: 3.672A pdb=" N LYS E2356 " --> pdb=" O PHE E2352 " (cutoff:3.500A) No H-bonds generated for 'chain 'E' and resid 2351 through 2356' Processing helix chain 'E' and resid 2357 through 2372 removed outlier: 3.512A pdb=" N MET E2361 " --> pdb=" O ALA E2357 " (cutoff:3.500A) Processing helix chain 'E' and resid 2392 through 2400 removed outlier: 3.655A pdb=" N TYR E2396 " --> pdb=" O CYS E2392 " (cutoff:3.500A) Processing helix chain 'E' and resid 2402 through 2418 Processing helix chain 'E' and resid 2424 through 2438 removed outlier: 3.638A pdb=" N LEU E2429 " --> pdb=" O VAL E2425 " (cutoff:3.500A) Processing helix chain 'E' and resid 2443 through 2471 removed outlier: 3.803A pdb=" N LEU E2447 " --> pdb=" O GLN E2443 " (cutoff:3.500A) Processing helix chain 'E' and resid 2534 through 2558 Processing helix chain 'E' and resid 2560 through 2567 Processing helix chain 'E' and resid 2588 through 2601 removed outlier: 3.623A pdb=" N ASP E2600 " --> pdb=" O PHE E2596 " (cutoff:3.500A) Processing helix chain 'E' and resid 2607 through 2626 removed outlier: 3.701A pdb=" N LYS E2611 " --> pdb=" O THR E2607 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 18 through 20 removed outlier: 3.978A pdb=" N GLN A 19 " --> pdb=" O ILE A 9 " (cutoff:3.500A) removed outlier: 3.802A pdb=" N ILE A 9 " --> pdb=" O GLN A 19 " (cutoff:3.500A) removed outlier: 6.923A pdb=" N VAL A 83 " --> pdb=" O LYS A 5 " (cutoff:3.500A) removed outlier: 4.530A pdb=" N ILE A 7 " --> pdb=" O SER A 81 " (cutoff:3.500A) removed outlier: 6.493A pdb=" N SER A 81 " --> pdb=" O ILE A 7 " (cutoff:3.500A) removed outlier: 4.245A pdb=" N ILE A 9 " --> pdb=" O PHE A 79 " (cutoff:3.500A) removed outlier: 6.467A pdb=" N PHE A 79 " --> pdb=" O ILE A 9 " (cutoff:3.500A) removed outlier: 4.425A pdb=" N GLU A 86 " --> pdb=" O GLU A 90 " (cutoff:3.500A) removed outlier: 4.616A pdb=" N GLU A 90 " --> pdb=" O GLU A 86 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 18 through 20 removed outlier: 3.978A pdb=" N GLN A 19 " --> pdb=" O ILE A 9 " (cutoff:3.500A) removed outlier: 3.802A pdb=" N ILE A 9 " --> pdb=" O GLN A 19 " (cutoff:3.500A) removed outlier: 6.923A pdb=" N VAL A 83 " --> pdb=" O LYS A 5 " (cutoff:3.500A) removed outlier: 4.530A pdb=" N ILE A 7 " --> pdb=" O SER A 81 " (cutoff:3.500A) removed outlier: 6.493A pdb=" N SER A 81 " --> pdb=" O ILE A 7 " (cutoff:3.500A) removed outlier: 4.245A pdb=" N ILE A 9 " --> pdb=" O PHE A 79 " (cutoff:3.500A) removed outlier: 6.467A pdb=" N PHE A 79 " --> pdb=" O ILE A 9 " (cutoff:3.500A) removed outlier: 3.883A pdb=" N VAL A 80 " --> pdb=" O ARG A 96 " (cutoff:3.500A) removed outlier: 3.672A pdb=" N LYS A 106 " --> pdb=" O ALA A 95 " (cutoff:3.500A) removed outlier: 4.440A pdb=" N TYR A 105 " --> pdb=" O VAL A 112 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 27 through 29 removed outlier: 6.638A pdb=" N THR A 28 " --> pdb=" O VAL A1186 " (cutoff:3.500A) removed outlier: 3.931A pdb=" N ILE A1185 " --> pdb=" O PHE A1153 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'A' and resid 1096 through 1100 Processing sheet with id=AA5, first strand: chain 'A' and resid 1199 through 1206 Processing sheet with id=AA6, first strand: chain 'B' and resid 16 through 17 removed outlier: 3.798A pdb=" N VAL B 78 " --> pdb=" O ARG B 100 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 28 through 30 removed outlier: 3.865A pdb=" N VAL B 139 " --> pdb=" O LEU B1141 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'B' and resid 1052 through 1057 Processing sheet with id=AA9, first strand: chain 'B' and resid 1188 through 1192 Processing sheet with id=AB1, first strand: chain 'C' and resid 579 through 580 1165 hydrogen bonds defined for protein. 3408 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 70 hydrogen bonds 140 hydrogen bond angles 0 basepair planarities 35 basepair parallelities 84 stacking parallelities Total time for adding SS restraints: 4.94 Time building geometry restraints manager: 2.35 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.21 - 1.33: 4746 1.33 - 1.45: 4470 1.45 - 1.58: 11132 1.58 - 1.70: 181 1.70 - 1.82: 179 Bond restraints: 20708 Sorted by residual: bond pdb=" N3B ANP B2000 " pdb=" PG ANP B2000 " ideal model delta sigma weight residual 1.795 1.618 0.177 2.00e-02 2.50e+03 7.81e+01 bond pdb=" N3B ANP A2000 " pdb=" PG ANP A2000 " ideal model delta sigma weight residual 1.795 1.627 0.168 2.00e-02 2.50e+03 7.07e+01 bond pdb=" O3A ANP B2000 " pdb=" PB ANP B2000 " ideal model delta sigma weight residual 1.700 1.591 0.109 2.00e-02 2.50e+03 2.98e+01 bond pdb=" O3A ANP A2000 " pdb=" PB ANP A2000 " ideal model delta sigma weight residual 1.700 1.601 0.099 2.00e-02 2.50e+03 2.45e+01 bond pdb=" C3' DT G 15 " pdb=" O3' DT G 15 " ideal model delta sigma weight residual 1.435 1.380 0.055 1.30e-02 5.92e+03 1.78e+01 ... (remaining 20703 not shown) Histogram of bond angle deviations from ideal: 0.00 - 3.05: 27744 3.05 - 6.11: 441 6.11 - 9.16: 51 9.16 - 12.21: 7 12.21 - 15.27: 3 Bond angle restraints: 28246 Sorted by residual: angle pdb=" O4' DA F 39 " pdb=" C1' DA F 39 " pdb=" N9 DA F 39 " ideal model delta sigma weight residual 108.00 102.88 5.12 7.00e-01 2.04e+00 5.36e+01 angle pdb=" O5' DT G 25 " pdb=" C5' DT G 25 " pdb=" C4' DT G 25 " ideal model delta sigma weight residual 109.40 114.41 -5.01 8.00e-01 1.56e+00 3.93e+01 angle pdb=" PB ANP B2000 " pdb=" N3B ANP B2000 " pdb=" PG ANP B2000 " ideal model delta sigma weight residual 126.95 111.83 15.12 3.00e+00 1.11e-01 2.54e+01 angle pdb=" O4' DA F 38 " pdb=" C1' DA F 38 " pdb=" N9 DA F 38 " ideal model delta sigma weight residual 108.00 104.60 3.40 7.00e-01 2.04e+00 2.36e+01 angle pdb=" O4' DA F 8 " pdb=" C1' DA F 8 " pdb=" N9 DA F 8 " ideal model delta sigma weight residual 108.00 104.71 3.29 7.00e-01 2.04e+00 2.21e+01 ... (remaining 28241 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 35.30: 11830 35.30 - 70.59: 706 70.59 - 105.89: 20 105.89 - 141.19: 0 141.19 - 176.49: 2 Dihedral angle restraints: 12558 sinusoidal: 5849 harmonic: 6709 Sorted by residual: dihedral pdb=" CA GLU B1111 " pdb=" C GLU B1111 " pdb=" N MET B1112 " pdb=" CA MET B1112 " ideal model delta harmonic sigma weight residual -180.00 -152.45 -27.55 0 5.00e+00 4.00e-02 3.04e+01 dihedral pdb=" CA VAL E2317 " pdb=" C VAL E2317 " pdb=" N GLN E2318 " pdb=" CA GLN E2318 " ideal model delta harmonic sigma weight residual 180.00 153.03 26.97 0 5.00e+00 4.00e-02 2.91e+01 dihedral pdb=" CA GLY E2351 " pdb=" C GLY E2351 " pdb=" N PHE E2352 " pdb=" CA PHE E2352 " ideal model delta harmonic sigma weight residual -180.00 -153.21 -26.79 0 5.00e+00 4.00e-02 2.87e+01 ... (remaining 12555 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.076: 2674 0.076 - 0.152: 472 0.152 - 0.228: 39 0.228 - 0.304: 8 0.304 - 0.380: 2 Chirality restraints: 3195 Sorted by residual: chirality pdb=" C3' ANP A2000 " pdb=" C2' ANP A2000 " pdb=" C4' ANP A2000 " pdb=" O3' ANP A2000 " both_signs ideal model delta sigma weight residual False -2.36 -2.74 0.38 2.00e-01 2.50e+01 3.60e+00 chirality pdb=" C3' ANP B2000 " pdb=" C2' ANP B2000 " pdb=" C4' ANP B2000 " pdb=" O3' ANP B2000 " both_signs ideal model delta sigma weight residual False -2.36 -2.71 0.35 2.00e-01 2.50e+01 3.04e+00 chirality pdb=" CB ILE B1098 " pdb=" CA ILE B1098 " pdb=" CG1 ILE B1098 " pdb=" CG2 ILE B1098 " both_signs ideal model delta sigma weight residual False 2.64 2.37 0.27 2.00e-01 2.50e+01 1.84e+00 ... (remaining 3192 not shown) Planarity restraints: 3284 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA GLU E1986 " -0.015 2.00e-02 2.50e+03 2.94e-02 8.64e+00 pdb=" C GLU E1986 " 0.051 2.00e-02 2.50e+03 pdb=" O GLU E1986 " -0.019 2.00e-02 2.50e+03 pdb=" N GLU E1987 " -0.017 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA ALA E2408 " 0.013 2.00e-02 2.50e+03 2.60e-02 6.77e+00 pdb=" C ALA E2408 " -0.045 2.00e-02 2.50e+03 pdb=" O ALA E2408 " 0.017 2.00e-02 2.50e+03 pdb=" N PHE E2409 " 0.015 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C ALA A1117 " -0.040 5.00e-02 4.00e+02 6.04e-02 5.83e+00 pdb=" N PRO A1118 " 0.104 5.00e-02 4.00e+02 pdb=" CA PRO A1118 " -0.031 5.00e-02 4.00e+02 pdb=" CD PRO A1118 " -0.034 5.00e-02 4.00e+02 ... (remaining 3281 not shown) Histogram of nonbonded interaction distances: 1.92 - 2.52: 146 2.52 - 3.11: 15154 3.11 - 3.71: 31071 3.71 - 4.30: 40894 4.30 - 4.90: 65934 Nonbonded interactions: 153199 Sorted by model distance: nonbonded pdb=" NE2 GLN B 141 " pdb="MG MG B2001 " model vdw 1.921 2.250 nonbonded pdb=" O3G ANP B2000 " pdb="MG MG B2001 " model vdw 1.924 2.170 nonbonded pdb=" OE1 GLN B 141 " pdb="MG MG B2001 " model vdw 2.059 2.170 nonbonded pdb=" NE2 GLN A 137 " pdb="MG MG A2001 " model vdw 2.160 2.250 nonbonded pdb=" O ILE E1315 " pdb=" OG1 THR E1319 " model vdw 2.184 3.040 ... (remaining 153194 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.00 ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.900 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.020 Extract box with map and model: 0.280 Check model and map are aligned: 0.060 Set scattering table: 0.050 Process input model: 20.370 Find NCS groups from input model: 0.090 Set up NCS constraints: 0.040 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.290 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 23.110 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8457 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.177 20708 Z= 0.431 Angle : 1.011 15.265 28246 Z= 0.615 Chirality : 0.058 0.380 3195 Planarity : 0.006 0.060 3284 Dihedral : 19.176 176.485 8182 Min Nonbonded Distance : 1.921 Molprobity Statistics. All-atom Clashscore : 6.66 Ramachandran Plot: Outliers : 0.04 % Allowed : 8.35 % Favored : 91.61 % Rotamer: Outliers : 0.15 % Allowed : 4.03 % Favored : 95.83 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z values with (uncertainties): Interpretation: poor |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores below are scaled independently, so they are not related in a simple way. whole: -2.35 (0.16), residues: 2240 helix: -0.67 (0.12), residues: 1382 sheet: -2.45 (0.39), residues: 157 loop : -3.01 (0.22), residues: 701 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.016 0.001 ARG E2407 TYR 0.026 0.003 TYR C 67 PHE 0.032 0.003 PHE E2161 TRP 0.016 0.002 TRP C 18 HIS 0.020 0.002 HIS B 66 Details of bonding type rmsd/Z covalent geometry : bond 0.00767 / 0.43 (20708) covalent geometry : angle 1.01146 / 0.61 (28246) hydrogen bonds : bond 0.13258 / 8.97 ( 1232) hydrogen bonds : angle 5.42272 / 3.84 ( 3548) *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4480 Ramachandran restraints generated. 2240 Oldfield, 0 Emsley, 2240 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4480 Ramachandran restraints generated. 2240 Oldfield, 0 Emsley, 2240 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 224 residues out of total 2037 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 221 time to evaluate : 0.754 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 1108 LEU cc_start: 0.9202 (tp) cc_final: 0.8823 (tp) REVERT: B 62 LEU cc_start: 0.9557 (tt) cc_final: 0.9314 (tp) REVERT: B 946 LEU cc_start: 0.8502 (mt) cc_final: 0.8271 (mt) REVERT: C 33 PHE cc_start: 0.8525 (t80) cc_final: 0.8022 (t80) REVERT: C 44 ILE cc_start: 0.8040 (mm) cc_final: 0.7742 (mm) REVERT: C 80 GLU cc_start: 0.8781 (tp30) cc_final: 0.8490 (tp30) REVERT: C 163 PHE cc_start: 0.8165 (m-80) cc_final: 0.7942 (m-80) REVERT: C 563 MET cc_start: 0.7166 (tpp) cc_final: 0.6725 (tpp) REVERT: C 583 GLU cc_start: 0.8317 (tp30) cc_final: 0.8062 (tp30) REVERT: E 1193 MET cc_start: 0.6705 (mtp) cc_final: 0.5866 (mtm) REVERT: E 1254 MET cc_start: 0.7620 (ptp) cc_final: 0.7093 (tpp) REVERT: E 1312 LEU cc_start: 0.9340 (mm) cc_final: 0.9065 (mm) REVERT: E 1859 ASP cc_start: 0.7551 (t0) cc_final: 0.7346 (t0) REVERT: E 2025 MET cc_start: 0.8667 (mpp) cc_final: 0.8081 (mtt) REVERT: E 2207 MET cc_start: 0.8982 (mmt) cc_final: 0.8533 (mmt) REVERT: E 2624 LEU cc_start: 0.8581 (mt) cc_final: 0.7916 (mt) outliers start: 3 outliers final: 1 residues processed: 224 average time/residue: 0.1462 time to fit residues: 50.8757 Evaluate side-chains 152 residues out of total 2037 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 1 poor density : 151 time to evaluate : 0.740 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain E residue 2424 ASP Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 235 random chunks: chunk 216 optimal weight: 0.9980 chunk 98 optimal weight: 10.0000 chunk 194 optimal weight: 3.9990 chunk 227 optimal weight: 50.0000 chunk 107 optimal weight: 7.9990 chunk 10 optimal weight: 10.0000 chunk 66 optimal weight: 9.9990 chunk 130 optimal weight: 8.9990 chunk 124 optimal weight: 0.8980 chunk 103 optimal weight: 0.7980 chunk 200 optimal weight: 0.0970 overall best weight: 1.3580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... A 66 HIS A 187 HIS A1114 ASN A1166 ASN B 134 ASN B 958 GLN B 973 ASN ** B1114 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 68 HIS C 158 ASN C 587 ASN E1343 HIS E1539 GLN E1567 GLN E1813 GLN E2086 GLN E2124 HIS E2153 HIS E2224 ASN E2236 ASN E2292 HIS E2294 GLN E2315 HIS E2372 ASN ** E2394 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** E2434 ASN E2551 GLN ** E2552 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 25 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3496 r_free = 0.3496 target = 0.094850 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 35)----------------| | r_work = 0.3064 r_free = 0.3064 target = 0.072018 restraints weight = 54527.876| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 34)----------------| | r_work = 0.3070 r_free = 0.3070 target = 0.072738 restraints weight = 37528.217| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 27)----------------| | r_work = 0.3078 r_free = 0.3078 target = 0.073160 restraints weight = 29225.679| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 4; iterations = 23)----------------| | r_work = 0.3084 r_free = 0.3084 target = 0.073466 restraints weight = 25767.220| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 5; iterations = 23)----------------| | r_work = 0.3086 r_free = 0.3086 target = 0.073544 restraints weight = 23569.587| |-----------------------------------------------------------------------------| r_work (final): 0.3075 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8422 moved from start: 0.1217 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.035 20708 Z= 0.163 Angle : 0.646 11.388 28246 Z= 0.349 Chirality : 0.041 0.229 3195 Planarity : 0.004 0.041 3284 Dihedral : 19.551 173.462 3489 Min Nonbonded Distance : 1.979 Molprobity Statistics. All-atom Clashscore : 7.71 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.87 % Favored : 95.13 % Rotamer: Outliers : 1.08 % Allowed : 13.16 % Favored : 85.76 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z values with (uncertainties): Interpretation: poor |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores below are scaled independently, so they are not related in a simple way. whole: -0.39 (0.18), residues: 2240 helix: 1.01 (0.13), residues: 1397 sheet: -1.67 (0.41), residues: 154 loop : -2.59 (0.22), residues: 689 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.004 0.000 ARG E2407 TYR 0.019 0.001 TYR B 137 PHE 0.019 0.001 PHE E1596 TRP 0.011 0.001 TRP C 18 HIS 0.009 0.001 HIS E2627 Details of bonding type rmsd/Z covalent geometry : bond 0.00317 / 0.16 (20708) covalent geometry : angle 0.64598 / 0.35 (28246) hydrogen bonds : bond 0.07339 / 4.95 ( 1232) hydrogen bonds : angle 4.54836 / 3.23 ( 3548) *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4480 Ramachandran restraints generated. 2240 Oldfield, 0 Emsley, 2240 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4480 Ramachandran restraints generated. 2240 Oldfield, 0 Emsley, 2240 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 201 residues out of total 2037 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 22 poor density : 179 time to evaluate : 0.742 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 175 MET cc_start: 0.7402 (tmm) cc_final: 0.7128 (tmm) REVERT: B 72 ARG cc_start: 0.6233 (mtp180) cc_final: 0.5902 (mtp180) REVERT: B 946 LEU cc_start: 0.8482 (mt) cc_final: 0.8146 (mt) REVERT: C 33 PHE cc_start: 0.8491 (t80) cc_final: 0.7893 (t80) REVERT: C 44 ILE cc_start: 0.8010 (mm) cc_final: 0.7724 (mm) REVERT: C 163 PHE cc_start: 0.8290 (m-80) cc_final: 0.8026 (m-80) REVERT: E 1193 MET cc_start: 0.6352 (mtp) cc_final: 0.5726 (mtm) REVERT: E 1251 MET cc_start: 0.7972 (mpp) cc_final: 0.7559 (mpp) REVERT: E 1254 MET cc_start: 0.7717 (ptp) cc_final: 0.7123 (tpp) REVERT: E 1312 LEU cc_start: 0.9327 (mm) cc_final: 0.9010 (mm) REVERT: E 1323 MET cc_start: 0.7988 (mtp) cc_final: 0.7770 (mtp) REVERT: E 2025 MET cc_start: 0.8612 (mpp) cc_final: 0.7992 (mtt) REVERT: E 2207 MET cc_start: 0.8951 (mmt) cc_final: 0.8705 (mmt) REVERT: E 2334 MET cc_start: 0.8269 (tpp) cc_final: 0.7918 (mtp) REVERT: E 2624 LEU cc_start: 0.8280 (mt) cc_final: 0.8036 (mt) outliers start: 22 outliers final: 11 residues processed: 191 average time/residue: 0.1355 time to fit residues: 42.1306 Evaluate side-chains 160 residues out of total 2037 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 11 poor density : 149 time to evaluate : 0.727 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 143 ILE Chi-restraints excluded: chain A residue 1154 VAL Chi-restraints excluded: chain B residue 147 MET Chi-restraints excluded: chain B residue 940 MET Chi-restraints excluded: chain B residue 960 LEU Chi-restraints excluded: chain B residue 1014 ILE Chi-restraints excluded: chain C residue 83 ILE Chi-restraints excluded: chain E residue 1688 THR Chi-restraints excluded: chain E residue 1812 MET Chi-restraints excluded: chain E residue 2015 ILE Chi-restraints excluded: chain E residue 2621 PHE Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 235 random chunks: chunk 218 optimal weight: 0.9990 chunk 100 optimal weight: 6.9990 chunk 115 optimal weight: 8.9990 chunk 211 optimal weight: 5.9990 chunk 54 optimal weight: 0.4980 chunk 77 optimal weight: 0.9990 chunk 126 optimal weight: 5.9990 chunk 167 optimal weight: 0.0970 chunk 206 optimal weight: 1.9990 chunk 128 optimal weight: 4.9990 chunk 118 optimal weight: 1.9990 overall best weight: 0.9184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... B 134 ASN C 613 GLN E1343 HIS E1977 ASN E2023 HIS E2202 GLN ** E2552 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.3521 r_free = 0.3521 target = 0.096257 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 39)----------------| | r_work = 0.3085 r_free = 0.3085 target = 0.073111 restraints weight = 54972.956| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 28)----------------| | r_work = 0.3100 r_free = 0.3100 target = 0.074182 restraints weight = 36852.727| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 27)----------------| | r_work = 0.3110 r_free = 0.3110 target = 0.074743 restraints weight = 26439.382| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 4; iterations = 23)----------------| | r_work = 0.3116 r_free = 0.3116 target = 0.075048 restraints weight = 23708.892| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 5; iterations = 26)----------------| | r_work = 0.3119 r_free = 0.3119 target = 0.075162 restraints weight = 21982.443| |-----------------------------------------------------------------------------| r_work (final): 0.3107 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8373 moved from start: 0.1669 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.031 20708 Z= 0.138 Angle : 0.587 9.241 28246 Z= 0.317 Chirality : 0.040 0.234 3195 Planarity : 0.003 0.041 3284 Dihedral : 19.448 178.340 3487 Min Nonbonded Distance : 1.992 Molprobity Statistics. All-atom Clashscore : 6.89 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.18 % Favored : 94.82 % Rotamer: Outliers : 2.06 % Allowed : 15.66 % Favored : 82.28 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z values with (uncertainties): Interpretation: poor |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores below are scaled independently, so they are not related in a simple way. whole: 0.52 (0.18), residues: 2240 helix: 1.74 (0.14), residues: 1393 sheet: -1.12 (0.41), residues: 155 loop : -2.32 (0.23), residues: 692 Max deviation from planes: Type MaxDev MeanDev LineInFile ARG 0.007 0.000 ARG B1110 TYR 0.027 0.001 TYR B 2 PHE 0.019 0.001 PHE E1596 TRP 0.008 0.001 TRP C 18 HIS 0.006 0.001 HIS E2627 Details of bonding type rmsd/Z covalent geometry : bond 0.00254 / 0.14 (20708) covalent geometry : angle 0.58684 / 0.32 (28246) hydrogen bonds : bond 0.06421 / 4.34 ( 1232) hydrogen bonds : angle 4.22595 / 3.00 ( 3548) *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 4480 Ramachandran restraints generated. 2240 Oldfield, 0 Emsley, 2240 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 4480 Ramachandran restraints generated. 2240 Oldfield, 0 Emsley, 2240 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 223 residues out of total 2037 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 42 poor density : 181 time to evaluate : 0.698 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 175 MET cc_start: 0.7363 (tmm) cc_final: 0.7097 (tmm) REVERT: B 946 LEU cc_start: 0.8416 (mt) cc_final: 0.8119 (mt) REVERT: B 1108 GLU cc_start: 0.6720 (mp0) cc_final: 0.6437 (mt-10) REVERT: C 33 PHE cc_start: 0.8442 (t80) cc_final: 0.7757 (t80) REVERT: C 163 PHE cc_start: 0.8310 (m-80) cc_final: 0.8012 (m-80) REVERT: E 1193 MET cc_start: 0.6213 (mtp) cc_final: 0.5671 (mtm) REVERT: E 1251 MET cc_start: 0.8088 (mpp) cc_final: 0.7722 (mpp) REVERT: E 1254 MET cc_start: 0.7737 (ptp) cc_final: 0.7202 (tpp) REVERT: E 1264 LYS cc_start: 0.8050 (ptpt) cc_final: 0.7706 (ptpt) REVERT: E 1312 LEU cc_start: 0.9352 (mm) cc_final: 0.9012 (mm) REVERT: E 1323 MET cc_start: 0.8412 (mtp) cc_final: 0.8190 (mtp) REVERT: E 1570 LEU cc_start: 0.8044 (mt) cc_final: 0.7771 (mt) REVERT: E 1858 LEU cc_start: 0.8299 (OUTLIER) cc_final: 0.7980 (mt) REVERT: E 2326 MET cc_start: 0.8591 (tpp) cc_final: 0.8100 (tpp) REVERT: E 2339 ASP cc_start: 0.8338 (OUTLIER) cc_final: 0.7864 (t70) REVERT: E 2361 MET cc_start: 0.8989 (mmm) cc_final: 0.8593 (tpp) outliers start: 42 outliers final: 22 residues processed: 210 average time/residue: 0.1129 time to fit residues: 39.9371 Evaluate side-chains 181 residues out of total 2037 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 24 poor density : 157 time to evaluate : 0.748 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 143 ILE Chi-restraints excluded: chain A residue 1154 VAL Chi-restraints excluded: chain A residue 1217 THR Chi-restraints excluded: chain B residue 29 VAL Chi-restraints excluded: chain B residue 147 MET Chi-restraints excluded: chain B residue 940 MET Chi-restraints excluded: chain B residue 960 LEU Chi-restraints excluded: chain B residue 1014 ILE Chi-restraints excluded: chain C residue 83 ILE Chi-restraints excluded: chain C residue 612 THR Chi-restraints excluded: chain E residue 1309 ASP Chi-restraints excluded: chain E residue 1320 SER Chi-restraints excluded: chain E residue 1343 HIS Chi-restraints excluded: chain E residue 1437 LEU Chi-restraints excluded: chain E residue 1747 THR Chi-restraints excluded: chain E residue 1858 LEU Chi-restraints excluded: chain E residue 1956 LEU Chi-restraints excluded: chain E residue 1977 ASN Chi-restraints excluded: chain E residue 2085 VAL Chi-restraints excluded: chain E residue 2201 ILE Chi-restraints excluded: chain E residue 2339 ASP Chi-restraints excluded: chain E residue 2344 GLU Chi-restraints excluded: chain E residue 2426 THR Chi-restraints excluded: chain E residue 2621 PHE Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/phenix/command_line/real_space_refine.py", line 9, in run_program(real_space_refine.Program) File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/iotbx/cli_parser.py", line 987, in run_program task.run() File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/phenix/programs/real_space_refine.py", line 206, in run self.rsr_result = rsr_wrappers.run_real_space_refine( ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/phenix/refinement/rsr/wrappers.py", line 52, in __init__ self._refined = macro_cycle_real_space.run( ^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/phenix/refinement/macro_cycle_real_space.py", line 306, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/phenix/refinement/macro_cycle_real_space.py", line 702, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/phenix/refinement/macro_cycle_real_space.py", line 755, in minimization_no_ncs ro = mmtbx.refinement.real_space.individual_sites.easy( ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 44, in __init__ self.weight = mmtbx.refinement.real_space.weight.run( ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/mmtbx/refinement/real_space/weight.py", line 84, in __init__ ro.refine( File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 362, in refine real_space_result = refinery( ^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 233, in __init__ refiner.refine( File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/mmtbx/refinement/real_space/individual_sites.py", line 100, in refine self.refined = maptbx.real_space_refinement_simple.lbfgs( ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/cctbx/maptbx/real_space_refinement_simple.py", line 168, in __init__ O.minimizer = scitbx.lbfgs.run( ^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/scitbx/lbfgs/__init__.py", line 270, in run return run_c_plus_plus( ^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/scitbx/lbfgs/__init__.py", line 131, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/cctbx/maptbx/real_space_refinement_simple.py", line 243, in compute_functional_and_gradients gr_e = O.geometry_restraints_manager.energies_sites( ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/cctbx/geometry_restraints/manager.py", line 1602, in energies_sites pair_proxies = self.pair_proxies( ^^^^^^^^^^^^^^^^^^ File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/cctbx/geometry_restraints/manager.py", line 1505, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer3/home/dcliebschner/04_cryoem/Phenix-dev-6108/lib/python3.11/site-packages/cctbx/geometry_restraints/manager.py", line 1386, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 56.0397 > 50: distance: 39 - 54: 5.052 distance: 43 - 65: 11.363 distance: 47 - 74: 23.770 distance: 50 - 54: 3.511 distance: 51 - 79: 7.856 distance: 54 - 55: 7.812 distance: 55 - 56: 11.013 distance: 55 - 58: 13.690 distance: 56 - 57: 10.037 distance: 56 - 65: 9.205 distance: 57 - 84: 22.840 distance: 58 - 59: 18.541 distance: 59 - 60: 17.663 distance: 59 - 61: 16.217 distance: 60 - 62: 14.890 distance: 61 - 63: 29.521 distance: 62 - 64: 18.779 distance: 63 - 64: 30.110 distance: 65 - 66: 4.891 distance: 66 - 67: 19.728 distance: 66 - 69: 7.942 distance: 67 - 68: 17.110 distance: 67 - 74: 10.628 distance: 68 - 95: 38.827 distance: 69 - 70: 11.920 distance: 70 - 71: 13.428 distance: 71 - 72: 14.242 distance: 71 - 73: 11.214 distance: 74 - 75: 29.399 distance: 75 - 76: 9.466 distance: 75 - 78: 24.698 distance: 76 - 77: 16.914 distance: 76 - 79: 16.851 distance: 77 - 104: 16.234 distance: 79 - 80: 9.223 distance: 80 - 81: 14.258 distance: 80 - 83: 22.321 distance: 81 - 82: 10.096 distance: 81 - 84: 7.965 distance: 82 - 115: 10.458 distance: 84 - 85: 9.699 distance: 85 - 86: 12.833 distance: 85 - 88: 15.874 distance: 86 - 87: 16.705 distance: 86 - 95: 8.910 distance: 87 - 120: 20.421 distance: 88 - 89: 18.055 distance: 89 - 90: 6.696 distance: 90 - 91: 8.344 distance: 91 - 92: 14.985 distance: 92 - 93: 7.982 distance: 92 - 94: 20.620 distance: 95 - 96: 15.897 distance: 96 - 97: 6.530 distance: 96 - 99: 8.189 distance: 97 - 98: 15.656 distance: 97 - 104: 9.920 distance: 98 - 129: 49.520 distance: 99 - 100: 19.546 distance: 100 - 101: 33.823 distance: 101 - 102: 30.428 distance: 102 - 103: 25.381 distance: 104 - 105: 17.559 distance: 105 - 106: 10.863 distance: 105 - 108: 15.529 distance: 106 - 107: 18.199 distance: 106 - 115: 20.449 distance: 107 - 138: 28.263 distance: 108 - 109: 19.845 distance: 109 - 110: 7.383 distance: 110 - 111: 17.157 distance: 111 - 112: 5.797 distance: 112 - 113: 15.713 distance: 112 - 114: 7.684 distance: 115 - 116: 10.882 distance: 116 - 117: 23.630 distance: 116 - 119: 4.032 distance: 117 - 118: 27.970 distance: 117 - 120: 24.905 distance: 118 - 143: 17.276 distance: 120 - 121: 19.594 distance: 121 - 122: 26.282 distance: 121 - 124: 40.631 distance: 122 - 123: 28.755 distance: 122 - 129: 23.056 distance: 124 - 125: 10.944 distance: 125 - 126: 38.405 distance: 126 - 127: 40.315 distance: 127 - 128: 56.040 distance: 129 - 130: 25.445 distance: 130 - 131: 20.922 distance: 130 - 133: 18.213 distance: 131 - 132: 14.497 distance: 131 - 138: 25.963 distance: 133 - 134: 14.386 distance: 134 - 135: 22.217 distance: 135 - 136: 39.727 distance: 136 - 137: 6.347 distance: 155 - 161: 3.015