Starting phenix.real_space_refine on Wed Jun 26 19:13:41 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hw4_35050/06_2024/8hw4_35050_updated.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hw4_35050/06_2024/8hw4_35050.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.52 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hw4_35050/06_2024/8hw4_35050.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hw4_35050/06_2024/8hw4_35050.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hw4_35050/06_2024/8hw4_35050_updated.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8hw4_35050/06_2024/8hw4_35050_updated.pdb" } resolution = 3.52 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.002 sd= 0.035 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 4 Type Number sf(0) Gaussians S 59 5.16 5 C 6997 2.51 5 N 1791 2.21 5 O 1899 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A ARG 57": "NH1" <-> "NH2" Residue "A PHE 106": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 153": "NH1" <-> "NH2" Residue "A PHE 169": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 196": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 217": "NH1" <-> "NH2" Residue "A TYR 231": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 232": "NH1" <-> "NH2" Residue "A ARG 248": "NH1" <-> "NH2" Residue "A ARG 261": "NH1" <-> "NH2" Residue "A ARG 391": "NH1" <-> "NH2" Residue "A ARG 402": "NH1" <-> "NH2" Residue "A ASP 416": "OD1" <-> "OD2" Residue "A ARG 487": "NH1" <-> "NH2" Residue "A GLU 493": "OE1" <-> "OE2" Residue "A TYR 504": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 518": "NH1" <-> "NH2" Residue "A ARG 526": "NH1" <-> "NH2" Residue "A GLU 610": "OE1" <-> "OE2" Residue "A GLU 707": "OE1" <-> "OE2" Residue "A GLU 734": "OE1" <-> "OE2" Residue "A ARG 756": "NH1" <-> "NH2" Residue "A PHE 772": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1011": "NH1" <-> "NH2" Residue "A PHE 1024": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1042": "NH1" <-> "NH2" Residue "A ARG 1054": "NH1" <-> "NH2" Residue "A ARG 1067": "NH1" <-> "NH2" Residue "A GLU 1081": "OE1" <-> "OE2" Residue "A ARG 1138": "NH1" <-> "NH2" Residue "A ARG 1144": "NH1" <-> "NH2" Residue "A ARG 1162": "NH1" <-> "NH2" Residue "A TYR 1164": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ARG 1166": "NH1" <-> "NH2" Residue "A ARG 1168": "NH1" <-> "NH2" Residue "A ARG 1193": "NH1" <-> "NH2" Residue "A PHE 1201": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 1281": "OE1" <-> "OE2" Residue "A GLU 1341": "OE1" <-> "OE2" Residue "A GLU 1346": "OE1" <-> "OE2" Residue "A PHE 1387": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 1400": "OE1" <-> "OE2" Residue "A ARG 1445": "NH1" <-> "NH2" Residue "A TYR 1490": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.05s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/modules/chem_data/mon_lib" Total number of atoms: 10746 Number of models: 1 Model: "" Number of chains: 1 Chain: "A" Number of atoms: 10746 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1366, 10746 Unusual residues: {'ZWY': 2} Classifications: {'peptide': 1364, 'undetermined': 2} Link IDs: {'PTRANS': 59, 'TRANS': 1304, None: 2} Not linked: pdbres="MET A1520 " pdbres="ZWY A1601 " Not linked: pdbres="ZWY A1601 " pdbres="ZWY A1602 " Chain breaks: 2 Time building chain proxies: 6.01, per 1000 atoms: 0.56 Number of scatterers: 10746 At special positions: 0 Unit cell: (74.9, 109.14, 153.01, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 4 Type Number sf(0) S 59 16.00 O 1899 8.00 N 1791 7.00 C 6997 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.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 Number of additional bonds: simple=1, symmetry=0 Coordination: Other bonds: Time building additional restraints: 4.28 Conformation dependent library (CDL) restraints added in 2.3 seconds 2716 Ramachandran restraints generated. 1358 Oldfield, 0 Emsley, 1358 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2560 Finding SS restraints... running ksdssp... Secondary structure from input PDB file: 68 helices and 5 sheets defined 56.3% alpha, 3.4% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 3.20 Creating SS restraints... Processing helix chain 'A' and resid 37 through 51 Proline residue: A 41 - end of helix Proline residue: A 50 - end of helix Processing helix chain 'A' and resid 53 through 58 removed outlier: 3.527A pdb=" N ARG A 57 " --> pdb=" O LEU A 53 " (cutoff:3.500A) removed outlier: 4.196A pdb=" N HIS A 58 " --> pdb=" O LEU A 54 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 53 through 58' Processing helix chain 'A' and resid 70 through 73 No H-bonds generated for 'chain 'A' and resid 70 through 73' Processing helix chain 'A' and resid 76 through 94 removed outlier: 4.239A pdb=" N ALA A 86 " --> pdb=" O CYS A 82 " (cutoff:3.500A) removed outlier: 3.502A pdb=" N ASP A 87 " --> pdb=" O VAL A 83 " (cutoff:3.500A) removed outlier: 4.033A pdb=" N LEU A 88 " --> pdb=" O SER A 84 " (cutoff:3.500A) removed outlier: 4.557A pdb=" N SER A 91 " --> pdb=" O ASP A 87 " (cutoff:3.500A) removed outlier: 4.139A pdb=" N GLY A 94 " --> pdb=" O TYR A 90 " (cutoff:3.500A) Processing helix chain 'A' and resid 103 through 127 removed outlier: 3.925A pdb=" N THR A 108 " --> pdb=" O VAL A 104 " (cutoff:3.500A) Proline residue: A 109 - end of helix removed outlier: 4.577A pdb=" N MET A 116 " --> pdb=" O VAL A 112 " (cutoff:3.500A) removed outlier: 3.552A pdb=" N LEU A 117 " --> pdb=" O GLY A 113 " (cutoff:3.500A) removed outlier: 3.951A pdb=" N LEU A 121 " --> pdb=" O LEU A 117 " (cutoff:3.500A) removed outlier: 4.190A pdb=" N GLN A 124 " --> pdb=" O THR A 120 " (cutoff:3.500A) Processing helix chain 'A' and resid 136 through 147 Processing helix chain 'A' and resid 149 through 162 removed outlier: 3.632A pdb=" N ILE A 156 " --> pdb=" O PHE A 152 " (cutoff:3.500A) removed outlier: 4.816A pdb=" N LEU A 157 " --> pdb=" O ARG A 153 " (cutoff:3.500A) removed outlier: 3.632A pdb=" N GLU A 162 " --> pdb=" O LEU A 158 " (cutoff:3.500A) Processing helix chain 'A' and resid 168 through 189 removed outlier: 3.502A pdb=" N ALA A 179 " --> pdb=" O TYR A 175 " (cutoff:3.500A) Processing helix chain 'A' and resid 214 through 218 Processing helix chain 'A' and resid 221 through 232 removed outlier: 3.710A pdb=" N LYS A 225 " --> pdb=" O TRP A 222 " (cutoff:3.500A) removed outlier: 3.648A pdb=" N ILE A 228 " --> pdb=" O LYS A 225 " (cutoff:3.500A) removed outlier: 3.840A pdb=" N GLY A 230 " --> pdb=" O ALA A 227 " (cutoff:3.500A) removed outlier: 4.077A pdb=" N ARG A 232 " --> pdb=" O TYR A 229 " (cutoff:3.500A) Processing helix chain 'A' and resid 237 through 239 No H-bonds generated for 'chain 'A' and resid 237 through 239' Processing helix chain 'A' and resid 254 through 262 removed outlier: 4.096A pdb=" N ALA A 259 " --> pdb=" O GLN A 255 " (cutoff:3.500A) Processing helix chain 'A' and resid 299 through 306 Processing helix chain 'A' and resid 315 through 339 removed outlier: 4.499A pdb=" N PHE A 326 " --> pdb=" O ASP A 322 " (cutoff:3.500A) removed outlier: 5.005A pdb=" N ILE A 327 " --> pdb=" O LEU A 323 " (cutoff:3.500A) removed outlier: 4.870A pdb=" N ASN A 328 " --> pdb=" O LEU A 324 " (cutoff:3.500A) Proline residue: A 329 - end of helix Processing helix chain 'A' and resid 346 through 357 removed outlier: 4.506A pdb=" N LEU A 351 " --> pdb=" O TRP A 348 " (cutoff:3.500A) removed outlier: 3.562A pdb=" N LEU A 355 " --> pdb=" O VAL A 352 " (cutoff:3.500A) Processing helix chain 'A' and resid 371 through 374 No H-bonds generated for 'chain 'A' and resid 371 through 374' Processing helix chain 'A' and resid 378 through 394 removed outlier: 3.811A pdb=" N MET A 386 " --> pdb=" O ARG A 382 " (cutoff:3.500A) Processing helix chain 'A' and resid 398 through 403 Processing helix chain 'A' and resid 406 through 412 removed outlier: 4.141A pdb=" N VAL A 410 " --> pdb=" O VAL A 406 " (cutoff:3.500A) removed outlier: 4.059A pdb=" N ASN A 411 " --> pdb=" O GLY A 407 " (cutoff:3.500A) removed outlier: 3.525A pdb=" N LEU A 412 " --> pdb=" O GLU A 408 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 406 through 412' Processing helix chain 'A' and resid 424 through 436 removed outlier: 4.273A pdb=" N ASN A 428 " --> pdb=" O ALA A 424 " (cutoff:3.500A) removed outlier: 4.722A pdb=" N LEU A 429 " --> pdb=" O PRO A 425 " (cutoff:3.500A) removed outlier: 4.150A pdb=" N LEU A 430 " --> pdb=" O PHE A 426 " (cutoff:3.500A) removed outlier: 4.815A pdb=" N TRP A 431 " --> pdb=" O LEU A 427 " (cutoff:3.500A) Proline residue: A 434 - end of helix Processing helix chain 'A' and resid 440 through 448 removed outlier: 3.523A pdb=" N LEU A 444 " --> pdb=" O ALA A 440 " (cutoff:3.500A) Processing helix chain 'A' and resid 450 through 464 Proline residue: A 464 - end of helix Processing helix chain 'A' and resid 468 through 472 Processing helix chain 'A' and resid 482 through 495 Processing helix chain 'A' and resid 498 through 503 Processing helix chain 'A' and resid 507 through 526 removed outlier: 3.510A pdb=" N LEU A 524 " --> pdb=" O GLY A 520 " (cutoff:3.500A) removed outlier: 3.526A pdb=" N LEU A 525 " --> pdb=" O GLU A 521 " (cutoff:3.500A) Processing helix chain 'A' and resid 530 through 532 No H-bonds generated for 'chain 'A' and resid 530 through 532' Processing helix chain 'A' and resid 542 through 557 Processing helix chain 'A' and resid 565 through 582 removed outlier: 4.160A pdb=" N VAL A 570 " --> pdb=" O GLU A 566 " (cutoff:3.500A) removed outlier: 4.021A pdb=" N ASN A 576 " --> pdb=" O VAL A 572 " (cutoff:3.500A) removed outlier: 3.877A pdb=" N ILE A 577 " --> pdb=" O SER A 573 " (cutoff:3.500A) removed outlier: 3.665A pdb=" N LEU A 580 " --> pdb=" O ASN A 576 " (cutoff:3.500A) Proline residue: A 581 - end of helix Processing helix chain 'A' and resid 584 through 606 removed outlier: 4.391A pdb=" N VAL A 597 " --> pdb=" O THR A 593 " (cutoff:3.500A) removed outlier: 3.933A pdb=" N SER A 598 " --> pdb=" O GLN A 594 " (cutoff:3.500A) removed outlier: 3.880A pdb=" N LYS A 600 " --> pdb=" O SER A 596 " (cutoff:3.500A) Processing helix chain 'A' and resid 668 through 675 removed outlier: 3.640A pdb=" N SER A 672 " --> pdb=" O SER A 668 " (cutoff:3.500A) removed outlier: 3.661A pdb=" N ALA A 673 " --> pdb=" O SER A 669 " (cutoff:3.500A) Processing helix chain 'A' and resid 705 through 710 Processing helix chain 'A' and resid 717 through 726 Processing helix chain 'A' and resid 730 through 735 Processing helix chain 'A' and resid 745 through 747 No H-bonds generated for 'chain 'A' and resid 745 through 747' Processing helix chain 'A' and resid 753 through 766 removed outlier: 3.768A pdb=" N GLN A 757 " --> pdb=" O GLY A 753 " (cutoff:3.500A) removed outlier: 3.755A pdb=" N ARG A 758 " --> pdb=" O GLY A 754 " (cutoff:3.500A) Processing helix chain 'A' and resid 777 through 780 removed outlier: 4.439A pdb=" N ALA A 780 " --> pdb=" O PRO A 777 " (cutoff:3.500A) No H-bonds generated for 'chain 'A' and resid 777 through 780' Processing helix chain 'A' and resid 783 through 792 Processing helix chain 'A' and resid 835 through 840 removed outlier: 4.010A pdb=" N GLN A 840 " --> pdb=" O PRO A 836 " (cutoff:3.500A) Processing helix chain 'A' and resid 844 through 849 removed outlier: 3.573A pdb=" N LEU A 849 " --> pdb=" O PHE A 845 " (cutoff:3.500A) Processing helix chain 'A' and resid 956 through 966 removed outlier: 3.664A pdb=" N TRP A 960 " --> pdb=" O LEU A 956 " (cutoff:3.500A) Processing helix chain 'A' and resid 971 through 996 removed outlier: 3.988A pdb=" N LEU A 977 " --> pdb=" O ALA A 973 " (cutoff:3.500A) Processing helix chain 'A' and resid 1006 through 1024 removed outlier: 3.823A pdb=" N SER A1009 " --> pdb=" O ASN A1006 " (cutoff:3.500A) removed outlier: 4.546A pdb=" N LEU A1012 " --> pdb=" O SER A1009 " (cutoff:3.500A) removed outlier: 5.221A pdb=" N ILE A1020 " --> pdb=" O ALA A1017 " (cutoff:3.500A) removed outlier: 4.815A pdb=" N LEU A1021 " --> pdb=" O LEU A1018 " (cutoff:3.500A) Processing helix chain 'A' and resid 1027 through 1052 removed outlier: 4.064A pdb=" N ALA A1032 " --> pdb=" O LEU A1028 " (cutoff:3.500A) removed outlier: 4.310A pdb=" N GLY A1036 " --> pdb=" O ALA A1032 " (cutoff:3.500A) removed outlier: 4.500A pdb=" N GLN A1039 " --> pdb=" O ALA A1035 " (cutoff:3.500A) removed outlier: 3.567A pdb=" N LEU A1044 " --> pdb=" O ALA A1040 " (cutoff:3.500A) removed outlier: 4.102A pdb=" N ASN A1051 " --> pdb=" O ALA A1047 " (cutoff:3.500A) Processing helix chain 'A' and resid 1057 through 1062 Processing helix chain 'A' and resid 1067 through 1070 No H-bonds generated for 'chain 'A' and resid 1067 through 1070' Processing helix chain 'A' and resid 1075 through 1081 removed outlier: 4.078A pdb=" N VAL A1079 " --> pdb=" O ASP A1075 " (cutoff:3.500A) Processing helix chain 'A' and resid 1083 through 1107 removed outlier: 3.617A pdb=" N PHE A1095 " --> pdb=" O LEU A1091 " (cutoff:3.500A) Processing helix chain 'A' and resid 1109 through 1141 removed outlier: 3.714A pdb=" N THR A1112 " --> pdb=" O PRO A1109 " (cutoff:3.500A) removed outlier: 3.556A pdb=" N VAL A1113 " --> pdb=" O LEU A1110 " (cutoff:3.500A) removed outlier: 4.229A pdb=" N LEU A1116 " --> pdb=" O VAL A1113 " (cutoff:3.500A) Proline residue: A1117 - end of helix removed outlier: 5.078A pdb=" N SER A1133 " --> pdb=" O ALA A1130 " (cutoff:3.500A) Processing helix chain 'A' and resid 1145 through 1150 Processing helix chain 'A' and resid 1154 through 1161 removed outlier: 5.360A pdb=" N SER A1159 " --> pdb=" O VAL A1155 " (cutoff:3.500A) removed outlier: 5.585A pdb=" N VAL A1160 " --> pdb=" O THR A1156 " (cutoff:3.500A) Processing helix chain 'A' and resid 1167 through 1204 removed outlier: 4.071A pdb=" N ILE A1173 " --> pdb=" O ASP A1169 " (cutoff:3.500A) removed outlier: 4.224A pdb=" N SER A1174 " --> pdb=" O PHE A1170 " (cutoff:3.500A) removed outlier: 3.831A pdb=" N ASN A1181 " --> pdb=" O LYS A1177 " (cutoff:3.500A) removed outlier: 4.710A pdb=" N SER A1184 " --> pdb=" O ALA A1180 " (cutoff:3.500A) Proline residue: A1187 - end of helix removed outlier: 3.823A pdb=" N VAL A1202 " --> pdb=" O GLY A1198 " (cutoff:3.500A) removed outlier: 3.786A pdb=" N GLY A1203 " --> pdb=" O VAL A1199 " (cutoff:3.500A) Processing helix chain 'A' and resid 1207 through 1215 Processing helix chain 'A' and resid 1223 through 1233 removed outlier: 4.000A pdb=" N LEU A1228 " --> pdb=" O GLY A1224 " (cutoff:3.500A) Processing helix chain 'A' and resid 1242 through 1247 Processing helix chain 'A' and resid 1255 through 1262 Processing helix chain 'A' and resid 1284 through 1286 No H-bonds generated for 'chain 'A' and resid 1284 through 1286' Processing helix chain 'A' and resid 1330 through 1337 removed outlier: 3.998A pdb=" N LEU A1334 " --> pdb=" O SER A1330 " (cutoff:3.500A) removed outlier: 3.751A pdb=" N PHE A1337 " --> pdb=" O THR A1333 " (cutoff:3.500A) Processing helix chain 'A' and resid 1354 through 1356 No H-bonds generated for 'chain 'A' and resid 1354 through 1356' Processing helix chain 'A' and resid 1359 through 1365 removed outlier: 3.821A pdb=" N SER A1364 " --> pdb=" O HIS A1360 " (cutoff:3.500A) removed outlier: 3.702A pdb=" N GLN A1365 " --> pdb=" O ASP A1361 " (cutoff:3.500A) Processing helix chain 'A' and resid 1380 through 1385 Processing helix chain 'A' and resid 1392 through 1401 Processing helix chain 'A' and resid 1405 through 1409 Processing helix chain 'A' and resid 1414 through 1416 No H-bonds generated for 'chain 'A' and resid 1414 through 1416' Processing helix chain 'A' and resid 1428 through 1442 removed outlier: 4.202A pdb=" N ARG A1442 " --> pdb=" O ARG A1438 " (cutoff:3.500A) Processing helix chain 'A' and resid 1458 through 1471 removed outlier: 3.619A pdb=" N GLN A1471 " --> pdb=" O THR A1467 " (cutoff:3.500A) Processing helix chain 'A' and resid 1507 through 1513 Processing helix chain 'A' and resid 1516 through 1519 No H-bonds generated for 'chain 'A' and resid 1516 through 1519' Processing sheet with id= A, first strand: chain 'A' and resid 820 through 823 Processing sheet with id= B, first strand: chain 'A' and resid 1289 through 1292 removed outlier: 3.683A pdb=" N ARG A1292 " --> pdb=" O GLU A1346 " (cutoff:3.500A) removed outlier: 3.705A pdb=" N GLU A1346 " --> pdb=" O ARG A1292 " (cutoff:3.500A) Processing sheet with id= C, first strand: chain 'A' and resid 1319 through 1323 removed outlier: 7.143A pdb=" N ARG A1492 " --> pdb=" O GLY A1320 " (cutoff:3.500A) removed outlier: 8.382A pdb=" N VAL A1322 " --> pdb=" O ARG A1492 " (cutoff:3.500A) removed outlier: 6.527A pdb=" N LEU A1494 " --> pdb=" O VAL A1322 " (cutoff:3.500A) removed outlier: 3.624A pdb=" N GLU A1503 " --> pdb=" O VAL A1495 " (cutoff:3.500A) removed outlier: 6.609A pdb=" N ASP A1497 " --> pdb=" O VAL A1501 " (cutoff:3.500A) removed outlier: 5.266A pdb=" N VAL A1501 " --> pdb=" O ASP A1497 " (cutoff:3.500A) Processing sheet with id= D, first strand: chain 'A' and resid 1446 through 1448 removed outlier: 5.995A pdb=" N THR A1476 " --> pdb=" O LEU A1447 " (cutoff:3.500A) No H-bonds generated for sheet with id= D Processing sheet with id= E, first strand: chain 'A' and resid 625 through 630 removed outlier: 7.085A pdb=" N VAL A 651 " --> pdb=" O ALA A 626 " (cutoff:3.500A) removed outlier: 4.441A pdb=" N THR A 628 " --> pdb=" O ILE A 649 " (cutoff:3.500A) removed outlier: 6.813A pdb=" N ILE A 649 " --> pdb=" O THR A 628 " (cutoff:3.500A) 419 hydrogen bonds defined for protein. 1059 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.69 Time building geometry restraints manager: 4.72 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: 2066 1.33 - 1.45: 2630 1.45 - 1.57: 6205 1.57 - 1.69: 7 1.69 - 1.81: 93 Bond restraints: 11001 Sorted by residual: bond pdb=" C15 ZWY A1602 " pdb=" C16 ZWY A1602 " ideal model delta sigma weight residual 1.549 1.277 0.272 2.00e-02 2.50e+03 1.85e+02 bond pdb=" C04 ZWY A1602 " pdb=" C05 ZWY A1602 " ideal model delta sigma weight residual 1.541 1.686 -0.145 2.00e-02 2.50e+03 5.28e+01 bond pdb=" C05 ZWY A1602 " pdb=" C06 ZWY A1602 " ideal model delta sigma weight residual 1.545 1.418 0.127 2.00e-02 2.50e+03 4.04e+01 bond pdb=" O20 ZWY A1602 " pdb=" S21 ZWY A1602 " ideal model delta sigma weight residual 1.628 1.740 -0.112 2.00e-02 2.50e+03 3.16e+01 bond pdb=" O20 ZWY A1601 " pdb=" S21 ZWY A1601 " ideal model delta sigma weight residual 1.628 1.740 -0.112 2.00e-02 2.50e+03 3.12e+01 ... (remaining 10996 not shown) Histogram of bond angle deviations from ideal: 99.10 - 106.09: 269 106.09 - 113.07: 6110 113.07 - 120.06: 3962 120.06 - 127.05: 4484 127.05 - 134.03: 153 Bond angle restraints: 14978 Sorted by residual: angle pdb=" N LEU A 365 " pdb=" CA LEU A 365 " pdb=" C LEU A 365 " ideal model delta sigma weight residual 114.56 106.76 7.80 1.27e+00 6.20e-01 3.77e+01 angle pdb=" N ILE A 366 " pdb=" CA ILE A 366 " pdb=" C ILE A 366 " ideal model delta sigma weight residual 112.29 106.54 5.75 9.40e-01 1.13e+00 3.74e+01 angle pdb=" N VAL A 75 " pdb=" CA VAL A 75 " pdb=" C VAL A 75 " ideal model delta sigma weight residual 112.96 107.67 5.29 1.00e+00 1.00e+00 2.79e+01 angle pdb=" C06 ZWY A1602 " pdb=" C07 ZWY A1602 " pdb=" C08 ZWY A1602 " ideal model delta sigma weight residual 120.63 104.90 15.73 3.00e+00 1.11e-01 2.75e+01 angle pdb=" C13 ZWY A1602 " pdb=" C14 ZWY A1602 " pdb=" C15 ZWY A1602 " ideal model delta sigma weight residual 122.95 108.50 14.45 3.00e+00 1.11e-01 2.32e+01 ... (remaining 14973 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 21.79: 6057 21.79 - 43.58: 494 43.58 - 65.36: 39 65.36 - 87.15: 17 87.15 - 108.94: 5 Dihedral angle restraints: 6612 sinusoidal: 2650 harmonic: 3962 Sorted by residual: dihedral pdb=" CD ARG A1245 " pdb=" NE ARG A1245 " pdb=" CZ ARG A1245 " pdb=" NH1 ARG A1245 " ideal model delta sinusoidal sigma weight residual 0.00 -83.40 83.40 1 1.00e+01 1.00e-02 8.50e+01 dihedral pdb=" CA PHE A 196 " pdb=" C PHE A 196 " pdb=" N PHE A 197 " pdb=" CA PHE A 197 " ideal model delta harmonic sigma weight residual -180.00 -150.03 -29.97 0 5.00e+00 4.00e-02 3.59e+01 dihedral pdb=" CA GLU A 164 " pdb=" C GLU A 164 " pdb=" N ILE A 165 " pdb=" CA ILE A 165 " ideal model delta harmonic sigma weight residual 180.00 -155.96 -24.04 0 5.00e+00 4.00e-02 2.31e+01 ... (remaining 6609 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.188: 1736 0.188 - 0.377: 8 0.377 - 0.565: 2 0.565 - 0.753: 2 0.753 - 0.941: 1 Chirality restraints: 1749 Sorted by residual: chirality pdb=" C15 ZWY A1602 " pdb=" C05 ZWY A1602 " pdb=" C14 ZWY A1602 " pdb=" C16 ZWY A1602 " both_signs ideal model delta sigma weight residual False 2.89 1.95 0.94 2.00e-01 2.50e+01 2.22e+01 chirality pdb=" C15 ZWY A1601 " pdb=" C05 ZWY A1601 " pdb=" C14 ZWY A1601 " pdb=" C16 ZWY A1601 " both_signs ideal model delta sigma weight residual False 2.89 2.17 0.72 2.00e-01 2.50e+01 1.28e+01 chirality pdb=" C05 ZWY A1602 " pdb=" C04 ZWY A1602 " pdb=" C06 ZWY A1602 " pdb=" C15 ZWY A1602 " both_signs ideal model delta sigma weight residual False -2.46 -1.86 -0.60 2.00e-01 2.50e+01 8.92e+00 ... (remaining 1746 not shown) Planarity restraints: 1864 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C12 ZWY A1601 " 0.194 2.00e-02 2.50e+03 1.73e-01 3.76e+02 pdb=" C13 ZWY A1601 " -0.196 2.00e-02 2.50e+03 pdb=" C14 ZWY A1601 " -0.192 2.00e-02 2.50e+03 pdb=" C15 ZWY A1601 " 0.001 2.00e-02 2.50e+03 pdb=" C19 ZWY A1601 " 0.194 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CD ARG A1245 " -1.092 9.50e-02 1.11e+02 4.89e-01 1.44e+02 pdb=" NE ARG A1245 " 0.066 2.00e-02 2.50e+03 pdb=" CZ ARG A1245 " 0.012 2.00e-02 2.50e+03 pdb=" NH1 ARG A1245 " -0.012 2.00e-02 2.50e+03 pdb=" NH2 ARG A1245 " -0.018 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C12 ZWY A1602 " -0.052 2.00e-02 2.50e+03 9.40e-02 1.10e+02 pdb=" C13 ZWY A1602 " 0.127 2.00e-02 2.50e+03 pdb=" C14 ZWY A1602 " -0.148 2.00e-02 2.50e+03 pdb=" C15 ZWY A1602 " 0.057 2.00e-02 2.50e+03 pdb=" C19 ZWY A1602 " 0.016 2.00e-02 2.50e+03 ... (remaining 1861 not shown) Histogram of nonbonded interaction distances: 2.14 - 2.69: 374 2.69 - 3.25: 11903 3.25 - 3.80: 17273 3.80 - 4.35: 22268 4.35 - 4.90: 35690 Nonbonded interactions: 87508 Sorted by model distance: nonbonded pdb=" OG1 THR A1479 " pdb=" OH TYR A1490 " model vdw 2.142 2.440 nonbonded pdb=" OE2 GLU A 617 " pdb=" ND1 HIS A 686 " model vdw 2.191 2.520 nonbonded pdb=" O GLY A1282 " pdb=" NH1 ARG A1286 " model vdw 2.236 2.520 nonbonded pdb=" O ILE A 396 " pdb=" OG1 THR A 397 " model vdw 2.242 2.440 nonbonded pdb=" OD2 ASP A 247 " pdb=" OH TYR A1186 " model vdw 2.247 2.440 ... (remaining 87503 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== 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.710 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: 3.810 Check model and map are aligned: 0.080 Set scattering table: 0.090 Process input model: 32.280 Find NCS groups from input model: 0.320 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.680 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 40.040 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7822 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.272 11001 Z= 0.501 Angle : 0.891 15.729 14978 Z= 0.461 Chirality : 0.063 0.941 1749 Planarity : 0.013 0.489 1864 Dihedral : 15.423 108.938 4052 Min Nonbonded Distance : 2.142 Molprobity Statistics. All-atom Clashscore : 14.67 Ramachandran Plot: Outliers : 0.00 % Allowed : 8.17 % Favored : 91.83 % Rotamer: Outliers : 0.09 % Allowed : 0.34 % Favored : 99.57 % 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: -2.38 (0.22), residues: 1358 helix: -1.57 (0.17), residues: 826 sheet: -0.76 (0.61), residues: 65 loop : -1.36 (0.29), residues: 467 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.017 0.002 TRP A1194 HIS 0.007 0.001 HIS A 369 PHE 0.030 0.002 PHE A1151 TYR 0.019 0.002 TYR A 207 ARG 0.003 0.000 ARG A1278 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2716 Ramachandran restraints generated. 1358 Oldfield, 0 Emsley, 1358 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2716 Ramachandran restraints generated. 1358 Oldfield, 0 Emsley, 1358 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 158 residues out of total 1167 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 157 time to evaluate : 1.088 Fit side-chains REVERT: A 361 MET cc_start: 0.7775 (tpt) cc_final: 0.7573 (tpt) REVERT: A 499 LYS cc_start: 0.8863 (ttpt) cc_final: 0.8624 (ttmt) REVERT: A 1115 ILE cc_start: 0.8961 (tp) cc_final: 0.8737 (tp) REVERT: A 1245 ARG cc_start: 0.7268 (OUTLIER) cc_final: 0.6802 (tpp-160) REVERT: A 1272 TRP cc_start: 0.7668 (m-10) cc_final: 0.7427 (m-10) REVERT: A 1432 GLN cc_start: 0.8535 (mp10) cc_final: 0.8328 (mp10) outliers start: 1 outliers final: 0 residues processed: 158 average time/residue: 0.1976 time to fit residues: 46.7257 Evaluate side-chains 134 residues out of total 1167 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 133 time to evaluate : 1.150 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 1245 ARG Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 136 random chunks: chunk 114 optimal weight: 9.9990 chunk 103 optimal weight: 0.6980 chunk 57 optimal weight: 20.0000 chunk 35 optimal weight: 0.8980 chunk 69 optimal weight: 10.0000 chunk 55 optimal weight: 8.9990 chunk 106 optimal weight: 0.0270 chunk 41 optimal weight: 0.1980 chunk 64 optimal weight: 0.4980 chunk 79 optimal weight: 4.9990 chunk 123 optimal weight: 0.6980 overall best weight: 0.4238 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 436 GLN ** A 576 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1365 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7711 moved from start: 0.1180 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 11001 Z= 0.193 Angle : 0.687 10.479 14978 Z= 0.341 Chirality : 0.044 0.228 1749 Planarity : 0.005 0.052 1864 Dihedral : 8.273 91.976 1593 Min Nonbonded Distance : 2.131 Molprobity Statistics. All-atom Clashscore : 11.67 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.19 % Favored : 93.81 % Rotamer: Outliers : 1.03 % Allowed : 9.08 % Favored : 89.89 % 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.98 (0.22), residues: 1358 helix: -1.19 (0.18), residues: 846 sheet: -0.41 (0.63), residues: 65 loop : -1.37 (0.29), residues: 447 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.001 TRP A1194 HIS 0.006 0.001 HIS A1403 PHE 0.019 0.001 PHE A1059 TYR 0.025 0.001 TYR A 720 ARG 0.002 0.000 ARG A1245 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2716 Ramachandran restraints generated. 1358 Oldfield, 0 Emsley, 1358 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2716 Ramachandran restraints generated. 1358 Oldfield, 0 Emsley, 1358 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 163 residues out of total 1167 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 12 poor density : 151 time to evaluate : 1.105 Fit side-chains revert: symmetry clash REVERT: A 580 LEU cc_start: 0.8675 (tp) cc_final: 0.8321 (tp) REVERT: A 678 MET cc_start: 0.7857 (tpt) cc_final: 0.7575 (mmt) REVERT: A 697 GLN cc_start: 0.8421 (mt0) cc_final: 0.7920 (pm20) REVERT: A 1031 MET cc_start: 0.7116 (mmt) cc_final: 0.6860 (mmt) REVERT: A 1052 LYS cc_start: 0.8270 (mmpt) cc_final: 0.7680 (mmtp) REVERT: A 1245 ARG cc_start: 0.7185 (OUTLIER) cc_final: 0.6816 (tpp-160) REVERT: A 1263 TYR cc_start: 0.8052 (m-10) cc_final: 0.7308 (m-10) REVERT: A 1337 PHE cc_start: 0.8629 (m-80) cc_final: 0.8426 (m-80) outliers start: 12 outliers final: 8 residues processed: 155 average time/residue: 0.2059 time to fit residues: 47.7516 Evaluate side-chains 136 residues out of total 1167 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 9 poor density : 127 time to evaluate : 1.130 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 415 VAL Chi-restraints excluded: chain A residue 453 LEU Chi-restraints excluded: chain A residue 471 VAL Chi-restraints excluded: chain A residue 629 ILE Chi-restraints excluded: chain A residue 704 THR Chi-restraints excluded: chain A residue 710 LEU Chi-restraints excluded: chain A residue 998 ASP Chi-restraints excluded: chain A residue 1114 VAL Chi-restraints excluded: chain A residue 1245 ARG Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 136 random chunks: chunk 68 optimal weight: 0.0980 chunk 38 optimal weight: 0.9990 chunk 102 optimal weight: 0.6980 chunk 84 optimal weight: 7.9990 chunk 34 optimal weight: 9.9990 chunk 123 optimal weight: 0.7980 chunk 133 optimal weight: 10.0000 chunk 110 optimal weight: 2.9990 chunk 122 optimal weight: 4.9990 chunk 42 optimal weight: 9.9990 chunk 99 optimal weight: 0.9990 overall best weight: 0.7184 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 436 GLN ** A 576 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1365 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7719 moved from start: 0.1505 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.034 11001 Z= 0.206 Angle : 0.650 9.582 14978 Z= 0.324 Chirality : 0.043 0.176 1749 Planarity : 0.005 0.056 1864 Dihedral : 7.774 88.287 1593 Min Nonbonded Distance : 2.136 Molprobity Statistics. All-atom Clashscore : 11.11 Ramachandran Plot: Outliers : 0.00 % Allowed : 5.45 % Favored : 94.55 % Rotamer: Outliers : 2.40 % Allowed : 12.17 % Favored : 85.43 % 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.71 (0.23), residues: 1358 helix: -0.93 (0.18), residues: 837 sheet: -0.06 (0.64), residues: 65 loop : -1.41 (0.29), residues: 456 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP A1194 HIS 0.004 0.001 HIS A1403 PHE 0.021 0.001 PHE A1059 TYR 0.011 0.001 TYR A 207 ARG 0.005 0.000 ARG A1483 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2716 Ramachandran restraints generated. 1358 Oldfield, 0 Emsley, 1358 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2716 Ramachandran restraints generated. 1358 Oldfield, 0 Emsley, 1358 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 176 residues out of total 1167 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 28 poor density : 148 time to evaluate : 1.073 Fit side-chains revert: symmetry clash revert: symmetry clash REVERT: A 215 LEU cc_start: 0.8752 (tt) cc_final: 0.8530 (mt) REVERT: A 580 LEU cc_start: 0.8716 (tp) cc_final: 0.8241 (tp) REVERT: A 678 MET cc_start: 0.7743 (OUTLIER) cc_final: 0.7518 (mmt) REVERT: A 741 GLN cc_start: 0.8249 (mm-40) cc_final: 0.7976 (mm-40) REVERT: A 1052 LYS cc_start: 0.8212 (mmpt) cc_final: 0.7751 (mmtp) outliers start: 28 outliers final: 17 residues processed: 159 average time/residue: 0.2048 time to fit residues: 48.4202 Evaluate side-chains 146 residues out of total 1167 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 18 poor density : 128 time to evaluate : 0.999 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 165 ILE Chi-restraints excluded: chain A residue 415 VAL Chi-restraints excluded: chain A residue 453 LEU Chi-restraints excluded: chain A residue 471 VAL Chi-restraints excluded: chain A residue 629 ILE Chi-restraints excluded: chain A residue 675 LEU Chi-restraints excluded: chain A residue 678 MET Chi-restraints excluded: chain A residue 704 THR Chi-restraints excluded: chain A residue 705 LEU Chi-restraints excluded: chain A residue 710 LEU Chi-restraints excluded: chain A residue 976 LEU Chi-restraints excluded: chain A residue 998 ASP Chi-restraints excluded: chain A residue 1113 VAL Chi-restraints excluded: chain A residue 1114 VAL Chi-restraints excluded: chain A residue 1200 GLU Chi-restraints excluded: chain A residue 1215 VAL Chi-restraints excluded: chain A residue 1380 LEU Chi-restraints excluded: chain A residue 1482 HIS Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 136 random chunks: chunk 122 optimal weight: 9.9990 chunk 93 optimal weight: 0.9990 chunk 64 optimal weight: 4.9990 chunk 13 optimal weight: 2.9990 chunk 59 optimal weight: 7.9990 chunk 83 optimal weight: 0.5980 chunk 124 optimal weight: 9.9990 chunk 131 optimal weight: 7.9990 chunk 117 optimal weight: 5.9990 chunk 35 optimal weight: 0.9990 chunk 109 optimal weight: 2.9990 overall best weight: 1.7188 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 576 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** A 810 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** ** A1365 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7797 moved from start: 0.1573 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.055 11001 Z= 0.331 Angle : 0.688 8.216 14978 Z= 0.344 Chirality : 0.045 0.201 1749 Planarity : 0.005 0.064 1864 Dihedral : 7.499 84.730 1589 Min Nonbonded Distance : 2.147 Molprobity Statistics. All-atom Clashscore : 13.28 Ramachandran Plot: Outliers : 0.00 % Allowed : 6.11 % Favored : 93.89 % Rotamer: Outliers : 3.51 % Allowed : 15.25 % Favored : 81.23 % 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.76 (0.22), residues: 1358 helix: -0.98 (0.18), residues: 837 sheet: 0.05 (0.65), residues: 65 loop : -1.43 (0.29), residues: 456 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.002 TRP A1194 HIS 0.004 0.001 HIS A1150 PHE 0.020 0.002 PHE A1059 TYR 0.015 0.001 TYR A 207 ARG 0.007 0.000 ARG A 756 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2716 Ramachandran restraints generated. 1358 Oldfield, 0 Emsley, 1358 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2716 Ramachandran restraints generated. 1358 Oldfield, 0 Emsley, 1358 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 179 residues out of total 1167 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 41 poor density : 138 time to evaluate : 1.087 Fit side-chains revert: symmetry clash revert: symmetry clash REVERT: A 580 LEU cc_start: 0.8774 (tp) cc_final: 0.8234 (tp) REVERT: A 678 MET cc_start: 0.7772 (OUTLIER) cc_final: 0.7480 (mmt) REVERT: A 741 GLN cc_start: 0.8205 (mm-40) cc_final: 0.7898 (mm-40) REVERT: A 757 GLN cc_start: 0.8868 (mp10) cc_final: 0.8267 (mp10) REVERT: A 1052 LYS cc_start: 0.8289 (mmpt) cc_final: 0.7699 (mmtp) outliers start: 41 outliers final: 27 residues processed: 158 average time/residue: 0.1793 time to fit residues: 43.4923 Evaluate side-chains 158 residues out of total 1167 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 28 poor density : 130 time to evaluate : 1.067 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 80 LEU Chi-restraints excluded: chain A residue 165 ILE Chi-restraints excluded: chain A residue 362 MET Chi-restraints excluded: chain A residue 415 VAL Chi-restraints excluded: chain A residue 453 LEU Chi-restraints excluded: chain A residue 459 MET Chi-restraints excluded: chain A residue 462 LEU Chi-restraints excluded: chain A residue 471 VAL Chi-restraints excluded: chain A residue 531 LEU Chi-restraints excluded: chain A residue 629 ILE Chi-restraints excluded: chain A residue 675 LEU Chi-restraints excluded: chain A residue 678 MET Chi-restraints excluded: chain A residue 704 THR Chi-restraints excluded: chain A residue 705 LEU Chi-restraints excluded: chain A residue 710 LEU Chi-restraints excluded: chain A residue 724 LEU Chi-restraints excluded: chain A residue 808 VAL Chi-restraints excluded: chain A residue 976 LEU Chi-restraints excluded: chain A residue 1086 VAL Chi-restraints excluded: chain A residue 1113 VAL Chi-restraints excluded: chain A residue 1114 VAL Chi-restraints excluded: chain A residue 1120 VAL Chi-restraints excluded: chain A residue 1139 LEU Chi-restraints excluded: chain A residue 1200 GLU Chi-restraints excluded: chain A residue 1215 VAL Chi-restraints excluded: chain A residue 1380 LEU Chi-restraints excluded: chain A residue 1482 HIS Chi-restraints excluded: chain A residue 1501 VAL Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/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-5370/modules/cctbx_project/iotbx/cli_parser.py", line 946, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/modules/phenix/phenix/programs/real_space_refine.py", line 200, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/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-5370/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 700, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 764, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/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-5370/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-5370/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-5370/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-5370/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-5370/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-5370/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/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-5370/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-5370/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1477, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1379, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5370/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1260, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 72.1301 > 50: distance: 21 - 25: 33.360 distance: 25 - 31: 55.841 distance: 26 - 29: 4.220 distance: 27 - 28: 55.734 distance: 27 - 32: 40.287 distance: 29 - 30: 41.088 distance: 30 - 31: 69.256 distance: 32 - 33: 29.561 distance: 33 - 34: 57.858 distance: 33 - 36: 58.817 distance: 34 - 35: 56.164 distance: 34 - 41: 57.062 distance: 35 - 69: 66.827 distance: 36 - 37: 38.599 distance: 38 - 39: 40.595 distance: 38 - 40: 39.451 distance: 41 - 42: 56.475 distance: 42 - 43: 40.319 distance: 42 - 45: 39.252 distance: 43 - 47: 40.664 distance: 44 - 77: 67.476 distance: 45 - 46: 70.004 distance: 47 - 48: 55.827 distance: 48 - 49: 56.889 distance: 48 - 51: 56.719 distance: 49 - 50: 55.499 distance: 51 - 52: 14.591 distance: 52 - 53: 58.204 distance: 52 - 54: 60.905 distance: 53 - 55: 55.080 distance: 54 - 56: 40.293 distance: 58 - 59: 36.576 distance: 59 - 60: 66.891 distance: 59 - 62: 5.736 distance: 60 - 61: 38.946 distance: 60 - 69: 38.747 distance: 62 - 63: 39.543 distance: 63 - 64: 62.948 distance: 63 - 65: 40.716 distance: 64 - 66: 11.447 distance: 65 - 67: 70.502 distance: 66 - 68: 39.271 distance: 69 - 70: 39.447 distance: 70 - 71: 40.738 distance: 70 - 73: 39.574 distance: 71 - 72: 55.819 distance: 71 - 77: 40.255 distance: 73 - 74: 40.207 distance: 74 - 75: 38.921 distance: 78 - 79: 56.974 distance: 78 - 81: 56.015 distance: 79 - 80: 69.129 distance: 84 - 85: 41.000 distance: 85 - 88: 36.355 distance: 86 - 87: 40.003 distance: 86 - 91: 6.420 distance: 88 - 89: 5.910 distance: 88 - 90: 62.271 distance: 91 - 92: 32.463 distance: 91 - 97: 49.514 distance: 92 - 93: 55.736 distance: 92 - 95: 56.947 distance: 93 - 94: 43.359 distance: 95 - 96: 50.022 distance: 96 - 97: 64.041 distance: 99 - 100: 38.814 distance: 99 - 102: 57.363 distance: 100 - 101: 39.770 distance: 100 - 104: 34.348