Designing molecular switches based on DNA-base mispairing

A unified computational view of DNA duplex, triplex, quadruplex and their donor-acceptor interactions

DNA can assume various structures as a result of interactions at atomic and molecular levels (e.g., hydrogen bonds, π–π stacking interactions, and electrostatic potentials), so understanding of the consequences of these interactions could guide development of ways to produce elaborate programmable DNA for applications in bio- and nanotechnology. We conducted advanced ab initio calculations to investigate nucleobase model structures by componentizing their donor-acceptor interactions. By unifying computational conditions, we compared the independent interactions of DNA duplexes, triplexes, and quadruplexes, which led us to evaluate a stability trend among Watson–Crick and Hoogsteen base pairing, stacking, and even ion binding. For a realistic solution-like environment, the influence of water molecules was carefully considered, and the potassium-ion preference of G-quadruplex was first analyzed at an ab initio level by considering both base-base and ion-water interactions. We devised new structure factors including hydrogen bond length, glycosidic vector angle, and twist angle, which were highly effective for comparison between computationally-predicted and experimentally-determined structures; we clarified the function of phosphate backbone during nucleobase ordering. The simulated tendency of net interaction energies agreed well with that of real world, and this agreement validates the potential of ab initio study to guide programming of complicated DNA constructs.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

Single-stranded DNA structural diversity: TAGGGT from monomers to dimers to tetramer formation

RATIONALE

DNA quadruplex structures have emerged as novel drug targets due to their role in preventing abnormal gene transcription and maintaining telomere stability. Trapped Ion Mobility Spectrometry-Mass Spectrometry (TIMS-MS), combined with theoretical modeling, is a powerful tool for studying the kinetic intermediates of DNA complexes formed in solution and interrogated in the gas phase after desolvation.

METHODS

A TAGGGT ssDNA sequence was purchased and studied in 10 mM ammonium acetate using nanospray electrospray ionization (nESI)-TIMS-MS in positive and negative ion mode. Collisional cross section (CCS) profiles were measured using internal calibration (Tune Mix). Theoretical structures were proposed based on molecular dynamics, charge location and geometry optimization for the most intense IMS bands based on the number of TAGGGT units, adduct form and charge states.

RESULTS

A distribution of monomeric, dimeric and tetrameric TAGGGT structures were formed in solution and separated in the gas phase based on their mobility and m/z value (e.g., [M + 2H]⁺² , [2M + 3H]⁺³ , [M - 2H]^-2 , [2M - 3H]^-3 , [4M + 4H]⁺⁴ , [4M + 3H + NH₄ ]⁺⁴ , [4M + 2H + 2NH₄ ]⁺⁴ and [4M + H + 3NH₄ ]⁺⁴ ). The high mobility resolution of the TIMS-MS analyzer permitted the observation of multiple CCS bands per molecular ion form. Comparison with theoretical candidate structures suggests that monomeric TAGGGT species are stabilized by A-T and G⁺ -G interactions, with the size of the conformer influenced by the proton location. In the case of the TAGGGT quadruplex, the protonated species displayed a broad CCS distribution, while six discrete conformers were stabilized by the presence of ammonium ions (n = 1-3).

CONCLUSIONS

This is the first observation of multiple conformations of TAGGGT complexes (n = 1, 2 and 4) in 10 mM ammonium acetate. Candidate structures with intramolecular interactions of the form of G⁺ -G and traditional A-T base pairing agreed with the experimental trends. Our results demonstrate the structural diversity of TAGGGT monomers, dimers and tetramers in the gas phase beyond the previously reported solution structure, using 10 mM ammonium acetate to replicate biological conditions.

Intramolecular hydrogen bond directed distribution of conformational populations in the derivatives of N′-benzylidenebenzohydrazide

The extensive NMR investigations reveal the presence of E-isomers in the derivative of N′-benzylidenebenzohydrazide. The different conformer populations are controlled by the strength of intramolecular hydrogen bonds.

Direct Selection Strategy for Isolating Aptamers with pH-Sensitive Binding Activity

An aptamer reagent that can switch its binding affinity in a pH-responsive manner would be highly valuable for many biomedical applications including imaging and drug delivery. Unfortunately, the discovery of such aptamers is difficult and only a few have been reported to date. Here we report the first experimental strategy for generating pH-responsive aptamers through direct selection. As an exemplar, we report streptavidin-binding aptamers that retain nanomolar affinity at pH 7.4 but exhibit a ∼100-fold decrease in affinity at pH 5.2. These aptamers were generated by incorporating a known streptavidin-binding DNA motif into an aptamer library and performing FACS-based screening at multiple pH conditions. Upon structural analysis, we found that one aptamer’s affinity-switching behavior is driven by a noncanonical G-A base-pair that controls its folding in a highly pH-dependent manner. We believe our strategy could be readily extended to other aptamer-target systems because it does not require a priori structural knowledge of the aptamer or the target.

On the misincorporation of nucleotides opposite mutagenic cyclic 1,N2-propanoguanine: A computational investigation

The misincorporation properties of exocyclic DNA adduct, cyclic 1,N²-propanoguanine with nucleobases have been investigated using DFT and DFT-D methods. Number of possible and stable mispairing conformations of cyclic 1,N²-propanoguanine with A,T,G and C have been considered for our investigation. The single point energy calculations have been carried out at the M06/6-311++G**, ωB97XD/6-311++G** and MP2/6-311++G** levels on corresponding optimized geometries. The reaction enthalpy values were employed at the M06/6-31 + G* and ωB97XD/6-31 + G* levels. The energies have been compared among the cyclic 1,N²-propanoguanine adduct with nucleobases to find the most stable conformer. The CPCM model was utilized on account of solvent phase and overall polarizability. The computed binding energies follow the order as CPr-Gua-G(2)(-23.2 kcal/mol) > CPr-Gua-C(1) (-16.1 kcal/mol) > CPr-Gua-A(2)(-10.6 kcal/mol) > CPr-Gua-T(2)(-9.6 kcal/mol) in the gas phase at M06 level, which indicates the guanine and cytosine are favorable for mispairing with the cyclic 1,N²-propanoguanine adduct. The obtained results using computational tools are in good agreement with the experimental observation.

pHAST (pH-Driven Aptamer Switch for Thrombin) Catch-and-Release of Target Protein

A pH-driven DNA nanomachine based on the human α-thrombin binding aptamer was designed for the specific catch-and-release of human α-thrombin at neutral and acidic pH, respectively. In neutral conditions, the thrombin aptamer component of the nanomachine is exposed and exists in the G-quadruplex conformation required to bind to the target protein. At slightly acidic pH, the polyadenine tail of the nanomachine becomes partially protonated and A+(anti)•G(syn) mispairing results in a conformational change, causing the target protein to be released. Förster resonance energy transfer (FRET) was used to monitor conformational switching over multiple pH cycles. Electrophoretic mobility shift assay (EMSA) and fluorescence anisotropy were used to show pH dependent protein binding and release by the nanomachine. This approach could be applied generally to existing G-rich aptamers to develop novel biosensors, theranostics, and nanoswitches.

The Role of Aromaticity, Hybridization, Electrostatics, and Covalency in Resonance-Assisted Hydrogen Bonds of Adenine-Thymine (AT) Base Pairs and Their Mimics

Hydrogen bonds play a crucial role in many biochemical processes and in supramolecular chemistry. In this study, we show quantum chemically that neither aromaticity nor other forms of π assistance are responsible for the enhanced stability of the hydrogen bonds in adenine-thymine (AT) DNA base pairs. This follows from extensive bonding analyses of AT and smaller analogs thereof, based on dispersion-corrected density functional theory (DFT). Removing the aromatic rings of either A or T has no effect on the Watson-Crick bond strength. Only when the smaller mimics become saturated, that is, when the hydrogen-bond acceptor and donor groups go from sp (2) to sp (3), does the stability of the resulting model complexes suddenly drop. Bonding analyses based on quantitative Kohn-Sham molecular orbital theory and corresponding energy decomposition analyses (EDA) show that the stronger hydrogen bonds in the unsaturated model complexes and in AT stem from stronger electrostatic interactions as well as enhanced donor-acceptor interactions in the σ-electron system, with the covalency being responsible for shortening the hydrogen bonds in these dimers.

Feasibility of occurrence of different types of protonated base pairs in RNA: a quantum chemical study

Protonated nucleobases have significant roles in facilitating catalytic functions of RNA, and in stabilizing different structural motifs. Reported pKa values of nucleobase protonation suggest that the population of neutral nucleobases is 10(3)-10(4) times higher than that of protonated nucleobases under physiological conditions (pH ∼ 7.4). Therefore, a molecular level understanding of various putative roles of protonated nucleobases cannot be achieved without addressing the question of how their occurrence propensities and stabilities are related to the free energy costs associated with the process of protonation under physiological conditions. With water as the proton donor, we use advanced QM methods to evaluate the site specific protonation propensities of nucleobases in terms of their associated free energy changes (ΔGprot). Quantitative follow up on the energetics of base pair formation and database search for evaluating their occurrence frequencies, reveal a lack of correlation between base pair stability and occurrence propensities on the one hand, and ease of protonation on the other. For example, although N7 protonated adenine (ΔGprot = 40.0 kcal mol(-1)) is found to participate in stable base pairing, base pairs involving N7 protonated guanine (ΔGprot = 36.8 kcal mol(-1)), on geometry optimization, converge to a minima where guanine transfers its extra proton to its partner base. Such observations, along with examples of weak base pairs involving N3 protonation of cytosine (ΔGprot = 37.0 kcal mol(-1)) are rationalized by analysing the protonation induced charge redistributions which are found to significantly influence, both positively and negatively, the hydrogen bonding potentials of different functional sites of individual nucleobases. Protonation induced charge redistribution is also found to strongly influence (i) the aromatic character of the rings of the participating bases and (ii) hydrogen bonding potential of the free edges of the protonated base pair. Comprehensive analysis of a non-redundant RNA crystal structure dataset further reveals that, while availability of stabilization possibilities determine the feasibility of occurrence of protonated bases, their occurrence context and specific functional roles are important factors determining their occurrence propensities.

Aggregation of nitroaniline in tetrahydrofuran through intriguing H-bond formation by sodium borohydride

The participation of sodium borohydride (NaBH4) in hydrogen bonding interactions and transient anion radical formation has been proved. Thus, the properties of NaBH4 are extended beyond the purview of its normal reducing capability and nucleophilic property. It is reported that ortho- and para-nitroanilines (NAs) form stable aggregates only in tetrahydrofuran (THF) in the presence of NaBH4 and unprecedented orange/red colorations are observed. The same recipe with nitrobenzene instead of nitroanilines (NAs) in the presence of NaBH4 evolves a transient rose red solution due to the formation of a highly fluorescent anion radical. Spectroscopic studies (UV-vis, fluorescence, RLS, Raman, NMR etc.) as well as theoretical calculations supplement the J-aggregate formation of NAs due to extensive hydrogen bonding. This is the first report where BH4(-) in THF has been shown to support such an aggregation process through H-bonding. It is further confirmed that stable intermolecular hydrogen bond-induced aggregation requires a geometrical match in both the nitro- and amino-functionalities attached to the phenyl ring with proper geometry. On the contrary, meta-nitroaniline remains as the odd man out and does not take part in such aggregation. Surprisingly, Au nanoparticles dismantle the J-aggregates of NA in THF. Explicit hydrogen bond formation in NA has been confirmed experimentally considering its promising applications in different fields including non-linear optics.

The role of N7 protonation of guanine in determining the structure, stability and function of RNA base pairs

Ab initio computations and bioinformatics studies reveal that stabilization of some important RNA structural motifs might involve N7 protonation of guanine.

Assessing the importance of proton transfer reactions in DNA

Although engineered by millions of years of evolution, the cellular machinery is not flawless, and errors regularly appear during DNA replication. The subsequent alteration of the stored genetic message results in a mutation and might be the starting point of important health disorders. The question therefore is what causes DNA mutations? All living organisms are constantly exposed to a number of external agents such as free radicals and to radiation, which may lead to induced mutations. There are also mutations happening without invoking the action of any exogenous element, the so-called spontaneous mutations. The former can be partially controlled by avoiding exposure to high-risk environments, while the latter are more intriguing because their origin is unclear and difficult to determine. As noted by Watson and Crick when they first discovered the DNA structure, the correct replication of DNA rests on the assumption that the base pairs remain in their most stable, canonical form. However, protons along the interbase hydrogen-bond network are not static entities. They can in fact interchange their positions in DNA bases through proton transfer (PT) reactions before strands unwind, giving rise to noncanonical structures defined as rare tautomers. The importance of these rare tautomers was also cleverly anticipated by Watson and Crick and some years later claimed by Löwdin to be a source of spontaneous mutations. In Watson and Crick's words: "It would be of interest to know the precise difference in free energy between the various tautomeric forms under physiological conditions." Unfortunately, rare tautomeric forms are very difficult to detect, so no direct and accurate free energy measure has been discerned. In contrast, theoretical chemistry is making good progress toward the quantification of PT reactions in DNA and their biological consequences. This Account touches upon the theoretical studies devoted to appraising the importance of rare tautomers as promoters of spontaneous mutations. We focus in particular on the crucial role played by the biological environment on DNA stability. It has now been demonstrated that valuable macroscopic predictions require not only highly accurate theories but also refined chemical models. Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations performed on short but complete DNA sequence fragments emerge in this context as the most adequate tools. In addition, these methods can be used to quantify the effect of different external agents on the PT tautomeric equilibria and, eventually, to conveniently handle them. This is the case for the possible alteration of the naturally observed mutation rate by exposure to intense electric fields. Theoretical predictions envision in this respect promising applications of ultrashort electric pulses in medicine to selectively modify the mutated/canonical ratio in DNA.

Design and Applications of Noncanonical DNA Base Pairs

While the Watson-Crick base pairs are known to stabilize the DNA double helix and play a vital role in storage/replication of genetic information, their replacement with non-Watson-Crick base pairs has recently been shown to have interesting practical applications. Nowadays, theoretical calculations are routinely performed on very complex systems to gain a better understanding of how molecules interact with each other. We not only bring together some of the basic concepts of how mispaired or unnatural nucleobases interact with each other but also look at how such an understanding influences the prediction of novel properties and development of new materials. We highlight the recent developments in this field of research. In this Perspective, we discuss the success of DFT methods, particularly, dispersion-corrected DFT, for applications such as pH-controlled molecular switching, electric-field-induced stacking of disk-like molecules with guanine quartets, and optical birefringence of alkali-metal-coordinated guanine quartets. The synergy between theoretical models and real applications is highlighted.

S. A, Pratik SM, Datta A, Nath R, Mandal S. Synthesis, structure, photocatalytic and magnetic properties of an oxo-bridged copper dimer

A new paramagnetic and semiconductor copper citrate dimer exhibits photocatalytic behaviour towards organic dye degradation.

Nucleic acid bases in anionic 2'-deoxyribonucleotides: a DFT/B3LYP study of structures, relative stability, and proton affinities

Protonation of nucleobases in anions of canonical 2'-deoxyribonucleotides has been investigated by the DFT computational study at the B3LYP/aug-cc-pvdz level of theory. It is demonstrated that the protonation leads to a significant decrease of conformational space of purine nucleotides while almost all conformers found for non-protonated molecules correspond to minima of the potential energy surface for protonated mdTMP and mdCMP. However, in all nucleotides, only one conformer is populated. This applies to all tautomers of protonated molecules except the mdTMP and mdCMP with the proton attached to the carbonyl group where a minor population of second conformer is observed. Protonation of nucleobase leads to significant elongation of the N-glycosidic bond. These findings agree well with suggestions that protonation of nucleobase is a first step in cleavage of the glycosidic bond. The oxygen atoms of both carbonyl groups of thymine and the N3 atom of the pyrimidine ring of cytosine, guanine, and adenine represent the most preferable sites for protonation of anions of 2'-deoxyrobonucleotides. The highest proton affinity is observed for the base in mdGMP and the lowest for the thymine moiety in mdTMP. It should be noted that calculated values of the proton affinities in anionic nucleotides are significantly higher (by 2-3 eV) than for nucleosides and neutral nucleotides. This allows assuming that the proton affinity of the base in DNA macromolecule may be tuned by changing the extent of shielding or neutralization of negative charge of the phosphate group.

Spectroscopic, electronic structure and natural bond analysis of 2-aminopyrimidine and 4-aminopyrazolo[3,4-d]pyrimidine: a comparative study

The FTIR and FT-Raman spectra of 2-aminopyrimidine (2-AP) and 4-aminopyrazolo[3,4-d]pyrimidine (4-APP) has been recorded in the region 4000-400 and 3500-100 cm(-1), respectively. The tautomeric stability, optimized geometry, frequency and intensity of the vibrational bands of 2-AP and 4-APP were obtained by the DFT level using 6-31G(d) and 6-31G(d,p) basis sets. The harmonic vibrational frequencies were calculated and the scaled values have been compared with experimental FTIR and FT-Raman spectra. A detailed interpretation of the infrared and Raman spectra of 2-AP and 4-APP are also reported based on total energy distribution (TED). The observed and the calculated frequencies are found to be in good agreement. The experimental spectra also coincide satisfactorily with those of theoretically simulated spectra. The (1)H and (13)C NMR spectra have been simulated using the gauge independent atomic orbital (GIAO) method. The theoretical UV-Vis spectrum of the compound using CIS method and the electronic properties, such as HOMO and LUMO energies, were performed by time-dependent DFT (TD-DFT) approach. The calculated HOMO and LUMO energies show that charge transfer occurs within molecule. The first order hyperpolarizability (β(0)) of these novel molecular system and related properties (β, α(0) and Δα) of 2-AP and 4-APP are calculated using DFT/6-31G(d) method on the finite-field approach. The Mulliken charges, the values of electric dipole moment (μ) of the molecule were computed using DFT calculations. The change in electron density (ED) in the σ(∗) antibonding orbitals and stabilization energies E(2) have been calculated by natural bond (NBO) analysis to give clear evidence of stabilization originating in the hyper conjugation of hydrogen-bonded interactions.

Cyclic AMP receptor protein-aequorin molecular switch for cyclic AMP

Molecular switches are designer molecules that combine the functionality of two individual proteins into one, capable of manifesting an "on/off" signal in response to a stimulus. These switches have unique properties and functionalities and thus, can be employed as nanosensors in a variety of applications. To that end, we have developed a bioluminescent molecular switch for cyclic AMP. Bioluminescence offers many advantages over fluorescence and other detection methods including the fact that there is essentially zero background signal in physiological fluids, allowing for more sensitive detection and monitoring. The switch was created by combining the properties of the cyclic AMP receptor protein (CRP), a transcriptional regulatory protein from E. Coli that binds selectively to cAMP with those of aequorin, a bioluminescent photoprotein native of the jellyfish Aequorea victoria . Genetic manipulation to split the genetic coding sequence of aequorin in two and genetically attach the fragments to the N and C termini of CRP resulted in a hybrid protein molecular switch. The conformational change experienced by CRP upon the binding of cyclic AMP is suspected to result in the observed loss of the bioluminescent signal from aequorin. The "on/off" bioluminescence can be modulated by cyclic AMP over a range of several orders of magnitude in a linear fashion in addition to the capacity to detect changes in cellular cyclic AMP of intact cells exposed to different external stimuli without the need to lyse the cells. We envision that the molecular switch could find applications in vitro as well as In Vivo cyclic AMP detection and/or imaging.