No abstract

Specialized DNA damage-bypass Y-family DNA polymerases contribute to cancer prevention by providing cellular tolerance to DNA damage that can lead to mutations and contribute to cancer progression by increasing genomic instability. Y-family polymerases can also bypass DNA adducts caused by chemotherapy agents. One of the four human Y-family DNA polymerases, DNA polymerase (pol) κ, has been shown to be specific for bypass of minor groove adducts and inhibited by major groove adducts. In addition, mutations in the gene encoding pol κ are associated with different types of cancers as well as with chemotherapy responses. We characterized nine variants of pol κ whose identity was inferred from cancer-associated single nucleotide polymorphisms for polymerization activity on undamaged and damaged DNA, their abilities to extend from mismatched or damaged base pairs at primer termini, and overall stability and dynamics. We find that these pol κ variants generally fall into three categories: similar activity to wild-type (WT) pol κ (L21F, I39T, P169T, F192C, and E292K), more active than WT pol κ (S423R), and less active than pol κ (R219I, R298H, and Y432S). Of these, only pol κ variants R298H and Y432S had markedly reduced thermal stability. Molecular dynamics (MD) simulations with undamaged DNA revealed that the active variant F192C and more active variant S423R with either correct or incorrect incoming nucleotide mimic WT pol κ with the correct incoming nucleotide, whereas the less active variants R219I, R298H, and Y432S with the correct incoming nucleotide mimic WT pol κ with the incorrect incoming nucleotide. Thus, the observations from MD simulations suggest a possible explanation for the observed experimental results that pol κ adopts specific active and inactive conformations that depend on both the protein variant and the identity of the DNA adduct.

Interactions in enzymes between catalytic and neighboring amino acids and how these interactions facilitate catalysis are examined. In examples from both natural and designed enzymes, it is shown that increases in catalytic rates may be achieved through elongation of the buffer range of the catalytic residues; such perturbations in the protonation equilibria are, in turn, achieved through enhanced coupling of the protonation equilibria of the active ionizable residues with those of other ionizable residues. The strongest coupling between protonation states for a pair of residues that deprotonate to form an anion (or a pair that accept a proton to form a cation) is achieved when the difference in the intrinsic pKas of the two residues is approximately within 1 pH unit. Thus, catalytic aspartates and glutamates are often coupled to nearby acidic residues. For an anion-forming residue coupled to a cation-forming residue, the elongated buffer range is achieved when the intrinsic pKa of the anion-forming residue is higher than the intrinsic pKa of the (conjugate acid of the) cation-forming residue. Therefore, the high pKa, anion-forming residues tyrosine and cysteine make good coupling partners for catalytic lysine residues. For the anion–cation pairs, the optimum difference in intrinsic pKas is a function of the energy of interaction between the residues. For the energy of interaction ε expressed in units of (ln 10)RT, the optimum difference in intrinsic pKas is within ∼1 pH unit of ε.

The process of DNA replication is carried out with high efficiency and accuracy by DNA polymerases. The replicative polymerase in E. coli is DNA Pol III, which is a complex of 10 different subunits that coordinates simultaneous replication on the leading and lagging strands. The 1160-residue Pol III alpha subunit is responsible for the polymerase activity and copies DNA accurately, making one error per 10 nucleotide incorporations. The goal of this research is to determine the residues that contribute to the activity of the polymerase subunit. Homology modeling and the computational methods of THEMATICS and POOL were used to predict functionally important amino acid residues through their computed chemical properties. Site-directed mutagenesis and biochemical assays were used to validate these predictions. Primer extension, steady-state single-nucleotide incorporation kinetics, and thermal denaturation assays were performed to understand the contribution of these residues to the function of the polymerase. This work shows that the top 15 residues predicted by POOL, a set that includes the three previously known catalytic aspartate residues, seven remote residues, plus five previously unexplored first-layer residues, are important for function. Six previously unidentified residues, R362, D405, K553, Y686, E688, and H760, are each essential to Pol III activity; three additional residues, Y340, R390, and K758, play important roles in activity.

The roles of local interactions in the laboratory evolution of a highly active, computationally designed retroaldolase (RA) are examined. Partial Order Optimum Likelihood (POOL) is used to identify catalytically important amino acid interactions in several RA95 enzyme variants. The series RA95.5, RA95.5-5, RA95.5-8, and RA95.5-8F, representing progress along an evolutionary trajectory with increasing activity, is examined. Computed measures of coupling between charged states of residues show that, as evolution proceeds and higher activities are achieved, electrostatic coupling between the biochemically active amino acids and other residues is increased. In silico residue scanning suggests multiple coupling partners for the catalytic lysine K83. The effects of two predicted partners, Y51 and E85, are tested using site-directed mutagenesis and kinetic analysis of the variants Y51F and E85Q. The Y51F variants show decreases in k cat relative to wild type, with the greatest losses observed for the more evolved constructs; they also exhibit significant decreases in k cat /K M across the series. Only modest decreases in k cat /K M are observed for the E85Q variants with little effect on k cat . Computed metrics of the degree of coupling between protonation states rise significantly as evolution proceeds and catalytic turnover rate increases. Specifically, the charge state of the catalytic lysine K83 becomes more strongly coupled to those of other amino acids as the enzyme evolves to a better catalyst.