Influence of local sequence context on damaged base conformation in human DNA polymerase ι: molecular dynamics studies of nucleotide incorporation opposite a benzo[a]pyrene-derived adenine lesion

Human DNA polymerase ι is a lesion bypass polymerase of the Y family, capable of incorporating nucleotides opposite a variety of lesions in both near error-free and error-prone bypass. With undamaged templating purines polymerase ι normally favors Hoogsteen base pairing. Polymerase ι can incorporate nucleotides opposite a benzo[a]pyrene-derived adenine lesion (dA*); while mainly error-free, the identity of misincorporated bases is influenced by local sequence context. We performed molecular modeling and molecular dynamics simulations to elucidate the structural basis for lesion bypass. Our results suggest that hydrogen bonds between the benzo[a]pyrenyl moiety and nearby bases limit the movement of the templating base to maintain the anti glycosidic bond conformation in the binary complex in a 5′-CAGA*TT-3′ sequence. This facilitates correct incorporation of dT via a Watson−Crick pair. In a 5′-TTTA*GA-3′ sequence the lesion does not form these hydrogen bonds, permitting dA* to rotate around the glycosidic bond to syn and incorporate dT via a Hoogsteen pair. With syn dA*, there is also an opportunity for increased misincorporation of dGTP. These results expand our understanding of the versatility and flexibility of polymerase ι and its lesion bypass functions in humans.     

NMR structure of duplex DNA containing the α‐OH‐PdG·dA base pair: A mutagenic intermediate of acrolein

Acrolein, a cell metabolic product and main component of cigarette smoke, reacts with DNA generating α‐OH‐PdG lesions, which have the ability to pair with dATP during replication thereby causing G to T transversions. We describe the solution structure of an 11‐mer DNA duplex containing the mutagenic α‐OH‐PdG·dA base pair intermediate, as determined by solution nuclear magnetic resonance (NMR) spectroscopy and retrained molecular dynamics (MD) simulations. The NMR data support a mostly regular right‐handed helix that is only perturbed at its center by the presence of the lesion. Undamaged residues of the duplex are in anti orientation, forming standard Watson‐Crick base pairs alignments. Duplication of proton signals at and near the damaged base pair reveals the presence of two enantiomeric duplexes, thus establishing the exocyclic nature of the lesion. The α‐OH‐PdG adduct assumes a syn conformation pairing to its partner dA base that is protonated at pH 6.6. The three‐dimensional structure obtained by restrained molecular dynamics simulations show hydrogen bond interactions that stabilize α‐OH‐PdG in a syn conformation and across the lesion containing base pair. We discuss the implications of the structures for the mutagenic bypass of acrolein lesions. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 391–401, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

DNA oligonucleotides with A, T, G or C opposite an abasic site: structure and dynamics

Abasic sites are common DNA lesions resulting from spontaneous depurination and excision of damaged nucleobases by DNA repair enzymes. However, the influence of the local sequence context on the structure of the abasic site and ultimately, its recognition and repair, remains elusive. In the present study, duplex DNAs with three different bases (G, C or T) opposite an abasic site have been synthesized in the same sequence context (5′-CCA AAG₆ XA₈C CGG G-3′, where X denotes the abasic site) and characterized by 2D NMR spectroscopy. Studies on a duplex DNA with an A opposite the abasic site in the same sequence has recently been reported [Chen,J., Dupradeau,F.-Y., Case,D.A., Turner,C.J. and Stubbe,J. (2007) Nuclear magnetic resonance structural studies and molecular modeling of duplex DNA containing normal and 4′-oxidized abasic sites. Biochemistry, 46, 3096–3107]. Molecular modeling based on NMR-derived distance and dihedral angle restraints and molecular dynamics calculations have been applied to determine structural models and conformational flexibility of each duplex. The results indicate that all four duplexes adopt an overall B-form conformation with each unpaired base stacked between adjacent bases intrahelically. The conformation around the abasic site is more perturbed when the base opposite to the lesion is a pyrimidine (C or T) than a purine (G or A). In both the former cases, the neighboring base pairs (G6-C21 and A8-T19) are closer to each other than those in B-form DNA. Molecular dynamics simulations reveal that transient H-bond interactions between the unpaired pyrimidine (C20 or T20) and the base 3′ to the abasic site play an important role in perturbing the local conformation. These results provide structural insight into the dynamics of abasic sites that are intrinsically modulated by the bases opposite the abasic site.     

Highly Mutagenic Exocyclic DNA Adducts Are Substrates for the Human Nucleotide Incision Repair Pathway

Background

Oxygen free radicals induce lipid peroxidation (LPO) that damages and breaks polyunsaturated fatty acids in cell membranes. LPO-derived aldehydes and hydroxyalkenals react with DNA leading to the formation of etheno(ε)-bases including 1,N ⁶-ethenoadenine (εA) and 3,N ⁴-ethenocytosine (εC). The εA and εC residues are highly mutagenic in mammalian cells and eliminated in the base excision repair (BER) pathway and/or by AlkB family proteins in the direct damage reversal process. BER initiated by DNA glycosylases is thought to be the major pathway for the removal of non-bulky endogenous base damage. Alternatively, in the nucleotide incision repair (NIR) pathway, the apurinic/apyrimidinic (AP) endonucleases can directly incise DNA duplex 5′ to a damaged base in a DNA glycosylase-independent manner.

Methodology/Principal Findings

Here we have characterized the substrate specificity of human major AP endonuclease 1, APE1, towards εA, εC, thymine glycol (Tg) and 7,8-dihydro-8-oxoguanine (8oxoG) residues when present in duplex DNA. APE1 cleaves oligonucleotide duplexes containing εA, εC and Tg, but not those containing 8oxoG. Activity depends strongly on sequence context. The apparent kinetic parameters of the reactions suggest that APE1 has a high affinity for DNA containing ε-bases but cleaves DNA duplexes at an extremely slow rate. Consistent with this observation, oligonucleotide duplexes containing an ε-base strongly inhibit AP site nicking activity of APE1 with IC₅₀ values in the range of 5–10 nM. MALDI-TOF MS analysis of the reaction products demonstrated that APE1-catalyzed cleavage of εA•T and εC•G duplexes generates, as expected, DNA fragments containing 5′-terminal ε-base residue.

Conclusions/Significance

The fact that ε-bases and Tg in duplex DNA are recognized and cleaved by APE1 in vitro, suggests that NIR may act as a backup pathway to BER to remove a large variety of genotoxic base lesions in human cells.     

NMR structure of a parallel-stranded DNA duplex at atomic resolution

DNA dodecamers have been designed with two cytosines on each end and intervening A and T stretches, such that the oligomers have fully complementary A:T base pairs when aligned in the parallel orientation. Spectroscopic (UV, CD and IR), NMR and molecular dynamics studies have shown that oligomers having the sequences d(CCATAATTTACC) and d(CCTATTAAATCC) form a parallel-stranded duplex when dissolved at 1:1 stoichiometry in aqueous solution. This is due to the C:C⁺ clamps on either end and extensive mismatches in the antiparallel orientation. The structure is stable at neutral and acidic pH. At higher temperatures, the duplex melts into single strands in a highly cooperative fashion. All adenine, cytosine and thymine nucleotides adopt the anti conformation with respect to the glycosidic bond. The A:T base pairs form reverse Watson–Crick base pairs. The duplex shows base stacking and NOEs between the base protons T(H6)/A(H8) and the sugar protons (H1′/H2′/H2″) of the preceding nucleotide, as has been observed in antiparallel duplexes. However, no NOEs are observed between base protons H2/H6/H8 of sequential nucleotides, though such NOEs are observed between T(CH₃) and A(H8). A three-dimensional structure of the parallel-stranded duplex at atomic resolution has been obtained using molecular dynamics simulations under NMR constraints. The simulated structures have torsional angles very similar to those found in B-DNA duplexes, but the base stacking and helicoid parameters are significantly different.     

Temperature and Electrolyte Optimization of the α-Hemolysin Latch Sensing Zone for Detection of Base Modification in Double-Stranded DNA

The latch region of the wild-type protein pore α-hemolysin (α-HL) constitutes a sensing zone for individual abasic sites (and furan analogs) in double-stranded DNA (dsDNA). The presence of an abasic site or furan within a DNA duplex, electrophoretically captured in the α-HL vestibule and positioned at the latch region, can be detected based on the current blockage prior to duplex unzipping. We investigated variations in blockage current as a function of temperature (12–35°C) and KCl concentration (0.15–1.0 M) to understand the origin of the current signature and to optimize conditions for identifying the base modification. In 1 M KCl solution, substitution of a furan for a cytosine base in the latch region results in an ∼8 kJ mol⁻¹ decrease in the activation energy for ion transport through the protein pore. This corresponds to a readily measured ∼2 pA increase in current at room temperature. Optimal resolution for detecting the presence of a furan in the latch region is achieved at lower KCl concentrations, where the noise in the measured blockage current is significantly lower. The noise associated with the blockage current also depends on the stability of the duplex (as measured from the melting temperature), where a greater noise in the measured blockage current is observed for less stable duplexes.     

Unzipping of A-Form DNA-RNA,A-Form DNA-PNA,and B-Form DNA-DNA in the α-Hemolysin Nanopore

Unzipping of double-stranded nucleic acids by an electric field applied across a wild-type α-hemolysin (αHL) nanopore provides structural information about different duplex forms. In this work, comparative studies on A-form DNA-RNA duplexes and B-form DNA-DNA duplexes with a single-stranded tail identified significant differences in the blockage current and the unzipping duration between the two helical forms. We observed that the B-form duplex blocks the channel 1.9 ± 0.2 pA more and unzips ∼15-fold more slowly than an A-form duplex at 120 mV. We developed a model to describe the dependence of duplex unzipping on structure. We demonstrate that the wider A-form duplex (d = 2.4 nm) is unable to enter the vestibule opening of αHL on the cis side, leading to unzipping outside of the nanopore with higher residual current and faster unzipping times. In contrast, the smaller B-form duplexes (d = 2.0 nm) enter the vestibule of αHL, resulting in decreased current blockages and slower unzipping. We investigated the effects of varying the length of the single-stranded overhang, and studied A-form DNA-PNA duplexes to provide additional support for the proposed model. This study identifies key differences between A- and B-form duplex unzipping that will be important in the design of future probe-based methods for detecting DNA or RNA.     

Genetic Control of Replication through N1-methyladenine in Human Cells

N1-methyl adenine (1-MeA) is formed in DNA by reaction with alkylating agents and naturally occurring methyl halides. The 1-MeA lesion impairs Watson-Crick base pairing and blocks normal DNA replication. Here we identify the translesion synthesis (TLS) DNA polymerases (Pols) required for replicating through 1-MeA in human cells and show that TLS through this lesion is mediated via three different pathways in which Pols ι and θ function in one pathway and Pols η and ζ, respectively, function in the other two pathways. Our biochemical studies indicate that in the Polι/Polθ pathway, Polι would carry out nucleotide insertion opposite 1-MeA from which Polθ would extend synthesis. In the Polη pathway, this Pol alone would function at both the nucleotide insertion and extension steps of TLS, and in the third pathway, Polζ would extend from the nucleotide inserted opposite 1-MeA by an as yet unidentified Pol. Whereas by pushing 1-MeA into the syn conformation and by forming Hoogsteen base pair with the T residue, Polι would carry out TLS opposite 1-MeA, the ability of Polη to replicate through 1-MeA suggests that despite its need for Watson-Crick hydrogen bonding, Polη can stabilize the adduct in its active site. Remarkably, even though Pols η and ι are quite error-prone at inserting nucleotides opposite 1-MeA, TLS opposite this lesion in human cells occurs in a highly error-free fashion. This suggests that the in vivo fidelity of TLS Pols is regulated by factors such as post-translational modifications, protein-protein interactions, and possibly others.     

Structural effect of the anticancer agent 6-thioguanine on duplex DNA

The incorporation of 6-thioguanine (S6G) into DNA is an essential step in the cytotoxic activity of thiopurines. However, the structural effects of this substitution on duplex DNA have not been fully characterized. Here, we present the solution structures of DNA duplexes containing S6G opposite thymine (S6G·T) and opposite cytosine (S6G·C), solved by high-resolution NMR spectroscopy and restrained molecular dynamics. The data indicate that both duplexes adopt right-handed helical conformations with all Watson–Crick hydrogen bonding in place. The S6G·T structures exhibit a wobble-type base pairing at the lesion site, with thymine shifted toward the major groove and S6G displaced toward the minor groove. Aside from the lesion site, the helices, including the flanking base pairs, are not highly perturbed by the presence of the lesion. Surprisingly, thermal dependence experiments suggest greater stability in the S6G-T mismatch than the S6G-C base pair.     

Solution structure of the DNA decamer duplex containing a 3'-T x T basepair of the cis-syn cyclobutane pyrimidine dimer: implication for the mutagenic property of the cis-syn dimer

The cis–syn dimer is the major DNA photoproduct produced by UV irradiation. In order to determine the origin of the mutagenic property of the cis–syn dimer, we used NMR restraints and molecular dynamics to determine the solution structure of a DNA decamer duplex containing a wobble pair between the 3′-T of the cis–syn dimer and the opposite T residue (CS/TA duplex). The solution structure of the CS/TA duplex revealed that the 3′-T·T base pair of the cis–syn dimer had base pair geometry that was significantly different from the canonical Watson–Crick base pair and caused destabilization and conformational distortion of its 3′-region. However, a 3′-T·A base pair at the cis–syn dimer within this related DNA decamer maintains the normal Watson–Crick base pair geometry and causes little distortion in the conformation of its 3′-side. Our results show that in spite of its stable hydrogen bonding, the insertion of a T residue opposite the 3′-T of the cis–syn dimer is inhibited by structural distortion caused by the 3′-T·T base pair. This may explain why the frequency of the 3′-T→A transversion, which is the major mutation produced by the cis–syn dimer, is only 4%.     

Impact of the C1' configuration of abasic sites on DNA duplex structure

Naturally occurring abasic sites in DNA exist as an equilibrium mixture of the aldehyde, the hydrated aldehyde, and the hemiacetal forms (dominant). The influence of the configuration of the C1' hydroxyl group of the hemiacetal form on duplex structure and abasic site repair has been examined using novel carbocyclic analogues. Both the alpha- and beta-forms of this novel abasic site were introduced into oligomeric DNA using the standard DMT-phosphoramidite approach in an automated solid-phase synthesizer. Solution structures of the d(CGTACXCATGC).d(GCATGAGTACG) duplex (where X is the alpha- or beta-anomer of the carbocyclic abasic site analogue) were determined by NMR spectroscopy and restrained molecular dynamics simulations. The structures were only minimally perturbed by the presence of either anomer of the abasic site. All residues adopted an anti conformation, and Watson-Crick alignments were observed on all base pairs of the duplexes. At the lesion site, the abasic residues and their partner adenines showed increased dynamic behavior but adopted intrahelical positions in the final refined structures. Incision of duplexes having the alpha- or beta-anomer of the carbocyclic abasic site by human AP endonuclease showed that the enzyme recognizes both configurations of the lesion and nicks the DNA backbone with similar efficiency. Our results challenge the suggestion that Ape1 is stereoselective and imply a plasticity at the active site of the enzyme for accommodating either anomer of the lesion.     

Dynamic Regulation of α-Actinin’s Calponin Homology Domains on F-Actin

α-Actinin is an essential actin cross-linker involved in cytoskeletal organization and dynamics. The molecular conformation of α-actinin’s actin-binding domain (ABD) regulates its association with actin and thus mutations in this domain can lead to severe pathogenic conditions. A point mutation at lysine 255 in human α-actinin-4 to glutamate increases the binding affinity resulting in stiffer cytoskeletal structures. The role of different ABD conformations and the effect of K255E mutation on ABD conformations remain elusive. To evaluate the impact of K255E mutation on ABD binding to actin we use all-atom molecular dynamics and free energy calculation methods and study the molecular mechanism of actin association in both wild-type α-actinin and in the K225E mutant. Our models illustrate that the strength of actin association is indeed sensitive to the ABD conformation, predict the effect of K255E mutation—based on simulations with the K237E mutant chicken α-actinin—and evaluate the mechanism of α-actinin binding to actin. Furthermore, our simulations showed that the calmodulin domain binding to the linker region was important for regulating the distance between actin and ABD. Our results provide valuable insights into the molecular details of this critical cellular phenomenon and further contribute to an understanding of cytoskeletal dynamics in health and disease.     

Structure of Human DNA Polymerase κ Inserting dATP Opposite an 8-OxoG DNA Lesion

Background

Oxygen-free radicals formed during normal aerobic cellular metabolism attack bases in DNA and 7,8-dihydro-8-oxoguanine (8-oxoG) is one of the major lesions formed. It is amongst the most mutagenic lesions in cells because of its dual coding potential, wherein 8-oxoG(syn) can pair with an A in addition to normal base pairing of 8-oxoG(anti) with a C. Human DNA polymerase κ (Polκ) is a member of the newly discovered Y-family of DNA polymerases that possess the ability to replicate through DNA lesions. To understand the basis of Polκ''s preference for insertion of an A opposite 8-oxoG lesion, we have solved the structure of Polκ in ternary complex with a template-primer presenting 8-oxoG in the active site and with dATP as the incoming nucleotide.

Methodology and Principal Findings

We show that the Polκ active site is well-adapted to accommodate 8-oxoG in the syn conformation. That is, the polymerase and the bound template-primer are almost identical in their conformations to that in the ternary complex with undamaged DNA. There is no steric hindrance to accommodating 8-oxoG in the syn conformation for Hoogsteen base-paring with incoming dATP.

Conclusions and Significance

The structure we present here is the first for a eukaryotic translesion synthesis (TLS) DNA polymerase with an 8-oxoG:A base pair in the active site. The structure shows why Polκ is more efficient at inserting an A opposite the 8-oxoG lesion than a C. The structure also provides a basis for why Polκ is more efficient at inserting an A opposite the lesion than other Y-family DNA polymerases.     

Error-prone lesion bypass by human DNA polymerase eta

DNA lesion bypass is an important cellular response to genomic damage during replication. Human DNA polymerase η (Polη), encoded by the Xeroderma pigmentosum variant (XPV) gene, is known for its activity of error-free translesion synthesis opposite a TT cis-syn cyclobutane dimer. Using purified human Polη, we have examined bypass activities of this polymerase opposite several other DNA lesions. Human Polη efficiently bypassed a template 8-oxoguanine, incorporating an A or a C opposite the lesion with similar efficiencies. Human Polη effectively bypassed a template abasic site, incorporating an A and less frequently a G opposite the lesion. Significant –1 deletion was also observed when the template base 5′ to the abasic site is a T. Human Polη partially bypassed a template (+)-trans-anti-benzo[a]pyrene-N²-dG and predominantly incorporated an A, less frequently a T, and least frequently a G or a C opposite the lesion. This specificity of nucleotide incorporation correlates well with the known mutation spectrum of (+)-trans-anti-benzo[a]pyrene-N²-dG lesion in mammalian cells. These results show that human Polη is capable of error-prone translesion DNA syntheses in vitro and suggest that Polη may bypass certain lesions with a mutagenic consequence in humans.     

Structural basis for uracil DNA glycosylase interaction with uracil: NMR study

Two dimensional (2D) NMR and molecular dynamics simulations have been used to determine the three dimensional (3D) structure of a hairpin DNA, d-CTA-GAGGATCC-TUTT-GGATCCT (22mer; abbreviated as U2-hairpin), which has uracil at the second position from the 5′ end of the tetraloop. The ¹H resonances of this hairpin have been assigned almost completely. NMR restrained molecular dynamics and energy minimization procedures have been used to describe the 3D structure of U2-hairpin. This study establishes that the stem of the hairpin adopts a right-handed B-DNA conformation, while the T₁₂ and T₁₅ nucleotides stack upon 3′ and 5′ ends of the stem, respectively. Further, T₁₄ stacks upon both T₁₂ and T₁₅. Though U₁₃ partially stacks upon T₁₄, no stacking interaction is observed between U₁₃ and T₁₂. All the individual nucleotide bases belonging to the stem and T₁₂ and T₁₅ of the loop adopt ‘anti’ conformation with respect to their sugar moiety, while the U₁₃ and T₁₄ of the loop are in ‘syn’ conformation. The turning phosphate in the loop is located between T₁₃ and T₁₄. This study and a concurrent NMR structural study on yet another hairpin DNA d-CTAGAGGAATAA-TTTU-GGATCCT (22mer; abbreviated as U4-hairpin), with uracil at the fourth position from the 5′ end of the tetraloop throw light upon various interactions which have been reported between Escherichia coli uracil DNA glycosylase (UDG) and uracil containing DNA. The of T₁₂ and α, β, γ, and ζ of U₁₃ and γ of T₁₄, which partially influence the local conformation of U₁₃ in U2-hairpin are all locked in ‘trans’ conformation. Such stretched out backbone conformation in the vicinity of U₁₃ could be the reason as to why the U2-hairpin is found to be the poor substrate for its interaction with UDG compared to the other substrates in which the uracil is at first, third and fourth positions of the tetraloop from its 5′ end, as reported earlier by Vinay and Varshney. This study shows that UDG actively promotes the flipping of uracil from a stacked conformation and rules out the possibility of UDG recognizing the flipped out uracil bases.     

The solution structure of [d(CGC)r(aaa)d(TTTGCG)](2): hybrid junctions flanked by DNA duplexes

The solution structure and hydration of the chimeric duplex [d(CGC)r(aaa)d(TTTGCG)]₂, in which the central hybrid segment is flanked by DNA duplexes at both ends, was determined using two-dimensional NMR, simulated annealing and restrained molecular dynamics. The solution structure of this chimeric duplex differs from the previously determined X-ray structure of the analogous B-DNA duplex [d(CGCAAATTTGCG)]₂ as well as NMR structure of the analogous A-RNA duplex [r(cgcaaauuugcg)]₂. Long-lived water molecules with correlation time τ_c longer than 0.3 ns were found close to the RNA adenine H2 and H1′ protons in the hybrid segment. A possible long-lived water molecule was also detected close to the methyl group of 7T in the RNA–DNA junction but not with the other two thymines (8T and 9T). This result correlates with the structural studies that only DNA residue 7T in the RNA–DNA junction adopts an O4′-endo sugar conformation, while the other DNA residues including 3C in the DNA–RNA junction, adopt C1′-exo or C2′-endo conformations. The exchange rates for RNA C2′-OH were found to be ~5–20 s^–1. This slow exchange rate may be due to the narrow minor groove width of [d(CGC)r(aaa)d(TTTGCG)]₂, which may trap the water molecules and restrict the dynamic motion of hydroxyl protons. The minor groove width of [d(CGC)r(aaa)d(TTTGCG)]₂ is wider than its B-DNA analog but narrower than that of the A-RNA analog. It was further confirmed by its titration with the minor groove binding drug distamycin. A possible 2:1 binding mode was found by the titration experiments, suggesting that this chimeric duplex contains a wider minor groove than its B-DNA analog but still narrow enough to hold two distamycin molecules. These distinct structural features and hydration patterns of this chimeric duplex provide a molecular basis for further understanding the structure and recognition of DNA·RNA hybrid and chimeric duplexes.     

Structural Basis for Promutagenicity of 8-Halogenated Guanine

8-Halogenated guanine (haloG), a major DNA adduct formed by reactive halogen species during inflammation, is a promutagenic lesion that promotes misincorporation of G opposite the lesion by various DNA polymerases. Currently, the structural basis for such misincorporation is unknown. To gain insights into the mechanism of misincorporation across haloG by polymerase, we determined seven x-ray structures of human DNA polymerase β (polβ) bound to DNA bearing 8-bromoguanine (BrG). We determined two pre-catalytic ternary complex structures of polβ with an incoming nonhydrolyzable dGTP or dCTP analog paired with templating BrG. We also determined five binary complex structures of polβ in complex with DNA containing BrG·C/T at post-insertion and post-extension sites. In the BrG·dGTP ternary structure, BrG adopts syn conformation and forms Hoogsteen base pairing with the incoming dGTP analog. In the BrG·dCTP ternary structure, BrG adopts anti conformation and forms Watson-Crick base pairing with the incoming dCTP analog. In addition, our polβ binary post-extension structures show Hoogsteen BrG·G base pair and Watson-Crick BrG·C base pair. Taken together, the first structures of haloG-containing DNA bound to a protein indicate that both BrG·G and BrG·C base pairs are accommodated in the active site of polβ. Our structures suggest that Hoogsteen-type base pairing between G and C8-modified G could be accommodated in the active site of a DNA polymerase, promoting G to C mutation.     

Structural Determinants of the Mechanical Stability of α-Catenin

α-Catenin plays a crucial role in cadherin-mediated adhesion by binding to β-catenin, F-actin, and vinculin, and its dysfunction is linked to a variety of cancers and developmental disorders. As a mechanotransducer in the cadherin complex at intercellular adhesions, mechanical and force-sensing properties of α-catenin are critical to its proper function. Biochemical data suggest that α-catenin adopts an autoinhibitory conformation, in the absence of junctional tension, and biophysical studies have shown that α-catenin is activated in a tension-dependent manner that in turn results in the recruitment of vinculin to strengthen the cadherin complex/F-actin linkage. However, the molecular switch mechanism from autoinhibited to the activated state remains unknown for α-catenin. Here, based on the results of an aggregate of 3 μs of molecular dynamics simulations, we have identified a dynamic salt-bridge network within the core M region of α-catenin that may be the structural determinant of the stability of the autoinhibitory conformation. According to our constant-force steered molecular dynamics simulations, the reorientation of the MII/MIII subdomains under force may constitute an initial step along the transition pathway. The simulations also suggest that the vinculin-binding domain (subdomain MI) is intrinsically much less stable than the other two subdomains in the M region (MII and MIII). Our findings reveal several key insights toward a complete understanding of the multistaged, force-induced conformational transition of α-catenin to the activated conformation.     

Role of helical constraints of the EBS1–IBS1 duplex of a group II intron on demarcation of the 5′ splice site

Recognition of the 5′ splice site by group II introns involves pairing between an exon binding sequence (EBS) 1 within the ID3 stem–loop of domain 1 and a complementary sequence at the 3′ end of exon 1 (IBS1). To identify the molecular basis for splice site definition of a group IIB ai5γ intron, we probed the solution structure of the ID3 stem–loop alone and upon binding of its IBS1 target by solution NMR. The ID3 stem was structured. The base of the ID3 loop was stacked but displayed a highly flexible EBS1 region. The flexibility of EBS1 appears to be a general feature of the ai5γ and the smaller Oceanobacillus iheyensis (O.i.) intron and may help in effective search of conformational space and prevent errors in splicing as a result of fortuitous base-pairing. Binding of IBS1 results in formation of a structured seven base pair duplex that terminates at the 5′ splice site in spite of the potential for additional A-U and G•U pairs. Comparison of these data with conformational features of EBS1–IBS1 duplexes extracted from published structures suggests that termination of the duplex and definition of the splice site are governed by constraints of the helical geometry within the ID3 loop. This feature and flexibility of the uncomplexed ID3 loop appear to be common for both the ai5γ and O.i. introns and may help to fine-tune elements of recognition in group II introns.     

Insights into structure, dynamics and hydration of locked nucleic acid (LNA) strand-based duplexes from molecular dynamics simulations

Locked nucleic acid (LNA) is a chemically modified nucleic acid with its sugar ring locked in an RNA-like (C3′-endo) conformation. LNAs show extraordinary thermal stabilities when hybridized with DNA, RNA or LNA itself. We performed molecular dynamics simulations on five isosequential duplexes (LNA–DNA, LNA–LNA, LNA–RNA, RNA–DNA and RNA–RNA) in order to characterize their structure, dynamics and hydration. Structurally, the LNA–DNA and LNA–RNA duplexes are found to be similar to regular RNA–DNA and RNA–RNA duplexes, whereas the LNA–LNA duplex is found to have its helix partly unwound and does not resemble RNA–RNA duplex in a number of properties. Duplexes with an LNA strand have on average longer interstrand phosphate distances compared to RNA–DNA and RNA–RNA duplexes. Furthermore, intrastrand phosphate distances in LNA strands are found to be shorter than in DNA and slightly shorter than in RNA. In case of induced sugar puckering, LNA is found to tune the sugar puckers in partner DNA strand toward C3′-endo conformations more efficiently than RNA. The LNA–LNA duplex has lesser backbone flexibility compared to the RNA–RNA duplex. Finally, LNA is less hydrated compared to DNA or RNA but is found to have a well-organized water structure.