Intermolecular CH···O/N H-bonds in the biologically important pairs of natural nucleobases: a thorough quantum-chemical study

Intermolecular hydrogen bonding delineates the stability of non-canonical adenine base pairs: a first-principles study

Non-canonical nucleobase pairs differ from canonical Watson-Crick (WC) pairs in their hydrogen bonding patterns. This study uses density functional theory with empirical dispersion correction to examine the stability and electronic properties of free adenine dimers stabilized by hydrogen bonds along the WC, Sugar (S), and Hoogsteen (H) edges. Dispersion correction is crucial for accurate interaction energy evaluation. The most stable adenine dimer is stabilized by N-H⋯N hydrogen bonds in gas and solvent phases. Binding energy decreases by ∼10.2 kcal mol-1 for dimers with both C-H⋯N and N-H⋯N bonds, increasing the donor-acceptor distance. However, with a sugar-phosphate backbone, dimers with C-H⋯N and N-H⋯N bonds have higher binding energy in an implicit solvent, emphasizing the role of C-H⋯N interactions in stability and nucleic acid folding dynamics. This study highlights noncovalent interactions, such as hydrogen bonding and π-π stacking, within adenine pairs with potential applications in biosensing and DNA-based self-assembly on nanomaterial interfaces.

IPPP-CLOPPA Analysis of the Hydrogen Bonds of the Adenine-Thymine Base Pair. Is C-H···O the Third Hydrogen Bond?

The IPPP-CLOPPA method is applied to investigate the feasibility of the C-H···O moiety as a third hydrogen bond in the adenine-thymine base pair, and the role of this intermolecular contact in its stability. For this purpose, an analysis of the interaction energy and the potential energy of the protons of the "conventional" intermolecular hydrogen bonds and the "unconventional" C-H···O contact is performed, in order to assess how much they contribute to the intermolecular stabilization energy of the base pair. On the same grounds, this study is complemented by the analysis of the molecular electric dipolar polarizability of the hydrogen bond moieties, in order to determine the information that this property can give about the electronic mechanisms that affect the stabilization of the hydrogen bonds, the influence of the π system on each one, and the cooperativity effects among them. The results obtained seem to confirm that the C-H···O moiety contributes to the stability of the adenine-thymine pair almost as much as the "conventional" hydrogen bonds do. Besides, this stabilization effect is strengthened by cooperativity between hydrogen bonds and, particularly and mostly, with the π system.

Profound importance of the conventional O-H⋯O hydrogen bond versus a considerable blue shift of the Csp2-H bond in complexes of substituted carbonyls and carboxyls

Using quantum chemical approaches, we investigated the conventional O-H⋯O and nonconventional Csp2-H⋯O hydrogen bonds between carboxylic acids and aldehydes in 21 stable complexes. The strength of complexes is determined by the conventional O-H⋯O bond together with the nonconventional Csp2-H⋯O hydrogen bond, in which the former one is 4-5 times as strong as the latter one. Proportional linear correlations of the interaction energy with both individual energies of the O-H⋯O and Csp2-H⋯O hydrogen bonds are proposed. Different impacts of electron-donating and electron-withdrawing groups in substituted formaldehyde and formic acid on characteristics of conventional and nonconventional hydrogen bonds, as well as the strength of both hydrogen bond types and complexes, are also evaluated. Following complexation, it is noteworthy that the largest blue shift of the Csp2-H stretching frequency in the Csp2-H⋯O bond up to 105.3 cm-1 in CH3CHO⋯FCOOH is due to a decisive role of the O-H⋯O hydrogen bond, which has been rarely reported in the literature. The obtained results show that the conventional O-H⋯O hydrogen bond plays a pivotal role in the significant blue shift of the Csp2-H stretching frequency in the nonconventional Csp2-H⋯O hydrogen bond. Remarkably, the considerable blue shift of the Csp2-H stretching frequency is found to be one H of C-H in formic acid substituted by the electron-withdrawing group and one H in formaldehyde substituted by the electron-donating group. In addition, the change in the Csp2-H stretching frequency following complexation is proportional to both changes of electron density in σ*(Csp2-H) and σ*(O-H) orbitals, in which a dominant role of σ*(O-H) versus σ*(Csp2-H) is observed.

Unveiling the Intermolecular Interactions between Drug 5-Fluorouracil and Watson-Crick/Hoogsteen Base Pairs: A Computational Analysis

The adsorption of 5-fluorouracil (5FU) on Watson-Crick (WC) base pairs and Hoogsteen (HT) base pairs has been studied using the dispersion-corrected density functional theory (DFT). The adsorption, binding energy, and thermochemistry for the drug 5FU on the WC and HT base pairs were determined. The most stable geometries were near planar geometry, and 5FU has a higher preference for WC than HT base pairs. The adsorption energies of 5FU on nucleobase pairs are consistently higher than pristine nucleobase pairs, indicating that nucleobase pair cleavage is less likely during the adsorption of the 5FU drug. The enthalpy change for the formation of 5FU-DNA base pairs is higher than that for the formation of 5FU-nucleobases and is enthalpy-driven. The E gap of AT base pairs is higher, suggesting that their chemical reactivity toward further reaction would be less than that of GC base pairs. The electron density difference (EDD) analysis shows a significant decrease in electron density in aromatic regions on the purine bases (adenine/guanine) compared to the pyrimidine bases. The MESP diagram of the stable 5FU-nucleobase pair complexes shows a directional interaction, with the positive regions in a molecule interacting with the negative region of other molecules. The atoms in molecule analysis show that the ρ(r) values of C=O···H-N are higher than those of N···H/N-H···O. The N···H intermolecular bonds between the base pair/drug and nucleobases are weak, closed shell interactions and are electrostatic in nature. The noncovalent interaction analysis shows that several new spikes are engendered along with an increase in their strength, which indicates that the H-bonding interactions are stronger and play a dominant role in stabilizing the complexes. Energy decomposition analysis shows that the drug-nucleobase pair complex has a marginal increase in the electrostatic contributions compared to nucleobase pair complexes.

Quantification of hydrogen bond energy based on equations using spectroscopic, structural, QTAIM-based, and NBO-based descriptors which calibrated by the molecular tailoring approach

CONTEXT

Hydrogen bonds critically influence the structure and properties of both organic molecules and biomolecules, as well as supramolecular assemblies. For this reason, the development and elaboration of methods for quantitative assessment of hydrogen bond energy is an urgent challenge. In this study, using a large series of hydroxycarbonyl aliphatic compounds with the O‒H···O = C intramolecular hydrogen bond, a bank of hydrogen bond descriptors was created, including spectroscopic, structural, QTAIM-based, and NBO-based parameters. It was shown that the O‒H vibration frequency, OH chemical shift as the spectroscopic descriptors, the O···H hydrogen bond length, O···O distance, and O‒H covalent bond length as the structural descriptors, the electron density and its Laplacian, electron potential energy density in the hydrogen bond critical point, the electron density at the ring critical point as the QTAIM-based descriptors change in a correlated manner. The same correlation is found in change of the charge transfer energy through a hydrogen bond, the occupancy of the O‒H bond antibonding orbital, the Wiberg indices of the O···H hydrogen bond, and the O‒H covalent bond, as well as the polarization of the O‒H bond, which are the NBO-based descriptors. It was also recognized that the specified descriptors from the spectroscopic, structural, QTAIM-based, and NBO-based categories are functionally related to the values of intramolecular hydrogen bond energy, quantified via the molecular tailoring approach. This allowed one to obtain a system of equations for quantitative estimation of intramolecular hydrogen bond energy based on the spectroscopic, structural, QTAIM, and NBO descriptors, which makes such quantification more dependable and reliable.

METHODS

To obtain the spectroscopic descriptors, the vibrational spectra and shielding constants were calculated using the GIAO method. Structural descriptors were obtained for the equilibrium geometry of molecules, calculated at the MP2(FC)/6-311 +  + (2d,2p) level using the Gaussian 09 program. The QTAIM-based descriptors were calculated using the AIMAll program within the framework of the quantum theory "Atoms in Molecules." The NBO-based descriptors were calculated using the NBO 3.1 program implemented into Gaussian 09. To quantify the energy of intramolecular hydrogen bonds, molecular fragmentation was used within the molecular tailoring approach.

C-H Groups as Donors in Hydrogen Bonds: A Historical Overview and Occurrence in Proteins and Nucleic Acids

Hydrogen bonds constitute a unique type of non-covalent interaction, with a critical role in biology. Until fairly recently, the canonical view held that these bonds occur between electronegative atoms, typically O and N, and that they are mostly electrostatic in nature. However, it is now understood that polarized C-H groups may also act as hydrogen bond donors in many systems, including biological macromolecules. First recognized from physical chemistry studies, C-H…X bonds were visualized with X-ray crystallography sixty years ago, although their true significance has only been recognized in the last few decades. This review traces the origins of the field and describes the occurrence and significance of the most important C-H…O bonds in proteins and nucleic acids.

Revealing the Reasons for Degeneration of Resonance-Assisted Hydrogen Bond on the Aromatic Platform: Calculations of Ortho-, Meta-, Para-Disubstituted Benzenes, and (Z)-(E)-Olefins

The energies of the O-H∙∙∙O=C intramolecular hydrogen bonds were compared quantitatively for the series of ortho-disubstituted benzenes and Z-isomers of olefins via a molecular tailoring approach. It was established that the hydrogen bond energy in the former series is significantly less than that in the latter one. The reason for lowering the hydrogen bond energy in the ortho-disubstituted benzenes compared to the Z-isomers of olefins is the decrease in the π-contribution to the total energy of the complex interaction, in which the hydrogen bond per se is enhanced by the resonance effect. By the example of the para- and meta-disubstituted benzenes, as well as E-isomers of olefins, it was explicitly shown that the aromatic ring is a much poorer conductor of the resonance effect compared to the double bond. The hydrogen bond in the ortho-disubstituted benzenes has a lower energy than a typical resonance-assisted hydrogen bond because the aromatic moiety cannot properly assist the hydrogen bond with a resonance effect. Thus, a hydrogen bond on an aromatic platform should fall into a special category, namely an aromaticity-assisted hydrogen bond, which is closer by nature to a simple hydrogen bond rather than to a resonance-assisted one.

Proline C-H Bonds as Loci for Proline Assembly via C-H/O Interactions

Proline residues within proteins lack a traditional hydrogen bond donor. However, the hydrogens of the proline ring are all sterically accessible, with polarized C-H bonds at Hα and Hδ that exhibit greater partial positive character and can be utilized as alternative sites for molecular recognition. C-H/O interactions, between proline C-H bonds and oxygen lone pairs, have been previously identified as modes of recognition within protein structures and for higher-order assembly of protein structures. In order to better understand intermolecular recognition of proline residues, a series of proline derivatives was synthesized, including 4R-hydroxyproline nitrobenzoate methyl ester, acylated on the proline nitrogen with bromoacetyl and glycolyl groups, and Boc-4S-(4-iodophenyl)hydroxyproline methyl amide. All three derivatives exhibited multiple close intermolecular C-H/O interactions in the crystallographic state, with H⋅⋅⋅O distances as close as 2.3 Å. These observed distances are well below the 2.72 Å sum of the van der Waals radii of H and O, and suggest that these interactions are particularly favorable. In order to generalize these results, we further analyzed the role of C-H/O interactions in all previously crystallized derivatives of these amino acids, and found that all 26 structures exhibited close intermolecular C-H/O interactions. Finally, we analyzed all proline residues in the Cambridge Structural Database of small-molecule crystal structures. We found that the majority of these structures exhibited intermolecular C-H/O interactions at proline C-H bonds, suggesting that C-H/O interactions are an inherent and important mode for recognition of and higher-order assembly at proline residues. Due to steric accessibility and multiple polarized C-H bonds, proline residues are uniquely positioned as sites for binding and recognition via C-H/O interactions.

RNABPDB: Molecular Modeling of RNA Structure-From Base Pair Analysis in Crystals to Structure Prediction

The stable three-dimensional structure of RNA is known to play several important biochemical roles, from post-transcriptional gene regulation to enzymatic action. These structures contain double-helical regions, which often have different types of non-canonical base pairs in addition to Watson-Crick base pairs. Hence, it is important to study their structures from experimentally obtained or even predicted ones, to understand their role, or to develop a drug against the potential targets. Molecular Modeling of RNA double helices containing non-canonical base pairs is a difficult process, particularly due to the unavailability of structural features of non-Watson-Crick base pairs. Here we show a composite web-server with an associated database that allows one to generate the structure of RNA double helix containing non-canonical base pairs using consensus parameters obtained from the database. The database classification is followed by an evaluation of the central tendency of the structural parameters as well as a quantitative estimation of interaction strengths. These parameters are used to construct three-dimensional structures of double helices composed of Watson-Crick and/or non-canonical base pairs. Our benchmark study to regenerate double-helical fragments of many experimentally derived RNA structures indicate very high accuracy. This composite server is expected to be highly useful in understanding functions of various pre-miRNA by modeling structures of the molecules and estimating binding efficiency. The database can be accessed from http://hdrnas.saha.ac.in/rnabpdb .

Guide to tuning the chalcone molecular properties based on the push-pull effect energy scale created via the molecular tailoring approach

Using the molecular tailoring approach, a total energy scale for the push-pull effect in the range from -40 to 100 kcal/mol is created for the wide series of neutral, charged and doubly charged compounds on the chalcone platform. Taking into account similar energy scale for hydrogen bonds, the strength of the push-pull effect is ranked in the seven categories, ranging from negative (anti-push-pull) to very weak and very strong push-pull effect. It is demonstrated that the molecular properties of chalcone can be tuned prior synthesis due to the created energy scale for the push-pull effect. The single bonds of the π-spacer in the chalcones are shortened, the double ones are lengthened, and the C=O bond vibrations are red shifted when the push-pull effect is enhanced along the energy scale. The HOMO and LUMO energies change systematically while the HOMO-LUMO energy gap narrows as the strength of the push-pull effect increases.

New ethionamide boosters and EthR2: structural and energetic analysis

Ethionamide (ETH) is a high-profile drug for the treatment of patients with multidrug-resistant Mycobacterium tuberculosis and, in order to produce its inhibitory effects, it needs to be bioactivated by monooxygenase EthA. This process is under the control of the transcriptional repressors EthR and EthR2, so that their inhibition results in the boosting of ethionamide activation. Herein, through crystallographic data and computer simulations, we calculated the interaction binding energies of four inhibitors with improved in vitro potency, namely BDM76060 (PDB ID: 6HS1), BDM72201 (PDB ID: 6HRX), BDM76150 (PDB ID: 6HS2) and BDM72719 (PDB ID: 6HRY), in complexes with the transcriptional repressor EthR2, using density functional theory (DFT) within the molecular fractionation with conjugated caps (MFCC) approach. It was observed that these ligands share the same binding site within a 10.0 Å radius of the EthR2 protein; however, their structural particularities have a significant impact on the global energies of systems. The BDM72201 and BDM72719 components are weakly attached to the binding site, while BDM76060 and BDM76150 components produce stronger bonds, corroborating with experimental studies demonstrating that BDM76060 and BDM76150 are more successful in producing inhibitory effects. BDM76060 and BDM76150 have many functional groups that increase the contact surface with the protein and attract a more significant number of amino acid residues, being able to produce polarities that generate stronger interactions. In the current scenario of a growing number of cases of bacterial resistance, the obtained data can be used to guide clinical trials of these inhibitors and other inhibitors that act on the alternative EthR2 pathway, focusing on improving the activity of ethionamide, its effectiveness, a reduction in the treatment time and exposure to cytotoxic effects.

Variable temperature FTIR spectra of polycrystalline purine nucleobases and estimating strengths of individual hydrogen bonds

In the first part of this work, we report the FTIR spectra of pure NH and isotopically substituted ND (10-15% D and 80-90% D) polycrystalline hypoxanthine, xanthine, adenine and guanine recorded in the 400-4000 cm^-1 range, as a function of temperature (10-300 K). We provide assignments of the stretching and out-of-plane bending amine (NH₂) and imine (NH) bands to the distinct H-bonds present in the crystal, based on the temperature sensitivity and isotopic exchange behavior. Empirical correlations between spectral and thermodynamic or structural parameters enabled us to estimate the energies and lengths of H-bonds in the studied nucleobase crystals and to correlate them with literature data. The empirical H-bonding energies are compared with H-bonding and stacking energies computed for hypoxanthine. In the second part, strategies for using the empirical correlations together with information extracted from quantum mechanical data (in particular from the Bader's quantum theory of atoms in molecules, QTAIM) for the evaluation of hydrogen bonding properties are discussed, and their advantages and drawbacks pointed out. The justification for a cooperative use of quantum-mechanical calculations with empirical spectra-energy correlations is discussed.

Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study

In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.

Strength of the [Z-I···Hal]- and [Z-Hal···I]- Halogen Bonds: Electron Density Properties and Halogen Bond Length as Estimators of Interaction Energy

Bond energy is the main characteristic of chemical bonds in general and of non-covalent interactions in particular. Simple methods of express estimates of the interaction energy, Eint, using relationships between Eint and a property which is easily accessible from experiment is of great importance for the characterization of non-covalent interactions. In this work, practically important relationships between Eint and electron density, its Laplacian, curvature, potential, kinetic, and total energy densities at the bond critical point as well as bond length were derived for the structures of the [Z-I···Hal]- and [Z-Hal···I]- types bearing halogen bonds and involving iodine as interacting atom(s) (totally 412 structures). The mean absolute deviations for the correlations found were 2.06-4.76 kcal/mol.

Carbonyl stretch of CH⋯O hydrogen-bonded methyl acetate in supercritical trifluoromethane

Infrared spectroscopy in the gas phase was used to study the formation reaction of the CH⋯O hydrogen bonding complex involving the CH group of trifluoromethane, as a hydrogen donor, and the carbonyl group of methyl acetate, as a hydrogen acceptor, under different (T, p) conditions. The hydrogen-bonded carbonyl stretch of the molecular pair was monitored in dilute mixtures of methyl acetate in trifluoromethane at near-critical temperatures, from gas- to liquid-like densities. In the gas region, it was possible to discriminate the carbonyl signal of the hydrogen-bonded complex from that of the free ester and have access to their relative concentration. The equilibrium constant of the hydrogen bonding reaction and the standard enthalpy and entropy changes in the process were determined using the spectroscopic data. CH⋯O bonding was favored by lowering temperature or pressurizing F₃CH in the mixture, remaining essentially no free carbonyl groups about the critical density. The carbonyl band of the hydrogen-bonded pair appeared as a single symmetric peak up to liquid-like densities, suggesting that the 1:1 methyl acetate-trifluoromethane complex has the most abundant stoichiometry. Spectral features as frequency shift and bandwidth of the hydrogen-bonded carbonyl were studied as a function of temperature and solvent-density. A bathochromic (red) vibrational shift was registered for the bound carbonyl band against density, with a sudden change in behavior in the near-critical region, while the width of this band remains mostly unresponsive.

A quantum chemical view of the interaction of RNA nucleobases and base pairs with the side chains of polar amino acids

Hydrogen bonding between amino acids and nucleobases is important for RNA-protein recognition. As a first step toward understanding the physicochemical features of these contacts, the present work employs density functional theory calculations to critically analyze the intrinsic structures and strength of all theoretically possible model hydrogen-bonded complexes involving RNA nucleobase edges and polar amino acid side chains. Our geometry optimizations uncover a number of unique complexes that involve variable hydrogen-bonding characteristics, including conventional donor-acceptor interactions, bifurcated interactions and single hydrogen-bonded contacts. Further, significant strength of these complexes in the gas phase (-27 kJ mol^-1 to -226 kJ mol^-1) and solvent phase (-19 kJ mol^-1 to -78 kJ mol^-1) points toward the ability of associated contacts to provide stability to RNA-protein complexes. More importantly, for the first time, our study uncovers the features of complexes involving protonated nucleobases, as well as those involving the weakly polar cysteine side chain, and thereby highlights their potential importance in biological processes that involve RNA-protein interactions. Additional analysis on select base pair-amino acid complexes uncovers the ability of amino acid side chain to simultaneously interact with both nucleobases of the base pair, and highlights the greater strength of such interactions compared to base-amino acid interactions. Overall, our analysis provides a basic physicochemical framework for understanding the molecular basis of nucleic acid-protein interactions. Further, our quantum chemical data can be used to design better algorithms for automated search of these contacts at the RNA-protein interface.Communicated by Ramaswamy H. Sarma.

Quantitative decomposition of resonance-assisted hydrogen bond energy in β-diketones into resonance and hydrogen bonding (π- and σ-) components using molecular tailoring and function-based approaches

Using the molecular tailoring and function-based approaches allows one to divide the energy of the O─H⋯O═C resonance-assisted hydrogen bond in a series of the β-diketones into resonance and hydrogen bonding components. The magnitude of the resonance component is assessed as about 6 kcal mol^-1 . This value increases by ca. 1 kcal mol^-1 on going from the weak to strong resonance-assisted hydrogen bonding. The magnitude of the hydrogen bonding component varies in the wide range from 2 to 20 kcal mol^-1 depending on the structure of the β-diketone in question.

Relationships between NMR shifts and interaction energies in biphenyls, alkanes, aza-alkanes, and oxa-alkanes with X─H… H─Y and X─H… Z (X, Y = C or N; Z = N or O) hydrogen bonding

Hydrogen-hydrogen C─H^… H─C bonding between the bay-area hydrogens in biphenyls, and more generally in congested alkanes, very strained polycyclic alkanes, and cis-2-butene, has been investigated by calculation of proton nuclear magnetic resonance (NMR) shifts and atom-atom interaction energies. Computed NMR shifts for all protons in the biphenyl derivatives correlate very well with experimental data, with zero intercept, unit slope, and a root mean square deviation of 0.06 ppm. For some congested alkanes, there is generally good agreement between computed values for a selected conformer and the experimental data, when it is available. In both cases, the shift of a given proton or pair of protons tends to increase with the corresponding interaction energy. Computed NMR shift differences for methylene protons in polycyclic alkanes, where one is involved in a very short contact ("in") and the other is not ("out"), show a rough correlation with the corresponding C─H^… H─C exchange energies. The "in" and "in,in" isomers of selected aza- and diaza-cycloalkanes, respectively, are X─H^… H─N hydrogen bonded, whereas the "out" and "in,out" isomers display X─H^… N hydrogen bonds (X = C or N). Oxa-alkanes and the "in" isomers of aza-oxa-alkanes are X─H^… O hydrogen bonded. There is a very good general correlation, including both N─H^… H─Y (Y = C or N) and N─H^… Z (Z = N or O) interactions, for NH proton shifts against the exchange energy. For "in" CH protons, the data for the different C─H^… H─Y and C─H^… Z interactions are much more dispersed and the overall shift/exchange energy correlation is less satisfactory.

Intramolecular tautomerization of the quercetin molecule due to the proton transfer: QM computational study

Quercetin molecule (3, 3', 4', 5, 7-pentahydroxyflavone, C15H10O7) is an important flavonoid compound of natural origin, consisting of two aromatic A and B rings linked through the C ring with endocyclic oxygen atom and five hydroxyl groups attached to the 3, 3', 4', 5 and 7 positions. This molecule is found in many foods and plants, and is known to have a wide range of therapeutic properties, like an anti-oxidant, anti-toxic, anti-inflammatory etc. In this study for the first time we have revealed and investigated the pathways of the tautomeric transformations for the most stable conformers of the isolated quercetin molecule (Brovarets' & Hovorun, 2019) via the intramolecular proton transfer. Energetic, structural, dynamical and polar characteristics of these transitions, in particular relative Gibbs free and electronic energies, characteristics of the intramolecular specific interactions-H-bonds and attractive van der Waals contacts, have been analysed in details. It was demonstrated that the most probable process among all investigated is the proton transfer from the O3H hydroxyl group of the C ring to the C2' carbon atom of the C2'H group of the B ring along the intramolecular O3H…C2' H-bond with the further formation of the C2'H2 group. It was established that the proton transfer from the hydroxyl groups to the carbon atoms of the neighboring CH groups is assisted at the transition states by the strong intramolecular HCH…O H-bond (~28.5 kcal∙mol-1). The least probable path of the proton transfer-from the C8H group to the endocyclic O1 oxygen atom-causes the decyclization of the C ring in some cases. It is shortly discussed the biological importance of the obtained results.

Structural properties and influence of solvent on the stability of telomeric four-stranded i-motif DNA

Repetitive cytosine rich i-motif forming sequences are abundant in the telomere, centromere and promoters of several oncogenes and in some instances are known to regulate transcription and gene expression. The in vivo existence of i-motif structures demands further insight into the factors affecting their formation and stability and development of better understanding of their gene regulatory functions. Most prior studies characterizing the conformational dynamics of i-motifs are based on i-motif forming synthetic constructs. Here, we present a systematic study on the stability and structural properties of biologically relevant i-motifs of telomeric and centromeric repeat fragments. Our results based on molecular dynamics simulations and quantum chemical calculations indicate that along with base pairing interactions within the i-motif core the overall folded conformation is associated with the stable C-HO sugar "zippers" in the narrow grooves and structured water molecules along the wide grooves. The stacked geometry of the hemi-protonated cytosine pairs within the i-motif core is mainly governed by the repulsive base stacking interaction. The loop sequence can affect the structural dynamics of the i-motif by altering the loop motion and backbone conformation. Overall this study provides microscopic insight into the i-motif structure that will be helpful to understand the structural aspect of mechanisms of gene regulation by i-motif DNA.

Novel Tautomerisation Mechanisms of the Biologically Important Conformers of the Reverse Löwdin, Hoogsteen, and Reverse Hoogsteen G*·C* DNA Base Pairs via Proton Transfer: A Quantum-Mechanical Survey

For the first time, in this study with the use of QM/QTAIM methods we have exhaustively investigated the tautomerization of the biologically-important conformers of the G*·C* DNA base pair-reverse Löwdin G*·C*(rWC), Hoogsteen G*'·C*(H), and reverse Hoogsteen G*'·C*(rH) DNA base pairs-via the single (SPT) or double (DPT) proton transfer along the neighboring intermolecular H-bonds. These tautomeric reactions finally lead to the formation of the novel G· C O 2 * (rWC), G N 2 * · C(rWC), G*'N2·C(rWC), G N 7 * · C(H), and G*'N7·C(rH) DNA base mispairs. Gibbs free energies of activation for these reactions are within the range 3.64-31.65 kcal·mol-1 in vacuum under normal conditions. All TSs are planar structures (Cs symmetry) with a single exception-the essentially non-planar transition state TSG*·C*(rWC)↔G+·C-(rWC) (C1 symmetry). Analysis of the kinetic parameters of the considered tautomerization reactions indicates that in reality only the reverse Hoogsteen G*'·C*(rH) base pair undergoes tautomerization. However, the population of its tautomerised state G*'N7·C(rH) amounts to an insignificant value-2.3·10-17. So, the G*·C*(rWC), G*'·C*(H), and G*'·C*(rH) base pairs possess a permanent tautomeric status, which does not depend on proton mobility along the neighboring H-bonds. The investigated tautomerization processes were analyzed in details by applying the author's unique methodology-sweeps of the main physical and chemical parameters along the intrinsic reaction coordinate (IRC). In general, the obtained data demonstrate the tautomeric mobility and diversity of the G*·C* DNA base pair.

Relationships between Interaction Energy and Electron Density Properties for Homo Halogen Bonds of the [(A)nY-X···X-Z(B)m] Type (X = Cl, Br, I)

Relationships between interaction energy (Eint) and electron density properties at the X···X bond critical point or the d(X···X) distance were established for the large set of structures [(A)nY-X···X-Z(B)m] bearing the halogen bonds Cl···Cl, Br···Br, and I···I (640 structures in total). The best estimator of Eint is the kinetic energy density (Gb), which reasonably approximates the whole set of the structures as -Eint = 0.128Gb2 - 0.82Gb + 1.66 (R2 = 0.91, mean absolute deviation 0.39 kcal/mol) and demonstrates low dispersion. The potential and kinetic energy densities, electron density, and the d(X···X) distance behave similarly as estimators of Eint for the individual series Cl···Cl, Br···Br, and I···I. A number of the Eint(property) correlations are recommended for the practical application in the express estimates of the strength of the homo-halogen bonds.

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.

Remarkable shifts of Csp2 -H and O-H stretching frequencies and stability of complexes of formic acid with formaldehydes and thioformaldehydes

Thirty-six stable complexes of formic acid with formaldehydes and thioformaldehydes were determined on the potential energy surface, in which the XCHO···HCOOH complexes are found to be more stable than the XCHS···HCOOH counterparts, with X = H, F, Cl, Br, CH₃ , NH₂ . All complexes are stabilized by hydrogen bonds, and their contribution to the total stabilization energy of the complexes increases in going from C-H···S to C-H···O to O-H···S and finally to O-H···O. Remarkably, a significant blueshift of C_sp2 -H bond by 81-96 cm^-1 in the C_sp2 -H···O hydrogen bond has hardly ever been reported, and a considerable redshift of O-H stretching frequency by 206-544 cm^-1 in the O-H···O/S hydrogen bonds is also predicted. The obtained results in our present work and previous literatures support that a distance contraction and a stretching frequency blueshift of C-H bond involving hydrogen bond depend mainly on its polarity and gas phase basicity of proton acceptor, besides the rearrangement of electron density due to complex formation. Markedly, we suggest the ratio of deprotonation enthalpy to proton affinity (R _c ) as an indicator to prospect for classification of hydrogen bonds. The symmetry adapted perturbation theory results show a larger role of attractive electrostatic term in XO-n as compared to that in XS-n and the electrostatic interaction is overwhelming dispersion or induction counterparts in stabilizing XO-n and XS-n, with n = 1, 2, 3. © 2019 Wiley Periodicals, Inc.

Estimating Strengths of Individual Hydrogen Bonds in RNA Base Pairs: Toward a Consensus between Different Computational Approaches

Noncoding RNA molecules are composed of a large variety of noncanonical base pairs that shape up their functionally competent folded structures. Each base pair is composed of at least two interbase hydrogen bonds (H-bonds). It is expected that the characteristic geometry and stability of different noncanonical base pairs are determined collectively by the properties of these interbase H-bonds. We have studied the ground-state electronic properties [using density functional theory (DFT) and DFT-D3-based methods] of all the 118 normal base pairs and 36 modified base pairs, belonging to 12 different geometric families (cis and trans of WW, WH, HH, WS, HS, and SS) that occur in a nonredundant set of high-resolution RNA crystal structures. Having addressed some of the limitations of the earlier approaches, we provide here a comprehensive compilation of the average energies of different types of interbase H-bonds (E _HB). We have also characterized each interbase H-bond using 13 different parameters that describe its geometry, charge distribution at its bond critical point (BCP), and n → σ*-type charge transfer from filled π orbitals of the H-bond acceptor to the empty antibonding orbital of the H-bond donor. On the basis of the extent of their linear correlation with the H-bonding energy, we have shortlisted five parameters to model linear equations for predicting E _HB values. They are (i) electron density at the BCP: ρ, (ii) its Laplacian: ∇²ρ, (iii) stabilization energy due to n → σ*-type charge transfer: E(2), (iv) donor-hydrogen distance, and (v) hydrogen-acceptor distance. We have performed single variable and multivariable linear regression analysis over the normal base pairs and have modeled sets of linear relationships between these five parameters and E _HB. Performance testing of our model over the set of modified base pairs shows promising results, at least for the moderately strong H-bonds.

Drug–DNA interaction, a joint DFT-D3/MD study on safranal as an anticancer and DNA nanostructure model

In this research, using a combination of quantum mechanics and molecular dynamic (MD) simulations, the interaction of safranal (2,6,6-trimethylcyclohexa-1,3-dien-1-carboxaldehyde) as an anti-cancer drug and Dickerson B-DNA was studied. MD simulations were executed for 35 ns in water. Binding energy analysis in three definite parts of the B-DNA and comparison between different contributions of the binding energy shows that the van der Waals energy part of the interaction is impressive among the standard molecular mechanic energy terms. On the basis of Gibbs energies, it is confirmed that the most important interactions in the safranal complex are related to the A–T and C–G rich regions, which is in agreement with the experimental data. Quantum theory of atoms in molecules and natural bond orbital analyses were applied. A diminution in the electronic chemical potential of the safranal–DNA complex in comparison with the isolated DNA, 0.026 and 0.022 au for the S1 region and 0.012 and 0.017 au for the S2 region, was obtained in the gas phase and water, respectively, which increases the complex stability. An enhancement in the electrophilicity character, during the complexation process, shows the electron charge flux between the safranal and DNA, especially in water. The strengths of the CH⋯O bonds at the center of safranal–DNA interaction were also evaluated. A mean value of 0.06 au for the electron density of the bond critical point of the H⋯O in the complex confirms the H-bond formation during the complexation.

Cooperativity effect of the ππ interaction between drug and DNA on intercalative binding induced by H-bonds: a QM/QTAIM investigation of the curcuminadenineH2O model system

In order to reveal the nature of intercalative binding of drug to DNA, the cooperativity effect of the ππ interaction was investigated in the curcuminadenineH2O model system by applying a combined QM and QTAIM computational approach. The H-bonds between the electron-donating group of curcumin and adenine induce the formation of the ππ stacking. The introduction of H2O weakens the H-bonding and ππ interactions, leading to an anti-cooperativity effect, as is confirmed by the AIM (atoms in molecules) and RDG (reduced density gradient) analysis. Thus, it can be inferred that the anti-cooperative effect is the main driving force for the intercalative binding of drug to DNA bases, which is in agreement with many experimental phenomena. Therefore, the designed DNA-targeted intercalating drugs should possess not only hydrophobic moieties, but also strong electron-donating groups bound to the DNA bases with H-bonds, which can slow the variation rates of the strengths of the H-bonding and ππ interactions between drug and DNA bases in the anti-cooperative process, leading to the intercalation formation. The enthalpy change is the major factor driving the positive thermodynamic cooperativity.

Extrapolation of phenolic compounds as multi-target agents against cancer and inflammation

Natural products acquire massive structural and chemical diversity, which cannot be coordinated by any synthetic libraries for small molecules and they are continuing to inspire novel discoveries in health sciences. We have performed the computational calculations for geometry optimization and prediction of electronic and structural properties of some plant phenolic compounds through Gaussian 09 program. Energies of molecular orbitals were computed, to mimic out the stabilities arising from charge delocalization and intramolecular interactions. This process indicated the eventual charge transfer within the molecules. The molecular docking and ADMET properties of these compounds with a novel anticancer (HER2) and anti-inflammatory (COX-2) targets revealed that two molecules were capable of inhibiting both the targets, and could be used as multi target inhibitors. Furthermore, molecular dynamics simulation studies were performed to elucidate the binding mechanism and the comparison of inhibitor's binding mode with diverse biological activities as anticancer and anti-inflammatory agents. A high-quality association was reported among quantum chemical, ADMET, docking, dynamics and MMGBSA results. Communicated By Ramaswamy H. Sarma.

Exploring the electron density localization in MoS2 nanoparticles using a localized-electron detector: Unraveling the origin of the one-dimensional metallic sites on MoS2 catalysts

The nature of the electron density localization in two MoS2 nanoclusters containing eight rows of Mo atoms, one with 100% sulphur coverage at the Mo edges (n8_100S) and the other with 50% coverage (n8_50S) was studied using a localized-electron detector function defined in the local moment representation. For n8_100S, pairs of neighboring S2 dimers cover the edges and the electron density localization function analysis shows the presence of a local triangular-shaped ring zone of highly delocalized electrons along these edges, which corresponds to a good metallic conductor zone. The optimized geometry analysis shows that the Mo-S2 bond length is much longer than that of the Mo-S bonds inside the cluster. The removal of one S atom from each sulphur dimer to create a cluster with 50% coverage produces a general compressive stress on the cluster optimized geometry, which shortens the Mo-S bond length, particularly at the edge. The electron density location function analysis shows that close to the cluster corners, a zone of highly delocalized electron zones with a characteristic semiconductor pattern and broken one-dimensional metallic ring was generated. These results suggest that the Mo-S2 bond elongation produced by the sulphur dimers is similar to a MoS2 monolayer under tensile strain and is the origin of the one-dimensional metallic sites at the Mo-edges. In general, the present findings show excellent agreement with the key features of the reported ambient pressure X-ray photoemission spectra and the corresponding simulated scanning tunneling microscopy images.

The intramolecular hydrogen bond as a unit of molecular electronics: Molecular switching controlled by overcrowded intramolecular three-centered hydrogen bond

The equilibrium geometry of the 2,5-bis-[2-(pyridin-2-yl)-vinyl]-1[Formula: see text]-pyrrole calculated at the MP2/6-311[Formula: see text]G([Formula: see text],[Formula: see text]) level of theory evidences the breaking of one of the components in the three-centered intramolecular hydrogen bond due to the steric strain. For this reason, the three-centered intramolecular hydrogen bonding turns out to be asymmetric interaction involving the major and minor components. However, the reversible switching between these components under an external impact is also possible. Two different stable states with unequal geometric and electronic structure are observed in the derivatives of the 2,5-bis-[2-(pyridin-2-yl)-vinyl]-1[Formula: see text]-pyrrole. These molecules represent novel molecular switches operating due to the pendulum-like transition between the nonequivalent two-centered components of the overcrowded three-centered intramolecular hydrogen bond. Implantation of hydrogen bond as a unit of the molecular scale device enhances potential of molecular electronics and could serve as a step towards the construction of artificial biological ensembles.