Starting phenix.real_space_refine (version: 1.21rc1) on Tue Oct 10 13:42:10 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h3m_34469/10_2023/8h3m_34469.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h3m_34469/10_2023/8h3m_34469.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.48 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h3m_34469/10_2023/8h3m_34469.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h3m_34469/10_2023/8h3m_34469.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h3m_34469/10_2023/8h3m_34469.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8h3m_34469/10_2023/8h3m_34469.pdb" } resolution = 2.48 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.001 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 88 5.16 5 C 13420 2.51 5 N 3458 2.21 5 O 3970 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A PHE 43": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 168": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 192": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 220": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A ASP 228": "OD1" <-> "OD2" Residue "A GLU 309": "OE1" <-> "OE2" Residue "A GLU 324": "OE1" <-> "OE2" Residue "A PHE 338": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 377": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A GLU 702": "OE1" <-> "OE2" Residue "A PHE 1103": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "A PHE 1121": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B ASP 53": "OD1" <-> "OD2" Residue "B PHE 65": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B GLU 281": "OE1" <-> "OE2" Residue "B PHE 329": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B TYR 612": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "B PHE 759": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 32": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 220": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C GLU 224": "OE1" <-> "OE2" Residue "C TYR 265": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 342": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 392": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C TYR 453": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C TYR 508": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 543": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 562": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 565": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C PHE 759": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "H PHE 100": "CD1" <-> "CD2" "CE1" <-> "CE2" Time to flip residues: 0.06s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-5097/modules/chem_data/mon_lib" Total number of atoms: 20936 Number of models: 1 Model: "" Number of chains: 8 Chain: "A" Number of atoms: 7039 Number of conformers: 1 Conformer: "" Number of residues, atoms: 910, 7039 Classifications: {'peptide': 910} Incomplete info: {'truncation_to_alanine': 23} Link IDs: {'PTRANS': 47, 'TRANS': 862} Chain breaks: 10 Unresolved non-hydrogen bonds: 69 Unresolved non-hydrogen angles: 82 Unresolved non-hydrogen dihedrals: 56 Unresolved non-hydrogen chiralities: 4 Planarities with less than four sites: {'GLU:plan': 4, 'ASN:plan1': 1, 'ASP:plan': 4} Unresolved non-hydrogen planarities: 31 Chain: "D" Number of atoms: 255 Number of conformers: 1 Conformer: "" Number of residues, atoms: 30, 255 Classifications: {'peptide': 30} Link IDs: {'PTRANS': 1, 'TRANS': 28} Chain breaks: 3 Chain: "B" Number of atoms: 5883 Number of conformers: 1 Conformer: "" Number of residues, atoms: 763, 5883 Classifications: {'peptide': 763} Incomplete info: {'truncation_to_alanine': 21} Link IDs: {'PTRANS': 40, 'TRANS': 722} Chain breaks: 13 Unresolved non-hydrogen bonds: 63 Unresolved non-hydrogen angles: 76 Unresolved non-hydrogen dihedrals: 50 Unresolved non-hydrogen chiralities: 4 Planarities with less than four sites: {'GLN:plan1': 1, 'GLU:plan': 3, 'ASN:plan1': 1, 'ASP:plan': 4} Unresolved non-hydrogen planarities: 31 Chain: "C" Number of atoms: 7213 Number of conformers: 1 Conformer: "" Number of residues, atoms: 930, 7213 Classifications: {'peptide': 930} Incomplete info: {'truncation_to_alanine': 27} Link IDs: {'PTRANS': 48, 'TRANS': 881} Chain breaks: 13 Unresolved non-hydrogen bonds: 90 Unresolved non-hydrogen angles: 111 Unresolved non-hydrogen dihedrals: 72 Unresolved non-hydrogen chiralities: 5 Planarities with less than four sites: {'GLN:plan1': 1, 'TYR:plan': 1, 'ASN:plan1': 1, 'ASP:plan': 6, 'GLU:plan': 4, 'ARG:plan': 1} Unresolved non-hydrogen planarities: 53 Chain: "H" Number of atoms: 238 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 238 Classifications: {'peptide': 29} Link IDs: {'PTRANS': 1, 'TRANS': 27} Chain breaks: 3 Chain: "A" Number of atoms: 126 Number of conformers: 1 Conformer: "" Number of residues, atoms: 9, 126 Unusual residues: {'NAG': 9} Classifications: {'undetermined': 9} Link IDs: {None: 8} Unresolved non-hydrogen bonds: 9 Unresolved non-hydrogen angles: 18 Unresolved non-hydrogen chiralities: 9 Chain: "B" Number of atoms: 84 Number of conformers: 1 Conformer: "" Number of residues, atoms: 6, 84 Unusual residues: {'NAG': 6} Classifications: {'undetermined': 6} Link IDs: {None: 5} Unresolved non-hydrogen bonds: 6 Unresolved non-hydrogen angles: 12 Unresolved non-hydrogen chiralities: 6 Chain: "C" Number of atoms: 98 Number of conformers: 1 Conformer: "" Number of residues, atoms: 7, 98 Unusual residues: {'NAG': 7} Classifications: {'undetermined': 7} Link IDs: {None: 6} Unresolved non-hydrogen bonds: 7 Unresolved non-hydrogen angles: 14 Unresolved non-hydrogen chiralities: 7 Time building chain proxies: 11.54, per 1000 atoms: 0.55 Number of scatterers: 20936 At special positions: 0 Unit cell: (131.6, 119.568, 168.448, 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 88 16.00 O 3970 8.00 N 3458 7.00 C 13420 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=30, symmetry=0 Simple disulfide: pdb=" SG CYS A 291 " - pdb=" SG CYS A 301 " distance=2.03 Simple disulfide: pdb=" SG CYS A 336 " - pdb=" SG CYS A 361 " distance=2.03 Simple disulfide: pdb=" SG CYS A 379 " - pdb=" SG CYS A 432 " distance=2.03 Simple disulfide: pdb=" SG CYS A 391 " - pdb=" SG CYS A 525 " distance=2.04 Simple disulfide: pdb=" SG CYS A 538 " - pdb=" SG CYS A 590 " distance=2.04 Simple disulfide: pdb=" SG CYS A 617 " - pdb=" SG CYS A 649 " distance=2.03 Simple disulfide: pdb=" SG CYS A 662 " - pdb=" SG CYS A 671 " distance=2.03 Simple disulfide: pdb=" SG CYS A 738 " - pdb=" SG CYS A 760 " distance=2.03 Simple disulfide: pdb=" SG CYS A 743 " - pdb=" SG CYS A 749 " distance=2.03 Simple disulfide: pdb=" SG CYS A1032 " - pdb=" SG CYS A1043 " distance=2.03 Simple disulfide: pdb=" SG CYS A1082 " - pdb=" SG CYS A1126 " distance=2.03 Simple disulfide: pdb=" SG CYS B 291 " - pdb=" SG CYS B 301 " distance=2.03 Simple disulfide: pdb=" SG CYS B 538 " - pdb=" SG CYS B 590 " distance=2.03 Simple disulfide: pdb=" SG CYS B 617 " - pdb=" SG CYS B 649 " distance=2.02 Simple disulfide: pdb=" SG CYS B 662 " - pdb=" SG CYS B 671 " distance=2.03 Simple disulfide: pdb=" SG CYS B 738 " - pdb=" SG CYS B 760 " distance=2.03 Simple disulfide: pdb=" SG CYS B 743 " - pdb=" SG CYS B 749 " distance=2.03 Simple disulfide: pdb=" SG CYS B1032 " - pdb=" SG CYS B1043 " distance=2.03 Simple disulfide: pdb=" SG CYS B1082 " - pdb=" SG CYS B1126 " distance=2.03 Simple disulfide: pdb=" SG CYS C 291 " - pdb=" SG CYS C 301 " distance=2.03 Simple disulfide: pdb=" SG CYS C 336 " - pdb=" SG CYS C 361 " distance=2.03 Simple disulfide: pdb=" SG CYS C 379 " - pdb=" SG CYS C 432 " distance=2.03 Simple disulfide: pdb=" SG CYS C 391 " - pdb=" SG CYS C 525 " distance=2.03 Simple disulfide: pdb=" SG CYS C 538 " - pdb=" SG CYS C 590 " distance=2.03 Simple disulfide: pdb=" SG CYS C 617 " - pdb=" SG CYS C 649 " distance=2.03 Simple disulfide: pdb=" SG CYS C 662 " - pdb=" SG CYS C 671 " distance=2.03 Simple disulfide: pdb=" SG CYS C 738 " - pdb=" SG CYS C 760 " distance=2.03 Simple disulfide: pdb=" SG CYS C 743 " - pdb=" SG CYS C 749 " distance=2.03 Simple disulfide: pdb=" SG CYS C1032 " - pdb=" SG CYS C1043 " distance=2.03 Simple disulfide: pdb=" SG CYS C1082 " - pdb=" SG CYS C1126 " 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=0, symmetry=0 Links applied NAG-ASN " NAG A1301 " - " ASN A1074 " " NAG A1302 " - " ASN A1134 " " NAG A1303 " - " ASN A1098 " " NAG A1304 " - " ASN A 616 " " NAG A1305 " - " ASN A 709 " " NAG A1306 " - " ASN A 801 " " NAG A1307 " - " ASN A 282 " " NAG A1308 " - " ASN A 343 " " NAG A1309 " - " ASN A 717 " " NAG B1301 " - " ASN B1074 " " NAG B1302 " - " ASN B1134 " " NAG B1303 " - " ASN B1098 " " NAG B1304 " - " ASN B 709 " " NAG B1305 " - " ASN B 801 " " NAG B1306 " - " ASN B 717 " " NAG C1301 " - " ASN C1074 " " NAG C1302 " - " ASN C1134 " " NAG C1303 " - " ASN C1098 " " NAG C1304 " - " ASN C 616 " " NAG C1305 " - " ASN C 709 " " NAG C1306 " - " ASN C 801 " " NAG C1307 " - " ASN C 717 " Time building additional restraints: 8.83 Conformation dependent library (CDL) restraints added in 3.0 seconds 5136 Ramachandran restraints generated. 2568 Oldfield, 0 Emsley, 2568 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 4998 Finding SS restraints... Secondary structure from input PDB file: 60 helices and 41 sheets defined 28.0% alpha, 27.2% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 2.54 Creating SS restraints... Processing helix chain 'A' and resid 294 through 304 Processing helix chain 'A' and resid 338 through 343 Processing helix chain 'A' and resid 364 through 368 Processing helix chain 'A' and resid 383 through 387 removed outlier: 3.747A pdb=" N LEU A 387 " --> pdb=" O PRO A 384 " (cutoff:3.500A) Processing helix chain 'A' and resid 502 through 505 Processing helix chain 'A' and resid 616 through 620 Processing helix chain 'A' and resid 737 through 743 Processing helix chain 'A' and resid 746 through 754 Processing helix chain 'A' and resid 758 through 783 removed outlier: 3.517A pdb=" N ASP A 775 " --> pdb=" O ALA A 771 " (cutoff:3.500A) Processing helix chain 'A' and resid 816 through 826 Processing helix chain 'A' and resid 849 through 856 removed outlier: 3.528A pdb=" N LYS A 856 " --> pdb=" O ALA A 852 " (cutoff:3.500A) Processing helix chain 'A' and resid 866 through 885 Processing helix chain 'A' and resid 886 through 890 Processing helix chain 'A' and resid 897 through 908 Processing helix chain 'A' and resid 912 through 919 removed outlier: 4.067A pdb=" N LEU A 916 " --> pdb=" O THR A 912 " (cutoff:3.500A) Processing helix chain 'A' and resid 919 through 940 removed outlier: 4.104A pdb=" N GLN A 935 " --> pdb=" O ILE A 931 " (cutoff:3.500A) removed outlier: 4.495A pdb=" N ASP A 936 " --> pdb=" O GLY A 932 " (cutoff:3.500A) Processing helix chain 'A' and resid 945 through 965 Processing helix chain 'A' and resid 966 through 968 No H-bonds generated for 'chain 'A' and resid 966 through 968' Processing helix chain 'A' and resid 976 through 982 Processing helix chain 'A' and resid 985 through 1033 removed outlier: 3.933A pdb=" N VAL A 991 " --> pdb=" O PRO A 987 " (cutoff:3.500A) Processing helix chain 'D' and resid 28 through 32 removed outlier: 3.509A pdb=" N TYR D 32 " --> pdb=" O PHE D 29 " (cutoff:3.500A) Processing helix chain 'B' and resid 294 through 304 Processing helix chain 'B' and resid 737 through 744 Processing helix chain 'B' and resid 746 through 754 Processing helix chain 'B' and resid 755 through 757 No H-bonds generated for 'chain 'B' and resid 755 through 757' Processing helix chain 'B' and resid 758 through 783 removed outlier: 3.513A pdb=" N ASP B 775 " --> pdb=" O ALA B 771 " (cutoff:3.500A) Processing helix chain 'B' and resid 816 through 826 removed outlier: 3.528A pdb=" N VAL B 826 " --> pdb=" O LEU B 822 " (cutoff:3.500A) Processing helix chain 'B' and resid 849 through 856 Processing helix chain 'B' and resid 866 through 885 Processing helix chain 'B' and resid 886 through 890 Processing helix chain 'B' and resid 897 through 909 Processing helix chain 'B' and resid 912 through 919 removed outlier: 4.076A pdb=" N LEU B 916 " --> pdb=" O THR B 912 " (cutoff:3.500A) Processing helix chain 'B' and resid 919 through 939 Processing helix chain 'B' and resid 945 through 965 Processing helix chain 'B' and resid 966 through 968 No H-bonds generated for 'chain 'B' and resid 966 through 968' Processing helix chain 'B' and resid 976 through 984 Processing helix chain 'B' and resid 985 through 1033 removed outlier: 4.067A pdb=" N VAL B 991 " --> pdb=" O PRO B 987 " (cutoff:3.500A) Processing helix chain 'C' and resid 294 through 304 Processing helix chain 'C' and resid 337 through 343 removed outlier: 3.608A pdb=" N VAL C 341 " --> pdb=" O PRO C 337 " (cutoff:3.500A) Processing helix chain 'C' and resid 349 through 353 Processing helix chain 'C' and resid 364 through 368 removed outlier: 3.687A pdb=" N VAL C 367 " --> pdb=" O ASP C 364 " (cutoff:3.500A) Processing helix chain 'C' and resid 383 through 389 removed outlier: 5.019A pdb=" N ASN C 388 " --> pdb=" O THR C 385 " (cutoff:3.500A) Processing helix chain 'C' and resid 403 through 410 Processing helix chain 'C' and resid 416 through 422 Processing helix chain 'C' and resid 438 through 443 removed outlier: 3.763A pdb=" N ASP C 442 " --> pdb=" O SER C 438 " (cutoff:3.500A) Processing helix chain 'C' and resid 616 through 620 removed outlier: 3.680A pdb=" N VAL C 620 " --> pdb=" O CYS C 617 " (cutoff:3.500A) Processing helix chain 'C' and resid 737 through 743 Processing helix chain 'C' and resid 746 through 754 removed outlier: 3.901A pdb=" N LEU C 754 " --> pdb=" O SER C 750 " (cutoff:3.500A) Processing helix chain 'C' and resid 757 through 783 removed outlier: 3.982A pdb=" N THR C 761 " --> pdb=" O GLY C 757 " (cutoff:3.500A) Processing helix chain 'C' and resid 816 through 826 Processing helix chain 'C' and resid 849 through 855 Processing helix chain 'C' and resid 866 through 885 Processing helix chain 'C' and resid 886 through 891 removed outlier: 3.953A pdb=" N ALA C 890 " --> pdb=" O TRP C 886 " (cutoff:3.500A) Processing helix chain 'C' and resid 897 through 909 Processing helix chain 'C' and resid 912 through 919 removed outlier: 4.148A pdb=" N LEU C 916 " --> pdb=" O THR C 912 " (cutoff:3.500A) Processing helix chain 'C' and resid 919 through 941 removed outlier: 3.736A pdb=" N GLN C 935 " --> pdb=" O ILE C 931 " (cutoff:3.500A) removed outlier: 3.799A pdb=" N ASP C 936 " --> pdb=" O GLY C 932 " (cutoff:3.500A) Processing helix chain 'C' and resid 945 through 965 Processing helix chain 'C' and resid 966 through 968 No H-bonds generated for 'chain 'C' and resid 966 through 968' Processing helix chain 'C' and resid 976 through 982 Processing helix chain 'C' and resid 985 through 1033 removed outlier: 3.752A pdb=" N VAL C 991 " --> pdb=" O PRO C 987 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 28 through 30 removed outlier: 8.127A pdb=" N ASN A 61 " --> pdb=" O TYR A 269 " (cutoff:3.500A) removed outlier: 6.346A pdb=" N TYR A 269 " --> pdb=" O ASN A 61 " (cutoff:3.500A) removed outlier: 6.273A pdb=" N PHE A 201 " --> pdb=" O ASP A 228 " (cutoff:3.500A) removed outlier: 4.260A pdb=" N ASP A 228 " --> pdb=" O PHE A 201 " (cutoff:3.500A) removed outlier: 6.614A pdb=" N ILE A 203 " --> pdb=" O LEU A 226 " (cutoff:3.500A) removed outlier: 6.800A pdb=" N HIS A 207 " --> pdb=" O ALA A 222 " (cutoff:3.500A) removed outlier: 8.287A pdb=" N ALA A 222 " --> pdb=" O HIS A 207 " (cutoff:3.500A) removed outlier: 7.143A pdb=" N VAL A 36 " --> pdb=" O LEU A 223 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 48 through 55 removed outlier: 3.842A pdb=" N ASP A 287 " --> pdb=" O LYS A 278 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 84 through 85 removed outlier: 3.772A pdb=" N GLY A 103 " --> pdb=" O LEU A 241 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'A' and resid 311 through 319 removed outlier: 7.193A pdb=" N VAL A 595 " --> pdb=" O THR A 315 " (cutoff:3.500A) removed outlier: 4.426A pdb=" N ASN A 317 " --> pdb=" O GLY A 593 " (cutoff:3.500A) removed outlier: 6.097A pdb=" N GLY A 593 " --> pdb=" O ASN A 317 " (cutoff:3.500A) removed outlier: 4.366A pdb=" N GLY A 594 " --> pdb=" O GLN A 613 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'A' and resid 325 through 328 removed outlier: 5.497A pdb=" N ASP A 574 " --> pdb=" O ILE A 587 " (cutoff:3.500A) removed outlier: 6.454A pdb=" N GLY A 566 " --> pdb=" O ASP A 574 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'A' and resid 354 through 358 removed outlier: 4.080A pdb=" N VAL A 395 " --> pdb=" O ILE A 358 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'A' and resid 391 through 392 Processing sheet with id=AA8, first strand: chain 'A' and resid 654 through 655 removed outlier: 5.957A pdb=" N GLU A 654 " --> pdb=" O ALA A 694 " (cutoff:3.500A) removed outlier: 8.419A pdb=" N THR A 696 " --> pdb=" O GLU A 654 " (cutoff:3.500A) removed outlier: 6.280A pdb=" N ALA A 672 " --> pdb=" O PRO A 665 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'A' and resid 702 through 704 Processing sheet with id=AB1, first strand: chain 'A' and resid 711 through 728 removed outlier: 6.664A pdb=" N SER A 711 " --> pdb=" O THR A1076 " (cutoff:3.500A) removed outlier: 6.953A pdb=" N THR A1076 " --> pdb=" O SER A 711 " (cutoff:3.500A) removed outlier: 6.615A pdb=" N ALA A 713 " --> pdb=" O ASN A1074 " (cutoff:3.500A) removed outlier: 6.777A pdb=" N ASN A1074 " --> pdb=" O ALA A 713 " (cutoff:3.500A) removed outlier: 6.195A pdb=" N GLU A1072 " --> pdb=" O PRO A 715 " (cutoff:3.500A) removed outlier: 6.971A pdb=" N GLY A1059 " --> pdb=" O SER A1055 " (cutoff:3.500A) removed outlier: 5.137A pdb=" N SER A1055 " --> pdb=" O GLY A1059 " (cutoff:3.500A) removed outlier: 6.661A pdb=" N VAL A1061 " --> pdb=" O PRO A1053 " (cutoff:3.500A) removed outlier: 6.488A pdb=" N LEU A1063 " --> pdb=" O SER A1051 " (cutoff:3.500A) removed outlier: 4.465A pdb=" N SER A1051 " --> pdb=" O LEU A1063 " (cutoff:3.500A) removed outlier: 6.923A pdb=" N VAL A1065 " --> pdb=" O LEU A1049 " (cutoff:3.500A) Processing sheet with id=AB2, first strand: chain 'A' and resid 711 through 728 removed outlier: 6.664A pdb=" N SER A 711 " --> pdb=" O THR A1076 " (cutoff:3.500A) removed outlier: 6.953A pdb=" N THR A1076 " --> pdb=" O SER A 711 " (cutoff:3.500A) removed outlier: 6.615A pdb=" N ALA A 713 " --> pdb=" O ASN A1074 " (cutoff:3.500A) removed outlier: 6.777A pdb=" N ASN A1074 " --> pdb=" O ALA A 713 " (cutoff:3.500A) removed outlier: 6.195A pdb=" N GLU A1072 " --> pdb=" O PRO A 715 " (cutoff:3.500A) removed outlier: 4.011A pdb=" N THR A1076 " --> pdb=" O SER A1097 " (cutoff:3.500A) removed outlier: 3.888A pdb=" N ALA A1078 " --> pdb=" O PHE A1095 " (cutoff:3.500A) removed outlier: 4.598A pdb=" N PHE A1095 " --> pdb=" O ALA A1078 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'A' and resid 733 through 736 removed outlier: 4.413A pdb=" N LYS A 733 " --> pdb=" O LEU A 861 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'A' and resid 788 through 790 removed outlier: 5.649A pdb=" N ILE A 788 " --> pdb=" O ASN C 703 " (cutoff:3.500A) No H-bonds generated for sheet with id=AB4 Processing sheet with id=AB5, first strand: chain 'A' and resid 1120 through 1125 removed outlier: 4.445A pdb=" N ALA A1087 " --> pdb=" O SER A1123 " (cutoff:3.500A) Processing sheet with id=AB6, first strand: chain 'B' and resid 28 through 30 removed outlier: 3.552A pdb=" N THR B 63 " --> pdb=" O VAL B 267 " (cutoff:3.500A) Processing sheet with id=AB7, first strand: chain 'B' and resid 47 through 51 removed outlier: 3.947A pdb=" N ASP B 287 " --> pdb=" O LYS B 278 " (cutoff:3.500A) Processing sheet with id=AB8, first strand: chain 'B' and resid 54 through 55 Processing sheet with id=AB9, first strand: chain 'B' and resid 170 through 171 removed outlier: 3.563A pdb=" N GLY B 107 " --> pdb=" O ARG B 237 " (cutoff:3.500A) Processing sheet with id=AC1, first strand: chain 'B' and resid 191 through 195 Processing sheet with id=AC2, first strand: chain 'B' and resid 311 through 318 removed outlier: 7.337A pdb=" N VAL B 595 " --> pdb=" O THR B 315 " (cutoff:3.500A) removed outlier: 4.641A pdb=" N ASN B 317 " --> pdb=" O GLY B 593 " (cutoff:3.500A) removed outlier: 6.373A pdb=" N GLY B 593 " --> pdb=" O ASN B 317 " (cutoff:3.500A) removed outlier: 4.041A pdb=" N GLY B 594 " --> pdb=" O GLN B 613 " (cutoff:3.500A) Processing sheet with id=AC3, first strand: chain 'B' and resid 538 through 543 Processing sheet with id=AC4, first strand: chain 'B' and resid 566 through 567 Processing sheet with id=AC5, first strand: chain 'B' and resid 654 through 655 removed outlier: 6.207A pdb=" N ALA B 672 " --> pdb=" O PRO B 665 " (cutoff:3.500A) Processing sheet with id=AC6, first strand: chain 'B' and resid 702 through 704 Processing sheet with id=AC7, first strand: chain 'B' and resid 711 through 728 removed outlier: 6.620A pdb=" N SER B 711 " --> pdb=" O THR B1076 " (cutoff:3.500A) removed outlier: 6.900A pdb=" N THR B1076 " --> pdb=" O SER B 711 " (cutoff:3.500A) removed outlier: 6.655A pdb=" N ALA B 713 " --> pdb=" O ASN B1074 " (cutoff:3.500A) removed outlier: 6.773A pdb=" N ASN B1074 " --> pdb=" O ALA B 713 " (cutoff:3.500A) removed outlier: 6.072A pdb=" N GLU B1072 " --> pdb=" O PRO B 715 " (cutoff:3.500A) removed outlier: 6.024A pdb=" N TYR B1067 " --> pdb=" O HIS B1048 " (cutoff:3.500A) removed outlier: 5.799A pdb=" N HIS B1048 " --> pdb=" O TYR B1067 " (cutoff:3.500A) Processing sheet with id=AC8, first strand: chain 'B' and resid 711 through 728 removed outlier: 6.620A pdb=" N SER B 711 " --> pdb=" O THR B1076 " (cutoff:3.500A) removed outlier: 6.900A pdb=" N THR B1076 " --> pdb=" O SER B 711 " (cutoff:3.500A) removed outlier: 6.655A pdb=" N ALA B 713 " --> pdb=" O ASN B1074 " (cutoff:3.500A) removed outlier: 6.773A pdb=" N ASN B1074 " --> pdb=" O ALA B 713 " (cutoff:3.500A) removed outlier: 6.072A pdb=" N GLU B1072 " --> pdb=" O PRO B 715 " (cutoff:3.500A) removed outlier: 3.696A pdb=" N THR B1076 " --> pdb=" O SER B1097 " (cutoff:3.500A) removed outlier: 3.830A pdb=" N ALA B1078 " --> pdb=" O PHE B1095 " (cutoff:3.500A) removed outlier: 4.641A pdb=" N PHE B1095 " --> pdb=" O ALA B1078 " (cutoff:3.500A) Processing sheet with id=AC9, first strand: chain 'B' and resid 733 through 736 removed outlier: 4.454A pdb=" N LYS B 733 " --> pdb=" O LEU B 861 " (cutoff:3.500A) Processing sheet with id=AD1, first strand: chain 'B' and resid 1120 through 1125 removed outlier: 4.660A pdb=" N ALA B1087 " --> pdb=" O SER B1123 " (cutoff:3.500A) Processing sheet with id=AD2, first strand: chain 'C' and resid 28 through 30 removed outlier: 8.154A pdb=" N ASN C 61 " --> pdb=" O TYR C 269 " (cutoff:3.500A) removed outlier: 6.423A pdb=" N TYR C 269 " --> pdb=" O ASN C 61 " (cutoff:3.500A) removed outlier: 6.555A pdb=" N PHE C 201 " --> pdb=" O ASP C 228 " (cutoff:3.500A) removed outlier: 4.388A pdb=" N ASP C 228 " --> pdb=" O PHE C 201 " (cutoff:3.500A) removed outlier: 6.655A pdb=" N ILE C 203 " --> pdb=" O LEU C 226 " (cutoff:3.500A) Processing sheet with id=AD3, first strand: chain 'C' and resid 47 through 55 removed outlier: 3.928A pdb=" N ASP C 287 " --> pdb=" O LYS C 278 " (cutoff:3.500A) Processing sheet with id=AD4, first strand: chain 'C' and resid 84 through 85 Processing sheet with id=AD5, first strand: chain 'C' and resid 311 through 319 removed outlier: 7.033A pdb=" N VAL C 595 " --> pdb=" O THR C 315 " (cutoff:3.500A) removed outlier: 4.515A pdb=" N ASN C 317 " --> pdb=" O GLY C 593 " (cutoff:3.500A) removed outlier: 6.289A pdb=" N GLY C 593 " --> pdb=" O ASN C 317 " (cutoff:3.500A) removed outlier: 4.205A pdb=" N GLY C 594 " --> pdb=" O GLN C 613 " (cutoff:3.500A) Processing sheet with id=AD6, first strand: chain 'C' and resid 324 through 325 removed outlier: 6.568A pdb=" N GLU C 324 " --> pdb=" O ASN C 540 " (cutoff:3.500A) removed outlier: 8.183A pdb=" N ASN C 542 " --> pdb=" O GLU C 324 " (cutoff:3.500A) removed outlier: 5.580A pdb=" N ASP C 574 " --> pdb=" O ILE C 587 " (cutoff:3.500A) Processing sheet with id=AD7, first strand: chain 'C' and resid 354 through 357 Processing sheet with id=AD8, first strand: chain 'C' and resid 361 through 362 removed outlier: 6.361A pdb=" N CYS C 361 " --> pdb=" O CYS C 525 " (cutoff:3.500A) Processing sheet with id=AD9, first strand: chain 'C' and resid 452 through 453 Processing sheet with id=AE1, first strand: chain 'C' and resid 654 through 655 removed outlier: 5.918A pdb=" N GLU C 654 " --> pdb=" O ALA C 694 " (cutoff:3.500A) removed outlier: 8.398A pdb=" N THR C 696 " --> pdb=" O GLU C 654 " (cutoff:3.500A) removed outlier: 6.469A pdb=" N ILE C 670 " --> pdb=" O ILE C 666 " (cutoff:3.500A) Processing sheet with id=AE2, first strand: chain 'C' and resid 711 through 728 removed outlier: 6.605A pdb=" N SER C 711 " --> pdb=" O THR C1076 " (cutoff:3.500A) removed outlier: 6.876A pdb=" N THR C1076 " --> pdb=" O SER C 711 " (cutoff:3.500A) removed outlier: 6.646A pdb=" N ALA C 713 " --> pdb=" O ASN C1074 " (cutoff:3.500A) removed outlier: 6.782A pdb=" N ASN C1074 " --> pdb=" O ALA C 713 " (cutoff:3.500A) removed outlier: 6.104A pdb=" N GLU C1072 " --> pdb=" O PRO C 715 " (cutoff:3.500A) removed outlier: 5.899A pdb=" N TYR C1067 " --> pdb=" O HIS C1048 " (cutoff:3.500A) removed outlier: 5.699A pdb=" N HIS C1048 " --> pdb=" O TYR C1067 " (cutoff:3.500A) Processing sheet with id=AE3, first strand: chain 'C' and resid 733 through 736 removed outlier: 4.462A pdb=" N LYS C 733 " --> pdb=" O LEU C 861 " (cutoff:3.500A) Processing sheet with id=AE4, first strand: chain 'C' and resid 1120 through 1125 removed outlier: 4.449A pdb=" N ALA C1087 " --> pdb=" O SER C1123 " (cutoff:3.500A) Processing sheet with id=AE5, first strand: chain 'C' and resid 1094 through 1097 961 hydrogen bonds defined for protein. 2592 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: 8.91 Time building geometry restraints manager: 9.38 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.35: 6624 1.35 - 1.48: 5598 1.48 - 1.60: 9056 1.60 - 1.73: 0 1.73 - 1.86: 113 Bond restraints: 21391 Sorted by residual: bond pdb=" N ASN A1108 " pdb=" CA ASN A1108 " ideal model delta sigma weight residual 1.456 1.492 -0.035 1.22e-02 6.72e+03 8.36e+00 bond pdb=" CG PRO C 426 " pdb=" CD PRO C 426 " ideal model delta sigma weight residual 1.503 1.412 0.091 3.40e-02 8.65e+02 7.15e+00 bond pdb=" N GLN A1106 " pdb=" CA GLN A1106 " ideal model delta sigma weight residual 1.454 1.487 -0.033 1.34e-02 5.57e+03 6.15e+00 bond pdb=" N ARG A1107 " pdb=" CA ARG A1107 " ideal model delta sigma weight residual 1.457 1.489 -0.032 1.33e-02 5.65e+03 5.75e+00 bond pdb=" CA ASN A 121 " pdb=" CB ASN A 121 " ideal model delta sigma weight residual 1.525 1.556 -0.031 1.29e-02 6.01e+03 5.74e+00 ... (remaining 21386 not shown) Histogram of bond angle deviations from ideal: 93.11 - 101.35: 28 101.35 - 109.59: 2303 109.59 - 117.82: 13791 117.82 - 126.06: 12716 126.06 - 134.30: 257 Bond angle restraints: 29095 Sorted by residual: angle pdb=" N PRO C 426 " pdb=" CD PRO C 426 " pdb=" CG PRO C 426 " ideal model delta sigma weight residual 103.20 93.11 10.09 1.50e+00 4.44e-01 4.52e+01 angle pdb=" CA PRO A 521 " pdb=" N PRO A 521 " pdb=" CD PRO A 521 " ideal model delta sigma weight residual 112.00 102.79 9.21 1.40e+00 5.10e-01 4.33e+01 angle pdb=" CA PRO C 426 " pdb=" CB PRO C 426 " pdb=" CG PRO C 426 " ideal model delta sigma weight residual 104.50 94.79 9.71 1.90e+00 2.77e-01 2.61e+01 angle pdb=" C GLU B 619 " pdb=" CA GLU B 619 " pdb=" CB GLU B 619 " ideal model delta sigma weight residual 109.29 115.76 -6.47 1.51e+00 4.39e-01 1.83e+01 angle pdb=" CA PRO C 426 " pdb=" N PRO C 426 " pdb=" CD PRO C 426 " ideal model delta sigma weight residual 112.00 106.01 5.99 1.40e+00 5.10e-01 1.83e+01 ... (remaining 29090 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.83: 11119 17.83 - 35.66: 1024 35.66 - 53.50: 276 53.50 - 71.33: 78 71.33 - 89.16: 32 Dihedral angle restraints: 12529 sinusoidal: 4801 harmonic: 7728 Sorted by residual: dihedral pdb=" CB CYS A1082 " pdb=" SG CYS A1082 " pdb=" SG CYS A1126 " pdb=" CB CYS A1126 " ideal model delta sinusoidal sigma weight residual -86.00 0.98 -86.98 1 1.00e+01 1.00e-02 9.09e+01 dihedral pdb=" CB CYS A 538 " pdb=" SG CYS A 538 " pdb=" SG CYS A 590 " pdb=" CB CYS A 590 " ideal model delta sinusoidal sigma weight residual -86.00 -16.94 -69.06 1 1.00e+01 1.00e-02 6.17e+01 dihedral pdb=" CB CYS B 617 " pdb=" SG CYS B 617 " pdb=" SG CYS B 649 " pdb=" CB CYS B 649 " ideal model delta sinusoidal sigma weight residual 93.00 30.18 62.82 1 1.00e+01 1.00e-02 5.21e+01 ... (remaining 12526 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.161: 3412 0.161 - 0.322: 1 0.322 - 0.483: 2 0.483 - 0.644: 0 0.644 - 0.805: 2 Chirality restraints: 3417 Sorted by residual: chirality pdb=" C1 NAG C1303 " pdb=" ND2 ASN C1098 " pdb=" C2 NAG C1303 " pdb=" O5 NAG C1303 " both_signs ideal model delta sigma weight residual False -2.40 -1.60 -0.80 2.00e-01 2.50e+01 1.62e+01 chirality pdb=" C1 NAG A1303 " pdb=" ND2 ASN A1098 " pdb=" C2 NAG A1303 " pdb=" O5 NAG A1303 " both_signs ideal model delta sigma weight residual False -2.40 -1.73 -0.67 2.00e-01 2.50e+01 1.14e+01 chirality pdb=" C1 NAG A1305 " pdb=" ND2 ASN A 709 " pdb=" C2 NAG A1305 " pdb=" O5 NAG A1305 " both_signs ideal model delta sigma weight residual False -2.40 -1.96 -0.44 2.00e-01 2.50e+01 4.94e+00 ... (remaining 3414 not shown) Planarity restraints: 3719 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C ALA A 520 " 0.089 5.00e-02 4.00e+02 1.28e-01 2.62e+01 pdb=" N PRO A 521 " -0.221 5.00e-02 4.00e+02 pdb=" CA PRO A 521 " 0.062 5.00e-02 4.00e+02 pdb=" CD PRO A 521 " 0.070 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CB ASN A1098 " 0.006 2.00e-02 2.50e+03 2.30e-02 6.61e+00 pdb=" CG ASN A1098 " -0.036 2.00e-02 2.50e+03 pdb=" OD1 ASN A1098 " 0.027 2.00e-02 2.50e+03 pdb=" ND2 ASN A1098 " -0.016 2.00e-02 2.50e+03 pdb=" C1 NAG A1303 " 0.019 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB ASN A 282 " 0.002 2.00e-02 2.50e+03 2.27e-02 6.45e+00 pdb=" CG ASN A 282 " 0.023 2.00e-02 2.50e+03 pdb=" OD1 ASN A 282 " -0.027 2.00e-02 2.50e+03 pdb=" ND2 ASN A 282 " 0.026 2.00e-02 2.50e+03 pdb=" C1 NAG A1307 " -0.025 2.00e-02 2.50e+03 ... (remaining 3716 not shown) Histogram of nonbonded interaction distances: 2.22 - 2.75: 2562 2.75 - 3.29: 19722 3.29 - 3.83: 35303 3.83 - 4.36: 42239 4.36 - 4.90: 72026 Nonbonded interactions: 171852 Sorted by model distance: nonbonded pdb=" OD1 ASP C 364 " pdb=" OG SER C 366 " model vdw 2.218 2.440 nonbonded pdb=" NE2 GLN C 804 " pdb=" OE1 GLN C 935 " model vdw 2.227 2.520 nonbonded pdb=" OD2 ASP C 578 " pdb=" OG1 THR C 581 " model vdw 2.233 2.440 nonbonded pdb=" NH1 ARG C 815 " pdb=" OD1 ASP C 820 " model vdw 2.235 2.520 nonbonded pdb=" OD1 ASP C 442 " pdb=" OH TYR C 451 " model vdw 2.269 2.440 ... (remaining 171847 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = (chain 'A' and (resid 27 through 108 or resid 116 through 130 or resid 168 throu \ gh 375 or resid 378 through 495 or resid 498 or resid 507 through 563 or (resid \ 564 and (name N or name CA or name C or name O or name CB )) or resid 565 throug \ h 570 or (resid 571 and (name N or name CA or name C or name O or name CB )) or \ resid 572 through 620 or resid 629 through 635 or resid 641 through 675 or resid \ 690 through 826 or resid 848 through 1106 or (resid 1107 and (name N or name CA \ or name C or name O or name CB )) or resid 1108 through 1138 or resid 1302 thro \ ugh 1308)) selection = (chain 'C' and (resid 27 through 404 or resid 429 through 443 or resid 451 or re \ sid 454 or resid 463 or resid 508 through 1138 or resid 1301 through 1307)) } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=0.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 9.030 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.020 Extract box with map and model: 27.360 Check model and map are aligned: 0.310 Set scattering table: 0.190 Process input model: 57.920 Find NCS groups from input model: 1.060 Set up NCS constraints: 0.100 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:12.220 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 108.210 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8483 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.091 21391 Z= 0.232 Angle : 0.646 12.783 29095 Z= 0.358 Chirality : 0.050 0.805 3417 Planarity : 0.005 0.128 3697 Dihedral : 16.022 89.161 7441 Min Nonbonded Distance : 2.218 Molprobity Statistics. All-atom Clashscore : 5.84 Ramachandran Plot: Outliers : 0.27 % Allowed : 3.58 % Favored : 96.14 % Rotamer: Outliers : 2.03 % Allowed : 16.77 % Favored : 81.20 % 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: -0.03 (0.16), residues: 2568 helix: 2.05 (0.21), residues: 614 sheet: 0.37 (0.20), residues: 582 loop : -1.21 (0.16), residues: 1372 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 5136 Ramachandran restraints generated. 2568 Oldfield, 0 Emsley, 2568 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 5136 Ramachandran restraints generated. 2568 Oldfield, 0 Emsley, 2568 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 242 residues out of total 2337 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 46 poor density : 196 time to evaluate : 2.239 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 46 outliers final: 38 residues processed: 234 average time/residue: 1.0710 time to fit residues: 294.2236 Evaluate side-chains 225 residues out of total 2337 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 38 poor density : 187 time to evaluate : 2.116 Switching outliers to nearest non-outliers outliers start: 38 outliers final: 34 residues processed: 4 average time/residue: 0.2217 time to fit residues: 4.5238 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 265 random chunks: chunk 223 optimal weight: 0.8980 chunk 200 optimal weight: 3.9990 chunk 111 optimal weight: 0.3980 chunk 68 optimal weight: 0.9990 chunk 135 optimal weight: 3.9990 chunk 107 optimal weight: 3.9990 chunk 207 optimal weight: 5.9990 chunk 80 optimal weight: 1.9990 chunk 126 optimal weight: 4.9990 chunk 154 optimal weight: 0.9990 chunk 240 optimal weight: 0.3980 overall best weight: 0.7384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 122 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 207 HIS A 394 ASN A 957 GLN ** B 30 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** B 762 GLN C 532 ASN C 641 ASN Total number of N/Q/H flips: 6 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8480 moved from start: 0.0474 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.062 21391 Z= 0.181 Angle : 0.552 7.236 29095 Z= 0.294 Chirality : 0.048 0.627 3417 Planarity : 0.004 0.061 3697 Dihedral : 4.494 30.759 2816 Min Nonbonded Distance : 2.203 Molprobity Statistics. All-atom Clashscore : 5.28 Ramachandran Plot: Outliers : 0.27 % Allowed : 3.86 % Favored : 95.87 % Rotamer: Outliers : 3.62 % Allowed : 15.36 % Favored : 81.02 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.04 % 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: 0.24 (0.16), residues: 2568 helix: 2.47 (0.20), residues: 610 sheet: 0.50 (0.20), residues: 579 loop : -1.12 (0.16), residues: 1379 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 5136 Ramachandran restraints generated. 2568 Oldfield, 0 Emsley, 2568 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 5136 Ramachandran restraints generated. 2568 Oldfield, 0 Emsley, 2568 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 272 residues out of total 2337 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 82 poor density : 190 time to evaluate : 2.231 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 82 outliers final: 50 residues processed: 249 average time/residue: 0.9651 time to fit residues: 285.7052 Evaluate side-chains 228 residues out of total 2337 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 50 poor density : 178 time to evaluate : 2.536 Switching outliers to nearest non-outliers outliers start: 50 outliers final: 35 residues processed: 15 average time/residue: 0.6423 time to fit residues: 15.4812 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 265 random chunks: chunk 133 optimal weight: 0.9980 chunk 74 optimal weight: 1.9990 chunk 200 optimal weight: 0.7980 chunk 163 optimal weight: 0.9980 chunk 66 optimal weight: 3.9990 chunk 241 optimal weight: 3.9990 chunk 260 optimal weight: 2.9990 chunk 214 optimal weight: 4.9990 chunk 239 optimal weight: 2.9990 chunk 82 optimal weight: 3.9990 chunk 193 optimal weight: 7.9990 overall best weight: 1.5584 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 122 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 394 ASN A 957 GLN B 30 ASN B 81 ASN B 703 ASN B 762 GLN B 965 GLN C 532 ASN C 641 ASN Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8511 moved from start: 0.0733 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.080 21391 Z= 0.292 Angle : 0.615 8.279 29095 Z= 0.326 Chirality : 0.050 0.670 3417 Planarity : 0.005 0.049 3697 Dihedral : 4.763 31.483 2816 Min Nonbonded Distance : 2.197 Molprobity Statistics. All-atom Clashscore : 5.88 Ramachandran Plot: Outliers : 0.27 % Allowed : 4.48 % Favored : 95.25 % Rotamer: Outliers : 3.88 % Allowed : 15.36 % Favored : 80.76 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.04 % 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: 0.11 (0.16), residues: 2568 helix: 2.37 (0.21), residues: 606 sheet: 0.45 (0.20), residues: 579 loop : -1.20 (0.15), residues: 1383 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 5136 Ramachandran restraints generated. 2568 Oldfield, 0 Emsley, 2568 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 5136 Ramachandran restraints generated. 2568 Oldfield, 0 Emsley, 2568 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 274 residues out of total 2337 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 88 poor density : 186 time to evaluate : 2.222 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 88 outliers final: 50 residues processed: 248 average time/residue: 1.0021 time to fit residues: 294.9522 Evaluate side-chains 227 residues out of total 2337 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 50 poor density : 177 time to evaluate : 2.489 Switching outliers to nearest non-outliers outliers start: 50 outliers final: 38 residues processed: 12 average time/residue: 0.3057 time to fit residues: 8.8657 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 265 random chunks: chunk 238 optimal weight: 2.9990 chunk 181 optimal weight: 3.9990 chunk 125 optimal weight: 0.5980 chunk 26 optimal weight: 5.9990 chunk 115 optimal weight: 4.9990 chunk 161 optimal weight: 2.9990 chunk 241 optimal weight: 0.6980 chunk 256 optimal weight: 0.6980 chunk 126 optimal weight: 1.9990 chunk 229 optimal weight: 0.8980 chunk 69 optimal weight: 0.7980 overall best weight: 0.7380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 81 ASN ** A 122 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** A 394 ASN A 957 GLN C 532 ASN C 641 ASN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8483 moved from start: 0.0787 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.072 21391 Z= 0.178 Angle : 0.541 8.084 29095 Z= 0.285 Chirality : 0.048 0.621 3417 Planarity : 0.004 0.044 3697 Dihedral : 4.455 29.996 2816 Min Nonbonded Distance : 2.226 Molprobity Statistics. All-atom Clashscore : 5.67 Ramachandran Plot: Outliers : 0.27 % Allowed : 3.74 % Favored : 95.99 % Rotamer: Outliers : 3.44 % Allowed : 16.33 % Favored : 80.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: 0.25 (0.16), residues: 2568 helix: 2.51 (0.20), residues: 618 sheet: 0.48 (0.21), residues: 569 loop : -1.13 (0.16), residues: 1381 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 5136 Ramachandran restraints generated. 2568 Oldfield, 0 Emsley, 2568 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 5136 Ramachandran restraints generated. 2568 Oldfield, 0 Emsley, 2568 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 268 residues out of total 2337 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 78 poor density : 190 time to evaluate : 2.300 Fit side-chains revert: symmetry clash revert: symmetry clash revert: symmetry clash revert: symmetry clash outliers start: 78 outliers final: 50 residues processed: 247 average time/residue: 0.9482 time to fit residues: 280.9641 Evaluate side-chains 230 residues out of total 2337 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 50 poor density : 180 time to evaluate : 2.287 Switching outliers to nearest non-outliers outliers start: 50 outliers final: 40 residues processed: 10 average time/residue: 0.5701 time to fit residues: 10.4777 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 265 random chunks: chunk 213 optimal weight: 0.8980 chunk 145 optimal weight: 0.6980 chunk 3 optimal weight: 3.9990 chunk 190 optimal weight: 0.8980 chunk 105 optimal weight: 5.9990 chunk 218 optimal weight: 3.9990 chunk 177 optimal weight: 5.9990 chunk 0 optimal weight: 2.9990 chunk 130 optimal weight: 1.9990 chunk 229 optimal weight: 1.9990 chunk 64 optimal weight: 2.9990 overall best weight: 1.2984 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Sorry: Reduce crashed with command 'molprobity.reduce -quiet -trim -'. Dumping stdin to file 'reduce_failure.pdb'. Return code: -15 Dumping stderr: