Starting phenix.real_space_refine on Sun Aug 4 14:54:48 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ky9_23072/08_2024/7ky9_23072.cif Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ky9_23072/08_2024/7ky9_23072.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=4.05 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ky9_23072/08_2024/7ky9_23072.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ky9_23072/08_2024/7ky9_23072.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ky9_23072/08_2024/7ky9_23072.cif" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ky9_23072/08_2024/7ky9_23072.cif" } resolution = 4.05 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } 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.002 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: 8 Type Number sf(0) Gaussians Al 1 5.89 5 P 2 5.49 5 Mg 2 5.21 5 S 53 5.16 5 C 6646 2.51 5 N 1690 2.21 5 O 1956 1.98 5 F 4 1.80 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A ARG 589": "NH1" <-> "NH2" Residue "A ARG 592": "NH1" <-> "NH2" Residue "A PHE 606": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 676": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 690": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 724": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 746": "NH1" <-> "NH2" Residue "A ARG 768": "NH1" <-> "NH2" Residue "A ASP 780": "OD1" <-> "OD2" Residue "A PHE 815": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 854": "NH1" <-> "NH2" Residue "A ARG 891": "NH1" <-> "NH2" Residue "A ARG 893": "NH1" <-> "NH2" Residue "A ARG 959": "NH1" <-> "NH2" Residue "A ARG 983": "NH1" <-> "NH2" Residue "A ARG 996": "NH1" <-> "NH2" Residue "A ARG 1009": "NH1" <-> "NH2" Residue "A ASP 1012": "OD1" <-> "OD2" Residue "A GLU 1023": "OE1" <-> "OE2" Residue "A PHE 1044": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1085": "NH1" <-> "NH2" Residue "A GLU 1099": "OE1" <-> "OE2" Residue "A GLU 1128": "OE1" <-> "OE2" Residue "A ARG 1130": "NH1" <-> "NH2" Residue "A ARG 1131": "NH1" <-> "NH2" Residue "A ARG 1147": "NH1" <-> "NH2" Residue "A GLU 1194": "OE1" <-> "OE2" Residue "A ARG 1196": "NH1" <-> "NH2" Residue "A ARG 1210": "NH1" <-> "NH2" Residue "A ARG 1214": "NH1" <-> "NH2" Residue "A GLU 1229": "OE1" <-> "OE2" Residue "A PHE 1254": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 1258": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 1268": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1300": "NH1" <-> "NH2" Residue "A ARG 1305": "NH1" <-> "NH2" Residue "A PHE 1328": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 1330": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A TYR 1332": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1351": "NH1" <-> "NH2" Residue "A PHE 1353": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 1373": "OE1" <-> "OE2" Residue "A TYR 1375": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1376": "NH1" <-> "NH2" Residue "A GLU 1407": "OE1" <-> "OE2" Residue "A TYR 1409": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1414": "NH1" <-> "NH2" Residue "A PHE 1431": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1436": "NH1" <-> "NH2" Residue "A PHE 1437": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1443": "NH1" <-> "NH2" Residue "A TYR 1447": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1455": "NH1" <-> "NH2" Residue "B ARG 51": "NH1" <-> "NH2" Residue "B ARG 52": "NH1" <-> "NH2" Residue "B GLU 55": "OE1" <-> "OE2" Residue "B ARG 62": "NH1" <-> "NH2" Residue "B ARG 73": "NH1" <-> "NH2" Residue "B PHE 86": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ASP 120": "OD1" <-> "OD2" Residue "B PHE 130": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 143": "NH1" <-> "NH2" Residue "B ARG 156": "NH1" <-> "NH2" Residue "B PHE 163": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 177": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 199": "NH1" <-> "NH2" Residue "B TYR 280": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 297": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B GLU 299": "OE1" <-> "OE2" Residue "B GLU 309": "OE1" <-> "OE2" Residue "B ARG 328": "NH1" <-> "NH2" Residue "B TYR 340": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 360": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 378": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 395": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ARG 401": "NH1" <-> "NH2" Time to flip residues: 0.04s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/chem_data/mon_lib" Total number of atoms: 10354 Number of models: 1 Model: "" Number of chains: 7 Chain: "A" Number of atoms: 7229 Number of conformers: 1 Conformer: "" Number of residues, atoms: 910, 7229 Classifications: {'peptide': 910} Link IDs: {'PTRANS': 25, 'TRANS': 884} Chain breaks: 2 Chain: "B" Number of atoms: 2959 Number of conformers: 1 Conformer: "" Number of residues, atoms: 365, 2959 Classifications: {'peptide': 365} Link IDs: {'PTRANS': 19, 'TRANS': 345} Chain: "C" Number of atoms: 51 Number of conformers: 1 Conformer: "" Number of residues, atoms: 4, 51 Unusual residues: {'MAN': 2, 'NAG': 2} Classifications: {'undetermined': 4} Link IDs: {None: 3} Unresolved non-hydrogen bonds: 3 Unresolved non-hydrogen angles: 6 Unresolved non-hydrogen dihedrals: 9 Unresolved non-hydrogen chiralities: 3 Chain: "D" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Chain: "E" Number of atoms: 28 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 28 Unusual residues: {'NAG': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Chain: "A" Number of atoms: 34 Number of conformers: 1 Conformer: "" Number of residues, atoms: 4, 34 Unusual residues: {' MG': 2, 'ADP': 1, 'ALF': 1} Classifications: {'undetermined': 4} Link IDs: {None: 3} Chain: "B" Number of atoms: 25 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 25 Unusual residues: {'MAN': 1, 'NAG': 1} Classifications: {'undetermined': 2} Link IDs: {None: 1} Unresolved non-hydrogen bonds: 2 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 6 Unresolved non-hydrogen chiralities: 2 Time building chain proxies: 6.46, per 1000 atoms: 0.62 Number of scatterers: 10354 At special positions: 0 Unit cell: (129.682, 114.814, 130.508, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 8 Type Number sf(0) S 53 16.00 P 2 15.00 Al 1 13.00 Mg 2 11.99 F 4 9.00 O 1956 8.00 N 1690 7.00 C 6646 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=2, symmetry=0 Simple disulfide: pdb=" SG CYS B 110 " - pdb=" SG CYS B 159 " distance=2.03 Simple disulfide: pdb=" SG CYS B 216 " - pdb=" SG CYS B 231 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 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=1, symmetry=0 Links applied ALPHA1-4 " NAG C 2 " - " MAN C 3 " BETA1-4 " NAG C 1 " - " NAG C 2 " " NAG D 1 " - " NAG D 2 " " NAG E 1 " - " NAG E 2 " NAG-ASN " NAG B1001 " - " ASN B 298 " " NAG C 1 " - " ASN B 240 " " NAG D 1 " - " ASN B 256 " " NAG E 1 " - " ASN B 332 " Number of additional bonds: simple=1, symmetry=0 Coordination: Other bonds: Time building additional restraints: 4.01 Conformation dependent library (CDL) restraints added in 1.9 seconds 2534 Ramachandran restraints generated. 1267 Oldfield, 0 Emsley, 1267 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2384 Finding SS restraints... Secondary structure from input PDB file: 52 helices and 7 sheets defined 51.5% alpha, 10.4% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.51 Creating SS restraints... Processing helix chain 'A' and resid 282 through 301 removed outlier: 3.600A pdb=" N SER A 301 " --> pdb=" O GLY A 297 " (cutoff:3.500A) Processing helix chain 'A' and resid 588 through 619 removed outlier: 3.682A pdb=" N ARG A 592 " --> pdb=" O SER A 588 " (cutoff:3.500A) removed outlier: 4.001A pdb=" N VAL A 615 " --> pdb=" O THR A 611 " (cutoff:3.500A) removed outlier: 3.949A pdb=" N GLY A 617 " --> pdb=" O GLY A 613 " (cutoff:3.500A) removed outlier: 4.134A pdb=" N TYR A 619 " --> pdb=" O VAL A 615 " (cutoff:3.500A) Processing helix chain 'A' and resid 640 through 654 Processing helix chain 'A' and resid 661 through 679 removed outlier: 5.005A pdb=" N GLU A 667 " --> pdb=" O TYR A 663 " (cutoff:3.500A) Processing helix chain 'A' and resid 680 through 684 Processing helix chain 'A' and resid 700 through 705 removed outlier: 3.694A pdb=" N GLY A 704 " --> pdb=" O ASP A 701 " (cutoff:3.500A) Processing helix chain 'A' and resid 739 through 749 Processing helix chain 'A' and resid 754 through 779 removed outlier: 3.511A pdb=" N SER A 779 " --> pdb=" O LEU A 775 " (cutoff:3.500A) Processing helix chain 'A' and resid 793 through 799 removed outlier: 3.873A pdb=" N VAL A 797 " --> pdb=" O SER A 793 " (cutoff:3.500A) removed outlier: 3.535A pdb=" N GLU A 798 " --> pdb=" O LYS A 794 " (cutoff:3.500A) Processing helix chain 'A' and resid 805 through 821 Processing helix chain 'A' and resid 843 through 857 Processing helix chain 'A' and resid 918 through 925 removed outlier: 4.200A pdb=" N ARG A 925 " --> pdb=" O VAL A 921 " (cutoff:3.500A) Processing helix chain 'A' and resid 934 through 949 Processing helix chain 'A' and resid 962 through 977 Processing helix chain 'A' and resid 983 through 995 Processing helix chain 'A' and resid 1013 through 1024 removed outlier: 3.796A pdb=" N SER A1017 " --> pdb=" O GLY A1013 " (cutoff:3.500A) Processing helix chain 'A' and resid 1036 through 1046 Processing helix chain 'A' and resid 1063 through 1069 removed outlier: 4.038A pdb=" N GLU A1067 " --> pdb=" O GLU A1063 " (cutoff:3.500A) Processing helix chain 'A' and resid 1071 through 1087 removed outlier: 3.930A pdb=" N VAL A1075 " --> pdb=" O ASP A1071 " (cutoff:3.500A) Processing helix chain 'A' and resid 1094 through 1104 Processing helix chain 'A' and resid 1116 through 1124 removed outlier: 3.782A pdb=" N LEU A1124 " --> pdb=" O LEU A1120 " (cutoff:3.500A) Processing helix chain 'A' and resid 1126 through 1137 Processing helix chain 'A' and resid 1149 through 1163 Processing helix chain 'A' and resid 1177 through 1183 Processing helix chain 'A' and resid 1193 through 1200 removed outlier: 3.792A pdb=" N GLN A1197 " --> pdb=" O GLU A1194 " (cutoff:3.500A) Processing helix chain 'A' and resid 1208 through 1210 No H-bonds generated for 'chain 'A' and resid 1208 through 1210' Processing helix chain 'A' and resid 1211 through 1218 removed outlier: 3.505A pdb=" N LEU A1217 " --> pdb=" O THR A1213 " (cutoff:3.500A) Processing helix chain 'A' and resid 1218 through 1249 Proline residue: A1232 - end of helix removed outlier: 3.712A pdb=" N VAL A1239 " --> pdb=" O PHE A1235 " (cutoff:3.500A) removed outlier: 4.323A pdb=" N ILE A1240 " --> pdb=" O TYR A1236 " (cutoff:3.500A) removed outlier: 5.091A pdb=" N PHE A1241 " --> pdb=" O LYS A1237 " (cutoff:3.500A) removed outlier: 3.870A pdb=" N THR A1242 " --> pdb=" O ASN A1238 " (cutoff:3.500A) removed outlier: 3.812A pdb=" N LEU A1245 " --> pdb=" O PHE A1241 " (cutoff:3.500A) Processing helix chain 'A' and resid 1250 through 1255 removed outlier: 4.622A pdb=" N ASN A1253 " --> pdb=" O ILE A1250 " (cutoff:3.500A) Processing helix chain 'A' and resid 1263 through 1271 removed outlier: 3.775A pdb=" N PHE A1267 " --> pdb=" O THR A1263 " (cutoff:3.500A) removed outlier: 3.945A pdb=" N TYR A1268 " --> pdb=" O TYR A1264 " (cutoff:3.500A) removed outlier: 5.668A pdb=" N ASN A1269 " --> pdb=" O LEU A1265 " (cutoff:3.500A) Processing helix chain 'A' and resid 1274 through 1283 removed outlier: 3.931A pdb=" N ILE A1278 " --> pdb=" O SER A1274 " (cutoff:3.500A) Processing helix chain 'A' and resid 1288 through 1295 Processing helix chain 'A' and resid 1296 through 1298 No H-bonds generated for 'chain 'A' and resid 1296 through 1298' Processing helix chain 'A' and resid 1299 through 1305 Processing helix chain 'A' and resid 1309 through 1334 removed outlier: 3.502A pdb=" N PHE A1313 " --> pdb=" O ASN A1309 " (cutoff:3.500A) Proline residue: A1331 - end of helix Processing helix chain 'A' and resid 1350 through 1374 Processing helix chain 'A' and resid 1380 through 1399 removed outlier: 3.622A pdb=" N PHE A1384 " --> pdb=" O PHE A1380 " (cutoff:3.500A) Processing helix chain 'A' and resid 1402 through 1406 Processing helix chain 'A' and resid 1410 through 1416 Processing helix chain 'A' and resid 1418 through 1447 removed outlier: 4.107A pdb=" N TRP A1422 " --> pdb=" O GLN A1418 " (cutoff:3.500A) removed outlier: 3.807A pdb=" N VAL A1427 " --> pdb=" O ALA A1423 " (cutoff:3.500A) Proline residue: A1435 - end of helix Processing helix chain 'A' and resid 1449 through 1461 removed outlier: 3.733A pdb=" N GLY A1461 " --> pdb=" O MET A1457 " (cutoff:3.500A) Processing helix chain 'B' and resid 56 through 61 Processing helix chain 'B' and resid 71 through 98 Proline residue: B 77 - end of helix Processing helix chain 'B' and resid 186 through 191 Processing helix chain 'B' and resid 204 through 210 Processing helix chain 'B' and resid 232 through 237 Processing helix chain 'B' and resid 267 through 273 Processing helix chain 'B' and resid 286 through 290 removed outlier: 3.879A pdb=" N GLU B 290 " --> pdb=" O PRO B 287 " (cutoff:3.500A) Processing helix chain 'B' and resid 309 through 314 removed outlier: 4.055A pdb=" N TRP B 314 " --> pdb=" O GLU B 310 " (cutoff:3.500A) Processing helix chain 'B' and resid 350 through 354 removed outlier: 3.683A pdb=" N PHE B 353 " --> pdb=" O VAL B 350 " (cutoff:3.500A) Processing helix chain 'B' and resid 372 through 399 Processing helix chain 'B' and resid 405 through 409 Processing sheet with id=AA1, first strand: chain 'A' and resid 1055 through 1057 removed outlier: 6.473A pdb=" N VAL A1113 " --> pdb=" O LEU A1144 " (cutoff:3.500A) removed outlier: 4.457A pdb=" N VAL A1143 " --> pdb=" O LEU A1028 " (cutoff:3.500A) removed outlier: 3.552A pdb=" N VAL A1030 " --> pdb=" O VAL A1143 " (cutoff:3.500A) removed outlier: 3.989A pdb=" N CYS A1145 " --> pdb=" O VAL A1030 " (cutoff:3.500A) removed outlier: 5.622A pdb=" N ILE A 709 " --> pdb=" O TRP A1029 " (cutoff:3.500A) removed outlier: 7.246A pdb=" N LEU A1031 " --> pdb=" O ILE A 709 " (cutoff:3.500A) removed outlier: 6.420A pdb=" N SER A 711 " --> pdb=" O LEU A1031 " (cutoff:3.500A) removed outlier: 6.904A pdb=" N TYR A 708 " --> pdb=" O LEU A1169 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 732 through 734 removed outlier: 5.520A pdb=" N VAL A 721 " --> pdb=" O ARG A1009 " (cutoff:3.500A) removed outlier: 3.732A pdb=" N ARG A1009 " --> pdb=" O VAL A 721 " (cutoff:3.500A) removed outlier: 3.775A pdb=" N LYS A 725 " --> pdb=" O ALA A1005 " (cutoff:3.500A) removed outlier: 7.484A pdb=" N ALA A1005 " --> pdb=" O LYS A 725 " (cutoff:3.500A) removed outlier: 5.757A pdb=" N CYS A 727 " --> pdb=" O GLY A1003 " (cutoff:3.500A) removed outlier: 7.515A pdb=" N GLY A1003 " --> pdb=" O CYS A 727 " (cutoff:3.500A) removed outlier: 6.253A pdb=" N ILE A 729 " --> pdb=" O LEU A1001 " (cutoff:3.500A) removed outlier: 10.236A pdb=" N LEU A1001 " --> pdb=" O ILE A 729 " (cutoff:3.500A) removed outlier: 6.848A pdb=" N LEU A 956 " --> pdb=" O LEU A1001 " (cutoff:3.500A) removed outlier: 4.490A pdb=" N GLY A1003 " --> pdb=" O LEU A 954 " (cutoff:3.500A) removed outlier: 6.878A pdb=" N LEU A 954 " --> pdb=" O GLY A1003 " (cutoff:3.500A) removed outlier: 4.953A pdb=" N ALA A1005 " --> pdb=" O ARG A 952 " (cutoff:3.500A) removed outlier: 6.461A pdb=" N ARG A 952 " --> pdb=" O ALA A1005 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 858 through 860 Processing sheet with id=AA4, first strand: chain 'A' and resid 1188 through 1190 Processing sheet with id=AA5, first strand: chain 'A' and resid 1377 through 1378 Processing sheet with id=AA6, first strand: chain 'B' and resid 126 through 129 removed outlier: 6.897A pdb=" N VAL B 173 " --> pdb=" O ILE B 329 " (cutoff:3.500A) removed outlier: 4.587A pdb=" N ILE B 329 " --> pdb=" O VAL B 173 " (cutoff:3.500A) removed outlier: 6.687A pdb=" N LEU B 175 " --> pdb=" O ILE B 327 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 119 through 120 removed outlier: 4.035A pdb=" N GLN B 341 " --> pdb=" O THR B 248 " (cutoff:3.500A) removed outlier: 4.101A pdb=" N THR B 248 " --> pdb=" O GLN B 341 " (cutoff:3.500A) removed outlier: 6.795A pdb=" N ASP B 343 " --> pdb=" O GLN B 246 " (cutoff:3.500A) removed outlier: 8.570A pdb=" N GLN B 246 " --> pdb=" O ASP B 343 " (cutoff:3.500A) 493 hydrogen bonds defined for protein. 1386 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: 3.36 Time building geometry restraints manager: 4.39 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.23 - 1.34: 3265 1.34 - 1.46: 2385 1.46 - 1.58: 4853 1.58 - 1.69: 3 1.69 - 1.81: 80 Bond restraints: 10586 Sorted by residual: bond pdb=" F1 ALF A1702 " pdb="AL ALF A1702 " ideal model delta sigma weight residual 1.684 1.785 -0.101 2.00e-02 2.50e+03 2.57e+01 bond pdb=" F2 ALF A1702 " pdb="AL ALF A1702 " ideal model delta sigma weight residual 1.685 1.785 -0.100 2.00e-02 2.50e+03 2.51e+01 bond pdb=" F4 ALF A1702 " pdb="AL ALF A1702 " ideal model delta sigma weight residual 1.686 1.785 -0.099 2.00e-02 2.50e+03 2.47e+01 bond pdb=" F3 ALF A1702 " pdb="AL ALF A1702 " ideal model delta sigma weight residual 1.685 1.784 -0.099 2.00e-02 2.50e+03 2.46e+01 bond pdb=" C5 MAN C 3 " pdb=" O5 MAN C 3 " ideal model delta sigma weight residual 1.418 1.461 -0.043 2.00e-02 2.50e+03 4.61e+00 ... (remaining 10581 not shown) Histogram of bond angle deviations from ideal: 89.89 - 107.90: 373 107.90 - 125.91: 13827 125.91 - 143.92: 161 143.92 - 161.93: 0 161.93 - 179.94: 2 Bond angle restraints: 14363 Sorted by residual: angle pdb=" F1 ALF A1702 " pdb="AL ALF A1702 " pdb=" F2 ALF A1702 " ideal model delta sigma weight residual 108.68 179.87 -71.19 3.00e+00 1.11e-01 5.63e+02 angle pdb=" F3 ALF A1702 " pdb="AL ALF A1702 " pdb=" F4 ALF A1702 " ideal model delta sigma weight residual 109.63 179.94 -70.31 3.00e+00 1.11e-01 5.49e+02 angle pdb=" F2 ALF A1702 " pdb="AL ALF A1702 " pdb=" F4 ALF A1702 " ideal model delta sigma weight residual 110.21 90.13 20.08 3.00e+00 1.11e-01 4.48e+01 angle pdb=" F2 ALF A1702 " pdb="AL ALF A1702 " pdb=" F3 ALF A1702 " ideal model delta sigma weight residual 109.59 89.89 19.70 3.00e+00 1.11e-01 4.31e+01 angle pdb=" F1 ALF A1702 " pdb="AL ALF A1702 " pdb=" F3 ALF A1702 " ideal model delta sigma weight residual 109.69 90.00 19.69 3.00e+00 1.11e-01 4.31e+01 ... (remaining 14358 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 28.79: 6163 28.79 - 57.58: 208 57.58 - 86.38: 18 86.38 - 115.17: 5 115.17 - 143.96: 2 Dihedral angle restraints: 6396 sinusoidal: 2679 harmonic: 3717 Sorted by residual: dihedral pdb=" CB CYS B 110 " pdb=" SG CYS B 110 " pdb=" SG CYS B 159 " pdb=" CB CYS B 159 " ideal model delta sinusoidal sigma weight residual -86.00 -157.60 71.60 1 1.00e+01 1.00e-02 6.57e+01 dihedral pdb=" O1B ADP A1701 " pdb=" O3A ADP A1701 " pdb=" PB ADP A1701 " pdb=" PA ADP A1701 " ideal model delta sinusoidal sigma weight residual -60.00 83.96 -143.96 1 2.00e+01 2.50e-03 4.34e+01 dihedral pdb=" O2A ADP A1701 " pdb=" O3A ADP A1701 " pdb=" PA ADP A1701 " pdb=" PB ADP A1701 " ideal model delta sinusoidal sigma weight residual -60.00 -150.44 90.44 1 2.00e+01 2.50e-03 2.42e+01 ... (remaining 6393 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.165: 1622 0.165 - 0.329: 2 0.329 - 0.494: 0 0.494 - 0.659: 0 0.659 - 0.824: 1 Chirality restraints: 1625 Sorted by residual: chirality pdb=" C1 NAG C 1 " pdb=" ND2 ASN B 240 " pdb=" C2 NAG C 1 " pdb=" O5 NAG C 1 " both_signs ideal model delta sigma weight residual False -2.40 -1.58 -0.82 2.00e-01 2.50e+01 1.70e+01 chirality pdb=" C1 MAN C 3 " pdb=" O4 NAG C 2 " pdb=" C2 MAN C 3 " pdb=" O5 MAN C 3 " both_signs ideal model delta sigma weight residual False 2.40 2.36 0.04 2.00e-02 2.50e+03 4.03e+00 chirality pdb=" C1 NAG D 2 " pdb=" O4 NAG D 1 " pdb=" C2 NAG D 2 " pdb=" O5 NAG D 2 " both_signs ideal model delta sigma weight residual False -2.40 -2.37 -0.03 2.00e-02 2.50e+03 2.09e+00 ... (remaining 1622 not shown) Planarity restraints: 1806 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB TRP A 648 " 0.023 2.00e-02 2.50e+03 2.22e-02 1.23e+01 pdb=" CG TRP A 648 " -0.060 2.00e-02 2.50e+03 pdb=" CD1 TRP A 648 " 0.028 2.00e-02 2.50e+03 pdb=" CD2 TRP A 648 " 0.002 2.00e-02 2.50e+03 pdb=" NE1 TRP A 648 " -0.004 2.00e-02 2.50e+03 pdb=" CE2 TRP A 648 " 0.002 2.00e-02 2.50e+03 pdb=" CE3 TRP A 648 " 0.004 2.00e-02 2.50e+03 pdb=" CZ2 TRP A 648 " 0.000 2.00e-02 2.50e+03 pdb=" CZ3 TRP A 648 " 0.004 2.00e-02 2.50e+03 pdb=" CH2 TRP A 648 " 0.001 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB ASN B 240 " -0.002 2.00e-02 2.50e+03 2.93e-02 1.07e+01 pdb=" CG ASN B 240 " 0.029 2.00e-02 2.50e+03 pdb=" OD1 ASN B 240 " -0.033 2.00e-02 2.50e+03 pdb=" ND2 ASN B 240 " 0.037 2.00e-02 2.50e+03 pdb=" C1 NAG C 1 " -0.031 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA ASN B 240 " 0.012 2.00e-02 2.50e+03 2.41e-02 5.80e+00 pdb=" C ASN B 240 " -0.042 2.00e-02 2.50e+03 pdb=" O ASN B 240 " 0.016 2.00e-02 2.50e+03 pdb=" N ASP B 241 " 0.014 2.00e-02 2.50e+03 ... (remaining 1803 not shown) Histogram of nonbonded interaction distances: 2.21 - 2.74: 1142 2.74 - 3.28: 11225 3.28 - 3.82: 17779 3.82 - 4.36: 20162 4.36 - 4.90: 33082 Nonbonded interactions: 83390 Sorted by model distance: nonbonded pdb=" OD1 ASP B 304 " pdb=" N ILE B 305 " model vdw 2.206 3.120 nonbonded pdb=" OE1 GLU A1193 " pdb=" F3 ALF A1702 " model vdw 2.216 2.990 nonbonded pdb=" OG1 THR A 718 " pdb=" O ARG A1009 " model vdw 2.218 3.040 nonbonded pdb=" O ARG B 188 " pdb=" OG SER B 192 " model vdw 2.224 3.040 nonbonded pdb=" CG2 THR A1342 " pdb=" O LEU A1346 " model vdw 2.244 3.460 ... (remaining 83385 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.00 Found NCS groups: ncs_group { reference = chain 'D' selection = chain 'E' } 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: 1.340 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.420 Check model and map are aligned: 0.070 Set scattering table: 0.090 Process input model: 31.370 Find NCS groups from input model: 0.130 Set up NCS constraints: 0.050 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:2.900 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 36.390 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7209 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.101 10586 Z= 0.244 Angle : 1.140 71.188 14363 Z= 0.474 Chirality : 0.049 0.824 1625 Planarity : 0.003 0.042 1802 Dihedral : 14.166 143.959 4006 Min Nonbonded Distance : 2.206 Molprobity Statistics. All-atom Clashscore : 15.79 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.29 % Favored : 91.71 % Rotamer: Outliers : 0.00 % Allowed : 0.63 % Favored : 99.37 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.25 (0.24), residues: 1267 helix: 0.13 (0.22), residues: 568 sheet: -0.22 (0.42), residues: 143 loop : -2.01 (0.28), residues: 556 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.060 0.002 TRP A 648 HIS 0.008 0.001 HIS B 363 PHE 0.027 0.002 PHE A1437 TYR 0.022 0.002 TYR B 79 ARG 0.019 0.001 ARG A1130 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2534 Ramachandran restraints generated. 1267 Oldfield, 0 Emsley, 1267 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2534 Ramachandran restraints generated. 1267 Oldfield, 0 Emsley, 1267 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 211 residues out of total 1113 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 211 time to evaluate : 1.170 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 771 MET cc_start: -0.0507 (tmm) cc_final: -0.0782 (tmm) REVERT: A 784 PHE cc_start: 0.8947 (p90) cc_final: 0.8497 (t80) REVERT: A 1130 ARG cc_start: 0.9604 (ptp90) cc_final: 0.9308 (ptp90) REVERT: A 1200 MET cc_start: 0.9229 (mpp) cc_final: 0.8903 (mpp) REVERT: A 1213 THR cc_start: 0.8914 (m) cc_final: 0.8662 (t) REVERT: A 1421 TYR cc_start: 0.7916 (t80) cc_final: 0.7499 (t80) REVERT: A 1457 MET cc_start: 0.9097 (ppp) cc_final: 0.8592 (ppp) REVERT: B 179 LEU cc_start: 0.8636 (mp) cc_final: 0.8370 (mp) REVERT: B 182 PHE cc_start: 0.8497 (t80) cc_final: 0.7809 (t80) outliers start: 0 outliers final: 0 residues processed: 211 average time/residue: 0.2006 time to fit residues: 61.3460 Evaluate side-chains 171 residues out of total 1113 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 171 time to evaluate : 1.124 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 127 random chunks: chunk 107 optimal weight: 0.0470 chunk 96 optimal weight: 6.9990 chunk 53 optimal weight: 4.9990 chunk 32 optimal weight: 7.9990 chunk 64 optimal weight: 5.9990 chunk 51 optimal weight: 6.9990 chunk 99 optimal weight: 7.9990 chunk 38 optimal weight: 0.8980 chunk 60 optimal weight: 6.9990 chunk 74 optimal weight: 0.9990 chunk 115 optimal weight: 5.9990 overall best weight: 2.5884 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 942 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1041 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1323 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1350 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** B 113 ASN ** B 215 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 412 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7307 moved from start: 0.1230 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.047 10586 Z= 0.234 Angle : 0.648 13.168 14363 Z= 0.325 Chirality : 0.044 0.165 1625 Planarity : 0.004 0.039 1802 Dihedral : 8.453 127.492 1606 Min Nonbonded Distance : 2.388 Molprobity Statistics. All-atom Clashscore : 15.45 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.89 % Favored : 92.11 % Rotamer: Outliers : 0.09 % Allowed : 7.37 % Favored : 92.54 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.22 (0.24), residues: 1267 helix: 0.09 (0.22), residues: 582 sheet: -0.26 (0.42), residues: 142 loop : -1.97 (0.28), residues: 543 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.022 0.002 TRP A 648 HIS 0.004 0.001 HIS A1104 PHE 0.017 0.002 PHE A1330 TYR 0.017 0.002 TYR A 678 ARG 0.006 0.001 ARG B 162 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2534 Ramachandran restraints generated. 1267 Oldfield, 0 Emsley, 1267 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2534 Ramachandran restraints generated. 1267 Oldfield, 0 Emsley, 1267 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 203 residues out of total 1113 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 202 time to evaluate : 1.156 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 784 PHE cc_start: 0.8989 (p90) cc_final: 0.8529 (t80) REVERT: A 791 PHE cc_start: 0.8230 (t80) cc_final: 0.7982 (t80) REVERT: A 1238 ASN cc_start: 0.8551 (t0) cc_final: 0.8213 (t0) REVERT: A 1303 ILE cc_start: 0.9370 (mm) cc_final: 0.9117 (mt) REVERT: A 1457 MET cc_start: 0.9153 (ppp) cc_final: 0.8897 (ppp) REVERT: B 168 ASP cc_start: 0.9519 (p0) cc_final: 0.9308 (t0) REVERT: B 169 MET cc_start: 0.9057 (mmp) cc_final: 0.8738 (tpp) REVERT: B 182 PHE cc_start: 0.8767 (t80) cc_final: 0.8023 (t80) REVERT: B 304 ASP cc_start: 0.7271 (t0) cc_final: 0.7047 (t0) REVERT: B 333 ASP cc_start: 0.8696 (m-30) cc_final: 0.8415 (m-30) outliers start: 1 outliers final: 0 residues processed: 202 average time/residue: 0.1953 time to fit residues: 57.7435 Evaluate side-chains 170 residues out of total 1113 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 170 time to evaluate : 1.196 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 127 random chunks: chunk 64 optimal weight: 7.9990 chunk 35 optimal weight: 5.9990 chunk 95 optimal weight: 6.9990 chunk 78 optimal weight: 2.9990 chunk 31 optimal weight: 7.9990 chunk 115 optimal weight: 1.9990 chunk 124 optimal weight: 6.9990 chunk 102 optimal weight: 5.9990 chunk 114 optimal weight: 3.9990 chunk 39 optimal weight: 20.0000 chunk 92 optimal weight: 9.9990 overall best weight: 4.1990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 942 HIS A1041 ASN ** A1323 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1350 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** B 215 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7412 moved from start: 0.2090 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.084 10586 Z= 0.345 Angle : 0.726 9.389 14363 Z= 0.374 Chirality : 0.047 0.170 1625 Planarity : 0.004 0.046 1802 Dihedral : 8.157 121.701 1606 Min Nonbonded Distance : 2.406 Molprobity Statistics. All-atom Clashscore : 20.12 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.00 % Favored : 91.00 % Rotamer: Outliers : 0.00 % Allowed : 5.75 % Favored : 94.25 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.39 (0.24), residues: 1267 helix: -0.11 (0.21), residues: 585 sheet: -0.48 (0.42), residues: 140 loop : -1.92 (0.28), residues: 542 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.018 0.002 TRP A1308 HIS 0.007 0.002 HIS A1350 PHE 0.025 0.002 PHE A1426 TYR 0.021 0.002 TYR A 678 ARG 0.008 0.001 ARG B 162 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2534 Ramachandran restraints generated. 1267 Oldfield, 0 Emsley, 1267 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2534 Ramachandran restraints generated. 1267 Oldfield, 0 Emsley, 1267 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 193 residues out of total 1113 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 193 time to evaluate : 1.211 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 784 PHE cc_start: 0.9025 (p90) cc_final: 0.8505 (t80) REVERT: A 1199 VAL cc_start: 0.8807 (m) cc_final: 0.8506 (m) REVERT: A 1255 ASP cc_start: 0.7590 (t0) cc_final: 0.7329 (p0) REVERT: A 1300 ARG cc_start: 0.8954 (mpt180) cc_final: 0.8611 (mmt-90) REVERT: A 1303 ILE cc_start: 0.9453 (mm) cc_final: 0.9239 (mt) REVERT: A 1457 MET cc_start: 0.9107 (ppp) cc_final: 0.8899 (ppp) REVERT: B 182 PHE cc_start: 0.8951 (t80) cc_final: 0.8278 (t80) REVERT: B 239 PHE cc_start: 0.8189 (t80) cc_final: 0.7799 (t80) REVERT: B 333 ASP cc_start: 0.8766 (m-30) cc_final: 0.8445 (m-30) REVERT: B 360 TYR cc_start: 0.9250 (t80) cc_final: 0.9043 (t80) outliers start: 0 outliers final: 0 residues processed: 193 average time/residue: 0.1988 time to fit residues: 56.3395 Evaluate side-chains 154 residues out of total 1113 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 154 time to evaluate : 1.164 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 127 random chunks: chunk 114 optimal weight: 0.4980 chunk 86 optimal weight: 2.9990 chunk 59 optimal weight: 0.0570 chunk 12 optimal weight: 8.9990 chunk 55 optimal weight: 5.9990 chunk 77 optimal weight: 2.9990 chunk 115 optimal weight: 2.9990 chunk 122 optimal weight: 1.9990 chunk 60 optimal weight: 5.9990 chunk 109 optimal weight: 0.8980 chunk 33 optimal weight: 9.9990 overall best weight: 1.2902 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A1041 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A1233 GLN ** A1323 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 125 HIS B 347 HIS Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7331 moved from start: 0.2137 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.041 10586 Z= 0.180 Angle : 0.629 10.016 14363 Z= 0.316 Chirality : 0.044 0.156 1625 Planarity : 0.004 0.043 1802 Dihedral : 7.558 118.324 1606 Min Nonbonded Distance : 2.433 Molprobity Statistics. All-atom Clashscore : 13.90 Ramachandran Plot: Outliers : 0.00 % Allowed : 7.89 % Favored : 92.11 % Rotamer: Outliers : 0.00 % Allowed : 4.22 % Favored : 95.78 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -1.22 (0.24), residues: 1267 helix: 0.06 (0.22), residues: 581 sheet: -0.46 (0.42), residues: 142 loop : -1.84 (0.28), residues: 544 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.018 0.001 TRP B 308 HIS 0.004 0.001 HIS B 363 PHE 0.029 0.002 PHE A1330 TYR 0.019 0.001 TYR A 678 ARG 0.003 0.000 ARG A1130 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2534 Ramachandran restraints generated. 1267 Oldfield, 0 Emsley, 1267 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2534 Ramachandran restraints generated. 1267 Oldfield, 0 Emsley, 1267 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 198 residues out of total 1113 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 198 time to evaluate : 1.036 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash REVERT: A 784 PHE cc_start: 0.9016 (p90) cc_final: 0.8492 (t80) REVERT: A 1199 VAL cc_start: 0.8638 (m) cc_final: 0.8323 (m) REVERT: A 1213 THR cc_start: 0.9122 (m) cc_final: 0.8837 (t) REVERT: A 1300 ARG cc_start: 0.8755 (mpt180) cc_final: 0.8415 (mmt-90) REVERT: A 1303 ILE cc_start: 0.9418 (mm) cc_final: 0.9178 (mt) REVERT: A 1457 MET cc_start: 0.9092 (ppp) cc_final: 0.8850 (ppp) REVERT: B 180 GLU cc_start: 0.8358 (tp30) cc_final: 0.8081 (mm-30) REVERT: B 182 PHE cc_start: 0.8807 (t80) cc_final: 0.7917 (t80) REVERT: B 239 PHE cc_start: 0.7796 (t80) cc_final: 0.7516 (t80) outliers start: 0 outliers final: 0 residues processed: 198 average time/residue: 0.2034 time to fit residues: 59.5743 Evaluate side-chains 163 residues out of total 1113 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 163 time to evaluate : 1.129 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 8, in run_program(real_space_refine.Program) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/iotbx/cli_parser.py", line 946, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/phenix/phenix/programs/real_space_refine.py", line 200, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 306, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 702, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 766, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 51, in __init__ gradients_method = gradients_method) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/mmtbx/refinement/real_space/weight.py", line 87, in __init__ rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 368, in refine rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 112, in refine states_collector = self.states_accumulator) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 174, in __init__ exception_handling_params=lbfgs_exception_handling_params) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 131, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 245, in compute_functional_and_gradients compute_gradients=True) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1479, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1381, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5409/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1262, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 68.2376 > 50: distance: 11 - 96: 28.247 distance: 56 - 117: 27.719 distance: 59 - 114: 35.216 distance: 68 - 100: 35.841 distance: 71 - 97: 16.258 distance: 86 - 91: 30.847 distance: 91 - 92: 12.280 distance: 91 - 241: 33.360 distance: 92 - 93: 36.903 distance: 92 - 95: 21.829 distance: 93 - 94: 17.018 distance: 93 - 97: 41.359 distance: 94 - 238: 31.995 distance: 95 - 96: 20.696 distance: 97 - 98: 27.599 distance: 98 - 99: 19.498 distance: 98 - 101: 17.017 distance: 99 - 100: 9.523 distance: 99 - 106: 23.538 distance: 101 - 102: 21.416 distance: 102 - 103: 19.352 distance: 103 - 104: 5.768 distance: 103 - 105: 22.057 distance: 106 - 107: 51.893 distance: 106 - 229: 29.867 distance: 107 - 108: 13.355 distance: 107 - 110: 13.991 distance: 108 - 109: 16.077 distance: 108 - 114: 43.333 distance: 109 - 226: 34.890 distance: 110 - 111: 33.304 distance: 110 - 112: 19.156 distance: 111 - 113: 28.467 distance: 115 - 116: 43.438 distance: 115 - 118: 39.809 distance: 116 - 117: 14.559 distance: 116 - 125: 39.425 distance: 118 - 119: 39.394 distance: 119 - 120: 5.454 distance: 120 - 121: 36.956 distance: 121 - 122: 20.949 distance: 122 - 123: 12.737 distance: 122 - 124: 18.968 distance: 125 - 126: 9.383 distance: 125 - 210: 28.437 distance: 126 - 127: 18.693 distance: 126 - 129: 40.547 distance: 127 - 128: 37.870 distance: 127 - 136: 50.705 distance: 128 - 207: 32.327 distance: 129 - 130: 64.994 distance: 130 - 131: 3.133 distance: 131 - 133: 44.236 distance: 133 - 135: 45.164 distance: 134 - 135: 35.668 distance: 137 - 138: 42.400 distance: 137 - 140: 39.357 distance: 138 - 139: 47.823 distance: 138 - 143: 39.785 distance: 143 - 144: 21.003 distance: 143 - 197: 26.995 distance: 144 - 145: 36.258 distance: 144 - 147: 57.022 distance: 145 - 146: 13.441 distance: 145 - 150: 41.958 distance: 146 - 194: 21.371 distance: 147 - 148: 42.099 distance: 147 - 149: 42.202 distance: 150 - 151: 41.086 distance: 150 - 156: 41.054 distance: 151 - 152: 47.316 distance: 151 - 154: 34.998 distance: 152 - 153: 38.744 distance: 152 - 157: 20.243 distance: 154 - 155: 7.781 distance: 155 - 156: 55.435 distance: 157 - 158: 17.508 distance: 158 - 159: 47.281 distance: 158 - 161: 37.936 distance: 159 - 160: 13.969 distance: 159 - 163: 26.580 distance: 164 - 165: 25.329 distance: 164 - 167: 39.561 distance: 165 - 166: 14.057 distance: 165 - 171: 28.864 distance: 167 - 168: 10.731 distance: 168 - 169: 15.051 distance: 168 - 170: 13.511