Stereoselective Single Step Cyclization to Give Belt-Functionalized Pillar[6]arenes.

Two macrocycles were synthesized through cyclization reactions of secondary benzylic alcohols, giving pillar[6]arenes with a methyl substituent at each belt position. These macrocycles form stereoselectively with only the rtctct isomer with alternating up and down orientations of the belt methyl groups definitively identified. Isolated yields were modest (7 and 9%), but the macrocycles are prepared in a single step from either a commercially available alcohol or a very readily prepared precursor. X-ray crystal structures of the macrocycles indicate they have a capsule-like structure, which is far from the conventional pillar shape. Density functional theory calculations reveal that the energy barrier required to obtain the pillar conformation is significantly higher for these belt-functionalized macrocycles than for conventional belt-unfunctionalized pillar[6]arenes.

Anatomy of noncovalent interactions between the nucleobases or ribose and π-containing amino acids in RNA-protein complexes.

A set of >300 nonredundant high-resolution RNA-protein complexes were rigorously searched for π-contacts between an amino acid side chain (W, H, F, Y, R, E and D) and an RNA nucleobase (denoted π-π interaction) or ribose moiety (denoted sugar-π). The resulting dataset of >1500 RNA-protein π-contacts were visually inspected and classified based on the interaction type, and amino acids and RNA components involved. More than 80% of structures searched contained at least one RNA-protein π-interaction, with π-π contacts making up 59% of the identified interactions. RNA-protein π-π and sugar-π contacts exhibit a range in the RNA and protein components involved, relative monomer orientations and quantum mechanically predicted binding energies. Interestingly, π-π and sugar-π interactions occur more frequently with RNA (4.8 contacts/structure) than DNA (2.6). Moreover, the maximum stability is greater for RNA-protein contacts than DNA-protein interactions. In addition to highlighting distinct differences between RNA and DNA-protein binding, this work has generated the largest dataset of RNA-protein π-interactions to date, thereby underscoring that RNA-protein π-contacts are ubiquitous in nature, and key to the stability and function of RNA-protein complexes.

The allosteric inhibition of glycine transporter 2 by bioactive lipid analgesics is controlled by penetration into a deep lipid cavity.

The role of lipids in modulating membrane protein function is an emerging and rapidly growing area of research. The rational design of lipids that target membrane proteins for the treatment of pathological conditions is a novel extension in this field and provides a step forward in our understanding of membrane transporters. Bioactive lipids show considerable promise as analgesics for the treatment of chronic pain and bind to a high-affinity allosteric-binding site on the human glycine transporter 2 (GlyT2 or SLC6A5). Here, we use a combination of medicinal chemistry, electrophysiology, and computational modeling to develop a rational structure-activity relationship for lipid inhibitors and demonstrate the key role of the lipid tail interactions for GlyT2 inhibition. Specifically, we examine how lipid inhibitor head group stereochemistry, tail length, and double-bond position promote enhanced inhibition. Overall, the l-stereoisomer is generally a better inhibitor than the d-stereoisomer, longer tail length correlates with greater potency, and the position of the double bond influences the activity of the inhibitor. We propose that the binding of the lipid inhibitor deep into the allosteric-binding pocket is critical for inhibition. Furthermore, this provides insight into the mechanism of inhibition of GlyT2 and highlights how lipids can modulate the activity of membrane proteins by binding to cavities between helices. The principles identified in this work have broader implications for the development of a larger class of compounds that could target SLC6 transporters for disease treatment.

Multiscale computational investigations of the translesion synthesis bypass of tobacco-derived DNA adducts: critical insights that complement experimental biochemical studies.

Among the numerous agents that damage DNA, tobacco products remain one of the most lethal and result in the most diverse set of DNA lesions. This perspective aims to provide an overview of computational work conducted to complement experimental biochemical studies on the mutagenicity of adducts derived from the most potent tobacco carcinogen, namely 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (nicotine-derived nitrosaminoketone or NNK). Lesions ranging from the smallest methylated thymine derivatives to the larger, flexible pyridyloxobutyl (POB) guanine adducts are considered. Insights are obtained from density functional theory (DFT) calculations and molecular dynamics (MD) simulations into the damaged nucleobase and nucleoside structures, the accommodation of the lesions in the active site of key human polymerases, the intrinsic base pairing potentials of the adducts, and dNTP incorporation opposite the lesions. Overall, the computational data provide atomic level information that can rationalize the differential mutagenic properties of tobacco-derived lesions and uncover important insights into the impact of adduct size, nucleobase, position, and chemical composition of the bulky moiety.

High resolution crystal structure of a KRAS promoter G-quadruplex reveals a dimer with extensive poly-A π-stacking interactions for small-molecule recognition.

Aberrant KRAS signaling is a driver of many cancers and yet remains an elusive target for drug therapy. The nuclease hypersensitive element of the KRAS promoter has been reported to form secondary DNA structures called G-quadruplexes (G4s) which may play important roles in regulating KRAS expression, and has spurred interest in structural elucidation studies of the KRAS G-quadruplexes. Here, we report the first high-resolution crystal structure (1.6 Å) of a KRAS G-quadruplex as a 5'-head-to-head dimer with extensive poly-A π-stacking interactions observed across the dimer. Molecular dynamics simulations confirmed that the poly-A π-stacking interactions are also maintained in the G4 monomers. Docking and molecular dynamics simulations with two G4 ligands that display high stabilization of the KRAS G4 indicated the poly-A loop was a binding site for these ligands in addition to the 5'-G-tetrad. Given sequence and structural variability in the loop regions provide the opportunity for small-molecule targeting of specific G4s, we envisage this high-resolution crystal structure for the KRAS G-quadruplex will aid in the rational design of ligands to selectively target KRAS.

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.

The role of plasmalogens, Forssman lipids, and sphingolipid hydroxylation in modulating the biophysical properties of the epithelial plasma membrane.

A coarse-grain model of the epithelial plasma membrane was developed from high-resolution lipidomic data and simulated using the MARTINI force field to characterize its biophysical properties. Plasmalogen lipids, Forssman glycosphingolipids, and hydroxylated Forssman glycosphingolipids and sphingomyelin were systematically added to determine their structural effects. Plasmalogen lipids have a minimal effect on the overall biophysical properties of the epithelial plasma membrane. In line with the hypothesized role of Forssman lipids in the epithelial apical membrane, the introduction of Forssman lipids initiates the formation of glycosphingolipid-rich nanoscale lipid domains, which also include phosphatidylethanolamine (PE), sphingomyelin (SM), and cholesterol (CHOL). This decreases the lateral diffusion in the extracellular leaflet, as well as the area per lipid of domain forming lipids, most notably PE. Finally, hydroxylation of the Forssman glycosphingolipids and sphingomyelin further modulates the lateral organization of the membrane. Through comparison to the previously studied average and neuronal plasma membranes, the impact of membrane lipid composition on membrane properties was characterized. Overall, this study furthers our understanding of the biophysical properties of complex membranes and the impact of lipid diversity in modulating membrane properties.

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.

Understanding the Link between Lipid Diversity and the Biophysical Properties of the Neuronal Plasma Membrane.

Cell membranes contain incredible diversity in the chemical structures of their individual lipid species and the ratios in which these lipids are combined to make membranes. Nevertheless, our current understanding of how each of these components affects the properties of the cell membrane remains elusive, in part due to the difficulties in studying the dynamics of membranes at high spatiotemporal resolution. In this work, we use coarse-grained molecular dynamics simulations to investigate how individual lipid species contribute to the biophysical properties of the neuronal plasma membrane. We progress through eight membranes of increasing chemical complexity, ranging from a simple POPC/CHOL membrane to a previously published neuronal plasma membrane [Ingólfsson, H. I., et al. (2017) Biophys. J. 113 (10), 2271-2280] containing 49 distinct lipid species. Our results show how subtle chemical changes can affect the properties of the membrane and highlight the lipid species that give the neuronal plasma membrane its unique biophysical properties. This work has potential far-reaching implications for furthering our understanding of cell membranes.

Computational insights into the mutagenicity of two tobacco-derived carcinogenic DNA lesions.

4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone is a potent carcinogen found in all tobacco products that leads to a variety of DNA lesions in cells, including O6-[4-oxo-4-(3-pyridyl)butyl]guanine (POB-G) and O6-[4-hydroxy-4-(3-pyridyl)butyl]guanine (PHB-G), which differ by only a single substituent in the bulky moiety. This work uses a multiscale computational approach to shed light on the intrinsic conformational and base-pairing preferences of POB-G and PHB-G, and the corresponding properties in DNA and the polymerase η active site. Our calculations reveal that both lesions form stable pairs with C and T, with the T pairs being the least distorted relative to canonical DNA. This rationalizes the experimentally reported mutational profile for POB-G and validates our computational model. The same approach predicts that PHB-G is more mutagenic than POB-G due to a difference in the bulky moiety hydrogen-bonding pattern, which increases the stability of the PHB-G:T pair. The mutagenicity of PHB-G is likely further increased by stabilization of an intercalated DNA conformation that is associated with deletion mutations. This work thereby uncovers structural explanations for the reported mutagenicity of POB-G, provides the first clues regarding the mutagenicity of PHB-G and complements a growing body of literature highlighting that subtle chemical changes can affect the biological outcomes of DNA adducts.

DNA base sequence effects on bulky lesion-induced conformational heterogeneity during DNA replication.

4-Aminobiphenyl (ABP) and its structure analog 2-aminofluorene (AF) are well-known carcinogens. In the present work, an unusual sequence effect in the 5'-CTTCTG1G2TCCTCATTC-3' DNA duplex is reported for ABP- and AF-modified G. Specifically, the ABP modification at G1 resulted in a mixture of 67% major groove B-type (B) and 33% stacked (S) conformers, while at the ABP modification at G2 exclusively resulted in the B-conformer. The AF modification at G1 and G2 lead to 25%:75% and 83%:17% B:S population ratios, respectively. These differences in preferred conformation are due to an interplay between stabilizing (hydrogen bonding and stacking that is enhanced by lesion planarity) and destabilizing (solvent exposure) forces at the lesion site. Furthermore, while the B-conformer is a thermodynamic stabilizer and the S-conformer is a destabilizer in duplex settings, the situation is reversed at the single strands/double strands (ss/ds) junction. Specifically, the twisted biphenyl is a better stacker at the ss/ds junction than the coplanar AF. Therefore, the ABP modification leads to a stronger strand binding affinity of the ss/ds junction than the AF modification. Overall, the current work provides conformational insights into the role of sequence and lesion effects in modulating DNA replication.

Computational Insight into the Differential Mutagenic Patterns of O-Methylthymine Lesions.

Differential mutagenic patterns were recently reported for O-methylated thymine lesions, which indicate that O4-methylthymine (O4-Me-T) frequently leads to G misinsertions, whereas O2-methylthymine (O2-Me-T) is primarily nonmutagenic. The reasons for these differences are unclear since both lesions similarly alter the Watson-Crick binding face of T. To rationalize these replication outcomes at a molecular level, this work uses density functional theory calculations and molecular dynamics simulations to probe the lesion base-pairing properties as well as lesion accommodation by human polymerase η (pol η) and post-extension DNA duplexes. O4-Me-T forms two strong hydrogen bonds with an opposing G in the active site of pol η, which rationalizes the observed lesion mutagenicity. Nevertheless, dATP insertion opposite O4-Me-T can proceed through water-mediated hydrogen bonding, which is similar to the pathway previously proposed for pol η bypass of abasic sites and other T alkylation lesions. In contrast, the position of O2-Me-T in the pol η active site is dynamic due to the presence of the aberrant methyl group on the minor groove side of DNA. In fact, the experimental replication outcomes can only be rationalized when the syn glycosidic orientation of O2-Me-T is considered, which stabilizes the pre-insertion complex by placing the damage in the polymerase open pocket on the major groove side of DNA. Although dATP insertion can occur opposite syn-O2-Me-T through a water-mediated pathway similar to O4-Me-T replication, rotation about the glycosidic bond precludes a stable pol η ternary complex corresponding to dGTP insertion, which correlates with the reported nonmutagenic bypass of O2-Me-T. In addition to providing structural insights into the differential mutagenicity of methylated T adducts, our data highlight an emerging theme in the literature for the replication of pyrimidine alkylation products in noncanonical glycosidic orientations and sets the stage for future work on the replication of other alkylated lesions by TLS polymerases.

DFT and MD Studies of Formaldehyde-Derived DNA Adducts: Molecular-Level Insights into the Differential Mispairing Potentials of the Adenine, Cytosine, and Guanine Lesions.

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone is a potent nicotine-based carcinogen that generates many DNA lesions, including the HOCH2-C, HOCH2-G, and HOCH2-A hydroxymethyl adducts. Despite all lesions containing an altered exocyclic amino group, which allows the hydroxymethyl group to be directed away from the Watson-Crick binding face, only the most persistent adenine adduct is mutagenic. As a first step toward understanding this differential mutagenicity, density functional theory (DFT) and molecular dynamics (MD) simulations were used to gain atomic-level structural details of these DNA damage products. DFT calculations reveal that all three lesions exhibit conformational diversity. However, regardless of the hydroxymethyl-nucleobase orientation, both DFT and MD simulations highlight that HOCH2-C and HOCH2-G form pairs with the canonical complementary base (G and C, respectively) that are structural and energetically preferred over mispairs. In contrast, depending on the hydroxymethyl-nucleobase orientation, the Watson-Crick HOCH2-A:T pair can become significantly destabilized relative to undamaged A:T. As a result, HOCH2-A mispairs with G, C, and A are energetically accessible and maintain key geometrical features of canonical DNA. Overall, our data directly correlate with the reported differential mutagenicity of the hydroxylmethyl lesions and will encourage future studies to further uncover the cellular impact of the most persistent adenine lesion.

Topology of RNA-protein nucleobase-amino acid π-π interactions and comparison to analogous DNA-protein π-π contacts.

The present work analyzed 120 high-resolution X-ray crystal structures and identified 335 RNA-protein π-interactions (154 nonredundant) between a nucleobase and aromatic (W, H, F, or Y) or acyclic (R, E, or D) π-containing amino acid. Each contact was critically analyzed (including using a visual inspection protocol) to determine the most prevalent composition, structure, and strength of π-interactions at RNA-protein interfaces. These contacts most commonly involve F and U, with U:F interactions comprising one-fifth of the total number of contacts found. Furthermore, the RNA and protein π-systems adopt many different relative orientations, although there is a preference for more parallel (stacked) arrangements. Due to the variation in structure, the strength of the intermolecular forces between the RNA and protein components (as determined from accurate quantum chemical calculations) exhibits a significant range, with most of the contacts providing significant stability to the associated RNA-protein complex (up to -65 kJ mol(-1)). Comparison to the analogous DNA-protein π-interactions emphasizes differences in RNA- and DNA-protein π-interactions at the molecular level, including the greater abundance of RNA contacts and the involvement of different nucleobase/amino acid residues. Overall, our results provide a clearer picture of the molecular basis of nucleic acid-protein binding and underscore the important role of these contacts in biology, including the significant contribution of π-π interactions to the stability of nucleic acid-protein complexes. Nevertheless, more work is still needed in this area in order to further appreciate the properties and roles of RNA nucleobase-amino acid π-interactions in nature.

Molecular Insights into the Translesion Synthesis of Benzyl-Guanine from Molecular Dynamics Simulations: Structural Evidence of Mutagenic and Nonmutagenic Replication.

DNA can be damaged by many compounds in our environment, and the resulting damaged DNA is commonly replicated by translesion synthesis (TLS) polymerases. Because the mechanism and efficiency of TLS are affected by the type of DNA damage, obtaining information for a variety of DNA adducts is critical. However, there is no structural information for the insertion of a dNTP opposite an O6-dG adduct, which is a particularly harmful class of DNA lesions. We used molecular dynamics (MD) simulations to investigate structural and energetic parameters that dictate preferred dNTP insertion opposite O6-benzyl-guanine (Bz-dG) by DNA polymerase IV, a prototypical TLS polymerase. Specifically, MD simulations were completed on all possible ternary insertion complexes and ternary -1 base deletion complexes with different Bz-dG conformations. Our data suggests that the purines are unlikely to be inserted opposite anti- or syn-Bz-dG, and dTTP is unlikely to be inserted opposite syn-Bz-dG, because of changes in the active site conformation, including critical hydrogen-bonding interactions and/or reaction-ready parameters compared to natural dG replication. In contrast, a preserved active site conformation suggests that dCTP can be inserted opposite either anti- or syn-Bz-dG and dTTP can be inserted opposite anti-Bz-dG. This is the first structural explanation for the experimentally observed preferential insertion of dCTP and misincorporation of dTTP opposite Bz-dG. Furthermore, we provide atomic level insight into why Bz-dG replication does not lead to deletion mutations, which is in contrast with the replication outcomes of other adducts. These findings provide a basis for understanding the replication of related O6-dG adducts.

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

Although translesion synthesis (TLS) polymerases play key roles in replicating DNA that contains nucleobase addition products (adducts), there are many unknowns about their function. The present work gains indispensable structural insights from molecular dynamics simulations on the replication of O6-benzyl-guanine (Bz-dG) prior to bond formation during dCTP insertion opposite the adduct by Dpo4. When combined with previous X-ray crystal structures of the Bz-dG extension complex, molecular details are now available for each stage during a single TLS replication cycle for this carcinogenic lesion. Our calculations illustrate that Bz-dG preferentially adopts an intercalated bulky moiety orientation in the Dpo4 preinsertion complex, which stabilizes the complex through Bz-dG interactions with the previously replicated 3'-base pair and positions the carcinogenic group in the dNTP binding site. Nevertheless, the maintained inherent flexibility of Bz-dG due to a stark lack of interactions with the polymerase or template DNA allows the bulky moiety to adopt a major groove position opposite an incoming dCTP in an orientation that is conducive for the experimentally observed nonmutagenic bypass. Comparison of Bz-dG and canonical dG replication clarifies that the experimentally observed decrease in dCTP binding affinity and replication efficiency upon adduct formation is likely caused by a combination of factors, including the required template nucleotide conformational change and destabilized template-dCTP hydrogen bonding. Although additional aspects of the replication process, such as the impact of the adduct on the nucleotidyl-transfer reaction, may also be important for fully rationalizing experimental replication data and must be considered in future work, the present contribution emphasizes the importance of considering the effect of DNA damage on the early stages of the TLS process.

DNA-protein π-interactions in nature: abundance, structure, composition and strength of contacts between aromatic amino acids and DNA nucleobases or deoxyribose sugar.

Four hundred twenty-eight high-resolution DNA-protein complexes were chosen for a bioinformatics study. Although 164 crystal structures (38% of those searched) contained no interactions, 574 discrete π-contacts between the aromatic amino acids and the DNA nucleobases or deoxyribose were identified using strict criteria, including visual inspection. The abundance and structure of the interactions were determined by unequivocally classifying the contacts as either π-π stacking, π-π T-shaped or sugar-π contacts. Three hundred forty-four nucleobase-amino acid π-π contacts (60% of all interactions identified) were identified in 175 of the crystal structures searched. Unprecedented in the literature, 230 DNA-protein sugar-π contacts (40% of all interactions identified) were identified in 137 crystal structures, which involve C-H···π and/or lone-pair···π interactions, contain any amino acid and can be classified according to sugar atoms involved. Both π-π and sugar-π interactions display a range of relative monomer orientations and therefore interaction energies (up to -50 (-70) kJ mol(-1) for neutral (charged) interactions as determined using quantum chemical calculations). In general, DNA-protein π-interactions are more prevalent than perhaps currently accepted and the role of such interactions in many biological processes may yet to be uncovered.

Computational Evaluation of Nucleotide Insertion Opposite Expanded and Widened DNA by the Translesion Synthesis Polymerase Dpo4.

Expanded (x) and widened (y) deoxyribose nucleic acids (DNA) have an extra benzene ring incorporated either horizontally (xDNA) or vertically (yDNA) between a natural pyrimidine base and the deoxyribose, or between the 5- and 6-membered rings of a natural purine. Far-reaching applications for (x,y)DNA include nucleic acid probes and extending the natural genetic code. Since modified nucleobases must encode information that can be passed to the next generation in order to be a useful extension of the genetic code, the ability of translesion (bypass) polymerases to replicate modified bases is an active area of research. The common model bypass polymerase DNA polymerase IV (Dpo4) has been previously shown to successfully replicate and extend past a single modified nucleobase on a template DNA strand. In the current study, molecular dynamics (MD) simulations are used to evaluate the accommodation of expanded/widened nucleobases in the Dpo4 active site, providing the first structural information on the replication of (x,y)DNA. Our results indicate that the Dpo4 catalytic (palm) domain is not significantly impacted by the (x,y)DNA bases. Instead, the template strand is displaced to accommodate the increased C1'-C1' base-pair distance. The structural insights unveiled in the present work not only increase our fundamental understanding of Dpo4 replication, but also reveal the process by which Dpo4 replicates (x,y)DNA, and thereby will contribute to the optimization of high fidelity and efficient polymerases for the replication of modified nucleobases.

Complex conformational heterogeneity of the highly flexible O6-benzyl-guanine DNA adduct.

The conformational preference of the O6-benzyl-guanine (BzG) adduct was computationally examined using nucleoside, nucleotide, and DNA models, which provided critical information about the potential mutagenic consequences and toxicity of the BzG adduct in our cells. Substantial conformational flexibility of the BzG moiety, including rotation of the bulky group with respect to the base and the internal conformation of the bulk moiety, is seen in the nucleoside and nucleotide models. This large conformational flexibility suggests the conformation adopted by BzG is dependent on the local environment of the BzG adduct. Upon incorporation of the adduct into the DNA helix, the BzG conformational flexibility is maintained. The range of BzG conformations adopted in DNA likely arises due to a combination of the long and flexible (-CH2-) linker, the small adduct size, and the lack of discrete interactions between the bulky moiety and G. Because of the conformational flexibility of the adduct, many DNA conformations are observed for BzG adducted DNA, including those not previously reported in the literature, and thus, a modified nomenclature for adducted DNA conformations is presented. Furthermore, the preferred conformation of BzG adducted DNA is greatly dependent on a number of factors, including the pairing nucleotide, the discrete interactions in the helix, and the solvation of the benzyl moiety. These factors in turn lead to a complicated mutagenic and toxic profile that may invoke pairing with natural C, mispairs, or deletion mutations, which is supported by previously reported experimental biochemical studies. Despite this complex mutagenic profile, pairing with C leads to the most stable helical structure, which is the first combined structural and energetic explanation for experimental studies reporting a higher rate of C incorporation than any other nucleobase upon BzG replication.

Standard role for a conserved aspartate or more direct involvement in deglycosylation? An ONIOM and MD investigation of adenine-DNA glycosylase.

8-Oxoguanine (OG) is one of the most frequently occurring forms of DNA damage and is particularly deleterious since it forms a stable Hoogsteen base pair with adenine (A). The repair of an OG:A mispair is initiated by adenine-DNA glycosylase (MutY), which hydrolyzes the sugar-nucleobase bond of the adenine residue before the lesion is processed by other proteins. MutY has been proposed to use a two-part chemical step involving protonation of the adenine nucleobase, followed by SN1 hydrolysis of the glycosidic bond. However, differences between a recent (fluorine recognition complex, denoted as the FLRC) crystal structure and the structure on which most mechanistic conclusions have been based to date (namely, the lesion recognition complex or LRC) raise questions regarding the mechanism used by MutY and the discrete role of various active-site residues. The present work uses both molecular dynamics (MD) and quantum mechanical (ONIOM) models to compare the active-site conformational dynamics in the two crystal structures, which suggests that only the understudied FLRC leads to a catalytically competent reactant. Indeed, all previous computational studies on MutY have been initiated from the LRC structure. Subsequently, for the first time, various mechanisms are examined with detailed ONIOM(M06-2X:PM6) reaction potential energy surfaces (PES) based on the FLRC structure, which significantly extends the mechanistic picture. Specifically, our work reveals that the reaction proceeds through a different route than the commonly accepted mechanism and the catalytic function of various active-site residues (Geobacillus stearothermophilus numbering). Specifically, contrary to proposals based on the LRC, E43 is determined to solely be involved in the initial adenine protonation step and not the deglycosylation reaction as the general base. Additionally, a novel catalytic role is proposed for Y126, whereby this residue plays a significant role in stabilizing the highly charged active site, primarily through interactions with E43. More importantly, D144 is found to explicitly catalyze the nucleobase dissociation step through partial nucleophilic attack. Although this is a more direct role than previously proposed for any other DNA glycosylase, comparison to previous work on other glycosylases justifies the larger contribution in the case of MutY and allows us to propose a unified role for the conserved Asp/Glu in the DNA glycosylases, as well as other enzymes that catalyze nucleotide deglycosylation reactions.