Starting phenix.real_space_refine (version: dev) on Fri May 13 14:09:52 2022 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ktq_23026/05_2022/7ktq_23026.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ktq_23026/05_2022/7ktq_23026.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.3 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ktq_23026/05_2022/7ktq_23026.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ktq_23026/05_2022/7ktq_23026.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ktq_23026/05_2022/7ktq_23026.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7ktq_23026/05_2022/7ktq_23026.pdb" } resolution = 3.3 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.000 sd= 0.001 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Convert atoms to be neutral Process input model Symmetric amino acids flipped None Time to flip residues: 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-4597/modules/chem_data/mon_lib" Total number of atoms: 12745 Number of models: 1 Model: "" Number of chains: 10 Chain: "A" Number of atoms: 817 Number of conformers: 1 Conformer: "" Number of residues, atoms: 99, 817 Classifications: {'peptide': 99} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 4, 'TRANS': 94} Chain: "B" Number of atoms: 627 Number of conformers: 1 Conformer: "" Number of residues, atoms: 79, 627 Classifications: {'peptide': 79} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 1, 'TRANS': 77} Chain: "C" Number of atoms: 809 Number of conformers: 1 Conformer: "" Number of residues, atoms: 105, 809 Classifications: {'peptide': 105} Link IDs: {'PTRANS': 5, 'TRANS': 99} Chain: "D" Number of atoms: 745 Number of conformers: 1 Conformer: "" Number of residues, atoms: 95, 745 Classifications: {'peptide': 95} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 2, 'TRANS': 92} Chain: "E" Number of atoms: 791 Number of conformers: 1 Conformer: "" Number of residues, atoms: 96, 791 Classifications: {'peptide': 96} Link IDs: {'PTRANS': 3, 'TRANS': 92} Chain: "F" Number of atoms: 627 Number of conformers: 1 Conformer: "" Number of residues, atoms: 79, 627 Classifications: {'peptide': 79} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 1, 'TRANS': 77} Chain: "G" Number of atoms: 800 Number of conformers: 1 Conformer: "" Number of residues, atoms: 104, 800 Classifications: {'peptide': 104} Link IDs: {'PTRANS': 5, 'TRANS': 98} Chain: "H" Number of atoms: 745 Number of conformers: 1 Conformer: "" Number of residues, atoms: 95, 745 Classifications: {'peptide': 95} Modifications used: {'COO': 1} Link IDs: {'PTRANS': 2, 'TRANS': 92} Chain: "I" Number of atoms: 3343 Number of conformers: 1 Conformer: "" Number of residues, atoms: 164, 3343 Classifications: {'DNA': 164} Link IDs: {'rna3p': 163} Chain: "J" Number of atoms: 3441 Number of conformers: 1 Conformer: "" Number of residues, atoms: 167, 3441 Classifications: {'DNA': 167} Link IDs: {'rna3p': 166} Time building chain proxies: 6.95, per 1000 atoms: 0.55 Number of scatterers: 12745 At special positions: 0 Unit cell: (76.95, 121.5, 147.15, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) S 12 16.00 P 331 15.00 O 3045 8.00 N 2384 7.00 C 6973 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.50 Amino acid : False - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 5.60 Conformation dependent library (CDL) restraints added in 945.9 milliseconds 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1408 Finding SS restraints... Secondary structure from input PDB file: 32 helices and 10 sheets defined 71.0% alpha, 3.2% beta 117 base pairs and 288 stacking pairs defined. Time for finding SS restraints: 5.11 Creating SS restraints... Processing helix chain 'A' and resid 44 through 57 Processing helix chain 'A' and resid 63 through 77 removed outlier: 3.709A pdb=" N PHE A 67 " --> pdb=" O ARG A 63 " (cutoff:3.500A) removed outlier: 3.884A pdb=" N ASP A 77 " --> pdb=" O GLU A 73 " (cutoff:3.500A) Processing helix chain 'A' and resid 85 through 114 Processing helix chain 'A' and resid 120 through 131 Processing helix chain 'B' and resid 30 through 41 Processing helix chain 'B' and resid 49 through 76 removed outlier: 3.551A pdb=" N GLU B 53 " --> pdb=" O LEU B 49 " (cutoff:3.500A) Processing helix chain 'B' and resid 82 through 93 removed outlier: 3.559A pdb=" N VAL B 86 " --> pdb=" O THR B 82 " (cutoff:3.500A) Processing helix chain 'C' and resid 16 through 21 Processing helix chain 'C' and resid 26 through 37 Processing helix chain 'C' and resid 45 through 73 removed outlier: 3.652A pdb=" N VAL C 49 " --> pdb=" O ALA C 45 " (cutoff:3.500A) Processing helix chain 'C' and resid 79 through 90 removed outlier: 3.710A pdb=" N ASP C 90 " --> pdb=" O ALA C 86 " (cutoff:3.500A) Processing helix chain 'C' and resid 90 through 97 Processing helix chain 'D' and resid 34 through 46 Processing helix chain 'D' and resid 52 through 81 Processing helix chain 'D' and resid 87 through 99 Processing helix chain 'D' and resid 101 through 120 Processing helix chain 'E' and resid 44 through 57 Processing helix chain 'E' and resid 63 through 77 removed outlier: 3.651A pdb=" N PHE E 67 " --> pdb=" O ARG E 63 " (cutoff:3.500A) removed outlier: 3.968A pdb=" N ASP E 77 " --> pdb=" O GLU E 73 " (cutoff:3.500A) Processing helix chain 'E' and resid 85 through 114 Processing helix chain 'E' and resid 120 through 132 Processing helix chain 'F' and resid 30 through 41 Processing helix chain 'F' and resid 49 through 76 removed outlier: 3.552A pdb=" N GLU F 53 " --> pdb=" O LEU F 49 " (cutoff:3.500A) Processing helix chain 'F' and resid 82 through 93 removed outlier: 3.559A pdb=" N VAL F 86 " --> pdb=" O THR F 82 " (cutoff:3.500A) Processing helix chain 'G' and resid 16 through 22 Processing helix chain 'G' and resid 26 through 37 Processing helix chain 'G' and resid 45 through 73 removed outlier: 3.674A pdb=" N VAL G 49 " --> pdb=" O ALA G 45 " (cutoff:3.500A) Processing helix chain 'G' and resid 79 through 90 removed outlier: 3.715A pdb=" N ASP G 90 " --> pdb=" O ALA G 86 " (cutoff:3.500A) Processing helix chain 'G' and resid 90 through 97 Processing helix chain 'H' and resid 34 through 46 Processing helix chain 'H' and resid 52 through 81 Processing helix chain 'H' and resid 87 through 99 Processing helix chain 'H' and resid 101 through 120 Processing sheet with id=AA1, first strand: chain 'A' and resid 83 through 84 removed outlier: 6.904A pdb=" N ARG A 83 " --> pdb=" O VAL B 81 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA1 Processing sheet with id=AA2, first strand: chain 'A' and resid 118 through 119 Processing sheet with id=AA3, first strand: chain 'B' and resid 96 through 98 removed outlier: 6.310A pdb=" N THR B 96 " --> pdb=" O THR G 101 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA3 Processing sheet with id=AA4, first strand: chain 'C' and resid 42 through 43 removed outlier: 7.173A pdb=" N ARG C 42 " --> pdb=" O ILE D 86 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA4 Processing sheet with id=AA5, first strand: chain 'C' and resid 77 through 78 Processing sheet with id=AA6, first strand: chain 'C' and resid 100 through 102 Processing sheet with id=AA7, first strand: chain 'E' and resid 83 through 84 removed outlier: 6.964A pdb=" N ARG E 83 " --> pdb=" O VAL F 81 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA7 Processing sheet with id=AA8, first strand: chain 'E' and resid 118 through 119 Processing sheet with id=AA9, first strand: chain 'G' and resid 42 through 43 removed outlier: 7.192A pdb=" N ARG G 42 " --> pdb=" O ILE H 86 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA9 Processing sheet with id=AB1, first strand: chain 'G' and resid 77 through 78 400 hydrogen bonds defined for protein. 1200 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 300 hydrogen bonds 596 hydrogen bond angles 0 basepair planarities 117 basepair parallelities 288 stacking parallelities Total time for adding SS restraints: 4.58 Time building geometry restraints manager: 5.91 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.21 - 1.33: 2203 1.33 - 1.45: 4715 1.45 - 1.57: 6049 1.57 - 1.69: 661 1.69 - 1.80: 22 Bond restraints: 13650 Sorted by residual: bond pdb=" O5' DA J 83 " pdb=" C5' DA J 83 " ideal model delta sigma weight residual 1.423 1.351 0.072 3.00e-02 1.11e+03 5.72e+00 bond pdb=" C3' DT J 106 " pdb=" O3' DT J 106 " ideal model delta sigma weight residual 1.422 1.487 -0.065 3.00e-02 1.11e+03 4.62e+00 bond pdb=" O5' DC I 64 " pdb=" C5' DC I 64 " ideal model delta sigma weight residual 1.423 1.360 0.063 3.00e-02 1.11e+03 4.41e+00 bond pdb=" O5' DG J 102 " pdb=" C5' DG J 102 " ideal model delta sigma weight residual 1.423 1.363 0.060 3.00e-02 1.11e+03 3.99e+00 bond pdb=" O5' DT J 17 " pdb=" C5' DT J 17 " ideal model delta sigma weight residual 1.423 1.366 0.057 3.00e-02 1.11e+03 3.56e+00 ... (remaining 13645 not shown) Histogram of bond angle deviations from ideal: 98.17 - 105.06: 1386 105.06 - 111.96: 7302 111.96 - 118.85: 3852 118.85 - 125.74: 6253 125.74 - 132.63: 1072 Bond angle restraints: 19865 Sorted by residual: angle pdb=" C ARG E 40 " pdb=" N TYR E 41 " pdb=" CA TYR E 41 " ideal model delta sigma weight residual 122.03 113.80 8.23 2.01e+00 2.48e-01 1.68e+01 angle pdb=" N ILE F 26 " pdb=" CA ILE F 26 " pdb=" C ILE F 26 " ideal model delta sigma weight residual 112.50 107.10 5.40 1.39e+00 5.18e-01 1.51e+01 angle pdb=" C4' DC J 82 " pdb=" C3' DC J 82 " pdb=" O3' DC J 82 " ideal model delta sigma weight residual 110.00 104.27 5.73 1.50e+00 4.44e-01 1.46e+01 angle pdb=" C3' DT J 106 " pdb=" O3' DT J 106 " pdb=" P DG J 107 " ideal model delta sigma weight residual 120.20 125.30 -5.10 1.50e+00 4.44e-01 1.16e+01 angle pdb=" N ILE B 26 " pdb=" CA ILE B 26 " pdb=" C ILE B 26 " ideal model delta sigma weight residual 112.50 107.99 4.51 1.39e+00 5.18e-01 1.05e+01 ... (remaining 19860 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 35.95: 5440 35.95 - 71.89: 1588 71.89 - 107.84: 21 107.84 - 143.79: 0 143.79 - 179.74: 5 Dihedral angle restraints: 7054 sinusoidal: 4872 harmonic: 2182 Sorted by residual: dihedral pdb=" C4' DC I 3 " pdb=" C3' DC I 3 " pdb=" O3' DC I 3 " pdb=" P DC I 4 " ideal model delta sinusoidal sigma weight residual -140.00 39.74 -179.74 1 3.50e+01 8.16e-04 1.57e+01 dihedral pdb=" C4' DG I 104 " pdb=" C3' DG I 104 " pdb=" O3' DG I 104 " pdb=" P DC I 105 " ideal model delta sinusoidal sigma weight residual 220.00 40.45 179.55 1 3.50e+01 8.16e-04 1.57e+01 dihedral pdb=" C4' DT J 106 " pdb=" C3' DT J 106 " pdb=" O3' DT J 106 " pdb=" P DG J 107 " ideal model delta sinusoidal sigma weight residual 220.00 42.30 177.70 1 3.50e+01 8.16e-04 1.57e+01 ... (remaining 7051 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.059: 1790 0.059 - 0.118: 418 0.118 - 0.177: 51 0.177 - 0.236: 1 0.236 - 0.295: 2 Chirality restraints: 2262 Sorted by residual: chirality pdb=" CG LEU G 115 " pdb=" CB LEU G 115 " pdb=" CD1 LEU G 115 " pdb=" CD2 LEU G 115 " both_signs ideal model delta sigma weight residual False -2.59 -2.29 -0.30 2.00e-01 2.50e+01 2.18e+00 chirality pdb=" C1' DG J 66 " pdb=" O4' DG J 66 " pdb=" C2' DG J 66 " pdb=" N9 DG J 66 " both_signs ideal model delta sigma weight residual False 2.42 2.18 0.24 2.00e-01 2.50e+01 1.40e+00 chirality pdb=" C4' DA J 83 " pdb=" C5' DA J 83 " pdb=" O4' DA J 83 " pdb=" C3' DA J 83 " both_signs ideal model delta sigma weight residual False -2.53 -2.31 -0.22 2.00e-01 2.50e+01 1.21e+00 ... (remaining 2259 not shown) Planarity restraints: 1368 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C1' DC I 106 " 0.005 2.00e-02 2.50e+03 3.37e-02 2.55e+01 pdb=" N1 DC I 106 " -0.025 2.00e-02 2.50e+03 pdb=" C2 DC I 106 " 0.083 2.00e-02 2.50e+03 pdb=" O2 DC I 106 " -0.045 2.00e-02 2.50e+03 pdb=" N3 DC I 106 " -0.009 2.00e-02 2.50e+03 pdb=" C4 DC I 106 " 0.016 2.00e-02 2.50e+03 pdb=" N4 DC I 106 " -0.014 2.00e-02 2.50e+03 pdb=" C5 DC I 106 " -0.004 2.00e-02 2.50e+03 pdb=" C6 DC I 106 " -0.007 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C1' DC J 87 " -0.043 2.00e-02 2.50e+03 2.87e-02 1.86e+01 pdb=" N1 DC J 87 " 0.068 2.00e-02 2.50e+03 pdb=" C2 DC J 87 " -0.027 2.00e-02 2.50e+03 pdb=" O2 DC J 87 " 0.013 2.00e-02 2.50e+03 pdb=" N3 DC J 87 " 0.002 2.00e-02 2.50e+03 pdb=" C4 DC J 87 " -0.005 2.00e-02 2.50e+03 pdb=" N4 DC J 87 " -0.007 2.00e-02 2.50e+03 pdb=" C5 DC J 87 " -0.005 2.00e-02 2.50e+03 pdb=" C6 DC J 87 " 0.004 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA HIS H 106 " -0.021 2.00e-02 2.50e+03 4.21e-02 1.77e+01 pdb=" C HIS H 106 " 0.073 2.00e-02 2.50e+03 pdb=" O HIS H 106 " -0.027 2.00e-02 2.50e+03 pdb=" N ALA H 107 " -0.025 2.00e-02 2.50e+03 ... (remaining 1365 not shown) Histogram of nonbonded interaction distances: 2.20 - 2.74: 1483 2.74 - 3.28: 10850 3.28 - 3.82: 24023 3.82 - 4.36: 29764 4.36 - 4.90: 41980 Nonbonded interactions: 108100 Sorted by model distance: nonbonded pdb=" OE1 GLU B 74 " pdb=" NH2 ARG D 96 " model vdw 2.201 2.520 nonbonded pdb=" OD2 ASP D 65 " pdb=" OH TYR F 98 " model vdw 2.248 2.440 nonbonded pdb=" C2' DA J 20 " pdb=" C7 DT J 21 " model vdw 2.274 3.860 nonbonded pdb=" OE1 GLU F 74 " pdb=" NH2 ARG H 96 " model vdw 2.295 2.520 nonbonded pdb=" O6 DG I 44 " pdb=" N4 DC J 124 " model vdw 2.313 2.520 ... (remaining 108095 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Found NCS groups: ncs_group { reference = (chain 'A' and resid 38 through 133) selection = chain 'E' } ncs_group { reference = chain 'B' selection = chain 'F' } ncs_group { reference = (chain 'C' and resid 14 through 117) selection = chain 'G' } ncs_group { reference = chain 'D' selection = chain 'H' } 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: 5 Type Number sf(0) Gaussians P 331 5.49 5 S 12 5.16 5 C 6973 2.51 5 N 2384 2.21 5 O 3045 1.98 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: 1.930 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.010 Extract box with map and model: 6.860 Check model and map are aligned: 0.210 Convert atoms to be neutral: 0.110 Process input model: 38.640 Find NCS groups from input model: 0.460 Set up NCS constraints: 0.050 Set refine NCS operators: 0.000 Set scattering table: 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.640 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 50.920 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8455 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.072 13650 Z= 0.456 Angle : 0.840 8.228 19865 Z= 0.523 Chirality : 0.050 0.295 2262 Planarity : 0.009 0.043 1368 Dihedral : 29.359 179.736 5646 Min Nonbonded Distance : 2.201 Molprobity Statistics. All-atom Clashscore : 6.42 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.49 % Favored : 98.51 % Rotamer Outliers : 0.16 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: -2.43 (0.24), residues: 736 helix: -1.35 (0.18), residues: 532 sheet: None (None), residues: 0 loop : -1.86 (0.34), residues: 204 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 189 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 188 time to evaluate : 0.979 Fit side-chains outliers start: 1 outliers final: 0 residues processed: 189 average time/residue: 0.4179 time to fit residues: 99.9594 Evaluate side-chains 126 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 126 time to evaluate : 1.402 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.3325 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 102 random chunks: chunk 86 optimal weight: 10.0000 chunk 77 optimal weight: 2.9990 chunk 42 optimal weight: 0.9990 chunk 26 optimal weight: 0.7980 chunk 52 optimal weight: 0.9990 chunk 41 optimal weight: 0.8980 chunk 79 optimal weight: 2.9990 chunk 30 optimal weight: 1.9990 chunk 48 optimal weight: 0.5980 chunk 59 optimal weight: 0.6980 chunk 92 optimal weight: 2.9990 overall best weight: 0.7982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 68 GLN B 25 ASN B 93 GLN D 81 ASN F 27 GLN G 38 ASN G 110 ASN Total number of N/Q/H flips: 7 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8467 moved from start: 0.1884 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.042 13650 Z= 0.215 Angle : 0.557 5.429 19865 Z= 0.335 Chirality : 0.033 0.165 2262 Planarity : 0.005 0.044 1368 Dihedral : 32.936 178.820 4164 Min Nonbonded Distance : 2.194 Molprobity Statistics. All-atom Clashscore : 5.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.36 % Favored : 98.64 % Rotamer Outliers : 1.76 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.12 (0.28), residues: 736 helix: 1.46 (0.21), residues: 540 sheet: None (None), residues: 0 loop : -1.25 (0.36), residues: 196 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 144 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 133 time to evaluate : 0.960 Fit side-chains outliers start: 11 outliers final: 7 residues processed: 137 average time/residue: 0.3209 time to fit residues: 60.3959 Evaluate side-chains 130 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 123 time to evaluate : 0.987 Switching outliers to nearest non-outliers outliers start: 7 outliers final: 0 residues processed: 7 average time/residue: 0.1083 time to fit residues: 2.6080 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 102 random chunks: chunk 51 optimal weight: 5.9990 chunk 28 optimal weight: 0.5980 chunk 77 optimal weight: 2.9990 chunk 63 optimal weight: 0.0570 chunk 25 optimal weight: 0.9980 chunk 92 optimal weight: 3.9990 chunk 100 optimal weight: 5.9990 chunk 82 optimal weight: 3.9990 chunk 31 optimal weight: 0.9990 chunk 74 optimal weight: 3.9990 chunk 91 optimal weight: 2.9990 overall best weight: 1.1302 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN F 93 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8518 moved from start: 0.2546 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.041 13650 Z= 0.257 Angle : 0.557 5.508 19865 Z= 0.334 Chirality : 0.034 0.140 2262 Planarity : 0.004 0.038 1368 Dihedral : 32.516 178.365 4164 Min Nonbonded Distance : 2.155 Molprobity Statistics. All-atom Clashscore : 5.45 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.04 % Favored : 97.96 % Rotamer Outliers : 0.48 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 2.13 (0.30), residues: 736 helix: 2.25 (0.21), residues: 542 sheet: None (None), residues: 0 loop : -1.07 (0.39), residues: 194 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 131 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 128 time to evaluate : 0.938 Fit side-chains outliers start: 3 outliers final: 0 residues processed: 129 average time/residue: 0.3204 time to fit residues: 55.0421 Evaluate side-chains 120 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 120 time to evaluate : 0.946 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 0.9951 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 102 random chunks: chunk 69 optimal weight: 0.9990 chunk 48 optimal weight: 0.9980 chunk 10 optimal weight: 0.9980 chunk 44 optimal weight: 0.5980 chunk 62 optimal weight: 0.9990 chunk 93 optimal weight: 0.9990 chunk 98 optimal weight: 3.9990 chunk 88 optimal weight: 4.9990 chunk 26 optimal weight: 5.9990 chunk 82 optimal weight: 3.9990 chunk 55 optimal weight: 0.6980 overall best weight: 0.8582 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN D 46 HIS F 93 GLN G 38 ASN G 73 ASN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8505 moved from start: 0.2725 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.042 13650 Z= 0.214 Angle : 0.527 5.381 19865 Z= 0.319 Chirality : 0.032 0.139 2262 Planarity : 0.004 0.037 1368 Dihedral : 32.415 178.413 4164 Min Nonbonded Distance : 2.145 Molprobity Statistics. All-atom Clashscore : 6.20 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.63 % Favored : 98.37 % Rotamer Outliers : 1.28 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 2.66 (0.30), residues: 736 helix: 2.66 (0.21), residues: 543 sheet: None (None), residues: 0 loop : -0.98 (0.39), residues: 193 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 136 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 128 time to evaluate : 1.040 Fit side-chains outliers start: 8 outliers final: 6 residues processed: 128 average time/residue: 0.2798 time to fit residues: 48.7995 Evaluate side-chains 127 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 121 time to evaluate : 0.847 Switching outliers to nearest non-outliers outliers start: 6 outliers final: 0 residues processed: 6 average time/residue: 0.1493 time to fit residues: 2.4259 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 102 random chunks: chunk 1 optimal weight: 1.9990 chunk 73 optimal weight: 5.9990 chunk 40 optimal weight: 0.5980 chunk 84 optimal weight: 6.9990 chunk 68 optimal weight: 0.9990 chunk 0 optimal weight: 3.9990 chunk 50 optimal weight: 0.5980 chunk 88 optimal weight: 3.9990 chunk 24 optimal weight: 1.9990 chunk 33 optimal weight: 0.0770 chunk 19 optimal weight: 0.5980 overall best weight: 0.5740 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN F 93 GLN G 38 ASN G 73 ASN H 81 ASN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8494 moved from start: 0.2774 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.042 13650 Z= 0.175 Angle : 0.518 6.109 19865 Z= 0.315 Chirality : 0.032 0.134 2262 Planarity : 0.004 0.036 1368 Dihedral : 32.394 178.369 4164 Min Nonbonded Distance : 2.144 Molprobity Statistics. All-atom Clashscore : 6.51 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.90 % Favored : 98.10 % Rotamer Outliers : 0.48 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 2.88 (0.31), residues: 736 helix: 2.85 (0.22), residues: 542 sheet: None (None), residues: 0 loop : -0.97 (0.39), residues: 194 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 131 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 128 time to evaluate : 0.868 Fit side-chains outliers start: 3 outliers final: 0 residues processed: 129 average time/residue: 0.3081 time to fit residues: 53.7267 Evaluate side-chains 124 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 124 time to evaluate : 0.875 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.0900 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 102 random chunks: chunk 57 optimal weight: 2.9990 chunk 24 optimal weight: 1.9990 chunk 98 optimal weight: 3.9990 chunk 81 optimal weight: 3.9990 chunk 45 optimal weight: 0.6980 chunk 8 optimal weight: 0.6980 chunk 32 optimal weight: 0.9980 chunk 51 optimal weight: 5.9990 chunk 95 optimal weight: 2.9990 chunk 11 optimal weight: 0.7980 chunk 56 optimal weight: 0.8980 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 93 GLN D 92 GLN F 93 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8500 moved from start: 0.2966 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.041 13650 Z= 0.205 Angle : 0.521 6.436 19865 Z= 0.315 Chirality : 0.032 0.132 2262 Planarity : 0.004 0.036 1368 Dihedral : 32.287 178.533 4164 Min Nonbonded Distance : 2.153 Molprobity Statistics. All-atom Clashscore : 6.20 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.04 % Favored : 97.96 % Rotamer Outliers : 1.28 % 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.03 (0.31), residues: 736 helix: 2.96 (0.22), residues: 542 sheet: None (None), residues: 0 loop : -0.93 (0.40), residues: 194 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 135 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 127 time to evaluate : 0.840 Fit side-chains outliers start: 8 outliers final: 4 residues processed: 128 average time/residue: 0.3046 time to fit residues: 52.6192 Evaluate side-chains 126 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 122 time to evaluate : 0.878 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.1154 time to fit residues: 1.9409 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 102 random chunks: chunk 72 optimal weight: 4.9990 chunk 55 optimal weight: 2.9990 chunk 83 optimal weight: 4.9990 chunk 98 optimal weight: 3.9990 chunk 61 optimal weight: 0.7980 chunk 59 optimal weight: 0.8980 chunk 45 optimal weight: 0.5980 chunk 60 optimal weight: 0.7980 chunk 39 optimal weight: 1.9990 chunk 58 optimal weight: 0.6980 chunk 29 optimal weight: 0.6980 overall best weight: 0.7180 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN F 93 GLN G 38 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8496 moved from start: 0.2988 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.044 13650 Z= 0.197 Angle : 0.518 5.569 19865 Z= 0.313 Chirality : 0.031 0.134 2262 Planarity : 0.004 0.035 1368 Dihedral : 32.273 178.540 4164 Min Nonbonded Distance : 2.151 Molprobity Statistics. All-atom Clashscore : 6.20 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.17 % Favored : 97.83 % Rotamer Outliers : 0.64 % 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.07 (0.31), residues: 736 helix: 3.00 (0.22), residues: 542 sheet: None (None), residues: 0 loop : -0.95 (0.41), residues: 194 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 133 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 129 time to evaluate : 0.962 Fit side-chains outliers start: 4 outliers final: 4 residues processed: 129 average time/residue: 0.2991 time to fit residues: 51.7842 Evaluate side-chains 128 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 124 time to evaluate : 0.967 Switching outliers to nearest non-outliers outliers start: 4 outliers final: 0 residues processed: 4 average time/residue: 0.2550 time to fit residues: 2.5628 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 102 random chunks: chunk 19 optimal weight: 0.6980 chunk 62 optimal weight: 0.6980 chunk 66 optimal weight: 4.9990 chunk 48 optimal weight: 0.5980 chunk 9 optimal weight: 1.9990 chunk 77 optimal weight: 2.9990 chunk 89 optimal weight: 4.9990 chunk 94 optimal weight: 2.9990 chunk 85 optimal weight: 5.9990 chunk 91 optimal weight: 2.9990 chunk 55 optimal weight: 0.5980 overall best weight: 0.9182 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN F 93 GLN G 38 ASN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8511 moved from start: 0.3237 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.042 13650 Z= 0.222 Angle : 0.527 5.371 19865 Z= 0.317 Chirality : 0.032 0.133 2262 Planarity : 0.004 0.036 1368 Dihedral : 32.232 178.961 4164 Min Nonbonded Distance : 2.154 Molprobity Statistics. All-atom Clashscore : 5.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.45 % Favored : 97.55 % Rotamer Outliers : 0.16 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 2.98 (0.31), residues: 736 helix: 2.90 (0.21), residues: 544 sheet: None (None), residues: 0 loop : -0.93 (0.41), residues: 192 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 129 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 128 time to evaluate : 0.972 Fit side-chains outliers start: 1 outliers final: 1 residues processed: 128 average time/residue: 0.3134 time to fit residues: 53.8522 Evaluate side-chains 128 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 127 time to evaluate : 0.938 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1063 time to fit residues: 1.2894 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 102 random chunks: chunk 39 optimal weight: 0.6980 chunk 71 optimal weight: 6.9990 chunk 28 optimal weight: 0.9990 chunk 82 optimal weight: 3.9990 chunk 86 optimal weight: 10.0000 chunk 91 optimal weight: 3.9990 chunk 60 optimal weight: 1.9990 chunk 96 optimal weight: 2.9990 chunk 59 optimal weight: 0.8980 chunk 45 optimal weight: 0.2980 chunk 67 optimal weight: 0.5980 overall best weight: 0.6982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN F 93 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8504 moved from start: 0.3252 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.041 13650 Z= 0.193 Angle : 0.520 5.777 19865 Z= 0.314 Chirality : 0.032 0.159 2262 Planarity : 0.004 0.035 1368 Dihedral : 32.220 178.945 4164 Min Nonbonded Distance : 2.062 Molprobity Statistics. All-atom Clashscore : 5.85 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.31 % Favored : 97.69 % Rotamer Outliers : 0.16 % 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.07 (0.31), residues: 736 helix: 2.98 (0.21), residues: 544 sheet: None (None), residues: 0 loop : -0.91 (0.41), residues: 192 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 128 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 127 time to evaluate : 0.967 Fit side-chains outliers start: 1 outliers final: 1 residues processed: 127 average time/residue: 0.3136 time to fit residues: 53.2942 Evaluate side-chains 127 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 126 time to evaluate : 0.917 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.1062 time to fit residues: 1.2686 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 102 random chunks: chunk 101 optimal weight: 4.9990 chunk 93 optimal weight: 2.9990 chunk 80 optimal weight: 2.9990 chunk 8 optimal weight: 0.5980 chunk 62 optimal weight: 2.9990 chunk 49 optimal weight: 1.9990 chunk 64 optimal weight: 0.8980 chunk 86 optimal weight: 10.0000 chunk 24 optimal weight: 0.1980 chunk 74 optimal weight: 3.9990 chunk 11 optimal weight: 0.9990 overall best weight: 0.9384 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN F 93 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8507 moved from start: 0.3387 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.042 13650 Z= 0.221 Angle : 0.522 5.815 19865 Z= 0.314 Chirality : 0.032 0.132 2262 Planarity : 0.004 0.035 1368 Dihedral : 32.130 179.222 4164 Min Nonbonded Distance : 2.144 Molprobity Statistics. All-atom Clashscore : 6.51 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.45 % Favored : 97.55 % Rotamer Outliers : 0.16 % 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.06 (0.31), residues: 736 helix: 2.98 (0.21), residues: 543 sheet: None (None), residues: 0 loop : -0.94 (0.41), residues: 193 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1472 Ramachandran restraints generated. 736 Oldfield, 0 Emsley, 736 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 130 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 129 time to evaluate : 0.924 Fit side-chains outliers start: 1 outliers final: 1 residues processed: 130 average time/residue: 0.3078 time to fit residues: 53.8016 Evaluate side-chains 125 residues out of total 625 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 1 poor density : 124 time to evaluate : 0.958 Switching outliers to nearest non-outliers outliers start: 1 outliers final: 0 residues processed: 1 average time/residue: 0.0956 time to fit residues: 1.2552 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 102 random chunks: chunk 22 optimal weight: 0.7980 chunk 81 optimal weight: 3.9990 chunk 33 optimal weight: 0.0970 chunk 83 optimal weight: 4.9990 chunk 10 optimal weight: 0.6980 chunk 14 optimal weight: 0.0870 chunk 71 optimal weight: 5.9990 chunk 4 optimal weight: 0.3980 chunk 58 optimal weight: 0.9980 chunk 92 optimal weight: 3.9990 chunk 54 optimal weight: 0.7980 overall best weight: 0.4156 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: B 93 GLN F 93 GLN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3852 r_free = 0.3852 target = 0.153681 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 56)----------------| | r_work = 0.3199 r_free = 0.3199 target = 0.104572 restraints weight = 17777.331| |-----------------------------------------------------------------------------| r_work (start): 0.3152 rms_B_bonded: 1.00 r_work: 0.3133 rms_B_bonded: 0.89 restraints_weight: 0.5000 r_work: 0.3118 rms_B_bonded: 0.99 restraints_weight: 0.2500 r_work: 0.3102 rms_B_bonded: 1.17 restraints_weight: 0.1250 r_work: 0.3083 rms_B_bonded: 1.42 restraints_weight: 0.0625 r_work: 0.3061 rms_B_bonded: 1.75 restraints_weight: 0.0312 r_work: 0.3034 rms_B_bonded: 2.17 restraints_weight: 0.0156 r_work: 0.3001 rms_B_bonded: 2.70 restraints_weight: 0.0078 r_work: 0.2961 rms_B_bonded: 3.37 restraints_weight: 0.0039 r_work: 0.2912 rms_B_bonded: 4.25 restraints_weight: 0.0020 r_work (final): 0.2912 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8697 moved from start: 0.3402 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.041 13650 Z= 0.169 Angle : 0.515 6.051 19865 Z= 0.311 Chirality : 0.031 0.168 2262 Planarity : 0.004 0.034 1368 Dihedral : 32.134 179.176 4164 Min Nonbonded Distance : 2.142 Molprobity Statistics. All-atom Clashscore : 6.51 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.31 % Favored : 97.69 % Rotamer Outliers : 0.16 % 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.18 (0.31), residues: 736 helix: 3.07 (0.22), residues: 543 sheet: None (None), residues: 0 loop : -0.91 (0.41), residues: 193 =============================================================================== Job complete usr+sys time: 1980.98 seconds wall clock time: 36 minutes 49.47 seconds (2209.47 seconds total)