Performance of various density functionals for the hydrogen bonds in DNA base pairs

Catalytic Behavior of Hydrogen Bonded Water in Oligomerization of Silicates

Even though much research has been done to demonstrate the oligomerization of zeolites and silicates, there has been almost no study that investigates the role of hydrogen bonds in these reaction pathways. This study demonstrates the catalytic activity of hydrogen bonds in the silicate oligomerization reaction pathway. The presence of hydrogen bonding has been shown to enhance the energetic favorability of the anionic-I mechanism. Catalysis is prevalent in the Si-OH rupture process of the reaction pathway. Simultaneously, the dependence of the activation barrier on the equatorial or apical nature of the cleaving hydroxyl group has also been shown. The preceding steps such as condensation and fluxional influence the strength of hydrogen bonds. An increase in hydrogen bond strength enhances its catalytic effect, leading to a higher extent of reduction in the activation barrier of the particular reaction step. Even though the quantum study focuses on the oligomerization of calcium silicate as a test case, it can be anticipated that such similar effects can be perceived in general for the oligomerization of silicates containing metallic ions in sol-gel chemistry and zeolite synthesis.

Quantum and classical effects in DNA point mutations: Watson-Crick tautomerism in AT and GC base pairs

Proton transfer along the hydrogen bonds of DNA can lead to the creation of short-lived, but biologically relevant point mutations that can further lead to gene mutation and, potentially, cancer. In this work, the energy landscape of the canonical A-T and G-C base pairs (standard, amino-keto) to tautomeric A*-T* and G*-C* (non-standard, imino-enol) Watson-Crick DNA base pairs is modelled with density functional theory and machine-learning nudge-elastic band methods. We calculate the energy barriers and tunnelling rates of hydrogen transfer between and within each base monomer (A, T, G and C). We show that the role of tunnelling in A-T tautomerisation is statistically unlikely due to the presence of a small reverse reaction barrier. On the contrary, the thermal populations of the G*-C* point mutation could be non-trivial and propagate through the replisome. For the direct intramolecular transfer, the reaction is hindered by a substantial energy barrier. However, our calculations indicate that tautomeric bases in their monomeric form have remarkably long lifetimes.

Predicted structure and selectivity of 3d transition metal complexes with glutamic N,N-bis(carboxymethyl) acid

Structure and selectivity of 3d transition metal complexes with glutamic N,N-bis(carboxymethyl) acid are analyzed and predicted from DFT calculations.

Effect of Alkali Metal Cations on Length and Strength of Hydrogen Bonds in DNA Base Pairs

For many years, non-covalently bonded complexes of nucleobases have attracted considerable interest. However, there is a lack of information about the nature of hydrogen bonding between nucleobases when the bonding is affected by metal coordination to one of the nucleobases, and how the individual hydrogen bonds and aromaticity of nucleobases respond to the presence of the metal cation. Here we report a DFT computational study of nucleobase pairs interacting with alkali metal cations. The metal cations contribute to the stabilization of the base pairs to varying degrees depending on their position. The energy decomposition analysis revealed that the nature of bonding between nucleobases does not change much upon metal coordination. The effect of the cations on individual hydrogen bonds were described by changes in VDD charges on frontier atoms, H-bond length, bond energy from NBO analysis, and the delocalization index from QTAIM calculations. The aromaticity changes were determined by a HOMA index.

Mutual Relations between Substituent Effect, Hydrogen Bonding, and Aromaticity in Adenine-Uracil and Adenine-Adenine Base Pairs

The electronic structure of substituted molecules is governed, to a significant extent, by the substituent effect (SE). In this paper, SEs in selected nucleic acid base pairs (Watson-Crick, Hoogsteen, adenine-adenine) are analyzed, with special emphasis on their influence on intramolecular interactions, aromaticity, and base pair hydrogen bonding. Quantum chemistry methods-DFT calculations, the natural bond orbital (NBO) approach, the Harmonic Oscillator Model of Aromaticity (HOMA) index, the charge of the substituent active region (cSAR) model, and the quantum theory of atoms in molecules (QTAIM)-are used to compare SEs acting on adenine moiety and H-bonds from various substitution positions. Comparisons of classical SEs in adenine with those observed in para- and meta-substituted benzenes allow for the better interpretation of the obtained results. Hydrogen bond stability and its other characteristics (e.g., covalency) can be significantly changed as a result of the SE, and its consequences are dependent on the substitution position. These changes allow us to investigate specific relations between H-bond parameters, leading to conclusions concerning the nature of hydrogen bonding in adenine dimers-e.g., H-bonds formed by five-membered ring nitrogen acceptor atoms have an inferior, less pronounced covalent nature as compared to those formed by six-membered ring nitrogen. The energies of individual H-bonds (obtained by the NBO method) are analyzed and compared to those predicted by the Espinosa-Molins-Lecomte (EML) model. Moreover, both SE and H-bonds can significantly affect the aromaticity of adenine rings; long-distance SEs on π-electron delocalization are also documented.

On the catalytic transfer hydrogenation of nitroarenes by a cubane-type Mo3S4 cluster hydride: disentangling the nature of the reaction mechanism

Cubane-type Mo₃S₄ cluster hydrides decorated with phosphine ligands are active catalysts for the transfer hydrogenation of nitroarenes to aniline derivatives in the presence of formic acid (HCOOH) and triethylamine (Et₃N). The process is highly selective and most of the cluster species involved in the catalytic cycle have been identified through reaction monitoring. Formation of a dihydrogen cluster intermediate has also been postulated based on previous kinetic and theoretical studies. However, the different steps involved in the transfer hydrogenation from the cluster to the nitroarene to finally produce aniline remain unclear. Herein, we report an in-depth computational investigation into this mechanism. Et₃N reduces the activation barrier associated with the formation of Mo-HHOOCH dihydrogen species. The global catalytic process is highly exergonic and occurs in three consecutive steps with nitrosobenzene and N-phenylhydroxylamine as reaction intermediates. Our computational findings explain how hydrogen is transferred from these Mo-HHOOCH dihydrogen adducts to nitrobenzene with the concomitant formation of nitrosobenzene and the formate substituted cluster. Then, a β-hydride elimination reaction accompanied by CO₂ release regenerates the cluster hydride. Two additional steps are needed for hydrogen transfer from the dihydrogen cluster to nitrosobenzene and N-phenylhydroxylamine to finally produce aniline. Our results show that the three metal centres in the Mo₃S₄ unit act independently, so the cluster can exist in up to ten different forms that are capable of opening a wide range of reaction paths. This behaviour reveals the outstanding catalytic possibilities of this kind of cluster complexes, which work as highly efficient catalytic machines.

Solvent effect on the intermolecular proton transfer of the Watson and Crick guanine-cytosine and adenine-thymine base pairs: a polarizable continuum model study

Herein we present our results on the study of the double proton transfer (DPT) mechanism in the adenine-thymine (AT) and guanine-cytosine (GC) base pairs, both in gas phase and in solution. The latter was modeled using the polarizable continuum method (PCM) in different solvents. According to our DFT calculations, the DPT may occur for both complexes in a stepwise mechanism in condensate phase. In gas phase only the GC base pair exhibits a concerted DPT mechanism. Using the Wigner's tunneling corrections to the transition state theory we demonstrate that such corrections are important for the prediction of the rate constants of both systems in gas and in condensate phase. We also show that (i) as the polarity of the medium decreases the equilibrium constant of the DPT reaction increases in both complexes, and (ii) that the equilibrium constant in the GC complex is four orders of magnitude larger than in AT. This observation suggests that the spontaneous mutations in DNA base pairs are more probable in GC than in AT.

Introducing DDEC6 atomic population analysis: part 4. Efficient parallel computation of net atomic charges, atomic spin moments, bond orders, and more

The DDEC6 method is one of the most accurate and broadly applicable atomic population analysis methods. It works for a broad range of periodic and non-periodic materials with no magnetism, collinear magnetism, and non-collinear magnetism irrespective of the basis set type. First, we show DDEC6 charge partitioning to assign net atomic charges corresponds to solving a series of 14 Lagrangians in order. Then, we provide flow diagrams for overall DDEC6 analysis, spin partitioning, and bond order calculations. We wrote an OpenMP parallelized Fortran code to provide efficient computations. We show that by storing large arrays as shared variables in cache line friendly order, memory requirements are independent of the number of parallel computing cores and false sharing is minimized. We show that both total memory required and the computational time scale linearly with increasing numbers of atoms in the unit cell. Using the presently chosen uniform grids, computational times of ∼9 to 94 seconds per atom were required to perform DDEC6 analysis on a single computing core in an Intel Xeon E5 multi-processor unit. Parallelization efficiencies were usually >50% for computations performed on 2 to 16 cores of a cache coherent node. As examples we study a B-DNA decamer, nickel metal, supercells of hexagonal ice crystals, six X@C₆₀ endohedral fullerene complexes, a water dimer, a Mn₁₂-acetate single molecule magnet exhibiting collinear magnetism, a Fe₄O₁₂N₄C₄₀H₅₂ single molecule magnet exhibiting non-collinear magnetism, and several spin states of an ozone molecule. Efficient parallel computation was achieved for systems containing as few as one and as many as >8000 atoms in a unit cell. We varied many calculation factors (e.g., grid spacing, code design, thread arrangement, etc.) and report their effects on calculation speed and precision. We make recommendations for excellent performance.

We parallelize the DDEC6 method to efficiently compute net atomic charges, atomic spin moments, and bond orders in diverse materials.

The Adhesion Mechanism of Marine Mussel Foot Protein: Adsorption of L-Dopa onα- and β-Cristobalite Silica Using Density Functional Theory

Marine mussels strongly adhere to various surfaces and endure their attachment under a variety of conditions. In order to understand the basic mechanism involved, we study the adsorption of L-dopa molecule on hydrophilic geminal and terminal isolated silanols of silica (001) surface. High content of modified amino acid L-dopa is found in the glue-like material secreted by the mussels through which it sticks to various surfaces under water. To understand the adsorption behavior, we have made use of periodic Density Functional Theory (DFT) study. The L-dopa molecule adheres to silica surfaces terminated with geminal and terminal silanols via its catechol part. In both cases, the adhesion is achieved through the formation of 4 H-bonds. A binding energy of 29.48 and 31.67 kcal/mol has been estimated, after the inclusion of dispersion energy, for geminal and terminal silanols of silica, respectively. These results suggest a relatively stronger adhesion of dopa molecule for surface with terminal isolated silanols.

B-DNA model systems in non-terran bio-solvents: implications for structure, stability and replication

We have computationally analyzed a comprehensive series of Watson–Crick and mismatched B-DNA base pairs, in the gas phase and in several solvents, including toluene, chloroform, ammonia, methanol and water, using dispersion-corrected density functional theory and implicit solvation.

Spinning around in Transition-Metal Chemistry

The great diversity and richness of transition metal chemistry, such as the features of an open d-shell, opened a way to numerous areas of scientific research and technological applications. Depending on the nature of the metal and its environment, there are often several energetically accessible spin states, and the progress in accurate theoretical treatment of this complicated phenomenon is presented in this Account. The spin state energetics of a transition metal complex can be predicted theoretically on the basis of density functional theory (DFT) or wave function based methodology, where DFT has advantages since it can be applied routinely to medium-to-large-sized molecules and spin-state consistent density functionals are now available. Additional factors such as the effect of the basis set, thermochemical contributions, solvation, relativity, and dispersion, have been investigated by many researchers, but challenges in unambiguous assignment of spin states still remain. The first DFT studies showed intrinsic spin-state preferences of hybrid functionals for high spin and early generalized gradient approximation functionals for low spin. Progress in the development of density functional approximations (DFAs) then led to a class of specially designed DFAs, such as OPBE, SSB-D, and S12g, and brought a very intriguing and fascinating observation that the spin states of transition metals and the S_N2 barriers of organic molecules are somehow intimately linked. Among the many noteworthy results that emerged from the search for the appropriate description of the complicated spin state preferences in transition metals, we mainly focused on the examination of the connection between the spin state and the structures or coordination modes of the transition metal complexes. Changes in spin states normally lead only to changes in the metal-ligand bond lengths, but to the best of our knowledge, the dapsox ligand showed the first example of a transition-metal complex where a change in spin state leads also to changes in the coordination, switching between pentagonal-bipyramidal and capped-octahedron. Moreover, we have summarized the results of the thorough study that corrected the experimental assignment of the nature of the recently synthesized Sc³⁺ adduct of [Fe^IV(O)(TMC)]²⁺ (TMC = 1,4,8,11-tetramethylcyclam) and firmly established that the Sc³⁺-capped iron-oxygen complex corresponds to high-spin Fe^III. Last, but not least, we have provided deeper insight and rationalization of the observation that unlike in metalloenzymes, where the Fe^IV-oxo is usually observed with high spin, biomimetic Fe^IV-oxo complexes typically have a intermediate spin state. Energy decomposition analyses on the trigonal-bypiramidal (TBP) and octahedral model systems with ammonia ligands have revealed that the interaction energy of the prepared metal ion in the intermediate spin state is much smaller for the TBP structure. This sheds light on the origin of the intermediate spin state of the biomimetic TBP Fe^IV-oxo complexes.

Investigating the adsorption mechanism of glycine in comparison with catechol on cristobalite surface using density functional theory for bio-adhesive materials

Amino acid proteins exist in Mussel's adhesive (mefp's) of which glycine has a significant amount. A density functional theory simulation study was performed in a belief that all the proteins in mefp's are responsible for the versatile adhesion.

An atlas of RNA base pairs involving modified nucleobases with optimal geometries and accurate energies

Posttranscriptional modifications greatly enhance the chemical information of RNA molecules, contributing to explain the diversity of their structures and functions. A significant fraction of RNA experimental structures available to date present modified nucleobases, with half of them being involved in H-bonding interactions with other bases, i.e. ‘modified base pairs’. Herein we present a systematic investigation of modified base pairs, in the context of experimental RNA structures. To this end, we first compiled an atlas of experimentally observed modified base pairs, for which we recorded occurrences and structural context. Then, for each base pair, we selected a representative for subsequent quantum mechanics calculations, to find out its optimal geometry and interaction energy. Our structural analyses show that most of the modified base pairs are non Watson–Crick like and are involved in RNA tertiary structure motifs. In addition, quantum mechanics calculations quantify and provide a rationale for the impact of the different modifications on the geometry and stability of the base pairs they participate in.

Binding energies of tyrosine kinase inhibitors: Error assessment of computational methods for imatinib and nilotinib binding

The binding energies of imatinib and nilotinib to tyrosine kinase have been determined by quantum mechanical (QM) computations, and compared with literature binding energy studies using molecular mechanics (MM). The potential errors in the computational methods include these critical factors. Errors in X-ray structures such as structural distortions and steric clashes give unrealistically high van der Waals energies, and erroneous binding energies.MM optimization gives a very different configuration to the QM optimization for nilotinib, whereas the imatinib ion gives similar configurations. Solvation energies are a major component of the overall binding energy. The QM based solvent model (PCM/SMD) gives different values from those used in the implicit PBSA solvent MM models. A major error in inhibitor—kinase binding lies in the non-polar solvation terms. Solvent transfer free energies and the required empirical solvent accessible surface area factors for nilotinib and imatinib ion to give the transfer free energies have been reverse calculated. These values differ from those used in the MM PBSA studies.An intertwined desolvation—conformational binding selectivity process is a balance of thermodynamic desolvation and intramolecular conformational kinetic control.The configurational entropies (TΔS) are minor error sources.

Theoretical modelling of epigenetically modified DNA sequences

We report herein a set of calculations designed to examine the effects of epigenetic modifications on the structure of DNA. The incorporation of methyl, hydroxymethyl, formyl and carboxy substituents at the 5-position of cytosine is shown to hardly affect the geometry of CG base pairs, but to result in rather larger changes to hydrogen-bond and stacking binding energies, as predicted by dispersion-corrected density functional theory (DFT) methods. The same modifications within double-stranded GCG and ACA trimers exhibit rather larger structural effects, when including the sugar-phosphate backbone as well as sodium counterions and implicit aqueous solvation. In particular, changes are observed in the buckle and propeller angles within base pairs and the slide and roll values of base pair steps, but these leave the overall helical shape of DNA essentially intact. The structures so obtained are useful as a benchmark of faster methods, including molecular mechanics (MM) and hybrid quantum mechanics/molecular mechanics (QM/MM) methods. We show that previously developed MM parameters satisfactorily reproduce the trimer structures, as do QM/MM calculations which treat bases with dispersion-corrected DFT and the sugar-phosphate backbone with AMBER. The latter are improved by inclusion of all six bases in the QM region, since a truncated model including only the central CG base pair in the QM region is considerably further from the DFT structure. This QM/MM method is then applied to a set of double-stranded DNA heptamers derived from a recent X-ray crystallographic study, whose size puts a DFT study beyond our current computational resources. These data show that still larger structural changes are observed than in base pairs or trimers, leading us to conclude that it is important to model epigenetic modifications within realistic molecular contexts.

Insights into the structure of intrastrand cross-link DNA lesion-containing oligonucleotides: G[8-5m]T and G[8-5]C from molecular dynamics simulations

Oxidatively generated complex DNA lesions occur more rarely than single-nucleotide defects, yet they play an important role in carcinogenesis and aging diseases because they have proved to be more mutagenic than simple lesions. Whereas their formation pathways are rather well understood, the field suffers from the absence of structural data that are crucial for interpreting the lack of repair. No experimental structures are available for oligonucleotides featuring such a lesion. Hence, the detailed structural basis of such damaged duplexes has remained elusive. We propose the use of explicit solvent molecular dynamics simulations to build up damaged oligonucleotides containing two intrastrand cross-link defects, namely, the guanine-thymine and guanine-cytosine defects. Each of these lesions, G[8-5m]T and G[8-5]C, is placed in the middle of a dodecameric sequence, which undergoes an important structural rearrangement that we monitor and analyze. In both duplexes, the structural evolution is dictated by the more favorable stacking of guanine G6, which aims to restore π-stacking with the 3' purine nucleobase. Subsequently, transient formation of hydrogen bonds with a strand shifting is observed. Our simulations are combined with density functional theory to rationalize the structural evolution. We report converging computational evidence that the G[8-5m]T- and G[8-5]C-containing structures evolve toward "abasic-like" duplexes, with a stabilization of the interstrand pairing noncovalent interactions. Meanwhile, both lesions restore B-helicity within tens of nanoseconds. The identification of plausible structures characterizes the last hydrogen abstraction step toward the formation of such defects as a non-innocent chemical reaction.

Modelling proton tunnelling in the adenine–thymine base pair

Environment-assisted quantum tunnelling is unlikely to play a significant role in adenine–thymine tautomerisation.

Experimental measurement and theory of substituent effects in π-hydrogen bonding: complexes of substituted phenols with benzene

IR spectroscopic experiments and theoretical DFT computations reveal the effects of aromatic substituents on π-hydrogen bonding between monosubstituted phenol derivatives and benzene. Simultaneous formation of two π-hydrogen bonds (red-shifting O-H···π and blue-shifting ortho-C-H···π) contribute to the stability of these complexes. The interaction of the acidic phenol O-H proton-donating group with the benzene π-system dominates the complex formation. The experimental shifts of O-H stretching frequencies for the different phenol complexes vary in the range 45-74 cm(-1). Strong effects on hydrogen-bonding energies and frequency shifts of electron-withdrawing aromatic substituents and very weak influence of electron-donating groups have been established. Experimental quantities and theoretical parameters are employed in rationalizing the properties of these complexes. The acidities of the proton-donating phenols describe quantitatively the hydrogen-bonding process. The results obtained provide clear evidence that, when the structural variations are in the proton-donating species, the substituent effects on π-hydrogen bonding follow classic mechanisms, comprising both resonance and direct through-space influences. The performance of three alternative DFT functionals (B3LYP, B97-D, and PBE0 combined with the 6-311++G(2df,2p) basis set) in predicting the O-H frequency shifts upon complexation is examined. For comparison, O-H frequency shifts for several complexes were also determined at MP2/6-31++G(d,p).

A fundamental understanding of catechol and water adsorption on a hydrophilic silica surface: exploring the underwater adhesion mechanism of mussels on an atomic scale

Mussels have a remarkable ability to bond to solid surfaces under water. From a microscopic perspective, the first step of this process is the adsorption of dopa molecules to the solid surface. In fact, it is the catechol part of the dopa molecule that is interacting with the surface. These molecules are able to make reversible bonds to a wide range of materials, even underwater. Previous experimental and theoretical efforts have produced only a limited understanding of the mechanism and quantitative details of the competitive adsorption of catechol and water on hydrophilic silica surfaces. In this work, we uncover the nature of this competitive absorption by atomic scale modeling of water and catechol adsorbed at the geminal (001) silica surface using density functional theory calculations. We find that catechol molecules displace preadsorbed water molecules and bond directly on the silica surface. Using molecular dynamics simulations, we observe this process in detail. We also calculate the interaction force as a function of distance, and observe a maximum of 0.5 nN of attraction. The catechol has a binding energy of 23 kcal/mol onto the silica surface with adsorbed water molecules.

Quantum chemical investigations of the water exchange mechanism on [Al(III)(H₂O)₅(L)]²⁺ as a function of the donor strength of the anionic L

Water exchange reactions of the complexes Al(H₂O)₅(L)]²⁺·H₂O for L = →OCN⁻, F⁻, CF₃⁻, →NC⁻, →CN⁻, Cl⁻, Br⁻, H⁻, SH⁻, OH⁻, →NCO⁻, →NCS⁻, →SCN⁻, CF₃CH₂⁻, CH₃⁻, Et⁻, i-Pr⁻ and t-but⁻, were studied by DFT calculations (B3LYP/6-311+G**). The reactions follow a dissociative (D) pathway for Al(H₂O)₅(L)]²⁺·H₂O to form the five-coordinate intermediate [Al(H₂O)₄L]²⁺·2H₂O for the weaker donor ligands of the series. On increasing the donor strength of L, the five-coordinate intermediate becomes significantly more stable than the reactant state. At this point there is a mechanistic changeover to an associative (A) pathway for [Al(H₂O)₄L]²⁺·2H₂O as reactant to form a six-coordinate intermediate Al(H₂O)₅(L)]²⁺·H₂O. For some of the anionic ligands L = → NC⁻, →CN⁻, Cl⁻, OH⁻ and →SCN⁻) the energy gap between the reactant and intermediate states is small, such that the water exchange mechanism lies in the boarder of dissociative and associative pathways. The water exchange process involves cis- and trans-orientated transition states to form the product state that is similar to the reactant state.

Higher order structural effects stabilizing the reverse Watson-Crick Guanine-Cytosine base pair in functional RNAs

The G:C reverse Watson-Crick (W:W trans) base pair, also known as Levitt base pair in the context of tRNAs, is a structurally and functionally important base pair that contributes to tertiary interactions joining distant domains in functional RNA molecules and also participates in metabolite binding in riboswitches. We previously indicated that the isolated G:C W:W trans base pair is a rather unstable geometry, and that dicationic metal binding to the Guanine base or posttranscriptional modification of the Guanine can increase its stability. Herein, we extend our survey and report on other H-bonding interactions that can increase the stability of this base pair. To this aim, we performed a bioinformatics search of the PDB to locate all the occurencies of G:C trans base pairs. Interestingly, 66% of the G:C trans base pairs in the PDB are engaged in additional H-bonding interactions with other bases, the RNA backbone or structured water molecules. High level quantum mechanical calculations on a data set of representative crystal structures were performed to shed light on the structural stability and energetics of the various crystallographic motifs. This analysis was extended to the binding of the preQ1 metabolite to a preQ1-II riboswitch.

B-DNA Structure and Stability as Function of Nucleic Acid Composition: Dispersion-Corrected DFT Study of Dinucleoside Monophosphate Single and Double Strands

We have computationally investigated the structure and stability of all 16 combinations of two out of the four natural DNA bases A, T, G and C in a di-2′-deoxyribonucleoside-monophosphate model DNA strand as well as in 10 double-strand model complexes thereof, using dispersion-corrected density functional theory (DFT-D). Optimized geometries with B-DNA conformation were obtained through the inclusion of implicit water solvent and, in the DNA models, of sodium counterions, to neutralize the negative charge of the phosphate groups. The results obtained allowed us to compare the relative stability of isomeric single and double strands. Moreover, the energy of the Watson–Crick pairing of complementary single strands to form double-helical structures was calculated. The latter furnished the following increasing stability trend of the double-helix formation energy: d(TpA)₂ <d(CpA)₂ <d(ApT)₂ <d(ApA)₂ <d(GpT)₂ <d(GpA)₂ <d(ApG)₂ <d(CpG)₂ <d(GpG)₂ <d(GpC)₂, where the energy differences between the last four dimers, d(ApG)₂, d(CpG)₂, d(GpG)₂ and d(GpC)₂, is within 4.0 kcal mol⁻¹, and the energy between the most and the least stable isomers is 13.4 kcal mol⁻¹. This trend shows that the formation energy essentially increases with the number of hydrogen bonds per base pair, that is two between A and T and three between G and C. Superimposed on this main trend are more subtle effects that depend on the order in which bases occur within a strand from the 5’- to the 3’-end.

Theoretical investigation of the coupling between hydrogen-atom transfer and stacking interaction in adenine-thymine dimers

Three different dimers of the adenine-thymine (A-T) base pair are studied to point out the changes of important properties (structure, atomic charge, energy and so on) induced by coupling between the movement of the atoms in the hydrogen bonds and the stacking interaction. The comparison of these results with those for the A-T monomer system explains the role of the stacking interaction in the hydrogen-atom transfer in this biologically important base pair. The results support the idea that this coupling depends on the exact dimer considered and is different for the N-N and N-O hydrogen bonds. In particular, the correlation between the hydrogen transfer and the stacking interaction is more relevant for the N-N bridge than for the N-O one. Also, the two different mechanisms of two-hydrogen transfer (step by step and concerted) can be modified by the stacking interaction between the base pairs.

Ab initio determination of coarse-grained interactions in double-stranded DNA

We derive the coarse-grained interactions between DNA nucleotides from ab initio total-energy calculations based on density functional theory (DFT). The interactions take into account base and sequence specificity, and are decomposed into physically distinct contributions that include hydrogen bonding, stacking interactions, backbone, and backbone-base interactions. The interaction energies of each contribution are calculated from DFT for a wide range of configurations and are fitted by simple analytical expressions for use in the coarse-grained model, which reduces each nucleotide into two sites. This model is not derived from experimental data, yet it successfully reproduces the stable B-DNA structure and gives good predictions for the persistence length. It may be used to realistically probe dynamics of DNA strands in various environments at the μs time scale and the μm length scale.

The origin of isotope-induced helical-sense bias in supramolecular polymers of benzene-1,3,5-tricarboxamides

The molecular origin of the isotope-induced diastereomeric enrichment in helical supramolecular polymers consisting of trialkylbenzene-1,3,5-tricarboxamides (BTAs) is studied using plane-wave DFT calculations. We demonstrate that the creation of a chiral center at the α-position of the alkyl chains of a BTA by H-D exchange leads to a small but notable preference for the formation of supramolecular hydrogen bonded structures with a particular helicity. The bias for one helical sense preference is caused by the orientation of the vibrational eigenmodes of the C-H and C-D stretching frequencies at the chiral center and by hyperconjugative destabilization of the anti C-H orbital.