Starting phenix.real_space_refine (version: dev) on Sun Dec 11 23:03:22 2022 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6h6b_0148/12_2022/6h6b_0148.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6h6b_0148/12_2022/6h6b_0148.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/6h6b_0148/12_2022/6h6b_0148.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6h6b_0148/12_2022/6h6b_0148.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6h6b_0148/12_2022/6h6b_0148.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/6h6b_0148/12_2022/6h6b_0148.pdb" } resolution = 3.4 write_initial_geo_file = False refinement { macro_cycles = 10 } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= 0.002 sd= 0.007 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 Convert atoms to be neutral Process input model Symmetric amino acids flipped Residue "A PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "B PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "C PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "D PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "E PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "F PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "G PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "H PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "I PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Residue "J PHE 94": "CD1" <-> "CD2" "CE1" <-> "CE2" "HD1" <-> "HD2" "HE1" <-> "HE2" Time to flip residues: 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4807/modules/chem_data/mon_lib" Total number of atoms: 8060 Number of models: 1 Model: "" Number of chains: 10 Chain: "A" Number of atoms: 806 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 806 Classifications: {'peptide': 58} Link IDs: {'TRANS': 57} Chain: "B" Number of atoms: 806 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 806 Classifications: {'peptide': 58} Link IDs: {'TRANS': 57} Chain: "C" Number of atoms: 806 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 806 Classifications: {'peptide': 58} Link IDs: {'TRANS': 57} Chain: "D" Number of atoms: 806 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 806 Classifications: {'peptide': 58} Link IDs: {'TRANS': 57} Chain: "E" Number of atoms: 806 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 806 Classifications: {'peptide': 58} Link IDs: {'TRANS': 57} Chain: "F" Number of atoms: 806 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 806 Classifications: {'peptide': 58} Link IDs: {'TRANS': 57} Chain: "G" Number of atoms: 806 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 806 Classifications: {'peptide': 58} Link IDs: {'TRANS': 57} Chain: "H" Number of atoms: 806 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 806 Classifications: {'peptide': 58} Link IDs: {'TRANS': 57} Chain: "I" Number of atoms: 806 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 806 Classifications: {'peptide': 58} Link IDs: {'TRANS': 57} Chain: "J" Number of atoms: 806 Number of conformers: 1 Conformer: "" Number of residues, atoms: 58, 806 Classifications: {'peptide': 58} Link IDs: {'TRANS': 57} Time building chain proxies: 3.95, per 1000 atoms: 0.49 Number of scatterers: 8060 At special positions: 0 Unit cell: (92.241, 90.579, 42.381, 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) O 800 8.00 N 680 7.00 C 2480 6.00 H 4100 1.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: 6.89 Conformation dependent library (CDL) restraints added in 728.3 milliseconds 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 980 Finding SS restraints... Secondary structure from input PDB file: 0 helices and 10 sheets defined 0.0% alpha, 36.2% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 0.34 Creating SS restraints... Processing sheet with id=AA1, first strand: chain 'C' and resid 46 through 56 removed outlier: 6.981A pdb=" N VAL B 48 " --> pdb=" O GLY C 47 " (cutoff:3.500A) removed outlier: 6.971A pdb=" N VAL C 49 " --> pdb=" O VAL B 48 " (cutoff:3.500A) removed outlier: 6.949A pdb=" N HIS B 50 " --> pdb=" O VAL C 49 " (cutoff:3.500A) removed outlier: 6.914A pdb=" N GLY C 51 " --> pdb=" O HIS B 50 " (cutoff:3.500A) removed outlier: 6.968A pdb=" N VAL A 48 " --> pdb=" O GLY B 47 " (cutoff:3.500A) removed outlier: 6.977A pdb=" N VAL B 49 " --> pdb=" O VAL A 48 " (cutoff:3.500A) removed outlier: 6.947A pdb=" N HIS A 50 " --> pdb=" O VAL B 49 " (cutoff:3.500A) removed outlier: 6.938A pdb=" N GLY B 51 " --> pdb=" O HIS A 50 " (cutoff:3.500A) removed outlier: 7.835A pdb=" N GLY A 47 " --> pdb=" O VAL H 48 " (cutoff:3.500A) removed outlier: 10.245A pdb=" N HIS H 50 " --> pdb=" O GLY A 47 " (cutoff:3.500A) removed outlier: 8.886A pdb=" N VAL A 49 " --> pdb=" O HIS H 50 " (cutoff:3.500A) removed outlier: 9.228A pdb=" N VAL H 52 " --> pdb=" O VAL A 49 " (cutoff:3.500A) removed outlier: 6.340A pdb=" N GLY A 51 " --> pdb=" O VAL H 52 " (cutoff:3.500A) removed outlier: 7.701A pdb=" N THR H 54 " --> pdb=" O GLY A 51 " (cutoff:3.500A) removed outlier: 6.350A pdb=" N ALA A 53 " --> pdb=" O THR H 54 " (cutoff:3.500A) removed outlier: 8.199A pdb=" N ALA H 56 " --> pdb=" O ALA A 53 " (cutoff:3.500A) removed outlier: 6.625A pdb=" N VAL A 55 " --> pdb=" O ALA H 56 " (cutoff:3.500A) removed outlier: 6.992A pdb=" N VAL G 48 " --> pdb=" O GLY H 47 " (cutoff:3.500A) removed outlier: 6.990A pdb=" N VAL H 49 " --> pdb=" O VAL G 48 " (cutoff:3.500A) removed outlier: 6.965A pdb=" N HIS G 50 " --> pdb=" O VAL H 49 " (cutoff:3.500A) removed outlier: 6.943A pdb=" N GLY H 51 " --> pdb=" O HIS G 50 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'C' and resid 63 through 66 removed outlier: 6.255A pdb=" N VAL B 63 " --> pdb=" O THR C 64 " (cutoff:3.500A) removed outlier: 8.547A pdb=" N VAL C 66 " --> pdb=" O VAL B 63 " (cutoff:3.500A) removed outlier: 6.932A pdb=" N ASN B 65 " --> pdb=" O VAL C 66 " (cutoff:3.500A) removed outlier: 6.258A pdb=" N VAL G 63 " --> pdb=" O THR H 64 " (cutoff:3.500A) removed outlier: 8.548A pdb=" N VAL H 66 " --> pdb=" O VAL G 63 " (cutoff:3.500A) removed outlier: 6.924A pdb=" N ASN G 65 " --> pdb=" O VAL H 66 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'C' and resid 70 through 71 removed outlier: 6.566A pdb=" N VAL B 70 " --> pdb=" O VAL C 71 " (cutoff:3.500A) removed outlier: 6.542A pdb=" N VAL A 70 " --> pdb=" O VAL B 71 " (cutoff:3.500A) removed outlier: 6.553A pdb=" N VAL G 70 " --> pdb=" O VAL H 71 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'C' and resid 74 through 77 removed outlier: 6.890A pdb=" N VAL B 77 " --> pdb=" O ALA C 76 " (cutoff:3.500A) removed outlier: 6.891A pdb=" N VAL A 77 " --> pdb=" O ALA B 76 " (cutoff:3.500A) removed outlier: 6.196A pdb=" N VAL A 74 " --> pdb=" O THR H 75 " (cutoff:3.500A) removed outlier: 6.478A pdb=" N VAL H 77 " --> pdb=" O VAL A 74 " (cutoff:3.500A) removed outlier: 7.649A pdb=" N ALA A 76 " --> pdb=" O VAL H 77 " (cutoff:3.500A) removed outlier: 6.903A pdb=" N VAL G 77 " --> pdb=" O ALA H 76 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'C' and resid 88 through 92 removed outlier: 6.965A pdb=" N ALA C 90 " --> pdb=" O ALA B 89 " (cutoff:3.500A) removed outlier: 6.901A pdb=" N ALA B 91 " --> pdb=" O ALA C 90 " (cutoff:3.500A) removed outlier: 6.996A pdb=" N THR C 92 " --> pdb=" O ALA B 91 " (cutoff:3.500A) removed outlier: 6.939A pdb=" N ALA B 90 " --> pdb=" O ALA A 89 " (cutoff:3.500A) removed outlier: 6.870A pdb=" N ALA A 91 " --> pdb=" O ALA B 90 " (cutoff:3.500A) removed outlier: 6.958A pdb=" N THR B 92 " --> pdb=" O ALA A 91 " (cutoff:3.500A) removed outlier: 8.462A pdb=" N ILE A 88 " --> pdb=" O ALA H 89 " (cutoff:3.500A) removed outlier: 7.977A pdb=" N ALA H 91 " --> pdb=" O ILE A 88 " (cutoff:3.500A) removed outlier: 9.217A pdb=" N ALA A 90 " --> pdb=" O ALA H 91 " (cutoff:3.500A) removed outlier: 6.958A pdb=" N ALA H 90 " --> pdb=" O ALA G 89 " (cutoff:3.500A) removed outlier: 6.889A pdb=" N ALA G 91 " --> pdb=" O ALA H 90 " (cutoff:3.500A) removed outlier: 6.979A pdb=" N THR H 92 " --> pdb=" O ALA G 91 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'F' and resid 46 through 56 removed outlier: 6.968A pdb=" N VAL D 48 " --> pdb=" O GLY F 47 " (cutoff:3.500A) removed outlier: 6.965A pdb=" N VAL F 49 " --> pdb=" O VAL D 48 " (cutoff:3.500A) removed outlier: 6.943A pdb=" N HIS D 50 " --> pdb=" O VAL F 49 " (cutoff:3.500A) removed outlier: 6.914A pdb=" N GLY F 51 " --> pdb=" O HIS D 50 " (cutoff:3.500A) removed outlier: 7.879A pdb=" N GLY D 47 " --> pdb=" O VAL E 48 " (cutoff:3.500A) removed outlier: 10.287A pdb=" N HIS E 50 " --> pdb=" O GLY D 47 " (cutoff:3.500A) removed outlier: 8.929A pdb=" N VAL D 49 " --> pdb=" O HIS E 50 " (cutoff:3.500A) removed outlier: 9.270A pdb=" N VAL E 52 " --> pdb=" O VAL D 49 " (cutoff:3.500A) removed outlier: 6.377A pdb=" N GLY D 51 " --> pdb=" O VAL E 52 " (cutoff:3.500A) removed outlier: 7.729A pdb=" N THR E 54 " --> pdb=" O GLY D 51 " (cutoff:3.500A) removed outlier: 6.377A pdb=" N ALA D 53 " --> pdb=" O THR E 54 " (cutoff:3.500A) removed outlier: 8.223A pdb=" N ALA E 56 " --> pdb=" O ALA D 53 " (cutoff:3.500A) removed outlier: 6.648A pdb=" N VAL D 55 " --> pdb=" O ALA E 56 " (cutoff:3.500A) removed outlier: 7.814A pdb=" N GLY E 47 " --> pdb=" O VAL J 48 " (cutoff:3.500A) removed outlier: 10.224A pdb=" N HIS J 50 " --> pdb=" O GLY E 47 " (cutoff:3.500A) removed outlier: 8.878A pdb=" N VAL E 49 " --> pdb=" O HIS J 50 " (cutoff:3.500A) removed outlier: 9.218A pdb=" N VAL J 52 " --> pdb=" O VAL E 49 " (cutoff:3.500A) removed outlier: 6.329A pdb=" N GLY E 51 " --> pdb=" O VAL J 52 " (cutoff:3.500A) removed outlier: 7.676A pdb=" N THR J 54 " --> pdb=" O GLY E 51 " (cutoff:3.500A) removed outlier: 6.330A pdb=" N ALA E 53 " --> pdb=" O THR J 54 " (cutoff:3.500A) removed outlier: 8.179A pdb=" N ALA J 56 " --> pdb=" O ALA E 53 " (cutoff:3.500A) removed outlier: 6.611A pdb=" N VAL E 55 " --> pdb=" O ALA J 56 " (cutoff:3.500A) removed outlier: 6.966A pdb=" N VAL I 48 " --> pdb=" O GLY J 47 " (cutoff:3.500A) removed outlier: 6.963A pdb=" N VAL J 49 " --> pdb=" O VAL I 48 " (cutoff:3.500A) removed outlier: 6.942A pdb=" N HIS I 50 " --> pdb=" O VAL J 49 " (cutoff:3.500A) removed outlier: 6.917A pdb=" N GLY J 51 " --> pdb=" O HIS I 50 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'F' and resid 63 through 66 removed outlier: 6.215A pdb=" N VAL D 63 " --> pdb=" O THR F 64 " (cutoff:3.500A) removed outlier: 8.507A pdb=" N VAL F 66 " --> pdb=" O VAL D 63 " (cutoff:3.500A) removed outlier: 6.889A pdb=" N ASN D 65 " --> pdb=" O VAL F 66 " (cutoff:3.500A) removed outlier: 6.610A pdb=" N THR D 64 " --> pdb=" O ASN E 65 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'F' and resid 70 through 71 removed outlier: 6.590A pdb=" N VAL D 70 " --> pdb=" O VAL F 71 " (cutoff:3.500A) removed outlier: 6.541A pdb=" N VAL I 70 " --> pdb=" O VAL J 71 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'F' and resid 74 through 77 removed outlier: 6.887A pdb=" N VAL D 77 " --> pdb=" O ALA F 76 " (cutoff:3.500A) removed outlier: 6.202A pdb=" N VAL D 74 " --> pdb=" O THR E 75 " (cutoff:3.500A) removed outlier: 6.469A pdb=" N VAL E 77 " --> pdb=" O VAL D 74 " (cutoff:3.500A) removed outlier: 7.646A pdb=" N ALA D 76 " --> pdb=" O VAL E 77 " (cutoff:3.500A) removed outlier: 6.212A pdb=" N VAL E 74 " --> pdb=" O THR J 75 " (cutoff:3.500A) removed outlier: 6.505A pdb=" N VAL J 77 " --> pdb=" O VAL E 74 " (cutoff:3.500A) removed outlier: 7.693A pdb=" N ALA E 76 " --> pdb=" O VAL J 77 " (cutoff:3.500A) removed outlier: 6.880A pdb=" N VAL I 77 " --> pdb=" O ALA J 76 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'F' and resid 88 through 92 removed outlier: 6.945A pdb=" N ALA F 90 " --> pdb=" O ALA D 89 " (cutoff:3.500A) removed outlier: 6.890A pdb=" N ALA D 91 " --> pdb=" O ALA F 90 " (cutoff:3.500A) removed outlier: 6.972A pdb=" N THR F 92 " --> pdb=" O ALA D 91 " (cutoff:3.500A) removed outlier: 8.494A pdb=" N ILE D 88 " --> pdb=" O ALA E 89 " (cutoff:3.500A) removed outlier: 7.987A pdb=" N ALA E 91 " --> pdb=" O ILE D 88 " (cutoff:3.500A) removed outlier: 9.230A pdb=" N ALA D 90 " --> pdb=" O ALA E 91 " (cutoff:3.500A) removed outlier: 8.461A pdb=" N ILE E 88 " --> pdb=" O ALA J 89 " (cutoff:3.500A) removed outlier: 7.993A pdb=" N ALA J 91 " --> pdb=" O ILE E 88 " (cutoff:3.500A) removed outlier: 9.228A pdb=" N ALA E 90 " --> pdb=" O ALA J 91 " (cutoff:3.500A) removed outlier: 6.942A pdb=" N ALA J 90 " --> pdb=" O ALA I 89 " (cutoff:3.500A) removed outlier: 6.872A pdb=" N ALA I 91 " --> pdb=" O ALA J 90 " (cutoff:3.500A) removed outlier: 6.965A pdb=" N THR J 92 " --> pdb=" O ALA I 91 " (cutoff:3.500A) 60 hydrogen bonds defined for protein. 180 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.52 Time building geometry restraints manager: 6.91 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 0.84 - 0.98: 4100 0.98 - 1.13: 0 1.13 - 1.27: 690 1.27 - 1.42: 800 1.42 - 1.57: 2490 Bond restraints: 8080 Sorted by residual: bond pdb=" CB VAL D 70 " pdb=" CG2 VAL D 70 " ideal model delta sigma weight residual 1.521 1.476 0.045 3.30e-02 9.18e+02 1.85e+00 bond pdb=" CB VAL I 70 " pdb=" CG2 VAL I 70 " ideal model delta sigma weight residual 1.521 1.476 0.045 3.30e-02 9.18e+02 1.84e+00 bond pdb=" CB VAL J 70 " pdb=" CG2 VAL J 70 " ideal model delta sigma weight residual 1.521 1.476 0.045 3.30e-02 9.18e+02 1.84e+00 bond pdb=" CB VAL C 70 " pdb=" CG2 VAL C 70 " ideal model delta sigma weight residual 1.521 1.476 0.045 3.30e-02 9.18e+02 1.84e+00 bond pdb=" CB VAL H 70 " pdb=" CG2 VAL H 70 " ideal model delta sigma weight residual 1.521 1.477 0.044 3.30e-02 9.18e+02 1.81e+00 ... (remaining 8075 not shown) Histogram of bond angle deviations from ideal: 83.12 - 92.71: 30 92.71 - 102.31: 0 102.31 - 111.90: 9868 111.90 - 121.50: 3005 121.50 - 131.09: 1757 Bond angle restraints: 14660 Sorted by residual: angle pdb=" CG2 VAL F 95 " pdb=" CB VAL F 95 " pdb=" HB VAL F 95 " ideal model delta sigma weight residual 108.00 83.12 24.88 3.00e+00 1.11e-01 6.88e+01 angle pdb=" CG2 VAL D 95 " pdb=" CB VAL D 95 " pdb=" HB VAL D 95 " ideal model delta sigma weight residual 108.00 83.12 24.88 3.00e+00 1.11e-01 6.88e+01 angle pdb=" CG2 VAL J 95 " pdb=" CB VAL J 95 " pdb=" HB VAL J 95 " ideal model delta sigma weight residual 108.00 83.18 24.82 3.00e+00 1.11e-01 6.84e+01 angle pdb=" CG2 VAL A 95 " pdb=" CB VAL A 95 " pdb=" HB VAL A 95 " ideal model delta sigma weight residual 108.00 83.19 24.81 3.00e+00 1.11e-01 6.84e+01 angle pdb=" CG2 VAL C 95 " pdb=" CB VAL C 95 " pdb=" HB VAL C 95 " ideal model delta sigma weight residual 108.00 83.19 24.81 3.00e+00 1.11e-01 6.84e+01 ... (remaining 14655 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 7.04: 2578 7.04 - 14.08: 365 14.08 - 21.13: 69 21.13 - 28.17: 38 28.17 - 35.21: 20 Dihedral angle restraints: 3070 sinusoidal: 1390 harmonic: 1680 Sorted by residual: dihedral pdb=" CA PHE B 94 " pdb=" C PHE B 94 " pdb=" N VAL B 95 " pdb=" CA VAL B 95 " ideal model delta harmonic sigma weight residual 180.00 154.59 25.41 0 5.00e+00 4.00e-02 2.58e+01 dihedral pdb=" CA PHE A 94 " pdb=" C PHE A 94 " pdb=" N VAL A 95 " pdb=" CA VAL A 95 " ideal model delta harmonic sigma weight residual 180.00 154.59 25.41 0 5.00e+00 4.00e-02 2.58e+01 dihedral pdb=" CA PHE I 94 " pdb=" C PHE I 94 " pdb=" N VAL I 95 " pdb=" CA VAL I 95 " ideal model delta harmonic sigma weight residual 180.00 154.60 25.40 0 5.00e+00 4.00e-02 2.58e+01 ... (remaining 3067 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.321: 710 0.321 - 0.642: 0 0.642 - 0.963: 0 0.963 - 1.284: 0 1.284 - 1.605: 10 Chirality restraints: 720 Sorted by residual: chirality pdb=" CB VAL F 95 " pdb=" CA VAL F 95 " pdb=" CG1 VAL F 95 " pdb=" CG2 VAL F 95 " both_signs ideal model delta sigma weight residual False -2.63 -1.02 -1.61 2.00e-01 2.50e+01 6.44e+01 chirality pdb=" CB VAL J 95 " pdb=" CA VAL J 95 " pdb=" CG1 VAL J 95 " pdb=" CG2 VAL J 95 " both_signs ideal model delta sigma weight residual False -2.63 -1.02 -1.60 2.00e-01 2.50e+01 6.44e+01 chirality pdb=" CB VAL E 95 " pdb=" CA VAL E 95 " pdb=" CG1 VAL E 95 " pdb=" CG2 VAL E 95 " both_signs ideal model delta sigma weight residual False -2.63 -1.02 -1.60 2.00e-01 2.50e+01 6.44e+01 ... (remaining 717 not shown) Planarity restraints: 1210 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C ALA A 76 " -0.006 2.00e-02 2.50e+03 1.07e-02 1.14e+00 pdb=" N VAL A 77 " 0.018 2.00e-02 2.50e+03 pdb=" CA VAL A 77 " -0.005 2.00e-02 2.50e+03 pdb=" H VAL A 77 " -0.008 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB HIS B 50 " -0.010 2.00e-02 2.50e+03 6.90e-03 1.07e+00 pdb=" CG HIS B 50 " 0.018 2.00e-02 2.50e+03 pdb=" ND1 HIS B 50 " -0.001 2.00e-02 2.50e+03 pdb=" CD2 HIS B 50 " -0.001 2.00e-02 2.50e+03 pdb=" CE1 HIS B 50 " 0.000 2.00e-02 2.50e+03 pdb=" NE2 HIS B 50 " -0.001 2.00e-02 2.50e+03 pdb=" HD1 HIS B 50 " 0.001 2.00e-02 2.50e+03 pdb=" HD2 HIS B 50 " -0.001 2.00e-02 2.50e+03 pdb=" HE1 HIS B 50 " -0.004 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB HIS I 50 " -0.010 2.00e-02 2.50e+03 6.87e-03 1.06e+00 pdb=" CG HIS I 50 " 0.018 2.00e-02 2.50e+03 pdb=" ND1 HIS I 50 " -0.001 2.00e-02 2.50e+03 pdb=" CD2 HIS I 50 " -0.002 2.00e-02 2.50e+03 pdb=" CE1 HIS I 50 " 0.000 2.00e-02 2.50e+03 pdb=" NE2 HIS I 50 " -0.001 2.00e-02 2.50e+03 pdb=" HD1 HIS I 50 " 0.000 2.00e-02 2.50e+03 pdb=" HD2 HIS I 50 " -0.001 2.00e-02 2.50e+03 pdb=" HE1 HIS I 50 " -0.004 2.00e-02 2.50e+03 ... (remaining 1207 not shown) Histogram of nonbonded interaction distances: 1.75 - 2.30: 1609 2.30 - 2.85: 18158 2.85 - 3.40: 16698 3.40 - 3.95: 25141 3.95 - 4.50: 36108 Nonbonded interactions: 97714 Sorted by model distance: nonbonded pdb=" OE1 GLU D 46 " pdb=" HZ1 LYS E 80 " model vdw 1.753 1.850 nonbonded pdb=" HB VAL B 95 " pdb="HG13 VAL B 95 " model vdw 1.764 1.952 nonbonded pdb=" HB VAL D 95 " pdb="HG13 VAL D 95 " model vdw 1.764 1.952 nonbonded pdb=" HB VAL C 95 " pdb="HG13 VAL C 95 " model vdw 1.766 1.952 nonbonded pdb=" HB VAL A 95 " pdb="HG13 VAL A 95 " model vdw 1.766 1.952 ... (remaining 97709 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' 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' } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Set scattering table Set to: electron Number of scattering types: 4 Type Number sf(0) Gaussians C 2480 2.51 5 N 680 2.21 5 O 800 1.98 5 H 4100 0.53 5 sf(0) = scattering factor at diffraction angle 0. 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: 0.820 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.130 Extract box with map and model: 5.070 Check model and map are aligned: 0.130 Convert atoms to be neutral: 0.070 Process input model: 29.000 Find NCS groups from input model: 0.420 Set up NCS constraints: 0.040 Set refine NCS operators: 0.000 Set scattering table: 0.010 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.900 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 38.590 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7351 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.008 0.045 3980 Z= 0.545 Angle : 1.101 10.230 5400 Z= 0.604 Chirality : 0.198 1.605 720 Planarity : 0.002 0.008 670 Dihedral : 8.131 35.210 1310 Min Nonbonded Distance : 2.530 Molprobity Statistics. All-atom Clashscore : 5.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 14.29 % Favored : 85.71 % Rotamer Outliers : 0.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.80 (0.21), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -3.65 (0.16), residues: 560 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 80 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 80 time to evaluate : 0.643 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 0 outliers final: 0 residues processed: 80 average time/residue: 0.2028 time to fit residues: 23.6419 Evaluate side-chains 54 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 54 time to evaluate : 0.720 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.7964 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 42 optimal weight: 0.8980 chunk 37 optimal weight: 0.8980 chunk 21 optimal weight: 5.9990 chunk 12 optimal weight: 0.5980 chunk 25 optimal weight: 10.0000 chunk 20 optimal weight: 9.9990 chunk 39 optimal weight: 0.9990 chunk 15 optimal weight: 0.9990 chunk 23 optimal weight: 4.9990 chunk 29 optimal weight: 1.9990 chunk 45 optimal weight: 8.9990 overall best weight: 0.8784 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 65 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7346 moved from start: 0.2309 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.034 3980 Z= 0.222 Angle : 0.805 11.090 5400 Z= 0.433 Chirality : 0.181 1.552 720 Planarity : 0.002 0.009 670 Dihedral : 5.761 17.623 580 Min Nonbonded Distance : 2.312 Molprobity Statistics. All-atom Clashscore : 3.84 Ramachandran Plot: Outliers : 0.00 % Allowed : 9.82 % Favored : 90.18 % Rotamer Outliers : 0.25 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -4.25 (0.23), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -3.23 (0.18), residues: 560 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 57 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 56 time to evaluate : 0.682 Fit side-chains revert: symmetry clash revert: symmetry clash outliers start: 1 outliers final: 0 residues processed: 56 average time/residue: 0.1780 time to fit residues: 14.8704 Evaluate side-chains 48 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 48 time to evaluate : 0.637 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.7467 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 25 optimal weight: 10.0000 chunk 14 optimal weight: 2.9990 chunk 37 optimal weight: 7.9990 chunk 30 optimal weight: 10.0000 chunk 12 optimal weight: 6.9990 chunk 45 optimal weight: 10.0000 chunk 49 optimal weight: 10.0000 chunk 40 optimal weight: 3.9990 chunk 15 optimal weight: 4.9990 chunk 36 optimal weight: 5.9990 chunk 44 optimal weight: 2.9990 overall best weight: 4.1990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 65 ASN H 65 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7543 moved from start: 0.3529 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.034 3980 Z= 0.351 Angle : 0.886 11.321 5400 Z= 0.481 Chirality : 0.187 1.591 720 Planarity : 0.002 0.010 670 Dihedral : 6.341 19.514 580 Min Nonbonded Distance : 2.083 Molprobity Statistics. All-atom Clashscore : 5.82 Ramachandran Plot: Outliers : 0.00 % Allowed : 15.54 % Favored : 84.46 % Rotamer Outliers : 1.50 % 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.04 (0.25), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -3.07 (0.19), residues: 560 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 64 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 58 time to evaluate : 0.663 Fit side-chains outliers start: 6 outliers final: 3 residues processed: 60 average time/residue: 0.1808 time to fit residues: 16.0432 Evaluate side-chains 49 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 46 time to evaluate : 0.631 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.0863 time to fit residues: 1.1985 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 34 optimal weight: 1.9990 chunk 23 optimal weight: 6.9990 chunk 5 optimal weight: 9.9990 chunk 21 optimal weight: 6.9990 chunk 30 optimal weight: 10.0000 chunk 45 optimal weight: 9.9990 chunk 48 optimal weight: 1.9990 chunk 43 optimal weight: 4.9990 chunk 13 optimal weight: 6.9990 chunk 40 optimal weight: 2.9990 chunk 27 optimal weight: 6.9990 overall best weight: 3.7990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: I 65 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7563 moved from start: 0.4160 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.030 3980 Z= 0.285 Angle : 0.850 10.972 5400 Z= 0.462 Chirality : 0.184 1.582 720 Planarity : 0.002 0.009 670 Dihedral : 6.116 20.901 580 Min Nonbonded Distance : 2.110 Molprobity Statistics. All-atom Clashscore : 5.69 Ramachandran Plot: Outliers : 0.00 % Allowed : 17.14 % Favored : 82.86 % Rotamer Outliers : 3.75 % 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.81 (0.27), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.90 (0.21), residues: 560 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 66 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 15 poor density : 51 time to evaluate : 0.633 Fit side-chains outliers start: 15 outliers final: 14 residues processed: 59 average time/residue: 0.1753 time to fit residues: 15.3049 Evaluate side-chains 59 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 45 time to evaluate : 0.625 Switching outliers to nearest non-outliers outliers start: 14 outliers final: 0 residues processed: 14 average time/residue: 0.0871 time to fit residues: 2.9181 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 0 optimal weight: 10.0000 chunk 35 optimal weight: 9.9990 chunk 19 optimal weight: 0.9990 chunk 41 optimal weight: 3.9990 chunk 33 optimal weight: 1.9990 chunk 24 optimal weight: 1.9990 chunk 43 optimal weight: 6.9990 chunk 12 optimal weight: 2.9990 chunk 16 optimal weight: 4.9990 chunk 9 optimal weight: 3.9990 chunk 28 optimal weight: 5.9990 overall best weight: 2.3990 ------------------------------------------------------------------------------- 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.7531 moved from start: 0.4480 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.030 3980 Z= 0.231 Angle : 0.799 11.253 5400 Z= 0.431 Chirality : 0.183 1.596 720 Planarity : 0.001 0.006 670 Dihedral : 5.585 17.408 580 Min Nonbonded Distance : 2.093 Molprobity Statistics. All-atom Clashscore : 4.83 Ramachandran Plot: Outliers : 0.00 % Allowed : 16.07 % Favored : 83.93 % Rotamer Outliers : 0.75 % 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.68 (0.28), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.80 (0.21), residues: 560 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 54 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 51 time to evaluate : 0.648 Fit side-chains outliers start: 3 outliers final: 2 residues processed: 52 average time/residue: 0.1846 time to fit residues: 14.3698 Evaluate side-chains 45 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 43 time to evaluate : 0.656 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0963 time to fit residues: 1.1097 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 11 optimal weight: 3.9990 chunk 48 optimal weight: 8.9990 chunk 40 optimal weight: 8.9990 chunk 22 optimal weight: 7.9990 chunk 4 optimal weight: 1.9990 chunk 16 optimal weight: 4.9990 chunk 25 optimal weight: 9.9990 chunk 46 optimal weight: 4.9990 chunk 5 optimal weight: 20.0000 chunk 27 optimal weight: 5.9990 chunk 35 optimal weight: 9.9990 overall best weight: 4.3990 ------------------------------------------------------------------------------- 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.7586 moved from start: 0.4870 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.030 3980 Z= 0.335 Angle : 0.870 11.062 5400 Z= 0.475 Chirality : 0.185 1.610 720 Planarity : 0.002 0.008 670 Dihedral : 6.031 21.788 580 Min Nonbonded Distance : 2.098 Molprobity Statistics. All-atom Clashscore : 6.44 Ramachandran Plot: Outliers : 0.00 % Allowed : 18.57 % Favored : 81.43 % Rotamer Outliers : 3.25 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.83 (0.28), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.91 (0.21), residues: 560 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 60 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 47 time to evaluate : 0.671 Fit side-chains outliers start: 13 outliers final: 11 residues processed: 52 average time/residue: 0.1821 time to fit residues: 14.2282 Evaluate side-chains 52 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 41 time to evaluate : 0.643 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 0 residues processed: 11 average time/residue: 0.0893 time to fit residues: 2.5410 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 27 optimal weight: 0.7980 chunk 40 optimal weight: 7.9990 chunk 48 optimal weight: 6.9990 chunk 30 optimal weight: 10.0000 chunk 29 optimal weight: 0.9990 chunk 22 optimal weight: 6.9990 chunk 19 optimal weight: 1.9990 chunk 28 optimal weight: 1.9990 chunk 14 optimal weight: 0.8980 chunk 9 optimal weight: 4.9990 chunk 32 optimal weight: 4.9990 overall best weight: 1.3386 ------------------------------------------------------------------------------- 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.7516 moved from start: 0.4993 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.031 3980 Z= 0.209 Angle : 0.773 11.037 5400 Z= 0.416 Chirality : 0.180 1.545 720 Planarity : 0.001 0.006 670 Dihedral : 5.240 15.066 580 Min Nonbonded Distance : 2.137 Molprobity Statistics. All-atom Clashscore : 4.08 Ramachandran Plot: Outliers : 0.00 % Allowed : 13.93 % Favored : 86.07 % Rotamer Outliers : 0.50 % 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.68 (0.28), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.79 (0.22), residues: 560 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 50 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 48 time to evaluate : 0.702 Fit side-chains revert: symmetry clash outliers start: 2 outliers final: 2 residues processed: 48 average time/residue: 0.1700 time to fit residues: 12.5072 Evaluate side-chains 45 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 43 time to evaluate : 0.624 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0914 time to fit residues: 1.0675 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 23 optimal weight: 7.9990 chunk 4 optimal weight: 10.0000 chunk 37 optimal weight: 5.9990 chunk 43 optimal weight: 8.9990 chunk 46 optimal weight: 6.9990 chunk 42 optimal weight: 3.9990 chunk 44 optimal weight: 2.9990 chunk 27 optimal weight: 6.9990 chunk 19 optimal weight: 3.9990 chunk 35 optimal weight: 9.9990 chunk 13 optimal weight: 5.9990 overall best weight: 4.5990 ------------------------------------------------------------------------------- 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.7583 moved from start: 0.5241 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.030 3980 Z= 0.343 Angle : 0.864 10.903 5400 Z= 0.472 Chirality : 0.184 1.549 720 Planarity : 0.002 0.006 670 Dihedral : 5.862 17.895 580 Min Nonbonded Distance : 2.118 Molprobity Statistics. All-atom Clashscore : 6.56 Ramachandran Plot: Outliers : 0.00 % Allowed : 17.86 % Favored : 82.14 % Rotamer Outliers : 1.75 % 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.84 (0.28), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.92 (0.21), residues: 560 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 48 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 41 time to evaluate : 0.640 Fit side-chains outliers start: 7 outliers final: 6 residues processed: 43 average time/residue: 0.1796 time to fit residues: 11.9125 Evaluate side-chains 46 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 40 time to evaluate : 0.597 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 0 residues processed: 6 average time/residue: 0.0982 time to fit residues: 1.7928 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 40 optimal weight: 8.9990 chunk 42 optimal weight: 5.9990 chunk 44 optimal weight: 8.9990 chunk 29 optimal weight: 6.9990 chunk 47 optimal weight: 2.9990 chunk 28 optimal weight: 4.9990 chunk 22 optimal weight: 2.9990 chunk 33 optimal weight: 0.9990 chunk 49 optimal weight: 0.9980 chunk 45 optimal weight: 0.9980 chunk 39 optimal weight: 0.6980 overall best weight: 1.3384 ------------------------------------------------------------------------------- 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.7511 moved from start: 0.5356 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.031 3980 Z= 0.208 Angle : 0.773 11.108 5400 Z= 0.416 Chirality : 0.181 1.545 720 Planarity : 0.001 0.006 670 Dihedral : 5.158 15.314 580 Min Nonbonded Distance : 2.151 Molprobity Statistics. All-atom Clashscore : 4.46 Ramachandran Plot: Outliers : 0.00 % Allowed : 14.29 % Favored : 85.71 % Rotamer Outliers : 0.50 % 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.74 (0.28), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.84 (0.21), residues: 560 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 42 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 40 time to evaluate : 0.709 Fit side-chains outliers start: 2 outliers final: 2 residues processed: 40 average time/residue: 0.1746 time to fit residues: 10.8682 Evaluate side-chains 40 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 38 time to evaluate : 0.674 Switching outliers to nearest non-outliers outliers start: 2 outliers final: 0 residues processed: 2 average time/residue: 0.0915 time to fit residues: 1.0721 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 4 optimal weight: 10.0000 chunk 30 optimal weight: 10.0000 chunk 24 optimal weight: 8.9990 chunk 31 optimal weight: 6.9990 chunk 42 optimal weight: 6.9990 chunk 12 optimal weight: 5.9990 chunk 36 optimal weight: 4.9990 chunk 5 optimal weight: 9.9990 chunk 11 optimal weight: 5.9990 chunk 39 optimal weight: 5.9990 chunk 16 optimal weight: 0.8980 overall best weight: 4.7788 ------------------------------------------------------------------------------- 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.7587 moved from start: 0.5485 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.031 3980 Z= 0.346 Angle : 0.862 10.844 5400 Z= 0.471 Chirality : 0.183 1.543 720 Planarity : 0.002 0.007 670 Dihedral : 5.857 18.795 580 Min Nonbonded Distance : 2.121 Molprobity Statistics. All-atom Clashscore : 6.81 Ramachandran Plot: Outliers : 0.00 % Allowed : 17.68 % Favored : 82.32 % Rotamer Outliers : 0.25 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -3.90 (0.27), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -2.96 (0.21), residues: 560 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1120 Ramachandran restraints generated. 560 Oldfield, 0 Emsley, 560 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 40 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 39 time to evaluate : 0.580 Fit side-chains outliers start: 1 outliers final: 1 residues processed: 40 average time/residue: 0.1794 time to fit residues: 10.8976 Evaluate side-chains 38 residues out of total 400 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 37 time to evaluate : 0.623 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1203 time to fit residues: 1.0465 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 50 random chunks: chunk 40 optimal weight: 8.9990 chunk 5 optimal weight: 10.0000 chunk 7 optimal weight: 6.9990 chunk 34 optimal weight: 5.9990 chunk 2 optimal weight: 3.9990 chunk 28 optimal weight: 5.9990 chunk 45 optimal weight: 6.9990 chunk 26 optimal weight: 6.9990 chunk 1 optimal weight: 3.9990 chunk 31 optimal weight: 0.6980 chunk 30 optimal weight: 10.0000 overall best weight: 4.1388 ------------------------------------------------------------------------------- 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.4783 r_free = 0.4783 target = 0.209318 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 36)----------------| | r_work = 0.4524 r_free = 0.4524 target = 0.181257 restraints weight = 13909.030| |-----------------------------------------------------------------------------| r_work (start): 0.4525 rms_B_bonded: 2.84 r_work: 0.4434 rms_B_bonded: 2.90 restraints_weight: 0.5000 r_work (final): 0.4434 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7765 moved from start: 0.5537 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.029 3980 Z= 0.314 Angle : 0.848 10.958 5400 Z= 0.462 Chirality : 0.183 1.559 720 Planarity : 0.002 0.006 670 Dihedral : 5.859 17.998 580 Min Nonbonded Distance : 2.136 Molprobity Statistics. All-atom Clashscore : 6.56 Ramachandran Plot: Outliers : 0.00 % Allowed : 19.64 % Favored : 80.36 % Rotamer Outliers : 0.50 % 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.04 (0.27), residues: 560 helix: None (None), residues: 0 sheet: None (None), residues: 0 loop : -3.07 (0.20), residues: 560 =============================================================================== Job complete usr+sys time: 2278.30 seconds wall clock time: 40 minutes 52.67 seconds (2452.67 seconds total)