Starting phenix.real_space_refine on Fri Jul 19 17:46:43 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rx5_24725/07_2024/7rx5_24725.cif Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rx5_24725/07_2024/7rx5_24725.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.4 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rx5_24725/07_2024/7rx5_24725.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rx5_24725/07_2024/7rx5_24725.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rx5_24725/07_2024/7rx5_24725.cif" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7rx5_24725/07_2024/7rx5_24725.cif" } resolution = 3.4 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 ------------------------------------------------------------------------------- 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.003 sd= 0.068 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: 3 Type Number sf(0) Gaussians C 5640 2.51 5 N 1720 2.21 5 O 1760 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped None Time to flip residues: 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5388/modules/chem_data/mon_lib" Total number of atoms: 9120 Number of models: 1 Model: "" Number of chains: 40 Chain: "B" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "C" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "D" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "E" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "F" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "G" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "H" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "I" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "J" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "K" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "L" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "M" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "N" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "O" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "P" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "Q" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "R" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "S" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "T" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "U" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "A" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "V" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "W" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "X" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "Y" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "Z" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "0" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "1" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "2" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "3" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "4" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "5" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "6" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "7" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "8" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "9" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "a" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "b" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "c" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Chain: "d" Number of atoms: 228 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 228 Classifications: {'peptide': 29} Link IDs: {'TRANS': 28} Time building chain proxies: 5.20, per 1000 atoms: 0.57 Number of scatterers: 9120 At special positions: 0 Unit cell: (89.64, 89.64, 105.84, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 3 Type Number sf(0) O 1760 8.00 N 1720 7.00 C 5640 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=0, symmetry=0 Time building additional restraints: 3.68 Conformation dependent library (CDL) restraints added in 1.6 seconds 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 2320 Finding SS restraints... Secondary structure from input PDB file: 40 helices and 0 sheets defined 100.0% alpha, 0.0% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.70 Creating SS restraints... Processing helix chain 'A' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA A 7 " --> pdb=" O GLU A 3 " (cutoff:3.500A) Processing helix chain 'B' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA B 7 " --> pdb=" O GLU B 3 " (cutoff:3.500A) Processing helix chain 'C' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA C 7 " --> pdb=" O GLU C 3 " (cutoff:3.500A) Processing helix chain 'D' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA D 7 " --> pdb=" O GLU D 3 " (cutoff:3.500A) Processing helix chain 'E' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA E 7 " --> pdb=" O GLU E 3 " (cutoff:3.500A) Processing helix chain 'F' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA F 7 " --> pdb=" O GLU F 3 " (cutoff:3.500A) Processing helix chain 'G' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA G 7 " --> pdb=" O GLU G 3 " (cutoff:3.500A) Processing helix chain 'H' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA H 7 " --> pdb=" O GLU H 3 " (cutoff:3.500A) Processing helix chain 'I' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA I 7 " --> pdb=" O GLU I 3 " (cutoff:3.500A) Processing helix chain 'J' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA J 7 " --> pdb=" O GLU J 3 " (cutoff:3.500A) Processing helix chain 'K' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA K 7 " --> pdb=" O GLU K 3 " (cutoff:3.500A) Processing helix chain 'L' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA L 7 " --> pdb=" O GLU L 3 " (cutoff:3.500A) Processing helix chain 'M' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA M 7 " --> pdb=" O GLU M 3 " (cutoff:3.500A) Processing helix chain 'N' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA N 7 " --> pdb=" O GLU N 3 " (cutoff:3.500A) Processing helix chain 'O' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA O 7 " --> pdb=" O GLU O 3 " (cutoff:3.500A) Processing helix chain 'P' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA P 7 " --> pdb=" O GLU P 3 " (cutoff:3.500A) Processing helix chain 'Q' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA Q 7 " --> pdb=" O GLU Q 3 " (cutoff:3.500A) Processing helix chain 'R' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA R 7 " --> pdb=" O GLU R 3 " (cutoff:3.500A) Processing helix chain 'S' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA S 7 " --> pdb=" O GLU S 3 " (cutoff:3.500A) Processing helix chain 'T' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA T 7 " --> pdb=" O GLU T 3 " (cutoff:3.500A) Processing helix chain 'U' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA U 7 " --> pdb=" O GLU U 3 " (cutoff:3.500A) Processing helix chain 'V' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA V 7 " --> pdb=" O GLU V 3 " (cutoff:3.500A) Processing helix chain 'W' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA W 7 " --> pdb=" O GLU W 3 " (cutoff:3.500A) Processing helix chain 'X' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA X 7 " --> pdb=" O GLU X 3 " (cutoff:3.500A) Processing helix chain 'Y' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA Y 7 " --> pdb=" O GLU Y 3 " (cutoff:3.500A) Processing helix chain 'Z' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA Z 7 " --> pdb=" O GLU Z 3 " (cutoff:3.500A) Processing helix chain '0' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA 0 7 " --> pdb=" O GLU 0 3 " (cutoff:3.500A) Processing helix chain '1' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA 1 7 " --> pdb=" O GLU 1 3 " (cutoff:3.500A) Processing helix chain '2' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA 2 7 " --> pdb=" O GLU 2 3 " (cutoff:3.500A) Processing helix chain '3' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA 3 7 " --> pdb=" O GLU 3 3 " (cutoff:3.500A) Processing helix chain '4' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA 4 7 " --> pdb=" O GLU 4 3 " (cutoff:3.500A) Processing helix chain '5' and resid 1 through 29 removed outlier: 3.912A pdb=" N ALA 5 7 " --> pdb=" O GLU 5 3 " (cutoff:3.500A) Processing helix chain '6' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA 6 7 " --> pdb=" O GLU 6 3 " (cutoff:3.500A) Processing helix chain '7' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA 7 7 " --> pdb=" O GLU 7 3 " (cutoff:3.500A) Processing helix chain '8' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA 8 7 " --> pdb=" O GLU 8 3 " (cutoff:3.500A) Processing helix chain '9' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA 9 7 " --> pdb=" O GLU 9 3 " (cutoff:3.500A) Processing helix chain 'a' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA a 7 " --> pdb=" O GLU a 3 " (cutoff:3.500A) Processing helix chain 'b' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA b 7 " --> pdb=" O GLU b 3 " (cutoff:3.500A) Processing helix chain 'c' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA c 7 " --> pdb=" O GLU c 3 " (cutoff:3.500A) Processing helix chain 'd' and resid 1 through 29 removed outlier: 3.911A pdb=" N ALA d 7 " --> pdb=" O GLU d 3 " (cutoff:3.500A) 960 hydrogen bonds defined for protein. 2880 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 0 hydrogen bonds 0 hydrogen bond angles 0 basepair planarities 0 basepair parallelities 0 stacking parallelities Total time for adding SS restraints: 2.65 Time building geometry restraints manager: 3.99 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.23 - 1.29: 1720 1.29 - 1.35: 1626 1.35 - 1.42: 374 1.42 - 1.48: 1240 1.48 - 1.54: 4200 Bond restraints: 9160 Sorted by residual: bond pdb=" N ILE B 18 " pdb=" CA ILE B 18 " ideal model delta sigma weight residual 1.461 1.466 -0.005 1.19e-02 7.06e+03 1.74e-01 bond pdb=" N ILE C 18 " pdb=" CA ILE C 18 " ideal model delta sigma weight residual 1.461 1.466 -0.005 1.19e-02 7.06e+03 1.74e-01 bond pdb=" N HIS F 22 " pdb=" CA HIS F 22 " ideal model delta sigma weight residual 1.459 1.464 -0.005 1.21e-02 6.83e+03 1.50e-01 bond pdb=" N HIS G 22 " pdb=" CA HIS G 22 " ideal model delta sigma weight residual 1.459 1.464 -0.005 1.21e-02 6.83e+03 1.50e-01 bond pdb=" N ILE N 18 " pdb=" CA ILE N 18 " ideal model delta sigma weight residual 1.461 1.465 -0.004 1.19e-02 7.06e+03 1.33e-01 ... (remaining 9155 not shown) Histogram of bond angle deviations from ideal: 106.11 - 111.13: 3564 111.13 - 116.15: 2328 116.15 - 121.18: 4508 121.18 - 126.20: 1840 126.20 - 131.22: 40 Bond angle restraints: 12280 Sorted by residual: angle pdb=" C ILE N 25 " pdb=" CA ILE N 25 " pdb=" CB ILE N 25 " ideal model delta sigma weight residual 111.97 110.97 1.00 1.28e+00 6.10e-01 6.07e-01 angle pdb=" C ILE O 25 " pdb=" CA ILE O 25 " pdb=" CB ILE O 25 " ideal model delta sigma weight residual 111.97 110.97 1.00 1.28e+00 6.10e-01 6.07e-01 angle pdb=" C ILE 8 25 " pdb=" CA ILE 8 25 " pdb=" CB ILE 8 25 " ideal model delta sigma weight residual 111.97 110.98 0.99 1.28e+00 6.10e-01 5.95e-01 angle pdb=" C ILE 9 25 " pdb=" CA ILE 9 25 " pdb=" CB ILE 9 25 " ideal model delta sigma weight residual 111.97 110.98 0.99 1.28e+00 6.10e-01 5.95e-01 angle pdb=" C ILE c 25 " pdb=" CA ILE c 25 " pdb=" CB ILE c 25 " ideal model delta sigma weight residual 111.97 110.99 0.98 1.28e+00 6.10e-01 5.82e-01 ... (remaining 12275 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 7.23: 4440 7.23 - 14.46: 680 14.46 - 21.68: 320 21.68 - 28.91: 200 28.91 - 36.14: 160 Dihedral angle restraints: 5800 sinusoidal: 2320 harmonic: 3480 Sorted by residual: dihedral pdb=" CB ARG 7 13 " pdb=" CG ARG 7 13 " pdb=" CD ARG 7 13 " pdb=" NE ARG 7 13 " ideal model delta sinusoidal sigma weight residual -180.00 -143.86 -36.14 3 1.50e+01 4.44e-03 6.24e+00 dihedral pdb=" CB ARG 6 13 " pdb=" CG ARG 6 13 " pdb=" CD ARG 6 13 " pdb=" NE ARG 6 13 " ideal model delta sinusoidal sigma weight residual -180.00 -143.86 -36.14 3 1.50e+01 4.44e-03 6.24e+00 dihedral pdb=" CB ARG B 13 " pdb=" CG ARG B 13 " pdb=" CD ARG B 13 " pdb=" NE ARG B 13 " ideal model delta sinusoidal sigma weight residual -180.00 -143.88 -36.12 3 1.50e+01 4.44e-03 6.23e+00 ... (remaining 5797 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.017: 782 0.017 - 0.035: 440 0.035 - 0.052: 58 0.052 - 0.070: 40 0.070 - 0.087: 80 Chirality restraints: 1400 Sorted by residual: chirality pdb=" CA ILE H 25 " pdb=" N ILE H 25 " pdb=" C ILE H 25 " pdb=" CB ILE H 25 " both_signs ideal model delta sigma weight residual False 2.43 2.52 -0.09 2.00e-01 2.50e+01 1.89e-01 chirality pdb=" CA ILE I 25 " pdb=" N ILE I 25 " pdb=" C ILE I 25 " pdb=" CB ILE I 25 " both_signs ideal model delta sigma weight residual False 2.43 2.52 -0.09 2.00e-01 2.50e+01 1.89e-01 chirality pdb=" CA ILE O 25 " pdb=" N ILE O 25 " pdb=" C ILE O 25 " pdb=" CB ILE O 25 " both_signs ideal model delta sigma weight residual False 2.43 2.52 -0.09 2.00e-01 2.50e+01 1.86e-01 ... (remaining 1397 not shown) Planarity restraints: 1640 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB HIS D 22 " 0.003 2.00e-02 2.50e+03 2.68e-03 1.07e-01 pdb=" CG HIS D 22 " -0.006 2.00e-02 2.50e+03 pdb=" ND1 HIS D 22 " 0.001 2.00e-02 2.50e+03 pdb=" CD2 HIS D 22 " 0.001 2.00e-02 2.50e+03 pdb=" CE1 HIS D 22 " 0.000 2.00e-02 2.50e+03 pdb=" NE2 HIS D 22 " 0.001 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB HIS E 22 " 0.003 2.00e-02 2.50e+03 2.68e-03 1.07e-01 pdb=" CG HIS E 22 " -0.006 2.00e-02 2.50e+03 pdb=" ND1 HIS E 22 " 0.001 2.00e-02 2.50e+03 pdb=" CD2 HIS E 22 " 0.001 2.00e-02 2.50e+03 pdb=" CE1 HIS E 22 " 0.000 2.00e-02 2.50e+03 pdb=" NE2 HIS E 22 " 0.001 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB HIS Y 22 " 0.002 2.00e-02 2.50e+03 2.65e-03 1.05e-01 pdb=" CG HIS Y 22 " -0.006 2.00e-02 2.50e+03 pdb=" ND1 HIS Y 22 " 0.001 2.00e-02 2.50e+03 pdb=" CD2 HIS Y 22 " 0.001 2.00e-02 2.50e+03 pdb=" CE1 HIS Y 22 " 0.000 2.00e-02 2.50e+03 pdb=" NE2 HIS Y 22 " 0.000 2.00e-02 2.50e+03 ... (remaining 1637 not shown) Histogram of nonbonded interaction distances: 2.29 - 2.81: 2424 2.81 - 3.33: 9756 3.33 - 3.86: 15358 3.86 - 4.38: 18174 4.38 - 4.90: 29490 Nonbonded interactions: 75202 Sorted by model distance: nonbonded pdb=" OE1 GLN Y 29 " pdb=" NH1 ARG 0 13 " model vdw 2.288 2.520 nonbonded pdb=" OE1 GLN Z 29 " pdb=" NH1 ARG 1 13 " model vdw 2.288 2.520 nonbonded pdb=" OE1 GLN S 29 " pdb=" NH1 ARG U 13 " model vdw 2.288 2.520 nonbonded pdb=" OE1 GLN R 29 " pdb=" NH1 ARG T 13 " model vdw 2.288 2.520 nonbonded pdb=" OE1 GLN U 29 " pdb=" NH1 ARG V 13 " model vdw 2.288 2.520 ... (remaining 75197 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 '0' selection = chain '1' selection = chain '2' selection = chain '3' selection = chain '4' selection = chain '5' selection = chain '6' selection = chain '7' selection = chain '8' selection = chain '9' selection = chain 'A' selection = chain 'B' selection = chain 'C' selection = chain 'D' selection = chain 'E' selection = chain 'F' selection = chain 'G' selection = chain 'H' selection = chain 'I' selection = chain 'J' selection = chain 'K' selection = chain 'L' selection = chain 'M' selection = chain 'N' selection = chain 'O' selection = chain 'P' selection = chain 'Q' selection = chain 'R' selection = chain 'S' selection = chain 'T' selection = chain 'U' selection = chain 'V' selection = chain 'W' selection = chain 'X' selection = chain 'Y' selection = chain 'Z' selection = chain 'a' selection = chain 'b' selection = chain 'c' selection = chain 'd' } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 1.880 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.010 Extract box with map and model: 0.330 Check model and map are aligned: 0.060 Set scattering table: 0.080 Process input model: 26.780 Find NCS groups from input model: 0.450 Set up NCS constraints: 0.070 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:2.650 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 32.310 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8156 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.002 0.008 9160 Z= 0.117 Angle : 0.246 1.318 12280 Z= 0.156 Chirality : 0.029 0.087 1400 Planarity : 0.001 0.006 1640 Dihedral : 11.982 36.141 3480 Min Nonbonded Distance : 2.288 Molprobity Statistics. All-atom Clashscore : 5.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % 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: 4.54 (0.24), residues: 1080 helix: 2.95 (0.15), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.002 HIS Y 22 TYR 0.002 0.001 TYR T 15 ARG 0.001 0.000 ARG P 17 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 80 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 80 time to evaluate : 0.970 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 80 average time/residue: 0.4211 time to fit residues: 40.9229 Evaluate side-chains 80 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 80 time to evaluate : 0.992 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 80 random chunks: chunk 67 optimal weight: 1.9990 chunk 60 optimal weight: 10.0000 chunk 33 optimal weight: 1.9990 chunk 20 optimal weight: 10.0000 chunk 40 optimal weight: 9.9990 chunk 32 optimal weight: 9.9990 chunk 62 optimal weight: 8.9990 chunk 24 optimal weight: 8.9990 chunk 38 optimal weight: 8.9990 chunk 46 optimal weight: 8.9990 chunk 72 optimal weight: 10.0000 overall best weight: 6.1990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8345 moved from start: 0.1918 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.014 0.061 9160 Z= 0.893 Angle : 0.724 4.026 12280 Z= 0.412 Chirality : 0.049 0.099 1400 Planarity : 0.007 0.055 1640 Dihedral : 3.778 7.306 1240 Min Nonbonded Distance : 2.175 Molprobity Statistics. All-atom Clashscore : 6.20 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 2.62 % Favored : 97.38 % 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.11 (0.24), residues: 1080 helix: 0.83 (0.15), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.009 0.004 HIS J 22 TYR 0.013 0.003 TYR b 15 ARG 0.009 0.002 ARG U 13 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 58 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 58 time to evaluate : 1.008 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 58 average time/residue: 0.5073 time to fit residues: 35.2085 Evaluate side-chains 58 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 58 time to evaluate : 1.069 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 80 random chunks: chunk 40 optimal weight: 5.9990 chunk 22 optimal weight: 8.9990 chunk 60 optimal weight: 9.9990 chunk 49 optimal weight: 0.6980 chunk 20 optimal weight: 7.9990 chunk 72 optimal weight: 10.0000 chunk 78 optimal weight: 9.9990 chunk 64 optimal weight: 10.0000 chunk 24 optimal weight: 8.9990 chunk 58 optimal weight: 10.0000 chunk 71 optimal weight: 2.9990 overall best weight: 5.3388 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8341 moved from start: 0.2041 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.012 0.051 9160 Z= 0.771 Angle : 0.637 3.424 12280 Z= 0.364 Chirality : 0.044 0.098 1400 Planarity : 0.006 0.046 1640 Dihedral : 3.829 7.605 1240 Min Nonbonded Distance : 2.158 Molprobity Statistics. All-atom Clashscore : 5.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 3.81 % Favored : 96.19 % 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.62 (0.24), residues: 1080 helix: 0.53 (0.15), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.008 0.004 HIS U 22 TYR 0.010 0.002 TYR O 15 ARG 0.007 0.001 ARG U 13 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 58 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 58 time to evaluate : 0.928 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 58 average time/residue: 0.4904 time to fit residues: 34.1492 Evaluate side-chains 58 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 58 time to evaluate : 0.979 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 80 random chunks: chunk 54 optimal weight: 5.9990 chunk 37 optimal weight: 2.9990 chunk 8 optimal weight: 9.9990 chunk 34 optimal weight: 4.9990 chunk 48 optimal weight: 10.0000 chunk 73 optimal weight: 0.9990 chunk 77 optimal weight: 0.6980 chunk 38 optimal weight: 9.9990 chunk 69 optimal weight: 0.6980 chunk 20 optimal weight: 9.9990 chunk 64 optimal weight: 9.9990 overall best weight: 2.0786 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8260 moved from start: 0.1655 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.023 9160 Z= 0.310 Angle : 0.357 1.926 12280 Z= 0.209 Chirality : 0.033 0.094 1400 Planarity : 0.003 0.022 1640 Dihedral : 3.518 7.367 1240 Min Nonbonded Distance : 2.209 Molprobity Statistics. All-atom Clashscore : 5.38 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 4.05 % Favored : 95.95 % 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.21 (0.24), residues: 1080 helix: 1.51 (0.15), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.005 0.002 HIS I 22 TYR 0.006 0.001 TYR b 15 ARG 0.003 0.001 ARG U 13 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 68 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 68 time to evaluate : 0.950 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 68 average time/residue: 0.4470 time to fit residues: 36.9601 Evaluate side-chains 68 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 68 time to evaluate : 0.933 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 80 random chunks: chunk 43 optimal weight: 0.6980 chunk 1 optimal weight: 3.9990 chunk 57 optimal weight: 0.9980 chunk 31 optimal weight: 0.7980 chunk 65 optimal weight: 1.9990 chunk 53 optimal weight: 2.9990 chunk 0 optimal weight: 10.0000 chunk 39 optimal weight: 2.9990 chunk 69 optimal weight: 1.9990 chunk 19 optimal weight: 0.6980 chunk 26 optimal weight: 3.9990 overall best weight: 1.0382 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8199 moved from start: 0.1400 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.015 9160 Z= 0.165 Angle : 0.276 2.845 12280 Z= 0.164 Chirality : 0.030 0.095 1400 Planarity : 0.002 0.013 1640 Dihedral : 3.376 6.651 1240 Min Nonbonded Distance : 2.237 Molprobity Statistics. All-atom Clashscore : 5.49 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 3.81 % Favored : 96.19 % 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: 3.66 (0.23), residues: 1080 helix: 2.41 (0.14), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.002 HIS R 22 TYR 0.003 0.001 TYR a 15 ARG 0.001 0.000 ARG 7 17 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 74 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 74 time to evaluate : 1.139 Fit side-chains REVERT: D 24 GLU cc_start: 0.7770 (tp30) cc_final: 0.7492 (tp30) REVERT: E 24 GLU cc_start: 0.7770 (tp30) cc_final: 0.7490 (tp30) REVERT: 8 24 GLU cc_start: 0.7715 (tp30) cc_final: 0.7498 (tp30) REVERT: 9 24 GLU cc_start: 0.7711 (tp30) cc_final: 0.7498 (tp30) outliers start: 0 outliers final: 0 residues processed: 74 average time/residue: 0.4359 time to fit residues: 39.2807 Evaluate side-chains 74 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 74 time to evaluate : 1.046 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 80 random chunks: chunk 69 optimal weight: 0.7980 chunk 15 optimal weight: 0.6980 chunk 45 optimal weight: 1.9990 chunk 19 optimal weight: 0.9980 chunk 77 optimal weight: 0.0980 chunk 64 optimal weight: 7.9990 chunk 35 optimal weight: 0.6980 chunk 6 optimal weight: 10.0000 chunk 25 optimal weight: 2.9990 chunk 40 optimal weight: 2.9990 chunk 74 optimal weight: 1.9990 overall best weight: 0.6580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8160 moved from start: 0.1330 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.022 9160 Z= 0.122 Angle : 0.261 4.079 12280 Z= 0.155 Chirality : 0.029 0.097 1400 Planarity : 0.001 0.009 1640 Dihedral : 3.272 6.504 1240 Min Nonbonded Distance : 2.255 Molprobity Statistics. All-atom Clashscore : 5.60 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 3.81 % Favored : 96.19 % 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: 4.60 (0.22), residues: 1080 helix: 2.98 (0.14), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.001 HIS b 22 TYR 0.003 0.000 TYR b 15 ARG 0.001 0.000 ARG 4 17 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 78 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 78 time to evaluate : 1.049 Fit side-chains REVERT: D 24 GLU cc_start: 0.7772 (tp30) cc_final: 0.7550 (tp30) REVERT: E 24 GLU cc_start: 0.7771 (tp30) cc_final: 0.7548 (tp30) outliers start: 0 outliers final: 0 residues processed: 78 average time/residue: 0.4199 time to fit residues: 40.0838 Evaluate side-chains 78 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 78 time to evaluate : 0.991 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 80 random chunks: chunk 8 optimal weight: 7.9990 chunk 44 optimal weight: 7.9990 chunk 56 optimal weight: 6.9990 chunk 43 optimal weight: 0.9990 chunk 65 optimal weight: 0.6980 chunk 77 optimal weight: 1.9990 chunk 48 optimal weight: 0.3980 chunk 47 optimal weight: 0.9980 chunk 35 optimal weight: 0.6980 chunk 30 optimal weight: 7.9990 chunk 46 optimal weight: 0.9990 overall best weight: 0.7582 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8165 moved from start: 0.1340 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.029 9160 Z= 0.139 Angle : 0.277 4.355 12280 Z= 0.163 Chirality : 0.029 0.101 1400 Planarity : 0.001 0.008 1640 Dihedral : 3.215 7.434 1240 Min Nonbonded Distance : 2.259 Molprobity Statistics. All-atom Clashscore : 6.09 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 4.05 % Favored : 95.95 % 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: 4.97 (0.22), residues: 1080 helix: 3.21 (0.14), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.001 HIS 1 22 TYR 0.002 0.000 TYR b 15 ARG 0.001 0.000 ARG 5 17 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 76 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 76 time to evaluate : 0.936 Fit side-chains REVERT: D 24 GLU cc_start: 0.7763 (tp30) cc_final: 0.7503 (tp30) REVERT: E 24 GLU cc_start: 0.7760 (tp30) cc_final: 0.7503 (tp30) outliers start: 0 outliers final: 0 residues processed: 76 average time/residue: 0.4045 time to fit residues: 37.5382 Evaluate side-chains 76 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 76 time to evaluate : 1.002 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 80 random chunks: chunk 23 optimal weight: 0.5980 chunk 15 optimal weight: 0.9990 chunk 14 optimal weight: 1.9990 chunk 49 optimal weight: 0.8980 chunk 52 optimal weight: 6.9990 chunk 38 optimal weight: 7.9990 chunk 7 optimal weight: 0.9990 chunk 60 optimal weight: 6.9990 chunk 70 optimal weight: 10.0000 chunk 73 optimal weight: 2.9990 chunk 67 optimal weight: 2.9990 overall best weight: 1.0986 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8190 moved from start: 0.1412 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.026 9160 Z= 0.178 Angle : 0.290 4.081 12280 Z= 0.169 Chirality : 0.030 0.102 1400 Planarity : 0.002 0.010 1640 Dihedral : 3.215 8.279 1240 Min Nonbonded Distance : 2.256 Molprobity Statistics. All-atom Clashscore : 6.09 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 4.05 % Favored : 95.95 % 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: 4.95 (0.22), residues: 1080 helix: 3.20 (0.14), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.002 HIS R 22 TYR 0.002 0.000 TYR b 15 ARG 0.001 0.000 ARG T 13 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 72 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 72 time to evaluate : 1.059 Fit side-chains REVERT: 6 27 LYS cc_start: 0.8436 (ttmm) cc_final: 0.8197 (tttp) REVERT: 7 27 LYS cc_start: 0.8432 (ttmm) cc_final: 0.8191 (tttp) outliers start: 0 outliers final: 0 residues processed: 72 average time/residue: 0.4377 time to fit residues: 38.1186 Evaluate side-chains 72 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 72 time to evaluate : 0.989 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 80 random chunks: chunk 71 optimal weight: 0.7980 chunk 73 optimal weight: 0.9980 chunk 43 optimal weight: 0.9990 chunk 31 optimal weight: 0.6980 chunk 56 optimal weight: 7.9990 chunk 22 optimal weight: 6.9990 chunk 64 optimal weight: 9.9990 chunk 67 optimal weight: 0.8980 chunk 47 optimal weight: 0.5980 chunk 75 optimal weight: 0.7980 chunk 46 optimal weight: 0.9980 overall best weight: 0.7580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8161 moved from start: 0.1374 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.025 9160 Z= 0.136 Angle : 0.273 3.893 12280 Z= 0.160 Chirality : 0.029 0.102 1400 Planarity : 0.001 0.009 1640 Dihedral : 3.217 8.773 1240 Min Nonbonded Distance : 2.261 Molprobity Statistics. All-atom Clashscore : 5.92 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 4.29 % Favored : 95.71 % 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: 5.19 (0.22), residues: 1080 helix: 3.35 (0.14), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.001 HIS b 22 TYR 0.002 0.000 TYR a 15 ARG 0.001 0.000 ARG 5 17 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 78 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 78 time to evaluate : 0.924 Fit side-chains REVERT: D 24 GLU cc_start: 0.7768 (tp30) cc_final: 0.7532 (tp30) REVERT: E 24 GLU cc_start: 0.7766 (tp30) cc_final: 0.7530 (tp30) outliers start: 0 outliers final: 0 residues processed: 78 average time/residue: 0.4405 time to fit residues: 41.7623 Evaluate side-chains 78 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 78 time to evaluate : 1.034 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 80 random chunks: chunk 36 optimal weight: 7.9990 chunk 52 optimal weight: 6.9990 chunk 79 optimal weight: 0.6980 chunk 73 optimal weight: 0.6980 chunk 63 optimal weight: 0.8980 chunk 6 optimal weight: 10.0000 chunk 49 optimal weight: 2.9990 chunk 38 optimal weight: 7.9990 chunk 50 optimal weight: 6.9990 chunk 67 optimal weight: 2.9990 chunk 19 optimal weight: 0.7980 overall best weight: 1.2182 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8195 moved from start: 0.1426 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.024 9160 Z= 0.193 Angle : 0.293 3.851 12280 Z= 0.170 Chirality : 0.031 0.101 1400 Planarity : 0.002 0.012 1640 Dihedral : 3.205 9.026 1240 Min Nonbonded Distance : 2.257 Molprobity Statistics. All-atom Clashscore : 5.92 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 4.29 % Favored : 95.71 % 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: 4.98 (0.22), residues: 1080 helix: 3.22 (0.14), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.002 HIS S 22 TYR 0.002 0.000 TYR b 15 ARG 0.001 0.000 ARG T 13 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2160 Ramachandran restraints generated. 1080 Oldfield, 0 Emsley, 1080 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 72 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 72 time to evaluate : 1.023 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 72 average time/residue: 0.4454 time to fit residues: 38.8101 Evaluate side-chains 72 residues out of total 840 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 72 time to evaluate : 1.063 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 80 random chunks: chunk 58 optimal weight: 10.0000 chunk 9 optimal weight: 1.9990 chunk 17 optimal weight: 0.5980 chunk 63 optimal weight: 0.7980 chunk 26 optimal weight: 6.9990 chunk 65 optimal weight: 0.7980 chunk 8 optimal weight: 6.9990 chunk 11 optimal weight: 0.6980 chunk 55 optimal weight: 0.8980 chunk 3 optimal weight: 2.9990 chunk 45 optimal weight: 0.9980 overall best weight: 0.7580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3528 r_free = 0.3528 target = 0.135426 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 50)----------------| | r_work = 0.3154 r_free = 0.3154 target = 0.108858 restraints weight = 8924.154| |-----------------------------------------------------------------------------| r_work (start): 0.3041 rms_B_bonded: 0.98 r_work: 0.2897 rms_B_bonded: 1.92 restraints_weight: 0.5000 r_work: 0.2725 rms_B_bonded: 3.34 restraints_weight: 0.2500 r_work (final): 0.2725 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8359 moved from start: 0.1376 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.002 0.024 9160 Z= 0.137 Angle : 0.273 3.720 12280 Z= 0.160 Chirality : 0.029 0.102 1400 Planarity : 0.001 0.009 1640 Dihedral : 3.212 9.099 1240 Min Nonbonded Distance : 2.263 Molprobity Statistics. All-atom Clashscore : 5.92 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Rotamer: Outliers : 0.00 % Allowed : 4.29 % Favored : 95.71 % 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: 5.27 (0.22), residues: 1080 helix: 3.40 (0.14), residues: 1080 sheet: None (None), residues: 0 loop : None (None), residues: 0 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.003 0.001 HIS b 22 TYR 0.002 0.000 TYR b 15 ARG 0.001 0.000 ARG 5 17 =============================================================================== Job complete usr+sys time: 1704.11 seconds wall clock time: 31 minutes 14.37 seconds (1874.37 seconds total)