Hypothesis driven single nucleotide polymorphism search (HyDn-SNP-S)

Effects of the Y432S Cancer-Associated Variant on the Reaction Mechanism of Human DNA Polymerase κ

Human DNA polymerases are vital for genetic information management. Their function involves catalyzing the synthesis of DNA strands with unparalleled accuracy, which ensures the fidelity and stability of the human genomic blueprint. Several disease-associated mutations and their functional impact on DNA polymerases have been reported. One particular polymerase, human DNA polymerase kappa (Pol κ), has been reported to be susceptible to several cancer-associated mutations. The Y432S mutation in Pol κ, associated with various cancers, is of interest due to its impact on polymerization activity and markedly reduced thermal stability. Here, we have used computational simulations to investigate the functional consequences of the Y432S using classical molecular dynamics (MD) and coupled quantum mechanics/molecular mechanics (QM/MM) methods. Our findings suggest that Y432S induces structural alterations in domains responsible for nucleotide addition and ternary complex stabilization while retaining structural features consistent with possible catalysis in the active site. Calculations of the minimum energy path associated with the reaction mechanism of the wild type (WT) and Y432S Pol κ indicate that, while both enzymes are catalytically competent (in terms of energetics and the active site's geometries), the cancer mutation results in an endoergic reaction and an increase in the catalytic barrier. Interactions with a third magnesium ion and environmental effects on nonbonded interactions, particularly involving key residues, contribute to the kinetic and thermodynamic distinctions between the WT and mutant during the catalytic reaction. The energetics and electronic findings suggest that active site residues favor the catalytic reaction with dCTP3- over dCTP4-.

Computational compensatory mutation discovery approach: Predicting a PARP1 variant rescue mutation

The prediction of protein mutations that affect function may be exploited for multiple uses. In the context of disease variants, the prediction of compensatory mutations that reestablish functional phenotypes could aid in the development of genetic therapies. In this work, we present an integrated approach that combines coevolutionary analysis and molecular dynamics (MD) simulations to discover functional compensatory mutations. This approach is employed to investigate possible rescue mutations of a poly(ADP-ribose) polymerase 1 (PARP1) variant, PARP1 V762A, associated with lung cancer and follicular lymphoma. MD simulations show PARP1 V762A exhibits noticeable changes in structural and dynamical behavior compared with wild-type (WT) PARP1. Our integrated approach predicts A755E as a possible compensatory mutation based on coevolutionary information, and molecular simulations indicate that the PARP1 A755E/V762A double mutant exhibits similar structural and dynamical behavior to WT PARP1. Our methodology can be broadly applied to a large number of systems where single-nucleotide polymorphisms have been identified as connected to disease and can shed light on the biophysical effects of such changes as well as provide a way to discover potential mutants that could restore WT-like functionality. This can, in turn, be further utilized in the design of molecular therapeutics that aim to mimic such compensatory effect.

Towards a comprehensive view of 8-oxo-7,8-dihydro-2'-deoxyguanosine: Highlighting the intertwined roles of DNA damage and epigenetics in genomic instability

8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG), a major product of DNA oxidation, is a pre-mutagenic lesion which is prone to mispair, if left unrepaired, with 2'-deoxyadenosine during DNA replication. While unrepaired or incompletely repaired 8-oxodG has classically been associated with genome instability and cancer, it has recently been reported to have a role in the epigenetic regulation of gene expression. Despite the growing collection of genome-wide 8-oxodG mapping studies that have been used to provide new insight on the functional nature of 8-oxodG within the genome, a comprehensive view that brings together the epigenetic and the mutagenic nature of the 8-oxodG is still lacking. To help address this gap, this review aims to provide (i) a description of the state-of-the-art knowledge on both the mutagenic and epigenetic roles of 8-oxodG; (ii) putative molecular models through which the 8-oxodG can cause genome instability; (iii) a possible molecular model on how 8-oxodG, acting as an epigenetic signal, could cause the translocations and deletions which are associated with cancer.

Biophysical and Mechanistic Models for Disease-Causing Protein Variants

The rapid decrease in DNA sequencing cost is revolutionizing medicine and science. In medicine, genome sequencing has revealed millions of missense variants that change protein sequences, yet we only understand the molecular and phenotypic consequences of a small fraction. Within protein science, high-throughput deep mutational scanning experiments enable us to probe thousands of variants in a single, multiplexed experiment. We review efforts that bring together these topics via experimental and computational approaches to determine the consequences of missense variants in proteins. We focus on the role of changes in protein stability as a driver for disease, and how experiments, biophysical models, and computation are providing a framework for understanding and predicting how changes in protein sequence affect cellular protein stability.

DNArCdb: A database of cancer biomarkers in DNA repair genes that includes variants related to multiple cancer phenotypes

Functioning DNA repair capabilities are vital for organisms to ensure that the biological information is preserved and correctly propagated. Disruptions in DNA repair pathways can result in the accumulation of DNA mutations, which may lead to onset of complex disease such as cancer. The discovery and characterization of cancer-related biomarkers may allow early diagnosis and targeted treatment, which could significantly contribute to the survival rates of cancer patients. To this end, we have applied a hypothesis driven bioinformatics approach to identify biomarkers related to 25 different DNA repair enzymes, in combination with structural analysis of six selected missense mutations of newly discovered SNPs that are associated with cancer phenotypes. Our search on 8 distinct cancer databases uncovered 43 missense SNPs that statistically significantly associated at least one phenotype. Moreover, nine of these missense SNPs are statistically significantly associated with two or more cancers. In addition, we have performed classical molecular dynamics to characterize the impact of rs10018786 on POLN, which results in the M310 L Pol ν variant, and rs3218784 on POLI, which results in the I236 M Pol ι. Our results suggest that both of these cancer-associated variants result in noticeable structural and dynamical changes compared with their respective wild-type proteins.

Characterization of Nine Cancer-Associated Variants in Human DNA Polymerase κ

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.

Not breathing is not an option: How to deal with oxidative DNA damage

Oxidative DNA damage constitutes a major threat to genetic integrity, and has thus been implicated in the pathogenesis of a wide variety of diseases, including cancer and neurodegeneration. 7,8-dihydro-8oxo-deoxyGuanine (8-oxo-G) is one of the best characterised oxidative DNA lesions, and it can give rise to point mutations due to its miscoding potential that instructs most DNA polymerases (Pols) to preferentially insert Adenine (A) opposite 8-oxo-G instead of the correct Cytosine (C). If uncorrected, A:8-oxo-G mispairs can give rise to C:G→A:T transversion mutations. Cells have evolved a variety of pathways to mitigate the mutational potential of 8-oxo-G that include i) mechanisms to avoid incorporation of oxidized nucleotides into DNA through nucleotide pool sanitisation enzymes (by MTH1, MTH2, MTH3 and NUDT5), ii) base excision repair (BER) of 8-oxo-G in DNA (involving MUTYH, OGG1, Pol λ, and other components of the BER machinery), and iii) faithful bypass of 8-oxo-G lesions during replication (using a switch between replicative Pols and Pol λ). In the following, the fate of 8-oxo-G in mammalian cells is reviewed in detail. The differential origins of 8-oxo-G in DNA and its consequences for genetic stability will be covered. This will be followed by a thorough discussion of the different mechanisms in place to cope with 8-oxo-G with an emphasis on Pol λ-mediated correct bypass of 8-oxo-G during MUTYH-initiated BER as well as replication across 8-oxo-G. Furthermore, the multitude of mechanisms in place to regulate key proteins involved in 8-oxo-G repair will be reviewed. Novel functions of 8-oxo-G as an epigenetic-like regulator and insights into the repair of 8-oxo-G within the cellular context will be touched upon. Finally, a discussion will outline the relevance of 8-oxo-G and the proteins involved in dealing with 8-oxo-G to human diseases with a special emphasis on cancer.

Computational Simulations of DNA Polymerases: Detailed Insights on Structure/Function/Mechanism from Native Proteins to Cancer Variants

Genetic information is vital in the cell cycle of DNA-based organisms. DNA polymerases (DNA Pols) are crucial players in transactions dealing with these processes. Therefore, the detailed understanding of the structure, function, and mechanism of these proteins has been the focus of significant effort. Computational simulations have been applied to investigate various facets of DNA polymerase structure and function. These simulations have provided significant insights over the years. This perspective presents the results of various computational studies that have been employed to research different aspects of DNA polymerases including detailed reaction mechanism investigation, mutagenicity of different metal cations, possible factors for fidelity synthesis, and discovery/functional characterization of cancer-related mutations on DNA polymerases.

Polymorphisms in POLG were associated with the prognosis and mtDNA content in hepatocellular carcinoma patients

BACKGROUND

POLG is a gene that codes for the catalytic subunit of the mitochondrial DNA polymerase, which is involved in the replication of mitochondrial DNA. Genetic variants in mitochondrial DNA polymerase-γ (POLG) have been associated with several malignancies. However, as an important metabolic tissue, association between genetic polymorphisms of POLG and the prognosis and mitochondrial DNA (mtDNA) content in hepatocellular carcinoma (HCC) remains unknown. Here we investigated the association between in POLG with the prognosis and mitochondrial DNA (mtDNA) content in hepatocellular carcinoma (HCC).

METHODS

Three nucleotide polymorphisms (SNPs) of rs1061316, rs2247233 and rs758130 in POLG were examined in 416 patients from two cohorts undergoing transcatheter arterial chemoembolization treatment. Leukocyte mtDNA content from 216 patients in cohort 2 was measured using a real-time PCR-based method. The association of SNPs with prognosis and of mtDNA content of patients was analyzed.

RESULTS

The rs758130 in POLG gene was significantly associated with the prognosis of patients in a dose-dependent manner. Moreover, GG genotype in rs1061316 showed significantly high mtDNA content, an indicator of better prognosis.

CONCLUSIONS

Our study for the first time demonstrates that rs1061316 and rs758130 in POLG is associated with the prognosis and leukocyte mtDNA content in HCC patients.

ALKBH7 Variant Related to Prostate Cancer Exhibits Altered Substrate Binding

The search for prostate cancer biomarkers has received increased attention and several DNA repair related enzymes have been linked to this dysfunction. Here we report a targeted search for single nucleotide polymorphisms (SNPs) and functional impact characterization of human ALKBH family dioxygenases related to prostate cancer. Our results uncovered a SNP of ALKBH7, rs7540, which is associated with prostate cancer disease in a statistically significantly manner in two separate cohorts, and maintained in African American men. Comparisons of molecular dynamics (MD) simulations on the wild-type and variant protein structures indicate that the resulting alteration in the enzyme induces a significant structural change that reduces ALKBH7’s ability to bind its cosubstrate. Experimental spectroscopy studies with purified proteins validate our MD predictions and corroborate the conclusion that this cancer-associated mutation affects productive cosubstrate binding in ALKBH7.

Improvements in personalized DNA sequencing have led to an increased interest in targeted biomarkers for therapeutic and diagnostic purposes. In this work, we report on a new biomarker for prostate cancer found through a targeted search for single nucleotide polymorphisms (SNPs) of the genes encoding human ALKBH family dioxygenases. Our results uncovered rs7540, which leads to a missense mutation in ALKBH7. Comparative molecular dynamics simulations on the wild type and SNP variant of the protein show that the mutation elicits a structural change that dramatically decreases ALKBH7’s affinity for its cosubstrate. This prediction is confirmed by experimental UV-Vis spectroscopy. Taken together, these results give important insights into a novel prostate-cancer related SNP and its impact on the structure and function of ALKBH7.

Estrogen Drives Cellular Transformation and Mutagenesis in Cells Expressing the Breast Cancer-Associated R438W DNA Polymerase Lambda Protein

Repair of DNA damage is critical for maintaining the genomic integrity of cells. DNA polymerase lambda (POLL/Pol λ) is suggested to function in base excision repair (BER) and nonhomologous end-joining (NHEJ), and is likely to play a role in damage tolerance at the replication fork. Here, using next-generation sequencing, it was discovered that the POLL rs3730477 single-nucleotide polymorphism (SNP) encoding R438W Pol λ was significantly enriched in the germlines of breast cancer patients. Expression of R438W Pol λ in human breast epithelial cells induces cellular transformation and chromosomal aberrations. The role of estrogen was assessed as it is commonly used in hormone replacement therapies and is a known breast cancer risk factor. Interestingly, the combination of estrogen treatment and the expression of the R438W Pol λ SNP drastically accelerated the rate of transformation. Estrogen exposure produces 8-oxoguanine lesions that persist in cells expressing R438W Pol λ compared with wild-type (WT) Pol λ-expressing cells. Unlike WT Pol λ, which performs error-free bypass of 8-oxoguanine lesions, expression of R438W Pol λ leads to an increase in mutagenesis and replicative stress in cells treated with estrogen. Together, these data suggest that individuals who carry the rs3730477 POLL germline variant have an increased risk of estrogen-associated breast cancer.

IMPLICATIONS

The Pol λ R438W mutation can serve as a biomarker to predict cancer risk and implicates that treatment with estrogen in individuals with this mutation may further increase their risk of breast cancer. Mol Cancer Res; 14(11); 1068-77. ©2016 AACR.