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.

The Impact of Antimicrobial Peptides on the Acinetobacter baumannii Inner Membrane Is Modulated by Lipid Polyunsaturation.

The Gram-negative pathogen Acinetobacter baumannii is a primary contributor to nosocomial multi-drug-resistant (MDR) infections. To combat the rise of MDR infections, novel features of A. baumannii need to be considered for the development of new treatment options. One such feature is the preferential scavenging of exogenous lipids, including host-derived polyunsaturated fatty acids (PUFAs), for membrane phospholipid synthesis. These alterations in membrane composition impact both the lipid chemistry and the membrane biophysical properties. In this work we examine how antimicrobial peptides (AMPs) interact with the inner membranes of A. baumannii in the presence and absence of polyunsaturated phospholipids. Using coarse-grained molecular dynamics simulations of complex A. baumannii inner membrane models derived from lipidomes of bacteria grown in the presence and absence of PUFAs, we examine the impact of the adsorption of four prototypical AMPs (CAMEL, LL-37, pexiganan, and magainin-2) on the membrane biophysical properties. Our simulations reveal that the impact of AMP adsorption on the membrane biophysical properties was dependent on both the membrane composition and the specific AMP involved. Both lipid headgroup charge and tail unsaturation played important roles in driving the interactions that occurred both within the membrane and between the membrane and AMPs. The changes to the membrane biophysical properties also showed a complex relationship with the AMP's physical properties, such as AMP charge, chain length, and charge-to-mass ratio. Cumulatively, this work highlights the importance of studying AMPs using a complex membrane environment and provides insights into the mechanistic action of AMPs in polyunsaturated lipid-rich bacterial membranes.

Membrane cholesterol regulates inhibition and substrate transport by the glycine transporter, GlyT2.

Membrane cholesterol binds to and modulates the function of various SLC6 neurotransmitter transporters, including stabilizing the outward-facing conformation of the dopamine and serotonin transporters. Here, we investigate how cholesterol binds to GlyT2 (SLC6A5), modulates glycine transport rate, and influences bioactive lipid inhibition of GlyT2. Bioactive lipid inhibitors are analgesics that bind to an allosteric site accessible from the extracellular solution when GlyT2 adopts an outward-facing conformation. Using molecular dynamics simulations, mutagenesis, and cholesterol depletion experiments, we show that bioactive lipid inhibition of glycine transport is modulated by the recruitment of membrane cholesterol to a binding site formed by transmembrane helices 1, 5, and 7. Recruitment involves cholesterol flipping from its membrane orientation, and insertion of the 3' hydroxyl group into the cholesterol binding cavity, close to the allosteric site. The synergy between cholesterol and allosteric inhibitors provides a novel mechanism of inhibition and a potential avenue for the development of potent GlyT2 inhibitors as alternative therapeutics for the treatment of neuropathic pain and therapeutics that target other SLC6 transporters.

Conserved FimH mutations in the global Escherichia coli ST131 multi-drug resistant lineage weaken interdomain interactions and alter adhesin function.

The binding of the type 1 fimbrial adhesin FimH to mannosylated receptors is allosterically regulated to enhance the fitness of uropathogenic Escherichia coli (UPEC) during urinary tract infection (UTI). Mutations in the two FimH domains (pilin and lectin) located outside the mannose binding pocket have been shown to influence mannose binding affinity, yet the details of the allostery mechanism are not fully elucidated. Here we characterised different FimH conformational states (termed low-affinity tense and high-affinity relaxed conformations) of natural FimH variants using molecular dynamics (MD) simulation techniques and report key structural dynamics differences between them. The clinically dominant FimH30 variant from the pandemic multidrug resistant E. coli ST131 lineage contains an R166H mutation that weakens FimH interdomain interactions and allows enhanced mannose interactions with pre-existing high-affinity relaxed conformations. When expressed in an isogenic ST131 strain background, FimH30 mediated high human cell adhesion and invasion, and enhanced biofilm formation over other variants. Collectively, our computational and experimental findings support a model of FimH protein allostery that is mediated by shifts in the pre-existing conformational equilibrium of FimH, additional to the sequential step-wise process of structural perturbations transmitted from one site to another within the protein. Importantly, it is the first study to shed light into how natural mutations in a clinically dominant FimH variant influence the protein's conformational landscape optimising its function for ST131 fitness at intestinal and extraintestinal niches.

Loss, Gain and Altered Function of GlyR α2 Subunit Mutations in Neurodevelopmental Disorders.

Glycine receptors (GlyRs) containing the α2 subunit govern cell fate, neuronal migration and synaptogenesis in the developing cortex and spinal cord. Rare missense variants and microdeletions in the X-linked GlyR α2 subunit gene (GLRA2) have been associated with human autism spectrum disorder (ASD), where they typically cause a loss-of-function via protein truncation, reduced cell-surface trafficking and/or reduced glycine sensitivity (e.g., GLRA2Δex8-9 and extracellular domain variants p.N109S and p.R126Q). However, the GlyR α2 missense variant p.R323L in the intracellular M3-M4 domain results in a gain-of-function characterized by slower synaptic decay times, longer duration active periods and increases in channel conductance. This study reports the functional characterization of four missense variants in GLRA2 associated with ASD or developmental disorders (p.V-22L, p.N38K, p.K213E, p.T269M) using a combination of bioinformatics, molecular dynamics simulations, cellular models of GlyR trafficking and electrophysiology in artificial synapses. The GlyR α2V-22L variant resulted in altered predicted signal peptide cleavage and a reduction in cell-surface expression, suggestive of a partial loss-of-function. Similarly, GlyR α2N38K homomers showed reduced cell-surface expression, a reduced affinity for glycine and a reduced magnitude of IPSCs in artificial synapses. By contrast, GlyR α2K213E homomers showed a slight reduction in cell-surface expression, but IPSCs were larger, with faster rise/decay times, suggesting a gain-of-function. Lastly, GlyR α2T269M homomers exhibited a high glycine sensitivity accompanied by a substantial leak current, suggestive of an altered function that could dramatically enhance glycinergic signaling. These results may explain the heterogeneity of clinical phenotypes associated with GLRA2 mutations and reveal that missense variants can result in a loss, gain or alteration of GlyR α2 function. In turn, these GlyR α2 missense variants are likely to either negatively or positively deregulate cortical progenitor homeostasis and neuronal migration in the developing brain, leading to changes in cognition, learning, and memory.

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.

Site of Cholesterol Oxidation Impacts Its Localization and Domain Formation in the Neuronal Plasma Membrane.

Cholesterol is integral to the structure of mammalian cell membranes. Oxidation of cholesterol alters how it behaves in the membrane and influences the membrane biophysical properties. Elevated levels of oxidized cholesterol are associated with neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, multiple sclerosis, and Huntington's disease. Previous work has investigated the impact of oxidized cholesterol in the context of simple model membrane systems. However, there is a growing body of literature that shows that complex membranes possessing physiological phospholipid distributions have different properties from those of binary or trinary model membranes. In the current work, the impact of oxidized cholesterol on the biophysical properties of a complex neuronal plasma membrane is investigated using coarse-grained Martini molecular dynamics simulations. Comparison of the native neuronal membrane to neuronal membranes containing 10% tail-oxidized or 10% head-oxidized cholesterol shows that the site of oxidization changes the behavior of the oxidized cholesterol in the membrane. Furthermore, species-specific domain formation is observed between each oxidized cholesterol and minor lipid classes. Although both tail-oxidized and head-oxidized cholesterols modulate the biophysical properties of the membrane, smaller changes are observed in the complex neuronal membrane than seen in the previous work on simple binary or trinary model membranes. This work highlights the presence of compensatory effects of lipid diversity in the complex neuronal membrane. Overall, this study improves our molecular-level understanding of the effects of oxidized cholesterol on the properties of neuronal tissue and emphasizes the importance of studying membranes with realistic lipid compositions.

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.

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.

Investigating the lipid fingerprint of SLC6 neurotransmitter transporters: a comparison of dDAT, hDAT, hSERT, and GlyT2.

The local lipid annulus, or "fingerprint", of four SLC6 transporters (dDAT, hDAT, hSERT, and GlyT2) embedded in a complex neuronal membrane were compared and characterised using molecular dynamics. Our analysis included the development of new tools to improve membrane leaflet detection and the analysis of leaflet-dependent properties. Overall, the lipid fingerprints of the four transporters are comprised of similar lipids when grouped by headgroup or tail saturation. The enrichment and depletion of specific lipids, including sites of cholesterol contacts, varies between transporters. The subtle differences in lipid fingerprints results in varying membrane biophysical properties near the transporter. Our results highlight that the lipid-fingerprint of SLC6 transporters in complex membranes is highly dependent on membrane composition. Our results further characterize how the presence and identity of membrane proteins affects the complex interplay of lipid-protein interactions, influencing the local lipid environment and membrane biophysical properties.

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.

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.

Identification of an allosteric binding site on the human glycine transporter, GlyT2, for bioactive lipid analgesics.

The treatment of chronic pain is poorly managed by current analgesics, and there is a need for new classes of drugs. We recently developed a series of bioactive lipids that inhibit the human glycine transporter GlyT2 (SLC6A5) and provide analgesia in animal models of pain. Here, we have used functional analysis of mutant transporters combined with molecular dynamics simulations of lipid-transporter interactions to understand how these bioactive lipids interact with GlyT2. This study identifies a novel extracellular allosteric modulator site formed by a crevice between transmembrane domains 5, 7, and 8, and extracellular loop 4 of GlyT2. Knowledge of this site could be exploited further in the development of drugs to treat pain, and to identify other allosteric modulators of the SLC6 family of transporters.

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.

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.

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.