Hydrogen bonding and stacking of DNA bases: a review of quantum-chemical ab initio studies

A density functional theory study on interactions in water-bridged dimeric complexes of lignin

Lignin is the main plant cell wall component responsible for recalcitrance in the process of lignocellulosic biomass conversion into biofuels. The recalcitrance and insolubility of lignin in different reaction media are due in part to the hydrogen bonds and π interactions that hold syringyl (S) and guaiacyl (G) units together and promote the formation of stable water-bridged dimeric complexes (WBDCs): S⋯G and S⋯S, in native lignin. The current understanding of how each type of interaction influences the stability of these complexes within lignin native cell walls is still limited. Here, we found by DFT calculations that hydrogen bonding is more dominant than π-stacking interaction between aromatic rings of WBDCs. Although there is a stronger interaction of hydrogen bonds between subunits and water and higher π-stacking interaction in the S⋯S complex compared to the S⋯G complex, the former complex is less thermodynamically stable than the latter due to the entropic contribution coming from the methoxy substituents in the S-unit. Our results demonstrate that the methoxylation degree of lignin units does not significantly influence the structural geometries of WBDCs; if anything, an enhanced dispersion interaction between ring aromatics results in quasi-sandwich geometries as found in "coiled" lignin structures in the xylem tissue of wood. In the same way as that with ionic liquids, polar solvents can dissolve S-lignin by favorable interactions with the aliphatic hydroxyl group in the α-position as the key site or the aromatic hydroxyl group as the secondary site.

Ab initio Modeling of Hydrogen Bonding of Remdesivir and Adenosine with Uridine

Remdesivir (RDV) emerged as an effective drug against the SARS-CoV-2 virus pandemic. One of the crucial steps in the mechanism of action of RDV is its incorporation into the growing RNA strand. RDV, an adenosine analogue, forms Watson-Crick (WC) type hydrogen bonds with uridine in the complementary strand and the strength of this interaction will control efficacy of RDV. While there is a plethora of structural and energetic information available about WC H-bonds in natural base pairs, the interaction of RDV with uridine has not been studied yet at the atomic level. In this article, we aim to bridge this gap, to understand RDV and its hydrogen bonding interactions, by employing density functional theory (DFT) at the M06-2X/cc-pVDZ level. The interaction energy, QTAIM analysis, NBO and SAPT2 are performed for RDV, adenosine, and their complex with uridine to gain insights into the nature of hydrogen bonding. The computations show that RDV has similar geometry, energetic, molecular orbitals, and aromaticity as adenosine, suggesting that RDV is an effective adenosine analogue. The important geometrical parameters, such as bond distances and red-shift in the stretching vibrational modes of adenosine, RDV and uridine identify two WC-type H-bonds. The relative strength of these two H-bonds is computed using QTAIM parameters and the computed hydrogen bond energy. Finally, the SAPT2 study is performed at the minima and at non-equilibrium base pair distances to understand the dominant intermolecular physical force. This study, based on a thorough analysis of a variety of computations, suggests that both adenosine and RDV have similar structure, energetic, and hydrogen bonding behaviour.

Human Translesion Synthesis Polymerases polκ and polη Perform Error-Free Replication across N2-dG Methyleugenol and Estragole DNA Adducts

The secondary metabolites of polypropanoids, methyleugenol (MEG), and estragole (EG), found in many herbs and spices, are commonly used as food flavoring agents and as ingredients in cosmetics. MEG and EG have been reported to cause hepatocarcinogenicity in rodents, human livers, and lung cells. The formation of N2-dG and N6-dA DNA adducts is primarily attributed to the carcinogenicity of these compounds. Therefore, these compounds have been classified as "possible human carcinogens" by the International Agency for Research on Cancer and "reasonably anticipated to be a human carcinogen" by the National Toxicology Program. Herein, we report the synthesis of the N2-MEG-dG and N2-EG-dG modified oligonucleotides to study the mutagenicity of these DNA adducts. Our studies show that N2-MEG-dG and N2-EG-dG could be bypassed by human translesion synthesis (TLS) polymerases hpolκ and hpolη in an error-free manner. The steady-state kinetics of dCTP incorporation by hpolκ across N2-MEG-dG and N2-EG-dG adducts show that the catalytic efficiencies (kcat/Km) were ∼2.5- and ∼4.4-fold higher, respectively, compared to the unmodified dG template. A full-length primer extension assay demonstrates that hpolκ exhibits better catalytic efficiency than hpolη. Molecular modeling and dynamics studies capturing pre-insertion, insertion, and post-insertion steps reveal the structural features associated with the efficient bypass of the N2-MEG-dG adduct by hpolκ and indicate the reorientation of the adduct in the active site allowing the successful insertion of the incoming nucleotide. Together, these results suggest that though hpolκ and hpolη perform error-free TLS across MEG and EG during DNA replication, the observed carcinogenicity of these adducts could be attributed to the involvement of other low fidelity polymerases.

Mechanism of Self-Healing Hydrogels and Application in Tissue Engineering

Self-healing hydrogels and traditional hydrogels both have three-dimensional polymeric networks that are capable of absorbing and retaining a large amount of water. Self-healing hydrogels can heal and restore damage automatically, and they can avoid premature failure of hydrogels caused by mechanical damage after implantation. The formation mechanism of self-healing hydrogels and the factors that hydrogels can load are various. Researchers can design hydrogels to meet the needs of different tissues through the diversity of hydrogels Therefore, it is necessary to summarize different self-healing mechanisms and different factors to achieve different functions. Here, we briefly reviewed the hydrogels designed by researchers in recent years according to the self-healing mechanism of water coagulation. Then, the factors for different functions of self-healing hydrogels in different tissues were statistically analyzed. We hope our work can provide effective support for researchers in the design process of self-healing hydrogel.

Nearest-Neighbor dsDNA Stability Analysis Using Alchemical Free-Energy Simulations

The thermodynamic stability of double-stranded (ds)DNA depends on its sequence. It is influenced by the base pairing and stacking with neighboring bases along DNA molecules. Semiempirical schemes are available that allow us to predict the thermodynamic stability of DNA sequences based on empirically derived nearest-neighbor contributions of base pairs formed in the context of all possible nearest-neighbor base pairs. Current molecular dynamics (MD) simulations allow one to simulate the dynamics of DNA molecules in good agreement with experimentally obtained structures and available data on conformational flexibility. However, the suitability of current force field methods to reproduce dsDNA stability and its sequence dependence has been much less well tested. We have employed alchemical free-energy simulations of whole base pair transversions in dsDNA and in unbound single-stranded partner molecules. Such transversions change the sequence context but not the nucleotide content or base pairing in dsDNA and allow a direct comparison with the empirical nearest-neighbor dsDNA stability model. For the alchemical free-energy changes in the unbound single-stranded (ss)DNA partner molecules, we tested different setups assuming either complete unstacking or unrestrained simulations with partial stacking in the unbound ssDNA. The free-energy simulations predicted nearest-neighbor effects of similar magnitude, as observed experimentally but showed overall limited correlation with experimental data. An inaccurate description of stacking interactions and other possible reasons such as the neglect of electronic polarization effects are discussed. The results indicate the need to improve the realistic description of stacking interactions in current molecular mechanic force fields.

Base-Pairing and Base-Stacking Contributions to Double-Stranded DNA Formation

Double-stranded (ds)DNA formation and dissociation are of fundamental biological importance. The negative DNA charge influences the dsDNA stability. However, the base pairing and the stacking between neighboring bases are responsible for the sequence-dependent stability of dsDNA. The stability of a dsDNA molecule can be estimated from empirical nearest-neighbor models based on contributions assigned to base-pair steps along the DNA and additional parameters because of DNA termini. In efforts to separate contributions, it has been concluded that base stacking dominates dsDNA stability, whereas base pairing contributes negligibly. Using a different model for dsDNA formation, we reanalyze dsDNA stability contributions and conclude that base stacking contributes already at the level of separate ssDNAs but that pairing contributions drive the dsDNA formation. The theoretical model also predicts that stability contributions of base-pair steps that contain only guanine/cytosine, mixed steps, and steps with only adenine/thymine follow the order 6:5:4, respectively, as expected based on the formed hydrogen bonds. The model is fully consistent with the available stacking data and the nearest-neighbor dsDNA parameters. It allows assigning a narrowly distributed value for the effective free energy contribution per formed hydrogen bond during dsDNA formation of -0.72 kcal·mol^-1 based entirely on the experimental data.

Label-Free Oligonucleotide-Based SPR Biosensor for the Detection of the Gene Mutation Causing Prothrombin-Related Thrombophilia

Prothrombin-related thrombophilia is a genetic disorder produced by a substitution of a single DNA base pair, replacing guanine with adenine, and is detected mainly by polymerase chain reaction (PCR). A suitable alternative that could detect the single point mutation without requiring sample amplification is the surface plasmon resonance (SPR) technique. SPR biosensors are of great interest: they offer a platform to monitor biomolecular interactions, are highly selective, and enable rapid analysis in real time. Oligonucleotide-based SPR biosensors can be used to differentiate complementary sequences from partially complementary or noncomplementary strands. In this work, a glass chip covered with an ultrathin (50 nm) gold film was modified with oligonucleotide strands complementary to the mutated or normal (nonmutated) DNA responsible for prothrombin-related thrombophilia, forming two detection platforms called mutated thrombophilia (MT) biosensor and normal thrombophilia (NT) biosensor. The results show that the hybridization response is obtained in 30 min, label free and with high reproducibility. The sensitivity obtained in both systems was approximately 4 ΔμRIU/nM. The dissociation constant and limits of detection calculated were 12.2 nM and 20 pM (3 fmol), respectively, for the MT biosensor, and 8.5 nM and 30 pM (4.5 fmol) for the NT biosensor. The two biosensors selectively recognize their complementary strand (mutated or normal) in buffer solution. In addition, each platform can be reused up to 24 times when the surface is regenerated with HCl. This work contributes to the design of the first SPR biosensor for the detection of prothrombin-related thrombophilia based on oligonucleotides with single point mutations, label-free and without the need to apply an amplification method.

Theoretical Studies of IR and NMR Spectral Changes Induced by Sigma-Hole Hydrogen, Halogen, Chalcogen, Pnicogen, and Tetrel Bonds in a Model Protein Environment

Various types of σ-hole bond complexes were formed with FX, HFY, H₂FZ, and H₃FT (X = Cl, Br, I; Y = S, Se, Te; Z = P, As, Sb; T = Si, Ge, Sn) as Lewis acid. In order to examine their interactions with a protein, N-methylacetamide (NMA), a model of the peptide linkage was used as the base. These noncovalent bonds were compared by computational means with H-bonds formed by NMA with XH molecules (X = F, Cl, Br, I). In all cases, the A-F bond, which lies opposite the base and is responsible for the σ-hole on the A atom (A refers to the bridging atom), elongates and its stretching frequency undergoes a shift to the red with a band intensification, much as what occurs for the X-H bond in a H-bond (HB). Unlike the NMR shielding decrease seen in the bridging proton of a H-bond, the shielding of the bridging A atom is increased. The spectroscopic changes within NMA are similar for H-bonds and the other noncovalent bonds. The C=O bond of the amide is lengthened and its stretching frequency red-shifted and intensified. The amide II band shifts to higher frequency and undergoes a small band weakening. The NMR shielding of the O atom directly involved in the bond rises, whereas the C and N atoms both undergo a shielding decrease. The frequency shifts of the amide I and II bands of the base as well as the shielding changes of the three pertinent NMA atoms correlate well with the strength of the noncovalent bond.

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.

Self-Healing Hydrogels: The Next Paradigm Shift in Tissue Engineering?

Given their durability and long-term stability, self-healable hydrogels have, in the past few years, emerged as promising replacements for the many brittle hydrogels currently being used in preclinical or clinical trials. To this end, the incompatibility between hydrogel toughness and rapid self-healing remains unaddressed, and therefore most of the self-healable hydrogels still face serious challenges within the dynamic and mechanically demanding environment of human organs/tissues. Furthermore, depending on the target tissue, the self-healing hydrogels must comply with a wide range of properties including electrical, biological, and mechanical. Notably, the incorporation of nanomaterials into double-network hydrogels is showing great promise as a feasible way to generate self-healable hydrogels with the above-mentioned attributes. Here, the recent progress in the development of multifunctional and self-healable hydrogels for various tissue engineering applications is discussed in detail. Their potential applications within the rapidly expanding areas of bioelectronic hydrogels, cyborganics, and soft robotics are further highlighted.

DFT investigation of role of N - H⋯O and N - H⋯π interactions in the stabilization of the hydrogen bonded complexes of anisole with aromatic amines

Theoretical investigations have been performed on hydrogen (H-) bonded complexes of two aromatic amines with anisole to investigate the effect of the methyl substituent on N – H⋯O and N – –H⋯π interactions. Natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM) analyses were done to elucidate the nature of H- bonding. In 1:1 complexes, the total interaction energy of N-methylaniline complex is higher than that of aniline complex. The existence of bond critical point between N–H of amine and oxygen of anisole confirms weak hydrogen bonding. The energy decomposition analysis showed the role of CT in stabilizing complexes.

Spectroscopic study of uracil, 1-methyluracil and 1-methyl-4-thiouracil: Hydrogen bond interactions in crystals and ab-initio molecular dynamics

Hydrogen bond networks in uracil, 1-methyluracil and 1-methyl-4-thiouracil were studied by ab initio molecular dynamics as well as analysis of the orbital interactions. The power spectra calculated by ab initio molecular dynamics for atoms involved in hydrogen bonds were analyzed. We calculated spectra by using anharmonic approximation based on the autocorrelation function of the atom positions obtained from the Born-Oppenheimer simulations. Our results show the differences between hydrogen bond networks in uracil and its methylated derivatives. The studied methylated derivatives, 1-methyluracil as well as 1-methyl-4-thiouracil, form dimeric structures in the crystal phase, while uracil does not form that kind of structures. The presence of sulfur atom instead oxygen atom reflects weakness of the hydrogen bonds that build dimers.

An Electron Density Source-Function Study of DNA Base Pairs in Their Neutral and Ionized Ground States†

The source function (SF) decomposes the electron density at any point into contributions from all other points in the molecule, complex, or crystal. The SF "illuminates" those regions in a molecule that most contribute to the electron density at a point of reference. When this point of reference is the bond critical point (BCP), a commonly used surrogate of chemical bonding, then the SF analysis at an atomic resolution within the framework of Bader's Quantum Theory of Atoms in Molecules returns the contribution of each atom in the system to the electron density at that BCP. The SF is used to locate the important regions that control the hydrogen bonds in both Watson-Crick (WC) DNA dimers (adenine:thymine (AT) and guanine:cytosine (GC)) which are studied in their neutral and their singly ionized (radical cationic and anionic) ground states. The atomic contributions to the electron density at the BCPs of the hydrogen bonds in the two dimers are found to be delocalized to various extents. Surprisingly, gaining or loosing an electron has similar net effects on some hydrogen bonds concealing subtle compensations traced to atomic sources contributions. Coarser levels of resolutions (groups, rings, and/or monomers-in-dimers) reveal that distant groups and rings often have non-negligible effects especially on the weaker hydrogen bonds such as the third weak CH⋅⋅⋅O hydrogen bond in AT. Interestingly, neither the purine nor the pyrimidine in the neutral or ionized forms dominate any given hydrogen bond despite that the former has more atoms that can act as source or sink for the density at its BCP. © 2018 Wiley Periodicals, Inc.

Occurrence and stability of lone pair-π stacking interactions between ribose and nucleobases in functional RNAs

The specific folding pattern and function of RNA molecules lies in various weak interactions, in addition to the strong base-base pairing and stacking. One of these relatively weak interactions, characterized by the stacking of the O4' atom of a ribose on top of the heterocycle ring of a nucleobase, has been known to occur but has largely been ignored in the description of RNA structures. We identified 2015 ribose-base stacking interactions in a high-resolution set of non-redundant RNA crystal structures. They are widespread in structured RNA molecules and are located in structural motifs other than regular stems. Over 50% of them involve an adenine, as we found ribose-adenine contacts to be recurring elements in A-minor motifs. Fewer than 50% of the interactions involve a ribose and a base of neighboring residues, while approximately 30% of them involve a ribose and a nucleobase at least four residues apart. Some of them establish inter-domain or inter-molecular contacts and often implicate functionally relevant nucleotides. In vacuo ribose-nucleobase stacking interaction energies were calculated by quantum mechanics methods. Finally, we found that lone pair-π stacking interactions also occur between ribose and aromatic amino acids in RNA-protein complexes.

An Extraordinarily Stable DNA Minidumbbell

The minidumbbell (MDB) is a new type of native DNA structure. At neutral pH, two TTTA or CCTG repeats can fold into the highly compact MDB with a melting temperature of ∼22 °C. Owing to the relatively low thermodynamic stability, MDBs have been proposed to be the structural intermediates that lead to efficient DNA repair escape and thus repeat expansions. In this study, we reveal that two CCTG repeats can also form an extraordinarily stable MDB with a melting temperature of ∼46 °C at pH 5.0. This unusual stability predominantly results from the formation of a three hydrogen bond C⁺·C mispair between the two minor groove cytosine residues. Due to the drastic stability change, the CCTG MDB, when combined with its complementary sequence, shows instant and complete structural conversions when the pH switches between 5.0 and 7.0, making the system serve as a simple and efficient pH-controlled molecular switch.

2,4-Ditellurouracil and its 5-fluoro derivative: Theoretical investigations of structural, energetics and ADME parameters

2,4-Ditellurouracil exhibits keto-enol tautomerism via different pathways resulting in seven tautomers. These pathways were studied in the gas phase using density functional theory method. The functionals used were BLYP, B3LYP and BHLYP and the basis sets were 6-311++G(d,p) for all atoms except that LanL2DZ ECP was used for tellurium atom only. The results indicate that the diketo form is more stable as observed for uracil and its sulfur and selenium analogues. The effect of introducing fluorine at position 5 was also investigated and the energy difference between the diketo and dienol forms is reduced. 2,4-Ditellurouracil and its 5-fluoro analogue are expected to exist exclusively as the diketo form due to the high interconversion energy barrier. We extended the investigation to predict ADME parameters of the most stable diketo and dienol tautomers in view of understanding their biological properties. This research enlightens keto-enol tautomerism of 2,4-ditellurouracil and its 5-fluoro derivative with additional insights to biological functions.

Cooperative vibrational properties of hydrogen bonds in Watson–Crick DNA base pairs

For the AT pair, Symst and Strech peaks further shift toward the red, giving the H-bonds an amplified effect (orange arrows).

The effect of sequence on the ionization of guanine in DNA

The accurate estimation of the ionization energies and understanding the nature of the ionized states of the nucleic acid bases (NABs) are crucial to the understanding of the DNA damage mechanism. The vertical ionization energy (VIE) of guanine is the lowest among the NABs and the ionization energies are strongly affected by the environment, such as solvation and characteristics of nearby NABs. Therefore, we investigate the sequence dependence of the VIEs of guanine in B-DNA. We use the equation of motion coupled cluster method for the estimation of ionization potential with single and double excitations (EOM-IP-CCSD) and density functional theory with dispersion corrected ωB97x-D for the estimation of VIEs. A significant amount of non-additivity or cooperativity, directly proportional to charge delocalization, is noticed in the change in VIE due to the interaction with the nearby NABs. While the change in VIE due to base pairing originates predominantly from charge-dipole interactions, stacking between base pairs is a more complicated balance of dispersion and charge-dipole interactions as well as stabilization due to the delocalization of the positive charge. The long range interactions are however dominated by 1/r(3) distance dependence which shows the major role played by charge-dipole interactions. The extent of localization of positive holes on guanine is also estimated for various sequences.

Hydrogen Bonding and Stability of Hybrid Organic-Inorganic Perovskites

In the past few years, the efficiency of solar cells based on hybrid organic-inorganic perovskites has exceeded the level needed for commercialization. However, existing perovskites solar cells (PSCs) suffer from several intrinsic instabilities, which prevent them from reaching industrial maturity, and stabilizing PSCs has become a critically important problem. Here we propose to stabilize PSCs chemically by strengthening the interactions between the organic cation and inorganic anion of the perovskite framework. In particular, we show that replacing the methylammonium cation with alternative protonated cations allows an increase in the stability of the perovskite by forming strong hydrogen bonds with the halide anions. This interaction also provides opportunities for tuning the electronic states near the bandgap. These mechanisms should have a universal character in different hybrid organic-inorganic framework materials that are widely used.

Aromatic Base Stacking in DNA: From ab initio Calculations to Molecular Dynamics Simulations

Abstract Aromatic stacking of nucleic acid bases is one of the key players in determining the structure and dynamics of nucleic acids. The arrangement of nucleic acid bases with extensive overlap of their aromatic rings gave rise to numerous often contradictory suggestions about the physical origins of stacking and the possible role of delocalized electrons in stacked aromatic π systems, leading to some confusion about the issue. The recent advance of computer hardware and software finally allowed the application of state of the art quantum-mechanical approaches with inclusion of electron correlation effects to study aromatic base stacking, now providing an ultimitate qualitative description of the phenomenon. Base stacking is determined by an interplay of the three most commonly encountered molecular interactions: dispersion attraction, electrostatic interaction, and short-range repulsion. Unusual (aromatic- stacking specific) energy contributions were in fact not evidenced and are not necessary to describe stacking. The currently used simple empirical potential form, relying on atom-centered constant point charges and Lennard-Jones van der Waals terms, is entirely able to reproduce the essential features of base stacking. Thus, we can conclude that base stacking is in principle one of the best described interactions in current molecular modeling and it allows to study base stacking in DNA using large-scale classical molecular dynamics simulations. Neglect of cooperativity of stacking appears to be the most serious approximation of the currently used force field form. This review summarizes recent developments in the field. It is written for an audience that is not necessarily expert in computational quantum chemistry and follows up on our previous contribution (Sponer et. al., J. Biomol. Struct. Dyn. 14, 117, (1997)). First, the applied methodology, its accuracy, and the physical nature of base stacking is briefly overviewed, including a comment on the accuracy of other molecular orbital methods and force fields. Then, base stacking is contrasted with hydrogen bonding, the other dominant force in nucleic acid structure. The sequence dependence and cooperativity of base stacking is commented on, and finally a brief introduction into recent progress in large-scale molecular dynamics simulations of nucleic acids is provided. Using four stranded DNA assemblies as an example, we demonstrate the efficacy of current molecular dynamics techniques that utilize refined and verified force fields in the study of stacking in nucleic acid molecules.

Hydrogen bonding and π–π interactions in imidazolium-chloride ionic liquid clusters

The importance of 1° and 2° hydrogen-bonding and anion–π⁺ interactions for ionic liquid structuring.

Influence of Magnetic Microparticles Isolation on Adenine Homonucleotides Structure

The electroactivity of purine and pyrimidine bases is the most important property of nucleic acids that is very useful for determining oligonucleotides using square wave voltammetry. This study was focused on the electrochemical behavior of adenine-containing oligonucleotides before and after their isolation using paramagnetic particles. Two peaks were detected-peak A related to the reduction of adenine base and another peak B involved in the interactions between individual adenine strands and contributes to the formation of various spatial structures. The influence of the number of adenine bases in the strand in the isolation process using paramagnetic particles was investigated too.

Computational Evaluation of Intermolecular Interactions of a Universal Base 3-Nitropyrrole in Stacked Dimers and DNA Duplexes

The stacking interactions between a universal base of 3-nitropyrrole (3NP) and four canonical nucleobases were studied by means of ab initio molecular orbital calculations. The stabilities of the complexes are comparable to those of the stacked dimers of canonical bases reported previously. The detailed analysis of the interaction energies revealed the importance of the dipole-dipole interaction included in the Hartree-Fock terms to determine the geometry dependence of the stacking energies. It was also clarified that the dispersion energies included in the electron-correlation terms were essential to obtain adequate stabilities. The contribution of the nitro group was evaluated by the comparative studies of pyrrole and 3NP. The increased molecular dipole moment and surface are expected to account for the enhancement of the stability of the stacked dimers containing 3NP. The force field parameters required for calculation of the molecular mechanics of 3NP were obtained for 3NP on the basis of these molecular orbital calculations. The energy-minimized structures obtained by the molecular mechanics calculations of 3NP accorded with those obtained by the molecular orbital calculations described above. A DNA duplex structure containing 3NP-A, 3NP-T, or 3NP-C was calculated by use of these force field parameters. In the case of 3NP-A, the computationally calculated structure was in good agreement with that previously determined by use of (1)H-NMR except for the orientation of the nitro group.

Intramolecular CH···O hydrogen bonds in the AI and BI DNA-like conformers of canonical nucleosides and their Watson-Crick pairs. Quantum chemical and AIM analysis

The aim of this work is to cast some light on the H-bonds in double-stranded DNA in its AI and BI forms. For this purpose, we have performed the MP2 and DFT quantum chemical calculations of the canonical nucleoside conformers, relative to the AI and BI DNA forms, and their Watson-Crick pairs, which were regarded as the simplest models of the double-stranded DNA. Based on the atoms-in-molecules analysis (AIM), five types of the CH···O hydrogen bonds, involving bases and sugar, were detected numerically from 1 to 3 per a conformer: C2'H···O5', C1'H···O2, C6H···O5', C8H···O5', and C6H···O4'. The energy values of H-bonds occupy the range of 2.3-5.6 kcal/mol, surely exceeding the kT value (0.62 kcal/mol). The nucleoside CH···O hydrogen bonds appeared to "survive" turns of bases against the sugar, sometimes in rather large ranges of the angle values, pertinent to certain conformations, which points out to the source of the DNA lability, necessary for the conformational adaptation in processes of its functioning. The calculation of the interactions in the dA·T nucleoside pair gives evidence, that additionally to the N6H···O4 and N1···N3H canonical H-bonds, between the bases adenine and thymine the third one (C2H···O2) is formed, which, though being rather weak (about 1 kcal/mol), satisfies the AIM criteria of H-bonding and may be classified as a true H-bond. The total energy of all the CH···O nontraditional intramolecular H-bonds in DNA nucleoside pairs appeared to be commensurable with the energy of H-bonds between the bases in Watson-Crick pairs, which implies their possible important role in the DNA shaping.

[Molecular docking: role of intermolecular contacts in formation of complexes of proteins with nucleotides and peptides]

Knowledge of 3D-structure of protein-ligand complex is a major prerequisite for understanding the functioning mechanism of cellular proteins and membrane receptors. This is also of a great help in rational drug design projects. In the present paper we briefly review the molecular docking approaches used to predict possible orientation of a ligand in the protein binding site. The recent trends to improve the accuracy and efficiency of docking algorithms are demonstrated with the results obtained in Laboratory of Biomolecular Modeling. Particular attention is paid to protein-ligand hydrophobic and stacking interactions responsible for molecular recognition of ligand fragments. Such type of interactions are not always adequately represented in scoring criteria of docking applications that leads to mismatch in 3D-structure complexes predictions. That is why further inquiry of methods to account for these interactions is now the area of active research.

Aromatic N versus aromatic F: bioisosterism discovered in RNA base pairing interactions leads to a novel class of universal base analogs

The thermodynamics of base pairing is of fundamental importance. Fluorinated base analogs are valuable tools for investigating pairing interactions. To understand the influence of direct base-base interactions in relation to the role of water, pairing free energies between natural nucleobases and fluorinated analogs are estimated by potential of mean force calculations. Compared to pairing of AU and GC, pairing involving fluorinated analogs is unfavorable by 0.5-1.0 kcal mol(-1). Decomposing the pairing free energies into enthalpic and entropic contributions reveals fundamental differences for Watson-Crick pairs compared to pairs involving fluorinated analogs. These differences originate from direct base-base interactions and contributions of water. Pairing free energies of fluorinated base analogs with natural bases are less unfavorable by 0.5-1.0 kcal mol(-1) compared to non-fluorinated analogs. This is attributed to stabilizing C-F(...)H-N dipolar interactions and stronger N(...)H-C hydrogen bonds, demonstrating direct and indirect influences of fluorine. 7-methyl-7H-purine and its 9-deaza analog (Z) have been suggested as members of a new class of non-fluorinated base analogs. Z is found to be the least destabilizing universal base in the context of RNA known to date. This is the first experimental evidence for nitrogen-containing heterocylces as bioisosteres of aromatic rings bearing fluorine atoms.