Conformational Flexibility of the Benzyl-Guanine Adduct in a Bypass Polymerase Active Site Permits Replication: Insights from Molecular Dynamics Simulations

Structural Rationalization for the Nonmutagenic and Mutagenic Bypass of the Tobacco-Derived O4-4-(3-Pyridyl)-4-oxobut-1-yl-thymine Lesion by Human Polymerase η: A Multiscale Computational Study

Tobacco-derived pyridyloxobutyl (POB) DNA adducts are unique due to the large size and flexibility of the alkyl chain connecting the pyridyl ring to the nucleobase. Recent experimental work suggests that the O4-4-(3-pyridyl)-4-oxobut-1-yl-T (O4-POB-T) lesion can undergo both nonmutagenic (dATP) and mutagenic (dGTP) insertion by the translesion synthesis (TLS) polymerase (pol) η in human cells. Interestingly, the mutagenic rate for O4-POB-T replication is reduced compared to that for the smaller O4-methylthymine (O4-Me-T) lesion, and O4-POB-T yields a different mutagenic profile than the O2-POB-T variant (dTTP insertion). The present work uses a combination of density functional theory calculations and molecular dynamics simulations to probe the impact of the size and flexibility of O4-POB-T on pol η replication outcomes. Due to changes in the Watson-Crick binding face upon damage of canonical T, O4-POB-T does not form favorable hydrogen-bonding interactions with A. Nevertheless, dATP is positioned for insertion in the pol η active site by a water chain to the template strand, which suggests a pol η replication pathway similar to that for abasic sites. Although a favorable O4-POB-T:G mispair forms in the pol η active site and DNA duplexes, the inherent dynamical nature of O4-POB-T periodically disrupts interstrand hydrogen bonding that would otherwise facilitate dGTP insertion and stabilize damaged DNA duplexes. In addition to explaining the origin of the experimentally reported pol η outcomes associated with O4-POB-T replication, comparison to structural data for the O4-Me-T and O2-POB-T adducts highlights an emerging common pathway for the nonmutagenic replication of thymine alkylated lesions by pol η, yet underscores the broader impacts of bulky moiety size, flexibility, and position on the associated mutagenic outcomes.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

The Effect of Cofactor Binding on the Conformational Plasticity of the Biological Receptors in Artificial Metalloenzymes: The Case Study of LmrR

The design of Artificial Metalloenzymes (ArMs), which result from the incorporation of organometallic cofactors into biological structures, has grown steadily in the last two decades and important new-to-Nature reactions have been reached. These type of exercises could greatly benefit from an understanding of the structural impact that the inclusion of organometallic moieties may have on the biological host. To date though, our understanding of this phenomenon is highly partial. This lack of knowledge is one of the elements that condition that first-generation ArMs generally display relatively poor catalytic profiles. In this work, we approach this matter by assessing the dynamics and stability of a series of ArMs resulting from the inclusion, via different anchoring strategies, of a variety of organometallic cofactors into the Lactococcal multidrug resistance regulator (LmrR) protein. To this aim, we coupled standard force field-based techniques such as Protein-Ligand Docking and Molecular Dynamics simulations with a variety of trajectory convergence analyses, capable of assessing both the stability and flexibility of the different systems under study upon the binding of cofactors. Together with the experimental evidence obtained in other studies, we provide an overview on how these changes can affect the catalytic outcomes obtained from the different ArMs. Fundamentally, our results show that the convergence analysis used in this work can assess how the inclusion of synthetic metallic cofactors in proteins can condition different structural modulations of their host. Those conformational modifications are key to the success of the desired catalytic activity and their proper identification can be wisely used to improve the quality and the rate of success of the ArMs.

Uncovering a unique approach for damaged DNA replication: A computational investigation of a mutagenic tobacco-derived thymine lesion

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone is a potent nicotine carcinogen that leads to many DNA lesions, the most persistent being the O2-[4-oxo-4-(3-pyridyl)butyl]thymine adduct (POB-T). Although the experimental mutagenic profile for the minor groove POB-T lesion has been previously reported, the findings are puzzling in terms of the human polymerases involved. Specifically, while pol κ typically replicates minor groove adducts, in vivo studies indicate pol η replicates POB-T despite being known for processing major groove adducts. Our multiscale modeling approach reveals that the canonical (anti) glycosidic orientation of POB-T can fit in the pol κ active site, but only a unique (syn) POB-T conformation is accommodated by pol η. These distinct binding orientations rationalize the differential in vitro mutagenic spectra based on the preferential stabilization of dGTP and dTTP opposite the lesion for pol κ and η, respectively. Overall, by uncovering the first evidence for the replication of a damaged pyrimidine in the syn glycosidic orientation, the current work provides the insight necessary to clarify a discrepancy in the DNA replication literature, expand the biological role of the critical human pol η, and understand the mutational signature in human cancers associated with tobacco exposure.