Starting phenix.real_space_refine on Tue Dec 31 15:45:25 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/7pkb_13467/12_2024/7pkb_13467.cif Found real_map, /net/cci-nas-00/data/ceres_data/7pkb_13467/12_2024/7pkb_13467.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.2 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/7pkb_13467/12_2024/7pkb_13467.map" default_real_map = "/net/cci-nas-00/data/ceres_data/7pkb_13467/12_2024/7pkb_13467.map" model { file = "/net/cci-nas-00/data/ceres_data/7pkb_13467/12_2024/7pkb_13467.cif" } default_model = "/net/cci-nas-00/data/ceres_data/7pkb_13467/12_2024/7pkb_13467.cif" } resolution = 3.2 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= 0.000 sd= 0.002 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 6 Type Number sf(0) Gaussians Ca 10 9.91 5 S 20 5.16 5 C 5290 2.51 5 N 1305 2.21 5 O 1550 1.98 5 H 7970 0.53 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 74 residue(s): 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5532/modules/chem_data/mon_lib" Total number of atoms: 16145 Number of models: 1 Model: "" Number of chains: 15 Chain: "E" Number of atoms: 3224 Number of conformers: 1 Conformer: "" Number of residues, atoms: 206, 3224 Classifications: {'peptide': 206} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PTRANS': 11, 'TRANS': 194} Chain: "D" Number of atoms: 3224 Number of conformers: 1 Conformer: "" Number of residues, atoms: 206, 3224 Classifications: {'peptide': 206} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PTRANS': 11, 'TRANS': 194} Chain: "C" Number of atoms: 3224 Number of conformers: 1 Conformer: "" Number of residues, atoms: 206, 3224 Classifications: {'peptide': 206} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PTRANS': 11, 'TRANS': 194} Chain: "B" Number of atoms: 3224 Number of conformers: 1 Conformer: "" Number of residues, atoms: 206, 3224 Classifications: {'peptide': 206} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PTRANS': 11, 'TRANS': 194} Chain: "A" Number of atoms: 3224 Number of conformers: 1 Conformer: "" Number of residues, atoms: 206, 3224 Classifications: {'peptide': 206} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PTRANS': 11, 'TRANS': 194} Chain: "E" Number of atoms: 2 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 2 Unusual residues: {' CA': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Chain: "D" Number of atoms: 2 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 2 Unusual residues: {' CA': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Chain: "C" Number of atoms: 2 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 2 Unusual residues: {' CA': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Chain: "B" Number of atoms: 2 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 2 Unusual residues: {' CA': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Chain: "A" Number of atoms: 2 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 2 Unusual residues: {' CA': 2} Classifications: {'undetermined': 2} Link IDs: {None: 1} Chain: "E" Number of atoms: 3 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 3 Classifications: {'water': 1} Chain: "D" Number of atoms: 3 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 3 Classifications: {'water': 1} Chain: "C" Number of atoms: 3 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 3 Classifications: {'water': 1} Chain: "B" Number of atoms: 3 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 3 Classifications: {'water': 1} Chain: "A" Number of atoms: 3 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 3 Classifications: {'water': 1} Time building chain proxies: 8.98, per 1000 atoms: 0.56 Number of scatterers: 16145 At special positions: 0 Unit cell: (113.696, 115.368, 53.504, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 6 Type Number sf(0) Ca 10 19.99 S 20 16.00 O 1550 8.00 N 1305 7.00 C 5290 6.00 H 7970 1.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=5, symmetry=0 Simple disulfide: pdb=" SG CYS E 36 " - pdb=" SG CYS E 97 " distance=2.03 Simple disulfide: pdb=" SG CYS D 36 " - pdb=" SG CYS D 97 " distance=2.04 Simple disulfide: pdb=" SG CYS C 36 " - pdb=" SG CYS C 97 " distance=2.03 Simple disulfide: pdb=" SG CYS B 36 " - pdb=" SG CYS B 97 " distance=2.04 Simple disulfide: pdb=" SG CYS A 36 " - pdb=" SG CYS A 97 " distance=2.03 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.00 Amino acid : True - 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.96 Conformation dependent library (CDL) restraints added in 1.2 seconds 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1890 Finding SS restraints... Secondary structure from input PDB file: 20 helices and 15 sheets defined 12.1% alpha, 37.9% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.78 Creating SS restraints... Processing helix chain 'E' and resid 42 through 46 removed outlier: 4.037A pdb=" N SER E 45 " --> pdb=" O GLU E 42 " (cutoff:3.500A) Processing helix chain 'E' and resid 147 through 151 Processing helix chain 'E' and resid 167 through 176 Processing helix chain 'E' and resid 186 through 190 removed outlier: 3.996A pdb=" N LEU E 190 " --> pdb=" O TRP E 187 " (cutoff:3.500A) Processing helix chain 'D' and resid 42 through 46 removed outlier: 4.019A pdb=" N SER D 45 " --> pdb=" O GLU D 42 " (cutoff:3.500A) Processing helix chain 'D' and resid 147 through 151 Processing helix chain 'D' and resid 167 through 176 Processing helix chain 'D' and resid 186 through 190 removed outlier: 4.073A pdb=" N LEU D 190 " --> pdb=" O TRP D 187 " (cutoff:3.500A) Processing helix chain 'C' and resid 42 through 46 removed outlier: 4.050A pdb=" N SER C 45 " --> pdb=" O GLU C 42 " (cutoff:3.500A) Processing helix chain 'C' and resid 147 through 151 Processing helix chain 'C' and resid 167 through 176 Processing helix chain 'C' and resid 186 through 190 removed outlier: 4.004A pdb=" N LEU C 190 " --> pdb=" O TRP C 187 " (cutoff:3.500A) Processing helix chain 'B' and resid 42 through 46 removed outlier: 4.040A pdb=" N SER B 45 " --> pdb=" O GLU B 42 " (cutoff:3.500A) Processing helix chain 'B' and resid 147 through 151 Processing helix chain 'B' and resid 167 through 176 Processing helix chain 'B' and resid 186 through 190 removed outlier: 3.987A pdb=" N LEU B 190 " --> pdb=" O TRP B 187 " (cutoff:3.500A) Processing helix chain 'A' and resid 42 through 46 removed outlier: 4.248A pdb=" N SER A 45 " --> pdb=" O GLU A 42 " (cutoff:3.500A) Processing helix chain 'A' and resid 147 through 151 Processing helix chain 'A' and resid 167 through 176 Processing helix chain 'A' and resid 186 through 190 removed outlier: 3.856A pdb=" N LEU A 190 " --> pdb=" O TRP A 187 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'E' and resid 114 through 115 Processing sheet with id=AA2, first strand: chain 'E' and resid 114 through 115 removed outlier: 3.513A pdb=" N ILE E 106 " --> pdb=" O GLU E 101 " (cutoff:3.500A) removed outlier: 6.561A pdb=" N THR E 34 " --> pdb=" O MET E 161 " (cutoff:3.500A) removed outlier: 4.358A pdb=" N MET E 161 " --> pdb=" O THR E 34 " (cutoff:3.500A) removed outlier: 6.785A pdb=" N CYS E 36 " --> pdb=" O VAL E 159 " (cutoff:3.500A) removed outlier: 4.079A pdb=" N VAL E 159 " --> pdb=" O CYS E 36 " (cutoff:3.500A) removed outlier: 4.768A pdb=" N ASN E 158 " --> pdb=" O LYS E 7 " (cutoff:3.500A) removed outlier: 10.009A pdb=" N LYS E 7 " --> pdb=" O ASN E 158 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'E' and resid 80 through 85 removed outlier: 3.507A pdb=" N TYR E 73 " --> pdb=" O PHE E 84 " (cutoff:3.500A) removed outlier: 6.962A pdb=" N TRP E 67 " --> pdb=" O SER E 50 " (cutoff:3.500A) removed outlier: 5.639A pdb=" N SER E 50 " --> pdb=" O TRP E 67 " (cutoff:3.500A) removed outlier: 3.996A pdb=" N ILE E 133 " --> pdb=" O LEU E 22 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'D' and resid 114 through 115 Processing sheet with id=AA5, first strand: chain 'D' and resid 114 through 115 removed outlier: 3.571A pdb=" N ILE D 106 " --> pdb=" O GLU D 101 " (cutoff:3.500A) removed outlier: 6.580A pdb=" N THR D 34 " --> pdb=" O MET D 161 " (cutoff:3.500A) removed outlier: 4.380A pdb=" N MET D 161 " --> pdb=" O THR D 34 " (cutoff:3.500A) removed outlier: 6.758A pdb=" N CYS D 36 " --> pdb=" O VAL D 159 " (cutoff:3.500A) removed outlier: 4.105A pdb=" N VAL D 159 " --> pdb=" O CYS D 36 " (cutoff:3.500A) removed outlier: 4.869A pdb=" N ASN D 158 " --> pdb=" O LYS D 7 " (cutoff:3.500A) removed outlier: 10.047A pdb=" N LYS D 7 " --> pdb=" O ASN D 158 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'D' and resid 80 through 85 removed outlier: 6.183A pdb=" N ASN D 61 " --> pdb=" O THR D 56 " (cutoff:3.500A) removed outlier: 7.410A pdb=" N THR D 56 " --> pdb=" O ASN D 61 " (cutoff:3.500A) removed outlier: 6.829A pdb=" N TRP D 67 " --> pdb=" O SER D 50 " (cutoff:3.500A) removed outlier: 5.596A pdb=" N SER D 50 " --> pdb=" O TRP D 67 " (cutoff:3.500A) removed outlier: 3.955A pdb=" N ILE D 133 " --> pdb=" O LEU D 22 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'C' and resid 114 through 115 Processing sheet with id=AA8, first strand: chain 'C' and resid 114 through 115 removed outlier: 3.525A pdb=" N ILE C 106 " --> pdb=" O GLU C 101 " (cutoff:3.500A) removed outlier: 6.588A pdb=" N THR C 34 " --> pdb=" O MET C 161 " (cutoff:3.500A) removed outlier: 4.400A pdb=" N MET C 161 " --> pdb=" O THR C 34 " (cutoff:3.500A) removed outlier: 6.811A pdb=" N CYS C 36 " --> pdb=" O VAL C 159 " (cutoff:3.500A) removed outlier: 4.137A pdb=" N VAL C 159 " --> pdb=" O CYS C 36 " (cutoff:3.500A) removed outlier: 4.831A pdb=" N ASN C 158 " --> pdb=" O LYS C 7 " (cutoff:3.500A) removed outlier: 10.073A pdb=" N LYS C 7 " --> pdb=" O ASN C 158 " (cutoff:3.500A) Processing sheet with id=AA9, first strand: chain 'C' and resid 80 through 85 removed outlier: 3.631A pdb=" N TYR C 73 " --> pdb=" O PHE C 84 " (cutoff:3.500A) removed outlier: 6.430A pdb=" N ASN C 61 " --> pdb=" O THR C 56 " (cutoff:3.500A) removed outlier: 7.418A pdb=" N THR C 56 " --> pdb=" O ASN C 61 " (cutoff:3.500A) removed outlier: 6.876A pdb=" N TRP C 67 " --> pdb=" O SER C 50 " (cutoff:3.500A) removed outlier: 5.560A pdb=" N SER C 50 " --> pdb=" O TRP C 67 " (cutoff:3.500A) removed outlier: 3.903A pdb=" N ILE C 133 " --> pdb=" O LEU C 22 " (cutoff:3.500A) Processing sheet with id=AB1, first strand: chain 'B' and resid 114 through 115 Processing sheet with id=AB2, first strand: chain 'B' and resid 114 through 115 removed outlier: 6.539A pdb=" N THR B 34 " --> pdb=" O MET B 161 " (cutoff:3.500A) removed outlier: 4.350A pdb=" N MET B 161 " --> pdb=" O THR B 34 " (cutoff:3.500A) removed outlier: 6.855A pdb=" N CYS B 36 " --> pdb=" O VAL B 159 " (cutoff:3.500A) removed outlier: 4.145A pdb=" N VAL B 159 " --> pdb=" O CYS B 36 " (cutoff:3.500A) removed outlier: 4.870A pdb=" N ASN B 158 " --> pdb=" O LYS B 7 " (cutoff:3.500A) removed outlier: 9.995A pdb=" N LYS B 7 " --> pdb=" O ASN B 158 " (cutoff:3.500A) Processing sheet with id=AB3, first strand: chain 'B' and resid 80 through 85 removed outlier: 3.520A pdb=" N TYR B 73 " --> pdb=" O PHE B 84 " (cutoff:3.500A) removed outlier: 6.204A pdb=" N ASN B 61 " --> pdb=" O THR B 56 " (cutoff:3.500A) removed outlier: 7.329A pdb=" N THR B 56 " --> pdb=" O ASN B 61 " (cutoff:3.500A) removed outlier: 6.878A pdb=" N TRP B 67 " --> pdb=" O SER B 50 " (cutoff:3.500A) removed outlier: 5.845A pdb=" N SER B 50 " --> pdb=" O TRP B 67 " (cutoff:3.500A) removed outlier: 3.863A pdb=" N ILE B 133 " --> pdb=" O LEU B 22 " (cutoff:3.500A) Processing sheet with id=AB4, first strand: chain 'A' and resid 114 through 115 Processing sheet with id=AB5, first strand: chain 'A' and resid 114 through 115 removed outlier: 3.551A pdb=" N ILE A 106 " --> pdb=" O GLU A 101 " (cutoff:3.500A) removed outlier: 6.606A pdb=" N THR A 34 " --> pdb=" O MET A 161 " (cutoff:3.500A) removed outlier: 4.393A pdb=" N MET A 161 " --> pdb=" O THR A 34 " (cutoff:3.500A) removed outlier: 6.756A pdb=" N CYS A 36 " --> pdb=" O VAL A 159 " (cutoff:3.500A) removed outlier: 4.115A pdb=" N VAL A 159 " --> pdb=" O CYS A 36 " (cutoff:3.500A) removed outlier: 4.900A pdb=" N ASN A 158 " --> pdb=" O LYS A 7 " (cutoff:3.500A) removed outlier: 10.127A pdb=" N LYS A 7 " --> pdb=" O ASN A 158 " (cutoff:3.500A) Processing sheet with id=AB6, first strand: chain 'A' and resid 80 through 85 removed outlier: 6.176A pdb=" N ASN A 61 " --> pdb=" O THR A 56 " (cutoff:3.500A) removed outlier: 7.334A pdb=" N THR A 56 " --> pdb=" O ASN A 61 " (cutoff:3.500A) removed outlier: 6.780A pdb=" N TRP A 67 " --> pdb=" O SER A 50 " (cutoff:3.500A) removed outlier: 5.785A pdb=" N SER A 50 " --> pdb=" O TRP A 67 " (cutoff:3.500A) removed outlier: 3.891A pdb=" N ILE A 133 " --> pdb=" O LEU A 22 " (cutoff:3.500A) 338 hydrogen bonds defined for protein. 810 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: 6.35 Time building geometry restraints manager: 4.73 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 0.84 - 1.03: 7970 1.03 - 1.22: 0 1.22 - 1.41: 3613 1.41 - 1.61: 4747 1.61 - 1.80: 30 Bond restraints: 16360 Sorted by residual: bond pdb=" CB THR D 98 " pdb=" CG2 THR D 98 " ideal model delta sigma weight residual 1.521 1.465 0.056 3.30e-02 9.18e+02 2.87e+00 bond pdb=" CB VAL D 94 " pdb=" CG1 VAL D 94 " ideal model delta sigma weight residual 1.521 1.470 0.051 3.30e-02 9.18e+02 2.36e+00 bond pdb=" CB VAL E 94 " pdb=" CG1 VAL E 94 " ideal model delta sigma weight residual 1.521 1.471 0.050 3.30e-02 9.18e+02 2.27e+00 bond pdb=" CB THR B 76 " pdb=" CG2 THR B 76 " ideal model delta sigma weight residual 1.521 1.472 0.049 3.30e-02 9.18e+02 2.16e+00 bond pdb=" CB VAL C 94 " pdb=" CG1 VAL C 94 " ideal model delta sigma weight residual 1.521 1.475 0.046 3.30e-02 9.18e+02 1.98e+00 ... (remaining 16355 not shown) Histogram of bond angle deviations from ideal: 0.00 - 5.63: 29479 5.63 - 11.27: 25 11.27 - 16.90: 1 16.90 - 22.53: 0 22.53 - 28.17: 30 Bond angle restraints: 29535 Sorted by residual: angle pdb=" CG1 VAL A 77 " pdb=" CB VAL A 77 " pdb=" HB VAL A 77 " ideal model delta sigma weight residual 108.00 79.83 28.17 3.00e+00 1.11e-01 8.81e+01 angle pdb=" CG2 VAL A 77 " pdb=" CB VAL A 77 " pdb=" HB VAL A 77 " ideal model delta sigma weight residual 108.00 80.17 27.83 3.00e+00 1.11e-01 8.61e+01 angle pdb=" CG2 ILE A 82 " pdb=" CB ILE A 82 " pdb=" HB ILE A 82 " ideal model delta sigma weight residual 109.00 81.65 27.35 3.00e+00 1.11e-01 8.31e+01 angle pdb=" CG1 VAL E 77 " pdb=" CB VAL E 77 " pdb=" HB VAL E 77 " ideal model delta sigma weight residual 108.00 80.70 27.30 3.00e+00 1.11e-01 8.28e+01 angle pdb=" CG2 THR A 173 " pdb=" CB THR A 173 " pdb=" HB THR A 173 " ideal model delta sigma weight residual 108.00 80.75 27.25 3.00e+00 1.11e-01 8.25e+01 ... (remaining 29530 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.99: 7037 17.99 - 35.98: 453 35.98 - 53.97: 179 53.97 - 71.96: 99 71.96 - 89.95: 7 Dihedral angle restraints: 7775 sinusoidal: 4100 harmonic: 3675 Sorted by residual: dihedral pdb=" CA GLU C 62 " pdb=" C GLU C 62 " pdb=" N ILE C 63 " pdb=" CA ILE C 63 " ideal model delta harmonic sigma weight residual 180.00 153.00 27.00 0 5.00e+00 4.00e-02 2.92e+01 dihedral pdb=" CA GLU B 62 " pdb=" C GLU B 62 " pdb=" N ILE B 63 " pdb=" CA ILE B 63 " ideal model delta harmonic sigma weight residual 180.00 153.26 26.74 0 5.00e+00 4.00e-02 2.86e+01 dihedral pdb=" CA GLU E 62 " pdb=" C GLU E 62 " pdb=" N ILE E 63 " pdb=" CA ILE E 63 " ideal model delta harmonic sigma weight residual 180.00 153.43 26.57 0 5.00e+00 4.00e-02 2.82e+01 ... (remaining 7772 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.349: 1225 0.349 - 0.698: 0 0.698 - 1.047: 0 1.047 - 1.395: 5 1.395 - 1.744: 5 Chirality restraints: 1235 Sorted by residual: chirality pdb=" CG LEU D 83 " pdb=" CB LEU D 83 " pdb=" CD1 LEU D 83 " pdb=" CD2 LEU D 83 " both_signs ideal model delta sigma weight residual False -2.59 -0.85 -1.74 2.00e-01 2.50e+01 7.61e+01 chirality pdb=" CG LEU C 83 " pdb=" CB LEU C 83 " pdb=" CD1 LEU C 83 " pdb=" CD2 LEU C 83 " both_signs ideal model delta sigma weight residual False -2.59 -0.90 -1.69 2.00e-01 2.50e+01 7.16e+01 chirality pdb=" CB VAL B 89 " pdb=" CA VAL B 89 " pdb=" CG1 VAL B 89 " pdb=" CG2 VAL B 89 " both_signs ideal model delta sigma weight residual False -2.63 -1.08 -1.55 2.00e-01 2.50e+01 5.98e+01 ... (remaining 1232 not shown) Planarity restraints: 2420 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB ASN E 183 " -0.032 2.00e-02 2.50e+03 4.02e-02 2.42e+01 pdb=" CG ASN E 183 " 0.025 2.00e-02 2.50e+03 pdb=" OD1 ASN E 183 " 0.002 2.00e-02 2.50e+03 pdb=" ND2 ASN E 183 " 0.065 2.00e-02 2.50e+03 pdb="HD21 ASN E 183 " 0.001 2.00e-02 2.50e+03 pdb="HD22 ASN E 183 " -0.061 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C VAL B 86 " 0.073 5.00e-02 4.00e+02 1.10e-01 1.93e+01 pdb=" N PRO B 87 " -0.190 5.00e-02 4.00e+02 pdb=" CA PRO B 87 " 0.058 5.00e-02 4.00e+02 pdb=" CD PRO B 87 " 0.059 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" CB ASN D 183 " 0.028 2.00e-02 2.50e+03 3.49e-02 1.82e+01 pdb=" CG ASN D 183 " -0.022 2.00e-02 2.50e+03 pdb=" OD1 ASN D 183 " -0.001 2.00e-02 2.50e+03 pdb=" ND2 ASN D 183 " -0.057 2.00e-02 2.50e+03 pdb="HD21 ASN D 183 " -0.001 2.00e-02 2.50e+03 pdb="HD22 ASN D 183 " 0.053 2.00e-02 2.50e+03 ... (remaining 2417 not shown) Histogram of nonbonded interaction distances: 1.53 - 2.14: 625 2.14 - 2.76: 30525 2.76 - 3.37: 45052 3.37 - 3.99: 63257 3.99 - 4.60: 96940 Nonbonded interactions: 236399 Sorted by model distance: nonbonded pdb=" HG1 THR B 41 " pdb=" O VAL B 153 " model vdw 1.531 2.450 nonbonded pdb=" HG1 THR C 41 " pdb=" O VAL C 153 " model vdw 1.541 2.450 nonbonded pdb=" OE1 GLU A 101 " pdb=" HG SER A 104 " model vdw 1.545 2.450 nonbonded pdb="HH22 ARG E 47 " pdb=" O SER E 151 " model vdw 1.556 2.450 nonbonded pdb=" HG1 THR A 41 " pdb=" O VAL A 153 " model vdw 1.557 2.450 ... (remaining 236394 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.00 Found NCS groups: ncs_group { reference = chain 'A' selection = chain 'B' selection = chain 'C' selection = chain 'D' selection = chain 'E' } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=1.00 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as 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: 3.620 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.010 Construct map_model_manager: 0.270 Extract box with map and model: 0.600 Check model and map are aligned: 0.120 Set scattering table: 0.150 Process input model: 37.140 Find NCS groups from input model: 0.460 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.010 Load rotamer database and sin/cos tables:2.830 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 45.280 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8400 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.056 8390 Z= 0.518 Angle : 0.903 14.543 11405 Z= 0.500 Chirality : 0.142 1.744 1235 Planarity : 0.006 0.110 1450 Dihedral : 13.281 89.946 2990 Min Nonbonded Distance : 1.950 Molprobity Statistics. All-atom Clashscore : 3.97 Ramachandran Plot: Outliers : 0.10 % Allowed : 5.10 % Favored : 94.80 % Rotamer: Outliers : 0.00 % Allowed : 0.44 % Favored : 99.56 % 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.98 (0.24), residues: 1020 helix: -1.59 (0.49), residues: 55 sheet: -0.19 (0.23), residues: 425 loop : -0.86 (0.25), residues: 540 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.028 0.006 TRP A 67 HIS 0.007 0.003 HIS D 95 PHE 0.018 0.003 PHE D 109 TYR 0.025 0.005 TYR D 54 ARG 0.063 0.012 ARG E 6 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 93 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 93 time to evaluate : 1.447 Fit side-chains REVERT: D 60 ASP cc_start: 0.7358 (p0) cc_final: 0.7086 (p0) REVERT: C 81 GLU cc_start: 0.7682 (tt0) cc_final: 0.7373 (tt0) REVERT: C 163 ASP cc_start: 0.7947 (p0) cc_final: 0.7715 (p0) REVERT: A 60 ASP cc_start: 0.7556 (p0) cc_final: 0.7051 (p0) outliers start: 0 outliers final: 0 residues processed: 93 average time/residue: 0.6111 time to fit residues: 72.6515 Evaluate side-chains 75 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 75 time to evaluate : 1.351 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 100 random chunks: chunk 84 optimal weight: 0.7980 chunk 75 optimal weight: 0.9980 chunk 42 optimal weight: 1.9990 chunk 25 optimal weight: 0.7980 chunk 51 optimal weight: 1.9990 chunk 40 optimal weight: 0.6980 chunk 78 optimal weight: 1.9990 chunk 30 optimal weight: 0.9990 chunk 47 optimal weight: 0.6980 chunk 58 optimal weight: 0.7980 chunk 90 optimal weight: 0.9990 overall best weight: 0.7580 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: E 1 GLN C 95 HIS B 95 HIS A 95 HIS Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8395 moved from start: 0.1278 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.036 8390 Z= 0.237 Angle : 0.724 15.242 11405 Z= 0.380 Chirality : 0.141 1.702 1235 Planarity : 0.004 0.060 1450 Dihedral : 6.208 22.592 1100 Min Nonbonded Distance : 2.156 Molprobity Statistics. All-atom Clashscore : 2.60 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.65 % Favored : 97.35 % Rotamer: Outliers : 0.33 % Allowed : 2.22 % Favored : 97.44 % 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.33 (0.25), residues: 1020 helix: 1.36 (0.60), residues: 55 sheet: 0.18 (0.23), residues: 425 loop : -0.67 (0.26), residues: 540 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.009 0.001 TRP D 162 HIS 0.004 0.002 HIS D 38 PHE 0.014 0.002 PHE E 66 TYR 0.021 0.001 TYR C 54 ARG 0.005 0.001 ARG A 188 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 97 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 94 time to evaluate : 1.184 Fit side-chains REVERT: D 60 ASP cc_start: 0.7406 (p0) cc_final: 0.7045 (p0) REVERT: C 81 GLU cc_start: 0.7783 (tt0) cc_final: 0.7508 (tt0) REVERT: A 60 ASP cc_start: 0.7583 (p0) cc_final: 0.7061 (p0) outliers start: 3 outliers final: 3 residues processed: 97 average time/residue: 0.6703 time to fit residues: 82.6835 Evaluate side-chains 79 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 76 time to evaluate : 1.378 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain D residue 134 ILE Chi-restraints excluded: chain B residue 134 ILE Chi-restraints excluded: chain A residue 134 ILE Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 100 random chunks: chunk 50 optimal weight: 1.9990 chunk 28 optimal weight: 0.9980 chunk 75 optimal weight: 1.9990 chunk 61 optimal weight: 0.9990 chunk 25 optimal weight: 1.9990 chunk 90 optimal weight: 0.4980 chunk 98 optimal weight: 1.9990 chunk 81 optimal weight: 1.9990 chunk 31 optimal weight: 0.9980 chunk 72 optimal weight: 0.7980 chunk 89 optimal weight: 0.8980 overall best weight: 0.8380 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: E 38 HIS E 95 HIS D 59 GLN D 95 HIS A 59 GLN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8407 moved from start: 0.1834 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.032 8390 Z= 0.250 Angle : 0.688 15.653 11405 Z= 0.360 Chirality : 0.142 1.701 1235 Planarity : 0.004 0.042 1450 Dihedral : 5.600 21.579 1100 Min Nonbonded Distance : 2.139 Molprobity Statistics. All-atom Clashscore : 1.98 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.33 % Favored : 96.67 % Rotamer: Outliers : 0.33 % Allowed : 3.89 % Favored : 95.78 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.30 (0.25), residues: 1020 helix: 1.60 (0.57), residues: 55 sheet: 0.14 (0.24), residues: 415 loop : -0.60 (0.27), residues: 550 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP D 162 HIS 0.004 0.002 HIS D 38 PHE 0.012 0.001 PHE A 109 TYR 0.020 0.001 TYR C 54 ARG 0.004 0.001 ARG B 188 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 86 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 83 time to evaluate : 1.352 Fit side-chains REVERT: D 81 GLU cc_start: 0.7702 (tt0) cc_final: 0.7312 (tt0) REVERT: C 81 GLU cc_start: 0.7678 (tt0) cc_final: 0.7295 (tt0) REVERT: A 60 ASP cc_start: 0.7556 (p0) cc_final: 0.7081 (p0) outliers start: 3 outliers final: 3 residues processed: 86 average time/residue: 0.6721 time to fit residues: 73.6729 Evaluate side-chains 79 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 76 time to evaluate : 1.457 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain D residue 134 ILE Chi-restraints excluded: chain B residue 134 ILE Chi-restraints excluded: chain A residue 134 ILE Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 100 random chunks: chunk 68 optimal weight: 0.9990 chunk 47 optimal weight: 1.9990 chunk 10 optimal weight: 1.9990 chunk 43 optimal weight: 0.9980 chunk 61 optimal weight: 1.9990 chunk 91 optimal weight: 0.9980 chunk 96 optimal weight: 0.3980 chunk 86 optimal weight: 1.9990 chunk 26 optimal weight: 1.9990 chunk 80 optimal weight: 2.9990 chunk 54 optimal weight: 2.9990 overall best weight: 1.0784 ------------------------------------------------------------------------------- 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.8422 moved from start: 0.2021 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.032 8390 Z= 0.288 Angle : 0.689 15.596 11405 Z= 0.360 Chirality : 0.142 1.701 1235 Planarity : 0.004 0.061 1450 Dihedral : 5.489 23.051 1100 Min Nonbonded Distance : 2.087 Molprobity Statistics. All-atom Clashscore : 2.11 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.73 % Favored : 96.27 % Rotamer: Outliers : 0.44 % Allowed : 5.00 % Favored : 94.56 % 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.50 (0.26), residues: 1020 helix: 1.47 (0.57), residues: 55 sheet: -0.03 (0.24), residues: 415 loop : -0.71 (0.27), residues: 550 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP D 162 HIS 0.007 0.002 HIS E 38 PHE 0.011 0.002 PHE A 109 TYR 0.020 0.001 TYR E 54 ARG 0.007 0.001 ARG A 188 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 80 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 4 poor density : 76 time to evaluate : 1.243 Fit side-chains REVERT: E 147 GLU cc_start: 0.7179 (mt-10) cc_final: 0.6926 (mt-10) REVERT: C 81 GLU cc_start: 0.7676 (tt0) cc_final: 0.7284 (tt0) REVERT: A 60 ASP cc_start: 0.7580 (p0) cc_final: 0.7048 (p0) outliers start: 4 outliers final: 3 residues processed: 80 average time/residue: 0.6706 time to fit residues: 67.7162 Evaluate side-chains 79 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 76 time to evaluate : 1.484 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain D residue 134 ILE Chi-restraints excluded: chain B residue 134 ILE Chi-restraints excluded: chain A residue 134 ILE Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 100 random chunks: chunk 1 optimal weight: 1.9990 chunk 71 optimal weight: 0.7980 chunk 39 optimal weight: 1.9990 chunk 82 optimal weight: 0.7980 chunk 66 optimal weight: 0.8980 chunk 0 optimal weight: 0.9980 chunk 49 optimal weight: 0.7980 chunk 86 optimal weight: 0.9980 chunk 24 optimal weight: 1.9990 chunk 32 optimal weight: 2.9990 chunk 87 optimal weight: 0.7980 overall best weight: 0.8180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 59 GLN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8413 moved from start: 0.2270 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.032 8390 Z= 0.230 Angle : 0.663 15.510 11405 Z= 0.344 Chirality : 0.141 1.687 1235 Planarity : 0.003 0.030 1450 Dihedral : 5.211 21.717 1100 Min Nonbonded Distance : 2.112 Molprobity Statistics. All-atom Clashscore : 2.11 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.02 % Favored : 95.98 % Rotamer: Outliers : 0.56 % Allowed : 5.56 % Favored : 93.89 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.50 (0.26), residues: 1020 helix: 1.57 (0.58), residues: 55 sheet: -0.03 (0.24), residues: 415 loop : -0.71 (0.27), residues: 550 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP D 162 HIS 0.006 0.002 HIS E 38 PHE 0.009 0.001 PHE D 39 TYR 0.016 0.001 TYR E 54 ARG 0.002 0.000 ARG E 116 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 83 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 5 poor density : 78 time to evaluate : 1.288 Fit side-chains REVERT: C 81 GLU cc_start: 0.7656 (tt0) cc_final: 0.7223 (tt0) REVERT: B 59 GLN cc_start: 0.8395 (tp40) cc_final: 0.8195 (tp40) REVERT: A 60 ASP cc_start: 0.7558 (p0) cc_final: 0.7093 (p0) outliers start: 5 outliers final: 3 residues processed: 83 average time/residue: 0.6796 time to fit residues: 71.4755 Evaluate side-chains 78 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 75 time to evaluate : 1.249 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain D residue 134 ILE Chi-restraints excluded: chain B residue 134 ILE Chi-restraints excluded: chain A residue 134 ILE Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 100 random chunks: chunk 19 optimal weight: 1.9990 chunk 56 optimal weight: 0.6980 chunk 23 optimal weight: 0.0980 chunk 96 optimal weight: 0.9980 chunk 80 optimal weight: 2.9990 chunk 44 optimal weight: 1.9990 chunk 8 optimal weight: 0.8980 chunk 32 optimal weight: 2.9990 chunk 50 optimal weight: 0.0980 chunk 93 optimal weight: 1.9990 chunk 10 optimal weight: 0.9990 overall best weight: 0.5580 ------------------------------------------------------------------------------- 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.8401 moved from start: 0.2456 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.032 8390 Z= 0.183 Angle : 0.643 15.380 11405 Z= 0.332 Chirality : 0.141 1.683 1235 Planarity : 0.003 0.027 1450 Dihedral : 4.945 19.107 1100 Min Nonbonded Distance : 2.148 Molprobity Statistics. All-atom Clashscore : 1.98 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.82 % Favored : 96.18 % Rotamer: Outliers : 0.33 % Allowed : 6.00 % Favored : 93.67 % 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.39 (0.26), residues: 1020 helix: 1.69 (0.59), residues: 55 sheet: 0.04 (0.25), residues: 415 loop : -0.66 (0.27), residues: 550 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP D 162 HIS 0.005 0.002 HIS E 38 PHE 0.010 0.001 PHE D 39 TYR 0.013 0.001 TYR E 54 ARG 0.002 0.000 ARG E 116 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 92 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 3 poor density : 89 time to evaluate : 1.392 Fit side-chains REVERT: E 147 GLU cc_start: 0.7190 (mt-10) cc_final: 0.6975 (mt-10) REVERT: D 16 ASP cc_start: 0.7762 (m-30) cc_final: 0.7310 (m-30) REVERT: A 60 ASP cc_start: 0.7549 (p0) cc_final: 0.7039 (p0) REVERT: A 112 ASP cc_start: 0.8776 (m-30) cc_final: 0.8504 (m-30) outliers start: 3 outliers final: 1 residues processed: 92 average time/residue: 0.6205 time to fit residues: 72.7437 Evaluate side-chains 81 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 1 poor density : 80 time to evaluate : 1.365 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 134 ILE Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 100 random chunks: chunk 55 optimal weight: 1.9990 chunk 70 optimal weight: 0.0770 chunk 54 optimal weight: 1.9990 chunk 81 optimal weight: 1.9990 chunk 96 optimal weight: 0.9980 chunk 60 optimal weight: 1.9990 chunk 58 optimal weight: 1.9990 chunk 44 optimal weight: 0.4980 chunk 59 optimal weight: 1.9990 chunk 38 optimal weight: 1.9990 chunk 57 optimal weight: 0.9990 overall best weight: 0.9142 ------------------------------------------------------------------------------- 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.8418 moved from start: 0.2500 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.031 8390 Z= 0.252 Angle : 0.658 15.602 11405 Z= 0.341 Chirality : 0.141 1.675 1235 Planarity : 0.003 0.034 1450 Dihedral : 5.039 20.764 1100 Min Nonbonded Distance : 2.110 Molprobity Statistics. All-atom Clashscore : 2.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.12 % Favored : 95.88 % Rotamer: Outliers : 0.22 % Allowed : 6.44 % Favored : 93.33 % 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.50 (0.26), residues: 1020 helix: 1.45 (0.59), residues: 55 sheet: -0.01 (0.25), residues: 415 loop : -0.73 (0.27), residues: 550 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP D 162 HIS 0.005 0.002 HIS E 38 PHE 0.010 0.001 PHE D 66 TYR 0.015 0.001 TYR E 54 ARG 0.002 0.000 ARG C 118 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 80 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 2 poor density : 78 time to evaluate : 1.381 Fit side-chains REVERT: E 147 GLU cc_start: 0.7212 (mt-10) cc_final: 0.7005 (mt-10) REVERT: D 16 ASP cc_start: 0.7746 (m-30) cc_final: 0.7283 (m-30) REVERT: A 60 ASP cc_start: 0.7510 (p0) cc_final: 0.7042 (p0) REVERT: A 112 ASP cc_start: 0.8798 (m-30) cc_final: 0.8570 (m-30) outliers start: 2 outliers final: 2 residues processed: 80 average time/residue: 0.6785 time to fit residues: 68.8860 Evaluate side-chains 75 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 2 poor density : 73 time to evaluate : 1.487 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain D residue 134 ILE Chi-restraints excluded: chain A residue 134 ILE Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 100 random chunks: chunk 29 optimal weight: 0.5980 chunk 18 optimal weight: 1.9990 chunk 61 optimal weight: 0.2980 chunk 65 optimal weight: 1.9990 chunk 47 optimal weight: 0.9990 chunk 8 optimal weight: 0.5980 chunk 75 optimal weight: 0.9990 chunk 87 optimal weight: 2.9990 chunk 92 optimal weight: 0.9990 chunk 84 optimal weight: 1.9990 chunk 89 optimal weight: 0.5980 overall best weight: 0.6182 ------------------------------------------------------------------------------- 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.8403 moved from start: 0.2617 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.032 8390 Z= 0.194 Angle : 0.640 15.468 11405 Z= 0.330 Chirality : 0.140 1.673 1235 Planarity : 0.003 0.052 1450 Dihedral : 4.878 19.238 1100 Min Nonbonded Distance : 2.130 Molprobity Statistics. All-atom Clashscore : 2.11 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.73 % Favored : 96.27 % Rotamer: Outliers : 0.11 % Allowed : 6.78 % Favored : 93.11 % 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.43 (0.26), residues: 1020 helix: 1.62 (0.59), residues: 55 sheet: 0.05 (0.25), residues: 415 loop : -0.70 (0.27), residues: 550 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP D 162 HIS 0.005 0.002 HIS E 38 PHE 0.011 0.001 PHE D 39 TYR 0.013 0.001 TYR E 54 ARG 0.005 0.000 ARG D 188 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 79 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 1 poor density : 78 time to evaluate : 1.335 Fit side-chains REVERT: D 16 ASP cc_start: 0.7749 (m-30) cc_final: 0.7289 (m-30) REVERT: C 81 GLU cc_start: 0.7548 (tt0) cc_final: 0.7263 (tt0) REVERT: A 60 ASP cc_start: 0.7476 (p0) cc_final: 0.7005 (p0) REVERT: A 112 ASP cc_start: 0.8790 (m-30) cc_final: 0.8564 (m-30) outliers start: 1 outliers final: 1 residues processed: 79 average time/residue: 0.7474 time to fit residues: 74.3197 Evaluate side-chains 78 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 1 poor density : 77 time to evaluate : 1.379 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 134 ILE Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 100 random chunks: chunk 92 optimal weight: 1.9990 chunk 54 optimal weight: 0.6980 chunk 39 optimal weight: 0.7980 chunk 70 optimal weight: 0.9980 chunk 27 optimal weight: 1.9990 chunk 81 optimal weight: 1.9990 chunk 84 optimal weight: 1.9990 chunk 89 optimal weight: 0.5980 chunk 58 optimal weight: 1.9990 chunk 94 optimal weight: 1.9990 chunk 57 optimal weight: 0.3980 overall best weight: 0.6980 ------------------------------------------------------------------------------- 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.8408 moved from start: 0.2655 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.031 8390 Z= 0.209 Angle : 0.643 15.533 11405 Z= 0.332 Chirality : 0.140 1.674 1235 Planarity : 0.003 0.048 1450 Dihedral : 4.854 19.124 1100 Min Nonbonded Distance : 2.123 Molprobity Statistics. All-atom Clashscore : 2.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.82 % Favored : 96.18 % Rotamer: Outliers : 0.11 % Allowed : 7.00 % Favored : 92.89 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -0.47 (0.26), residues: 1020 helix: 1.56 (0.59), residues: 55 sheet: 0.02 (0.25), residues: 415 loop : -0.73 (0.27), residues: 550 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP D 162 HIS 0.005 0.002 HIS E 38 PHE 0.010 0.001 PHE D 39 TYR 0.014 0.001 TYR E 54 ARG 0.005 0.000 ARG D 188 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 78 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 1 poor density : 77 time to evaluate : 1.309 Fit side-chains REVERT: D 16 ASP cc_start: 0.7739 (m-30) cc_final: 0.7371 (m-30) REVERT: C 81 GLU cc_start: 0.7531 (tt0) cc_final: 0.7231 (tt0) REVERT: A 60 ASP cc_start: 0.7512 (p0) cc_final: 0.7043 (p0) REVERT: A 112 ASP cc_start: 0.8795 (m-30) cc_final: 0.8564 (m-30) outliers start: 1 outliers final: 1 residues processed: 78 average time/residue: 0.6527 time to fit residues: 64.8457 Evaluate side-chains 77 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 1 poor density : 76 time to evaluate : 1.557 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 134 ILE Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 100 random chunks: chunk 45 optimal weight: 0.9990 chunk 66 optimal weight: 1.9990 chunk 99 optimal weight: 1.9990 chunk 91 optimal weight: 0.1980 chunk 79 optimal weight: 0.9990 chunk 8 optimal weight: 0.9980 chunk 61 optimal weight: 0.5980 chunk 48 optimal weight: 0.9990 chunk 63 optimal weight: 0.5980 chunk 84 optimal weight: 1.9990 chunk 24 optimal weight: 0.9980 overall best weight: 0.6780 ------------------------------------------------------------------------------- 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.8405 moved from start: 0.2713 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.032 8390 Z= 0.202 Angle : 0.640 15.533 11405 Z= 0.330 Chirality : 0.140 1.672 1235 Planarity : 0.003 0.054 1450 Dihedral : 4.820 19.203 1100 Min Nonbonded Distance : 2.126 Molprobity Statistics. All-atom Clashscore : 2.05 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.43 % Favored : 96.57 % Rotamer: Outliers : 0.11 % Allowed : 7.22 % Favored : 92.67 % 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.47 (0.26), residues: 1020 helix: 1.62 (0.59), residues: 55 sheet: 0.03 (0.25), residues: 415 loop : -0.74 (0.27), residues: 550 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP D 162 HIS 0.005 0.002 HIS E 38 PHE 0.010 0.001 PHE D 39 TYR 0.013 0.001 TYR E 54 ARG 0.005 0.000 ARG D 188 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 2040 Ramachandran restraints generated. 1020 Oldfield, 0 Emsley, 1020 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 80 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 1 poor density : 79 time to evaluate : 1.312 Fit side-chains REVERT: D 16 ASP cc_start: 0.7719 (m-30) cc_final: 0.7345 (m-30) REVERT: C 81 GLU cc_start: 0.7519 (tt0) cc_final: 0.7222 (tt0) REVERT: A 60 ASP cc_start: 0.7495 (p0) cc_final: 0.7031 (p0) REVERT: A 112 ASP cc_start: 0.8789 (m-30) cc_final: 0.8564 (m-30) outliers start: 1 outliers final: 1 residues processed: 80 average time/residue: 0.6433 time to fit residues: 66.3742 Evaluate side-chains 79 residues out of total 900 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 1 poor density : 78 time to evaluate : 1.728 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 134 ILE Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 100 random chunks: chunk 73 optimal weight: 0.8980 chunk 11 optimal weight: 0.7980 chunk 22 optimal weight: 0.7980 chunk 79 optimal weight: 0.8980 chunk 33 optimal weight: 0.9990 chunk 81 optimal weight: 0.5980 chunk 10 optimal weight: 0.5980 chunk 14 optimal weight: 2.9990 chunk 69 optimal weight: 0.8980 chunk 4 optimal weight: 1.9990 chunk 57 optimal weight: 0.2980 overall best weight: 0.6180 ------------------------------------------------------------------------------- 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.2838 r_free = 0.2838 target = 0.079674 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 34)----------------| | r_work = 0.2565 r_free = 0.2565 target = 0.064320 restraints weight = 38704.461| |-----------------------------------------------------------------------------| r_work (start): 0.2599 rms_B_bonded: 2.53 r_work: 0.2481 rms_B_bonded: 2.52 restraints_weight: 0.5000 r_work: 0.2337 rms_B_bonded: 4.14 restraints_weight: 0.2500 r_work (final): 0.2337 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8628 moved from start: 0.2759 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.032 8390 Z= 0.192 Angle : 0.637 15.516 11405 Z= 0.328 Chirality : 0.140 1.672 1235 Planarity : 0.003 0.053 1450 Dihedral : 4.756 19.229 1100 Min Nonbonded Distance : 2.132 Molprobity Statistics. All-atom Clashscore : 1.98 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.63 % Favored : 96.37 % Rotamer: Outliers : 0.22 % Allowed : 7.33 % Favored : 92.44 % 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.44 (0.26), residues: 1020 helix: 1.67 (0.59), residues: 55 sheet: 0.06 (0.25), residues: 415 loop : -0.74 (0.27), residues: 550 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.010 0.001 TRP C 162 HIS 0.005 0.002 HIS E 38 PHE 0.011 0.001 PHE D 39 TYR 0.012 0.001 TYR E 54 ARG 0.006 0.000 ARG D 188 Origin is already at (0, 0, 0), no shifts will be applied =============================================================================== Job complete usr+sys time: 4167.13 seconds wall clock time: 75 minutes 30.33 seconds (4530.33 seconds total)