Starting phenix.real_space_refine (version: 1.21rc1) on Wed Apr 26 11:34:33 2023 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mqe_23944/04_2023/7mqe_23944.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mqe_23944/04_2023/7mqe_23944.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.69 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mqe_23944/04_2023/7mqe_23944.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mqe_23944/04_2023/7mqe_23944.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mqe_23944/04_2023/7mqe_23944.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7mqe_23944/04_2023/7mqe_23944.pdb" } resolution = 3.69 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.001 sd= 0.012 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Anisotropic ADP refinement not supported. Converting 10872 atoms to isotropic. 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: 5 Type Number sf(0) Gaussians P 16 5.49 5 S 48 5.16 5 C 7752 2.51 5 N 2257 2.21 5 O 2090 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "A GLU 124": "OE1" <-> "OE2" Residue "C PHE 78": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "C GLU 124": "OE1" <-> "OE2" Residue "C GLU 167": "OE1" <-> "OE2" Residue "M PHE 78": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "M GLU 124": "OE1" <-> "OE2" Residue "M TYR 151": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "O PHE 78": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "O GLU 167": "OE1" <-> "OE2" Time to flip residues: 0.04s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21rc1-4938/modules/chem_data/mon_lib" Total number of atoms: 12163 Number of models: 1 Model: "" Number of chains: 16 Chain: "A" Number of atoms: 1505 Number of conformers: 1 Conformer: "" Number of residues, atoms: 205, 1505 Classifications: {'peptide': 205} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 42} Link IDs: {'PTRANS': 5, 'TRANS': 199} Unresolved non-hydrogen bonds: 128 Unresolved non-hydrogen angles: 159 Unresolved non-hydrogen dihedrals: 97 Unresolved non-hydrogen chiralities: 6 Planarities with less than four sites: {'GLU:plan': 10, 'ASN:plan1': 2, 'ARG:plan': 2, 'ASP:plan': 11} Unresolved non-hydrogen planarities: 89 Chain: "C" Number of atoms: 1482 Number of conformers: 1 Conformer: "" Number of residues, atoms: 205, 1482 Classifications: {'peptide': 205} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 50} Link IDs: {'PTRANS': 5, 'TRANS': 199} Unresolved non-hydrogen bonds: 151 Unresolved non-hydrogen angles: 189 Unresolved non-hydrogen dihedrals: 113 Unresolved non-hydrogen chiralities: 8 Planarities with less than four sites: {'GLN:plan1': 2, 'GLU:plan': 8, 'ASN:plan1': 4, 'ARG:plan': 3, 'ASP:plan': 13} Unresolved non-hydrogen planarities: 106 Chain: "E" Number of atoms: 1448 Number of conformers: 1 Conformer: "" Number of residues, atoms: 201, 1448 Classifications: {'peptide': 201} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 47} Link IDs: {'PTRANS': 5, 'TRANS': 195} Chain breaks: 1 Unresolved non-hydrogen bonds: 149 Unresolved non-hydrogen angles: 187 Unresolved non-hydrogen dihedrals: 111 Unresolved non-hydrogen chiralities: 9 Planarities with less than four sites: {'GLU:plan': 10, 'ARG:plan': 4, 'ASN:plan1': 4, 'ASP:plan': 11} Unresolved non-hydrogen planarities: 105 Chain: "G" Number of atoms: 1483 Number of conformers: 1 Conformer: "" Number of residues, atoms: 205, 1483 Classifications: {'peptide': 205} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 49} Link IDs: {'PTRANS': 5, 'TRANS': 199} Unresolved non-hydrogen bonds: 150 Unresolved non-hydrogen angles: 187 Unresolved non-hydrogen dihedrals: 113 Unresolved non-hydrogen chiralities: 7 Planarities with less than four sites: {'GLN:plan1': 2, 'GLU:plan': 9, 'ASN:plan1': 3, 'ARG:plan': 3, 'ASP:plan': 13} Unresolved non-hydrogen planarities: 107 Chain: "I" Number of atoms: 1445 Number of conformers: 1 Conformer: "" Number of residues, atoms: 201, 1445 Classifications: {'peptide': 201} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 48} Link IDs: {'PTRANS': 5, 'TRANS': 195} Chain breaks: 1 Unresolved non-hydrogen bonds: 152 Unresolved non-hydrogen angles: 191 Unresolved non-hydrogen dihedrals: 113 Unresolved non-hydrogen chiralities: 9 Planarities with less than four sites: {'GLU:plan': 10, 'ARG:plan': 4, 'ASN:plan1': 3, 'ASP:plan': 13} Unresolved non-hydrogen planarities: 108 Chain: "K" Number of atoms: 1467 Number of conformers: 1 Conformer: "" Number of residues, atoms: 202, 1467 Classifications: {'peptide': 202} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 44} Link IDs: {'PTRANS': 5, 'TRANS': 196} Chain breaks: 1 Unresolved non-hydrogen bonds: 139 Unresolved non-hydrogen angles: 173 Unresolved non-hydrogen dihedrals: 105 Unresolved non-hydrogen chiralities: 6 Planarities with less than four sites: {'GLU:plan': 10, 'ASN:plan1': 3, 'ARG:plan': 3, 'ASP:plan': 12} Unresolved non-hydrogen planarities: 100 Chain: "M" Number of atoms: 1488 Number of conformers: 1 Conformer: "" Number of residues, atoms: 205, 1488 Classifications: {'peptide': 205} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 46} Link IDs: {'PTRANS': 5, 'TRANS': 199} Unresolved non-hydrogen bonds: 145 Unresolved non-hydrogen angles: 181 Unresolved non-hydrogen dihedrals: 109 Unresolved non-hydrogen chiralities: 6 Planarities with less than four sites: {'GLN:plan1': 3, 'ARG:plan': 2, 'ASN:plan1': 3, 'GLU:plan': 9, 'ASP:plan': 13} Unresolved non-hydrogen planarities: 106 Chain: "O" Number of atoms: 1477 Number of conformers: 1 Conformer: "" Number of residues, atoms: 205, 1477 Classifications: {'peptide': 205} Modifications used: {'COO': 1} Incomplete info: {'truncation_to_alanine': 50} Link IDs: {'PTRANS': 5, 'TRANS': 199} Unresolved non-hydrogen bonds: 156 Unresolved non-hydrogen angles: 198 Unresolved non-hydrogen dihedrals: 114 Unresolved non-hydrogen chiralities: 10 Planarities with less than four sites: {'GLN:plan1': 2, 'GLU:plan': 9, 'ASN:plan1': 5, 'ARG:plan': 3, 'ASP:plan': 13} Unresolved non-hydrogen planarities: 113 Chain: "B" Number of atoms: 46 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 46 Classifications: {'RNA': 2} Modifications used: {'rna3p_pur': 2} Link IDs: {'rna3p': 1} Chain: "D" Number of atoms: 46 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 46 Classifications: {'RNA': 2} Modifications used: {'rna3p_pur': 2} Link IDs: {'rna3p': 1} Chain: "F" Number of atoms: 46 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 46 Classifications: {'RNA': 2} Modifications used: {'rna3p_pur': 2} Link IDs: {'rna3p': 1} Chain: "H" Number of atoms: 46 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 46 Classifications: {'RNA': 2} Modifications used: {'rna3p_pur': 2} Link IDs: {'rna3p': 1} Chain: "J" Number of atoms: 46 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 46 Classifications: {'RNA': 2} Modifications used: {'rna3p_pur': 2} Link IDs: {'rna3p': 1} Chain: "L" Number of atoms: 46 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 46 Classifications: {'RNA': 2} Modifications used: {'rna3p_pur': 2} Link IDs: {'rna3p': 1} Chain: "N" Number of atoms: 46 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 46 Classifications: {'RNA': 2} Modifications used: {'rna3p_pur': 2} Link IDs: {'rna3p': 1} Chain: "P" Number of atoms: 46 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 46 Classifications: {'RNA': 2} Modifications used: {'rna3p_pur': 2} Link IDs: {'rna3p': 1} Time building chain proxies: 7.54, per 1000 atoms: 0.62 Number of scatterers: 12163 At special positions: 0 Unit cell: (107.88, 117.8, 105.4, 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 48 16.00 P 16 15.00 O 2090 8.00 N 2257 7.00 C 7752 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.00 Conformation dependent library (CDL) restraints added in 2.2 seconds 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 3098 Finding SS restraints... Secondary structure from input PDB file: 56 helices and 8 sheets defined 45.5% alpha, 13.1% beta 0 base pairs and 0 stacking pairs defined. Time for finding SS restraints: 1.08 Creating SS restraints... Processing helix chain 'A' and resid 62 through 71 Processing helix chain 'A' and resid 80 through 93 removed outlier: 3.838A pdb=" N ASP A 84 " --> pdb=" O PHE A 80 " (cutoff:3.500A) removed outlier: 3.657A pdb=" N ALA A 86 " --> pdb=" O ARG A 82 " (cutoff:3.500A) Processing helix chain 'A' and resid 102 through 110 Processing helix chain 'A' and resid 117 through 127 Processing helix chain 'A' and resid 145 through 155 removed outlier: 3.739A pdb=" N ASP A 155 " --> pdb=" O TYR A 151 " (cutoff:3.500A) Processing helix chain 'A' and resid 158 through 173 Processing helix chain 'A' and resid 175 through 196 Proline residue: A 187 - end of helix Processing helix chain 'C' and resid 62 through 71 Processing helix chain 'C' and resid 80 through 93 removed outlier: 3.844A pdb=" N ASP C 84 " --> pdb=" O PHE C 80 " (cutoff:3.500A) removed outlier: 3.618A pdb=" N ALA C 86 " --> pdb=" O ARG C 82 " (cutoff:3.500A) Processing helix chain 'C' and resid 102 through 110 Processing helix chain 'C' and resid 117 through 127 Processing helix chain 'C' and resid 145 through 155 removed outlier: 3.622A pdb=" N ASP C 155 " --> pdb=" O TYR C 151 " (cutoff:3.500A) Processing helix chain 'C' and resid 158 through 173 Processing helix chain 'C' and resid 175 through 196 Proline residue: C 187 - end of helix Processing helix chain 'E' and resid 64 through 71 Processing helix chain 'E' and resid 80 through 93 removed outlier: 3.859A pdb=" N ASP E 84 " --> pdb=" O PHE E 80 " (cutoff:3.500A) removed outlier: 3.656A pdb=" N ALA E 86 " --> pdb=" O ARG E 82 " (cutoff:3.500A) Processing helix chain 'E' and resid 102 through 110 Processing helix chain 'E' and resid 117 through 127 Processing helix chain 'E' and resid 145 through 155 Processing helix chain 'E' and resid 158 through 173 Processing helix chain 'E' and resid 175 through 196 Proline residue: E 187 - end of helix Processing helix chain 'G' and resid 62 through 71 Processing helix chain 'G' and resid 80 through 93 removed outlier: 3.921A pdb=" N ASP G 84 " --> pdb=" O PHE G 80 " (cutoff:3.500A) removed outlier: 3.724A pdb=" N ALA G 86 " --> pdb=" O ARG G 82 " (cutoff:3.500A) Processing helix chain 'G' and resid 102 through 110 Processing helix chain 'G' and resid 117 through 127 Processing helix chain 'G' and resid 145 through 155 Processing helix chain 'G' and resid 158 through 173 Processing helix chain 'G' and resid 175 through 196 Proline residue: G 187 - end of helix Processing helix chain 'I' and resid 62 through 70 Processing helix chain 'I' and resid 80 through 93 removed outlier: 3.924A pdb=" N ASP I 84 " --> pdb=" O PHE I 80 " (cutoff:3.500A) removed outlier: 3.685A pdb=" N ALA I 86 " --> pdb=" O ARG I 82 " (cutoff:3.500A) Processing helix chain 'I' and resid 102 through 110 Processing helix chain 'I' and resid 117 through 127 Processing helix chain 'I' and resid 145 through 155 Processing helix chain 'I' and resid 158 through 173 Processing helix chain 'I' and resid 175 through 196 Proline residue: I 187 - end of helix Processing helix chain 'K' and resid 62 through 71 Processing helix chain 'K' and resid 80 through 93 removed outlier: 3.912A pdb=" N ASP K 84 " --> pdb=" O PHE K 80 " (cutoff:3.500A) removed outlier: 3.705A pdb=" N ALA K 86 " --> pdb=" O ARG K 82 " (cutoff:3.500A) Processing helix chain 'K' and resid 102 through 110 Processing helix chain 'K' and resid 117 through 127 Processing helix chain 'K' and resid 145 through 155 removed outlier: 3.666A pdb=" N ASP K 155 " --> pdb=" O TYR K 151 " (cutoff:3.500A) Processing helix chain 'K' and resid 158 through 173 Processing helix chain 'K' and resid 175 through 196 Proline residue: K 187 - end of helix Processing helix chain 'M' and resid 62 through 71 Processing helix chain 'M' and resid 80 through 93 removed outlier: 3.827A pdb=" N ASP M 84 " --> pdb=" O PHE M 80 " (cutoff:3.500A) removed outlier: 3.656A pdb=" N ALA M 86 " --> pdb=" O ARG M 82 " (cutoff:3.500A) Processing helix chain 'M' and resid 102 through 110 Processing helix chain 'M' and resid 117 through 127 Processing helix chain 'M' and resid 145 through 155 removed outlier: 3.581A pdb=" N ASP M 155 " --> pdb=" O TYR M 151 " (cutoff:3.500A) Processing helix chain 'M' and resid 158 through 173 Processing helix chain 'M' and resid 175 through 196 Proline residue: M 187 - end of helix Processing helix chain 'O' and resid 62 through 71 Processing helix chain 'O' and resid 80 through 93 removed outlier: 3.834A pdb=" N ASP O 84 " --> pdb=" O PHE O 80 " (cutoff:3.500A) removed outlier: 3.583A pdb=" N ALA O 86 " --> pdb=" O ARG O 82 " (cutoff:3.500A) Processing helix chain 'O' and resid 102 through 110 Processing helix chain 'O' and resid 117 through 127 Processing helix chain 'O' and resid 145 through 155 removed outlier: 3.534A pdb=" N ASP O 155 " --> pdb=" O TYR O 151 " (cutoff:3.500A) Processing helix chain 'O' and resid 158 through 173 Processing helix chain 'O' and resid 175 through 196 Proline residue: O 187 - end of helix Processing sheet with id=AA1, first strand: chain 'A' and resid 4 through 8 removed outlier: 6.267A pdb=" N ARG A 5 " --> pdb=" O VAL A 52 " (cutoff:3.500A) removed outlier: 7.446A pdb=" N GLN A 54 " --> pdb=" O ARG A 5 " (cutoff:3.500A) removed outlier: 6.746A pdb=" N HIS A 7 " --> pdb=" O GLN A 54 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'C' and resid 4 through 8 removed outlier: 6.084A pdb=" N ARG C 5 " --> pdb=" O VAL C 52 " (cutoff:3.500A) removed outlier: 7.281A pdb=" N GLN C 54 " --> pdb=" O ARG C 5 " (cutoff:3.500A) removed outlier: 6.605A pdb=" N HIS C 7 " --> pdb=" O GLN C 54 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'E' and resid 4 through 8 removed outlier: 6.276A pdb=" N ARG E 5 " --> pdb=" O VAL E 52 " (cutoff:3.500A) removed outlier: 7.388A pdb=" N GLN E 54 " --> pdb=" O ARG E 5 " (cutoff:3.500A) removed outlier: 6.689A pdb=" N HIS E 7 " --> pdb=" O GLN E 54 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'G' and resid 4 through 8 removed outlier: 6.240A pdb=" N ARG G 5 " --> pdb=" O VAL G 52 " (cutoff:3.500A) removed outlier: 7.308A pdb=" N GLN G 54 " --> pdb=" O ARG G 5 " (cutoff:3.500A) removed outlier: 6.668A pdb=" N HIS G 7 " --> pdb=" O GLN G 54 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'I' and resid 3 through 8 removed outlier: 6.458A pdb=" N GLU I 3 " --> pdb=" O ALA I 50 " (cutoff:3.500A) removed outlier: 7.531A pdb=" N VAL I 52 " --> pdb=" O GLU I 3 " (cutoff:3.500A) removed outlier: 6.195A pdb=" N ARG I 5 " --> pdb=" O VAL I 52 " (cutoff:3.500A) removed outlier: 7.310A pdb=" N GLN I 54 " --> pdb=" O ARG I 5 " (cutoff:3.500A) removed outlier: 6.578A pdb=" N HIS I 7 " --> pdb=" O GLN I 54 " (cutoff:3.500A) removed outlier: 6.662A pdb=" N VAL I 22 " --> pdb=" O ILE I 77 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'K' and resid 4 through 8 removed outlier: 6.216A pdb=" N ARG K 5 " --> pdb=" O VAL K 52 " (cutoff:3.500A) removed outlier: 7.340A pdb=" N GLN K 54 " --> pdb=" O ARG K 5 " (cutoff:3.500A) removed outlier: 6.726A pdb=" N HIS K 7 " --> pdb=" O GLN K 54 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'M' and resid 3 through 8 removed outlier: 6.330A pdb=" N GLU M 3 " --> pdb=" O ALA M 50 " (cutoff:3.500A) removed outlier: 7.388A pdb=" N VAL M 52 " --> pdb=" O GLU M 3 " (cutoff:3.500A) removed outlier: 6.115A pdb=" N ARG M 5 " --> pdb=" O VAL M 52 " (cutoff:3.500A) removed outlier: 7.334A pdb=" N GLN M 54 " --> pdb=" O ARG M 5 " (cutoff:3.500A) removed outlier: 6.727A pdb=" N HIS M 7 " --> pdb=" O GLN M 54 " (cutoff:3.500A) Processing sheet with id=AA8, first strand: chain 'O' and resid 3 through 8 removed outlier: 6.412A pdb=" N GLU O 3 " --> pdb=" O ALA O 50 " (cutoff:3.500A) removed outlier: 7.415A pdb=" N VAL O 52 " --> pdb=" O GLU O 3 " (cutoff:3.500A) removed outlier: 5.998A pdb=" N ARG O 5 " --> pdb=" O VAL O 52 " (cutoff:3.500A) removed outlier: 7.261A pdb=" N GLN O 54 " --> pdb=" O ARG O 5 " (cutoff:3.500A) removed outlier: 6.660A pdb=" N HIS O 7 " --> pdb=" O GLN O 54 " (cutoff:3.500A) 589 hydrogen bonds defined for protein. 1767 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: 3.98 Time building geometry restraints manager: 5.39 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.34: 3987 1.34 - 1.46: 2396 1.46 - 1.58: 5953 1.58 - 1.69: 24 1.69 - 1.81: 64 Bond restraints: 12424 Sorted by residual: bond pdb=" C ARG C 67 " pdb=" N LEU C 68 " ideal model delta sigma weight residual 1.335 1.288 0.047 1.36e-02 5.41e+03 1.18e+01 bond pdb=" N VAL O 41 " pdb=" CA VAL O 41 " ideal model delta sigma weight residual 1.458 1.494 -0.035 1.20e-02 6.94e+03 8.71e+00 bond pdb=" N VAL G 128 " pdb=" CA VAL G 128 " ideal model delta sigma weight residual 1.457 1.495 -0.038 1.32e-02 5.74e+03 8.30e+00 bond pdb=" N VAL C 40 " pdb=" CA VAL C 40 " ideal model delta sigma weight residual 1.459 1.493 -0.034 1.19e-02 7.06e+03 8.14e+00 bond pdb=" N VAL C 41 " pdb=" CA VAL C 41 " ideal model delta sigma weight residual 1.459 1.493 -0.034 1.20e-02 6.94e+03 8.08e+00 ... (remaining 12419 not shown) Histogram of bond angle deviations from ideal: 98.31 - 105.46: 222 105.46 - 112.60: 6748 112.60 - 119.75: 4117 119.75 - 126.90: 5697 126.90 - 134.04: 204 Bond angle restraints: 16988 Sorted by residual: angle pdb=" N GLN A 135 " pdb=" CA GLN A 135 " pdb=" C GLN A 135 " ideal model delta sigma weight residual 112.54 107.67 4.87 1.22e+00 6.72e-01 1.59e+01 angle pdb=" CA VAL G 128 " pdb=" C VAL G 128 " pdb=" O VAL G 128 " ideal model delta sigma weight residual 121.59 117.83 3.76 1.03e+00 9.43e-01 1.33e+01 angle pdb=" CA LEU C 43 " pdb=" C LEU C 43 " pdb=" O LEU C 43 " ideal model delta sigma weight residual 121.68 117.50 4.18 1.18e+00 7.18e-01 1.25e+01 angle pdb=" N GLN G 135 " pdb=" CA GLN G 135 " pdb=" C GLN G 135 " ideal model delta sigma weight residual 112.89 108.53 4.36 1.24e+00 6.50e-01 1.23e+01 angle pdb=" CA LEU O 43 " pdb=" C LEU O 43 " pdb=" O LEU O 43 " ideal model delta sigma weight residual 121.46 117.55 3.91 1.17e+00 7.31e-01 1.12e+01 ... (remaining 16983 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 18.08: 6844 18.08 - 36.16: 331 36.16 - 54.24: 67 54.24 - 72.33: 25 72.33 - 90.41: 19 Dihedral angle restraints: 7286 sinusoidal: 2538 harmonic: 4748 Sorted by residual: dihedral pdb=" CA GLN M 134 " pdb=" C GLN M 134 " pdb=" N GLN M 135 " pdb=" CA GLN M 135 " ideal model delta harmonic sigma weight residual 180.00 154.05 25.95 0 5.00e+00 4.00e-02 2.69e+01 dihedral pdb=" CA PHE O 78 " pdb=" C PHE O 78 " pdb=" N HIS O 79 " pdb=" CA HIS O 79 " ideal model delta harmonic sigma weight residual 180.00 157.61 22.39 0 5.00e+00 4.00e-02 2.01e+01 dihedral pdb=" CA PHE I 78 " pdb=" C PHE I 78 " pdb=" N HIS I 79 " pdb=" CA HIS I 79 " ideal model delta harmonic sigma weight residual 180.00 157.67 22.33 0 5.00e+00 4.00e-02 1.99e+01 ... (remaining 7283 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.055: 1384 0.055 - 0.110: 540 0.110 - 0.164: 121 0.164 - 0.219: 12 0.219 - 0.274: 7 Chirality restraints: 2064 Sorted by residual: chirality pdb=" CA LEU C 38 " pdb=" N LEU C 38 " pdb=" C LEU C 38 " pdb=" CB LEU C 38 " both_signs ideal model delta sigma weight residual False 2.51 2.78 -0.27 2.00e-01 2.50e+01 1.87e+00 chirality pdb=" CA LEU O 38 " pdb=" N LEU O 38 " pdb=" C LEU O 38 " pdb=" CB LEU O 38 " both_signs ideal model delta sigma weight residual False 2.51 2.77 -0.26 2.00e-01 2.50e+01 1.75e+00 chirality pdb=" CA VAL O 41 " pdb=" N VAL O 41 " pdb=" C VAL O 41 " pdb=" CB VAL O 41 " both_signs ideal model delta sigma weight residual False 2.44 2.69 -0.25 2.00e-01 2.50e+01 1.57e+00 ... (remaining 2061 not shown) Planarity restraints: 2095 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CA ALA M 153 " -0.010 2.00e-02 2.50e+03 1.98e-02 3.91e+00 pdb=" C ALA M 153 " 0.034 2.00e-02 2.50e+03 pdb=" O ALA M 153 " -0.013 2.00e-02 2.50e+03 pdb=" N SER M 154 " -0.012 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA ARG C 37 " 0.010 2.00e-02 2.50e+03 1.92e-02 3.68e+00 pdb=" C ARG C 37 " -0.033 2.00e-02 2.50e+03 pdb=" O ARG C 37 " 0.013 2.00e-02 2.50e+03 pdb=" N LEU C 38 " 0.011 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CA ARG O 37 " 0.009 2.00e-02 2.50e+03 1.89e-02 3.57e+00 pdb=" C ARG O 37 " -0.033 2.00e-02 2.50e+03 pdb=" O ARG O 37 " 0.012 2.00e-02 2.50e+03 pdb=" N LEU O 38 " 0.011 2.00e-02 2.50e+03 ... (remaining 2092 not shown) Histogram of nonbonded interaction distances: 2.31 - 2.83: 3527 2.83 - 3.34: 10699 3.34 - 3.86: 19724 3.86 - 4.38: 22813 4.38 - 4.90: 40034 Nonbonded interactions: 96797 Sorted by model distance: nonbonded pdb=" NH1 ARG I 82 " pdb=" OD2 ASP I 193 " model vdw 2.308 2.520 nonbonded pdb=" OH TYR C 112 " pdb=" OE2 GLU C 124 " model vdw 2.327 2.440 nonbonded pdb=" OG1 THR K 26 " pdb=" OD1 ASP K 84 " model vdw 2.428 2.440 nonbonded pdb=" OG1 THR G 26 " pdb=" OD1 ASP G 84 " model vdw 2.435 2.440 nonbonded pdb=" O LYS I 103 " pdb=" OG SER I 106 " model vdw 2.442 2.440 ... (remaining 96792 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 2 through 5 and (name N or name CA or name C or name O or \ name CB )) or resid 6 through 25 or (resid 26 and (name N or name CA or name C \ or name O or name CB )) or resid 27 or (resid 28 and (name N or name CA or name \ C or name O or name CB )) or resid 29 through 46 or (resid 47 and (name N or nam \ e CA or name C or name O or name CB )) or resid 48 through 56 or (resid 57 and ( \ name N or name CA or name C or name O or name CB )) or resid 58 through 63 or (r \ esid 64 and (name N or name CA or name C or name O or name CB )) or resid 65 thr \ ough 66 or (resid 67 and (name N or name CA or name C or name O or name CB )) or \ resid 68 through 105 or (resid 106 and (name N or name CA or name C or name O o \ r name CB )) or resid 107 through 114 or (resid 115 and (name N or name CA or na \ me C or name O or name CB )) or resid 116 through 118 or (resid 119 through 120 \ and (name N or name CA or name C or name O or name CB )) or resid 121 through 12 \ 3 or (resid 124 and (name N or name CA or name C or name O or name CB )) or resi \ d 125 through 126 or (resid 127 through 131 and (name N or name CA or name C or \ name O or name CB )) or resid 136 through 148 or (resid 149 through 150 and (nam \ e N or name CA or name C or name O or name CB )) or resid 151 through 163 or (re \ sid 164 through 165 and (name N or name CA or name C or name O or name CB )) or \ resid 166 through 176 or (resid 177 and (name N or name CA or name C or name O o \ r name CB )) or resid 178 through 200 or (resid 201 and (name N or name CA or na \ me C or name O or name CB )) or resid 202 through 206)) selection = (chain 'C' and ((resid 2 through 5 and (name N or name CA or name C or name O or \ name CB )) or resid 6 through 63 or (resid 64 and (name N or name CA or name C \ or name O or name CB )) or resid 65 through 66 or (resid 67 and (name N or name \ CA or name C or name O or name CB )) or resid 68 through 70 or (resid 71 through \ 73 and (name N or name CA or name C or name O or name CB )) or resid 74 through \ 93 or (resid 94 and (name N or name CA or name C or name O or name CB )) or res \ id 95 through 119 or (resid 120 and (name N or name CA or name C or name O or na \ me CB )) or resid 121 through 123 or (resid 124 and (name N or name CA or name C \ or name O or name CB )) or resid 125 through 131 or resid 136 through 148 or (r \ esid 149 through 150 and (name N or name CA or name C or name O or name CB )) or \ resid 151 through 164 or (resid 165 and (name N or name CA or name C or name O \ or name CB )) or resid 166 or (resid 167 through 168 and (name N or name CA or n \ ame C or name O or name CB )) or resid 169 through 176 or (resid 177 and (name N \ or name CA or name C or name O or name CB )) or resid 178 through 179 or (resid \ 180 through 181 and (name N or name CA or name C or name O or name CB )) or res \ id 182 through 206)) selection = (chain 'E' and ((resid 2 through 5 and (name N or name CA or name C or name O or \ name CB )) or resid 6 through 15 or (resid 16 and (name N or name CA or name C \ or name O or name CB )) or resid 17 through 25 or (resid 26 and (name N or name \ CA or name C or name O or name CB )) or resid 27 or (resid 28 and (name N or nam \ e CA or name C or name O or name CB )) or resid 29 through 56 or (resid 57 and ( \ name N or name CA or name C or name O or name CB )) or resid 58 through 60 or (r \ esid 61 through 62 and (name N or name CA or name C or name O or name CB )) or r \ esid 63 through 66 or (resid 67 and (name N or name CA or name C or name O or na \ me CB )) or resid 68 through 105 or (resid 106 and (name N or name CA or name C \ or name O or name CB )) or resid 107 through 154 or (resid 155 and (name N or na \ me CA or name C or name O or name CB )) or resid 156 through 163 or (resid 164 t \ hrough 165 and (name N or name CA or name C or name O or name CB )) or resid 166 \ through 172 or (resid 173 through 174 and (name N or name CA or name C or name \ O or name CB )) or resid 175 through 176 or (resid 177 and (name N or name CA or \ name C or name O or name CB )) or resid 178 through 179 or (resid 180 through 1 \ 81 and (name N or name CA or name C or name O or name CB )) or resid 182 through \ 200 or (resid 201 and (name N or name CA or name C or name O or name CB )) or r \ esid 202 through 206)) selection = (chain 'G' and ((resid 2 through 5 and (name N or name CA or name C or name O or \ name CB )) or resid 6 through 25 or (resid 26 and (name N or name CA or name C \ or name O or name CB )) or resid 27 through 56 or (resid 57 and (name N or name \ CA or name C or name O or name CB )) or resid 58 through 63 or (resid 64 and (na \ me N or name CA or name C or name O or name CB )) or resid 65 through 105 or (re \ sid 106 and (name N or name CA or name C or name O or name CB )) or resid 107 th \ rough 114 or (resid 115 and (name N or name CA or name C or name O or name CB )) \ or resid 116 through 127 or (resid 128 through 131 and (name N or name CA or na \ me C or name O or name CB )) or resid 136 through 148 or (resid 149 through 150 \ and (name N or name CA or name C or name O or name CB )) or resid 151 through 15 \ 4 or (resid 155 and (name N or name CA or name C or name O or name CB )) or resi \ d 156 through 164 or (resid 165 and (name N or name CA or name C or name O or na \ me CB )) or resid 166 or (resid 167 through 168 and (name N or name CA or name C \ or name O or name CB )) or resid 169 through 172 or (resid 173 through 174 and \ (name N or name CA or name C or name O or name CB )) or resid 175 through 206)) selection = (chain 'I' and ((resid 2 through 5 and (name N or name CA or name C or name O or \ name CB )) or resid 6 through 15 or (resid 16 and (name N or name CA or name C \ or name O or name CB )) or resid 17 through 27 or (resid 28 and (name N or name \ CA or name C or name O or name CB )) or resid 29 through 56 or (resid 57 and (na \ me N or name CA or name C or name O or name CB )) or resid 58 through 66 or (res \ id 67 and (name N or name CA or name C or name O or name CB )) or resid 68 throu \ gh 69 or (resid 70 through 73 and (name N or name CA or name C or name O or name \ CB )) or resid 74 through 93 or (resid 94 and (name N or name CA or name C or n \ ame O or name CB )) or resid 95 through 105 or (resid 106 and (name N or name CA \ or name C or name O or name CB )) or resid 107 through 119 or (resid 120 and (n \ ame N or name CA or name C or name O or name CB )) or resid 121 through 128 or ( \ resid 129 through 138 and (name N or name CA or name C or name O or name CB )) o \ r resid 139 through 154 or (resid 155 and (name N or name CA or name C or name O \ or name CB )) or resid 156 through 176 or (resid 177 and (name N or name CA or \ name C or name O or name CB )) or resid 178 through 179 or (resid 180 through 18 \ 1 and (name N or name CA or name C or name O or name CB )) or resid 182 through \ 206)) selection = (chain 'K' and ((resid 2 through 5 and (name N or name CA or name C or name O or \ name CB )) or resid 6 through 25 or (resid 26 and (name N or name CA or name C \ or name O or name CB )) or resid 27 or (resid 28 and (name N or name CA or name \ C or name O or name CB )) or resid 29 through 56 or (resid 57 and (name N or nam \ e CA or name C or name O or name CB )) or resid 58 through 60 or (resid 61 throu \ gh 62 and (name N or name CA or name C or name O or name CB )) or resid 63 or (r \ esid 64 and (name N or name CA or name C or name O or name CB )) or resid 65 thr \ ough 105 or (resid 106 and (name N or name CA or name C or name O or name CB )) \ or resid 107 through 114 or (resid 115 and (name N or name CA or name C or name \ O or name CB )) or resid 116 through 118 or (resid 119 through 120 and (name N o \ r name CA or name C or name O or name CB )) or resid 121 through 127 or (resid 1 \ 28 through 131 and (name N or name CA or name C or name O or name CB )) or (resi \ d 136 through 138 and (name N or name CA or name C or name O or name CB )) or re \ sid 139 through 148 or (resid 149 through 150 and (name N or name CA or name C o \ r name O or name CB )) or resid 151 through 163 or (resid 164 through 165 and (n \ ame N or name CA or name C or name O or name CB )) or resid 166 through 172 or ( \ resid 173 through 174 and (name N or name CA or name C or name O or name CB )) o \ r resid 175 through 176 or (resid 177 and (name N or name CA or name C or name O \ or name CB )) or resid 178 through 206)) selection = (chain 'M' and ((resid 2 through 5 and (name N or name CA or name C or name O or \ name CB )) or resid 6 through 15 or (resid 16 and (name N or name CA or name C \ or name O or name CB )) or resid 17 through 25 or (resid 26 and (name N or name \ CA or name C or name O or name CB )) or resid 27 or (resid 28 and (name N or nam \ e CA or name C or name O or name CB )) or resid 29 through 56 or (resid 57 and ( \ name N or name CA or name C or name O or name CB )) or resid 58 through 63 or (r \ esid 64 and (name N or name CA or name C or name O or name CB )) or resid 65 thr \ ough 66 or (resid 67 and (name N or name CA or name C or name O or name CB )) or \ resid 68 through 69 or (resid 70 through 73 and (name N or name CA or name C or \ name O or name CB )) or resid 74 through 93 or (resid 94 and (name N or name CA \ or name C or name O or name CB )) or resid 95 through 105 or (resid 106 and (na \ me N or name CA or name C or name O or name CB )) or resid 107 through 114 or (r \ esid 115 and (name N or name CA or name C or name O or name CB )) or resid 116 t \ hrough 123 or (resid 124 and (name N or name CA or name C or name O or name CB ) \ ) or resid 125 through 127 or (resid 128 through 131 and (name N or name CA or n \ ame C or name O or name CB )) or resid 136 through 148 or (resid 149 through 150 \ and (name N or name CA or name C or name O or name CB )) or resid 151 through 1 \ 64 or (resid 165 and (name N or name CA or name C or name O or name CB )) or res \ id 166 through 176 or (resid 177 and (name N or name CA or name C or name O or n \ ame CB )) or resid 178 through 179 or (resid 180 through 181 and (name N or name \ CA or name C or name O or name CB )) or resid 182 through 200 or (resid 201 and \ (name N or name CA or name C or name O or name CB )) or resid 202 through 206)) \ selection = (chain 'O' and ((resid 2 through 5 and (name N or name CA or name C or name O or \ name CB )) or resid 6 through 56 or (resid 57 and (name N or name CA or name C \ or name O or name CB )) or resid 58 through 60 or (resid 61 through 62 and (name \ N or name CA or name C or name O or name CB )) or resid 63 through 66 or (resid \ 67 and (name N or name CA or name C or name O or name CB )) or resid 68 through \ 69 or (resid 70 through 73 and (name N or name CA or name C or name O or name C \ B )) or resid 74 through 105 or (resid 106 and (name N or name CA or name C or n \ ame O or name CB )) or resid 107 through 118 or (resid 119 through 120 and (name \ N or name CA or name C or name O or name CB )) or resid 121 through 123 or (res \ id 124 and (name N or name CA or name C or name O or name CB )) or resid 125 thr \ ough 131 or resid 136 through 148 or (resid 149 through 150 and (name N or name \ CA or name C or name O or name CB )) or resid 151 through 164 or (resid 165 and \ (name N or name CA or name C or name O or name CB )) or resid 166 or (resid 167 \ through 168 and (name N or name CA or name C or name O or name CB )) or resid 16 \ 9 through 176 or (resid 177 and (name N or name CA or name C or name O or name C \ B )) or resid 178 through 179 or (resid 180 through 181 and (name N or name CA o \ r name C or name O or name CB )) or resid 182 through 206)) } ncs_group { reference = chain 'B' selection = chain 'D' selection = chain 'F' selection = chain 'H' selection = chain 'J' selection = chain 'L' selection = chain 'N' selection = chain 'P' } 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 group one per residue Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 0.100 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.140 Construct map_model_manager: 0.040 Extract box with map and model: 1.800 Check model and map are aligned: 0.200 Set scattering table: 0.120 Process input model: 35.300 Find NCS groups from input model: 0.970 Set up NCS constraints: 0.090 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:2.880 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 41.640 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7760 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.071 12424 Z= 0.545 Angle : 0.840 7.763 16988 Z= 0.510 Chirality : 0.059 0.274 2064 Planarity : 0.005 0.032 2095 Dihedral : 13.453 90.408 4188 Min Nonbonded Distance : 2.308 Molprobity Statistics. All-atom Clashscore : 1.58 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.74 % Favored : 98.26 % Rotamer Outliers : 0.77 % 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.67 (0.20), residues: 1607 helix: 0.34 (0.18), residues: 726 sheet: 0.22 (0.24), residues: 253 loop : 0.98 (0.27), residues: 628 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 195 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 8 poor density : 187 time to evaluate : 1.499 Fit side-chains outliers start: 8 outliers final: 0 residues processed: 195 average time/residue: 0.2362 time to fit residues: 66.7189 Evaluate side-chains 156 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 156 time to evaluate : 1.413 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 1.8300 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 160 random chunks: chunk 135 optimal weight: 8.9990 chunk 121 optimal weight: 10.0000 chunk 67 optimal weight: 0.9980 chunk 41 optimal weight: 10.0000 chunk 81 optimal weight: 9.9990 chunk 64 optimal weight: 5.9990 chunk 125 optimal weight: 1.9990 chunk 48 optimal weight: 1.9990 chunk 76 optimal weight: 0.7980 chunk 93 optimal weight: 0.3980 chunk 145 optimal weight: 0.2980 overall best weight: 0.8982 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 8 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 54 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** G 8 GLN ** G 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** G 205 HIS ** I 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** K 8 GLN ** K 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** M 8 GLN M 16 ASN Total number of N/Q/H flips: 5 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7697 moved from start: 0.1620 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.017 12424 Z= 0.128 Angle : 0.446 7.470 16988 Z= 0.239 Chirality : 0.043 0.130 2064 Planarity : 0.003 0.022 2095 Dihedral : 5.937 43.523 1922 Min Nonbonded Distance : 2.309 Molprobity Statistics. All-atom Clashscore : 2.43 Ramachandran Plot: Outliers : 0.00 % Allowed : 0.62 % Favored : 99.38 % Rotamer Outliers : 0.68 % 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.03 (0.22), residues: 1607 helix: 1.90 (0.20), residues: 730 sheet: 1.10 (0.24), residues: 253 loop : 0.96 (0.28), residues: 624 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 183 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 7 poor density : 176 time to evaluate : 1.343 Fit side-chains outliers start: 7 outliers final: 3 residues processed: 179 average time/residue: 0.2252 time to fit residues: 59.3401 Evaluate side-chains 146 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 143 time to evaluate : 1.463 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1142 time to fit residues: 2.4969 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 160 random chunks: chunk 80 optimal weight: 20.0000 chunk 45 optimal weight: 6.9990 chunk 120 optimal weight: 10.0000 chunk 98 optimal weight: 0.7980 chunk 40 optimal weight: 5.9990 chunk 145 optimal weight: 0.9990 chunk 157 optimal weight: 10.0000 chunk 129 optimal weight: 8.9990 chunk 144 optimal weight: 1.9990 chunk 49 optimal weight: 0.8980 chunk 116 optimal weight: 1.9990 overall best weight: 1.3386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 8 GLN A 42 GLN E 90 HIS ** G 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** I 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** K 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** O 8 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7729 moved from start: 0.1995 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.002 0.018 12424 Z= 0.153 Angle : 0.425 7.663 16988 Z= 0.227 Chirality : 0.043 0.131 2064 Planarity : 0.003 0.023 2095 Dihedral : 5.411 40.041 1922 Min Nonbonded Distance : 2.247 Molprobity Statistics. All-atom Clashscore : 2.30 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.31 % Favored : 98.69 % Rotamer Outliers : 1.25 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 2.02 (0.22), residues: 1607 helix: 1.91 (0.20), residues: 731 sheet: 1.43 (0.27), residues: 257 loop : 0.75 (0.28), residues: 619 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 175 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 162 time to evaluate : 1.548 Fit side-chains outliers start: 13 outliers final: 11 residues processed: 169 average time/residue: 0.2215 time to fit residues: 56.5327 Evaluate side-chains 156 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 11 poor density : 145 time to evaluate : 1.526 Switching outliers to nearest non-outliers outliers start: 11 outliers final: 0 residues processed: 11 average time/residue: 0.1226 time to fit residues: 4.7238 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 160 random chunks: chunk 143 optimal weight: 8.9990 chunk 109 optimal weight: 10.0000 chunk 75 optimal weight: 7.9990 chunk 16 optimal weight: 4.9990 chunk 69 optimal weight: 7.9990 chunk 97 optimal weight: 0.9980 chunk 146 optimal weight: 2.9990 chunk 154 optimal weight: 2.9990 chunk 76 optimal weight: 0.9980 chunk 138 optimal weight: 6.9990 chunk 41 optimal weight: 0.0010 overall best weight: 1.5990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** G 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** I 32 GLN ** I 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** I 90 HIS ** K 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** O 8 GLN O 54 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7740 moved from start: 0.2221 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.020 12424 Z= 0.171 Angle : 0.419 8.512 16988 Z= 0.224 Chirality : 0.043 0.127 2064 Planarity : 0.003 0.023 2095 Dihedral : 5.277 38.206 1922 Min Nonbonded Distance : 2.260 Molprobity Statistics. All-atom Clashscore : 2.73 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.80 % Favored : 98.20 % Rotamer Outliers : 0.29 % 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.90 (0.22), residues: 1607 helix: 1.79 (0.20), residues: 731 sheet: 1.63 (0.28), residues: 257 loop : 0.59 (0.28), residues: 619 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 151 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 148 time to evaluate : 1.331 Fit side-chains outliers start: 3 outliers final: 3 residues processed: 149 average time/residue: 0.2020 time to fit residues: 46.2912 Evaluate side-chains 145 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 142 time to evaluate : 1.338 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1191 time to fit residues: 2.4302 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 160 random chunks: chunk 128 optimal weight: 7.9990 chunk 87 optimal weight: 9.9990 chunk 2 optimal weight: 10.0000 chunk 115 optimal weight: 4.9990 chunk 63 optimal weight: 1.9990 chunk 131 optimal weight: 8.9990 chunk 106 optimal weight: 8.9990 chunk 0 optimal weight: 10.0000 chunk 78 optimal weight: 5.9990 chunk 138 optimal weight: 3.9990 chunk 39 optimal weight: 0.3980 overall best weight: 3.4788 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 54 GLN ** C 135 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** C 205 HIS ** G 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** I 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** K 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** M 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** O 8 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7802 moved from start: 0.2319 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.005 0.032 12424 Z= 0.327 Angle : 0.491 8.225 16988 Z= 0.260 Chirality : 0.045 0.129 2064 Planarity : 0.003 0.026 2095 Dihedral : 5.626 37.177 1922 Min Nonbonded Distance : 2.266 Molprobity Statistics. All-atom Clashscore : 2.55 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.68 % Favored : 97.32 % Rotamer Outliers : 1.35 % 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.53 (0.22), residues: 1607 helix: 1.41 (0.20), residues: 731 sheet: 1.50 (0.29), residues: 257 loop : 0.47 (0.28), residues: 619 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 155 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 14 poor density : 141 time to evaluate : 1.501 Fit side-chains outliers start: 14 outliers final: 13 residues processed: 146 average time/residue: 0.2043 time to fit residues: 46.5319 Evaluate side-chains 146 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 13 poor density : 133 time to evaluate : 1.434 Switching outliers to nearest non-outliers outliers start: 13 outliers final: 0 residues processed: 13 average time/residue: 0.1127 time to fit residues: 4.6545 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 160 random chunks: chunk 52 optimal weight: 8.9990 chunk 139 optimal weight: 1.9990 chunk 30 optimal weight: 3.9990 chunk 90 optimal weight: 4.9990 chunk 38 optimal weight: 1.9990 chunk 154 optimal weight: 9.9990 chunk 128 optimal weight: 4.9990 chunk 71 optimal weight: 9.9990 chunk 12 optimal weight: 4.9990 chunk 51 optimal weight: 10.0000 chunk 81 optimal weight: 0.6980 overall best weight: 2.7388 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 135 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** G 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** I 59 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** K 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** M 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** O 8 GLN O 54 GLN O 135 GLN Total number of N/Q/H flips: 3 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7791 moved from start: 0.2458 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.029 12424 Z= 0.265 Angle : 0.457 8.466 16988 Z= 0.242 Chirality : 0.044 0.126 2064 Planarity : 0.003 0.025 2095 Dihedral : 5.563 38.479 1922 Min Nonbonded Distance : 2.273 Molprobity Statistics. All-atom Clashscore : 2.90 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.55 % Favored : 97.45 % Rotamer Outliers : 0.29 % 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.46 (0.22), residues: 1607 helix: 1.38 (0.20), residues: 731 sheet: 1.55 (0.29), residues: 257 loop : 0.36 (0.28), residues: 619 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 145 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 142 time to evaluate : 1.361 Fit side-chains outliers start: 3 outliers final: 3 residues processed: 143 average time/residue: 0.2141 time to fit residues: 46.6828 Evaluate side-chains 137 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 134 time to evaluate : 1.321 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1204 time to fit residues: 2.4713 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 160 random chunks: chunk 149 optimal weight: 0.0050 chunk 17 optimal weight: 5.9990 chunk 88 optimal weight: 10.0000 chunk 113 optimal weight: 0.1980 chunk 87 optimal weight: 0.9990 chunk 130 optimal weight: 6.9990 chunk 86 optimal weight: 1.9990 chunk 154 optimal weight: 6.9990 chunk 96 optimal weight: 4.9990 chunk 94 optimal weight: 9.9990 chunk 71 optimal weight: 9.9990 overall best weight: 1.6400 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 135 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** E 59 GLN ** G 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** I 59 GLN ** K 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** M 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7762 moved from start: 0.2596 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.023 12424 Z= 0.177 Angle : 0.417 8.723 16988 Z= 0.221 Chirality : 0.043 0.124 2064 Planarity : 0.003 0.024 2095 Dihedral : 5.303 37.986 1922 Min Nonbonded Distance : 2.278 Molprobity Statistics. All-atom Clashscore : 2.77 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.30 % Favored : 97.70 % Rotamer Outliers : 0.39 % 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.64 (0.22), residues: 1607 helix: 1.54 (0.20), residues: 731 sheet: 1.75 (0.29), residues: 257 loop : 0.35 (0.28), residues: 619 *********************** REFINEMENT MACRO_CYCLE 7 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 142 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 138 time to evaluate : 1.435 Fit side-chains outliers start: 4 outliers final: 3 residues processed: 139 average time/residue: 0.2109 time to fit residues: 45.0940 Evaluate side-chains 139 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 136 time to evaluate : 1.378 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1187 time to fit residues: 2.4963 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=2.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 160 random chunks: chunk 95 optimal weight: 8.9990 chunk 61 optimal weight: 9.9990 chunk 92 optimal weight: 2.9990 chunk 46 optimal weight: 10.0000 chunk 30 optimal weight: 1.9990 chunk 29 optimal weight: 7.9990 chunk 98 optimal weight: 7.9990 chunk 105 optimal weight: 0.2980 chunk 76 optimal weight: 7.9990 chunk 14 optimal weight: 0.0050 chunk 121 optimal weight: 8.9990 overall best weight: 2.6600 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 135 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** G 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** K 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** O 54 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7787 moved from start: 0.2604 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.029 12424 Z= 0.260 Angle : 0.459 9.093 16988 Z= 0.242 Chirality : 0.044 0.129 2064 Planarity : 0.003 0.025 2095 Dihedral : 5.476 38.353 1922 Min Nonbonded Distance : 2.265 Molprobity Statistics. All-atom Clashscore : 2.90 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.11 % Favored : 96.89 % Rotamer Outliers : 0.58 % 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.43 (0.22), residues: 1607 helix: 1.39 (0.20), residues: 731 sheet: 1.62 (0.29), residues: 257 loop : 0.25 (0.28), residues: 619 *********************** REFINEMENT MACRO_CYCLE 8 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 146 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 6 poor density : 140 time to evaluate : 1.393 Fit side-chains outliers start: 6 outliers final: 3 residues processed: 142 average time/residue: 0.2202 time to fit residues: 47.3782 Evaluate side-chains 139 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 136 time to evaluate : 1.625 Switching outliers to nearest non-outliers outliers start: 3 outliers final: 0 residues processed: 3 average time/residue: 0.1109 time to fit residues: 2.6226 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 160 random chunks: chunk 140 optimal weight: 9.9990 chunk 147 optimal weight: 1.9990 chunk 134 optimal weight: 6.9990 chunk 143 optimal weight: 9.9990 chunk 86 optimal weight: 9.9990 chunk 62 optimal weight: 5.9990 chunk 112 optimal weight: 6.9990 chunk 44 optimal weight: 4.9990 chunk 129 optimal weight: 9.9990 chunk 135 optimal weight: 0.0770 chunk 94 optimal weight: 0.9980 overall best weight: 2.8144 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 134 GLN ** C 135 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** G 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** K 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** M 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** O 54 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** O 191 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7792 moved from start: 0.2671 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.004 0.029 12424 Z= 0.272 Angle : 0.467 8.837 16988 Z= 0.246 Chirality : 0.044 0.129 2064 Planarity : 0.003 0.024 2095 Dihedral : 5.573 38.890 1922 Min Nonbonded Distance : 2.260 Molprobity Statistics. All-atom Clashscore : 2.94 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.05 % Favored : 96.95 % 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: 1.31 (0.22), residues: 1607 helix: 1.33 (0.20), residues: 731 sheet: 1.51 (0.29), residues: 257 loop : 0.14 (0.27), residues: 619 *********************** REFINEMENT MACRO_CYCLE 9 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 146 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 141 time to evaluate : 1.444 Fit side-chains outliers start: 5 outliers final: 5 residues processed: 142 average time/residue: 0.2159 time to fit residues: 46.9280 Evaluate side-chains 144 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 139 time to evaluate : 1.544 Switching outliers to nearest non-outliers outliers start: 5 outliers final: 0 residues processed: 5 average time/residue: 0.1249 time to fit residues: 3.0820 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 160 random chunks: chunk 151 optimal weight: 10.0000 chunk 92 optimal weight: 0.3980 chunk 72 optimal weight: 5.9990 chunk 105 optimal weight: 3.9990 chunk 159 optimal weight: 2.9990 chunk 146 optimal weight: 5.9990 chunk 126 optimal weight: 0.9990 chunk 13 optimal weight: 3.9990 chunk 98 optimal weight: 0.7980 chunk 77 optimal weight: 6.9990 chunk 100 optimal weight: 8.9990 overall best weight: 1.8386 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 135 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** G 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** K 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** M 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** O 54 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7769 moved from start: 0.2734 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.003 0.025 12424 Z= 0.194 Angle : 0.435 9.171 16988 Z= 0.229 Chirality : 0.043 0.126 2064 Planarity : 0.003 0.024 2095 Dihedral : 5.354 38.120 1922 Min Nonbonded Distance : 2.269 Molprobity Statistics. All-atom Clashscore : 3.07 Ramachandran Plot: Outliers : 0.00 % Allowed : 2.74 % Favored : 97.26 % Rotamer Outliers : 0.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 1.45 (0.22), residues: 1607 helix: 1.47 (0.20), residues: 731 sheet: 1.66 (0.30), residues: 257 loop : 0.15 (0.28), residues: 619 ********************** REFINEMENT MACRO_CYCLE 10 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 3214 Ramachandran restraints generated. 1607 Oldfield, 0 Emsley, 1607 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 139 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 139 time to evaluate : 1.531 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 139 average time/residue: 0.2193 time to fit residues: 46.9686 Evaluate side-chains 137 residues out of total 1413 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 137 time to evaluate : 1.556 Switching outliers to nearest non-outliers outliers start: 0 outliers final: 0 residues processed: 0 time to fit residues: 2.0589 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=1.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 160 random chunks: chunk 135 optimal weight: 10.0000 chunk 38 optimal weight: 1.9990 chunk 117 optimal weight: 7.9990 chunk 18 optimal weight: 1.9990 chunk 35 optimal weight: 5.9990 chunk 127 optimal weight: 4.9990 chunk 53 optimal weight: 9.9990 chunk 130 optimal weight: 4.9990 chunk 16 optimal weight: 6.9990 chunk 23 optimal weight: 9.9990 chunk 111 optimal weight: 8.9990 overall best weight: 3.9990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** C 135 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** E 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** G 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** K 42 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** ** O 54 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3162 r_free = 0.3162 target = 0.109157 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 49)----------------| | r_work = 0.2964 r_free = 0.2964 target = 0.095126 restraints weight = 16683.198| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 25)----------------| | r_work = 0.2985 r_free = 0.2985 target = 0.096579 restraints weight = 12731.583| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 43)----------------| | r_work = 0.3006 r_free = 0.3006 target = 0.097945 restraints weight = 10478.613| |-----------------------------------------------------------------------------| r_work (final): 0.2992 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8100 moved from start: 0.2731 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd. rmsZ for bonds and angles. Bond : 0.006 0.035 12424 Z= 0.372 Angle : 0.523 9.354 16988 Z= 0.273 Chirality : 0.046 0.134 2064 Planarity : 0.003 0.025 2095 Dihedral : 5.781 39.548 1922 Min Nonbonded Distance : 2.248 Molprobity Statistics. All-atom Clashscore : 3.07 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.61 % Favored : 96.39 % Rotamer Outliers : 0.29 % 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.06 (0.22), residues: 1607 helix: 1.17 (0.20), residues: 730 sheet: 1.29 (0.30), residues: 257 loop : 0.00 (0.27), residues: 620 =============================================================================== Job complete usr+sys time: 1886.06 seconds wall clock time: 35 minutes 27.39 seconds (2127.39 seconds total)