CH/π Interactions in DNA and Proteins. A Theoretical Study

Influence of conventional hydrogen bonds in the intercalation of phenanthroline derivatives with DNA: The important role of the sugar and phosphate backbone

The influence of hydrogen bonds in model intercalated systems between guanine-cytosine and adenine-thymine DNA base pairs (bps) was analyzed with the popular intercalator 1,10-phenanthroline (phen) and derivatives obtained by substitution with OH and NH2 groups in positions 4 and 7. Semiempirical and Density Functional Theory (DFT) methods were used both including dispersion effects: PM6-DH2, M06-2X and B3LYP-D3 along with the recently developed near linear-scaling coupled cluster method DLPNO-CCSD(T) for benchmark calculations. Our results given by QTAIM and non-covalent interaction analysis confirmed the existence of hydrogen bonds created by OH and NH2 . The trends in the energy decomposition analysis for the interaction energy, ΔEint , showed that the ΔEelstat contributions are equal or even a little bit higher than the values for ΔEdisp . Such important ΔEelstat attractive contribution comes mainly from the conventional hydrogen bonds formed by OH and NH2 functional groups with DNA not only with bps but specially with the sugar and phosphate backbone. This behavior is very different from that of phen and other classical intercalators that cannot form conventional hydrogen bonds, where the ΔEdisp is the most important attractive contribution to the ΔEint . The inclusion of explicit water molecules in molecular dynamics simulations showed, as a general trend, that the hydrogen bonds with the bps disappear during the simulations but those with the sugar and phosphate backbone remain in time, which highlights the important role of the sugar and phosphate backbone in the stabilization of these systems.

Greener Strategy for Lupanine Purification from Lupin Bean Wastewaters Using a Molecularly Imprinted Polymer

Lupanine is an alkaloid used in the pharma industry as a building block or precursor in the synthesis of sparteine and also explored for drug synthesis in the pharma industry as a chiral selector. This alkaloid is found in lupin bean processing wastewaters originated from the debittering process to make these beans edible. In this work, a computational chemistry approach was taken to design molecularly imprinted polymers (MIPs) selecting itaconic acid, a biobased building block, as a functional monomer that can provide higher affinities for lupanine. MIP-1 was prepared using lupanine as the template, itaconic acid as a functional monomer, and ethylene glycol dimethacrylate as a cross-linker by bulk polymerization. Lupanine was concentrated from lupin bean wastewater by nanofiltration, extracted with ethyl acetate, and purified using the synthesized MIP. MIP-1 was able to selectively recognize lupanine and improve the purity of lupanine from 78 to 88%, with 82% recovery of the alkaloid. These results show the potential application of this strategy to render the industrial process more sustainable.

From groove binding to intercalation: unravelling the weak interactions and other factors modulating the modes of interaction between methylated phenanthroline-based drugs and duplex DNA

Several antitumor drugs base their cytotoxicity on their capacity to intercalate between base pairs of DNA. Nevertheless, it has been established that the mechanism of intercalation of drugs in DNA starts with the prior groove binding mode of interaction of the drug with DNA. Sometimes, for some kind of flat small molecules, groove binding does not produce any cytotoxic effect and the fast transition of such flat small molecules to the cytotoxic intercalation mode is desirable. This is the case of methylated phenanthroline (phen) derivatives, where, changes in the substitution in the position and number of methyl groups determine their capability as cytotoxic compounds and, therefore, it is a way for the modulation of cytotoxic effects. In this work, we studied this modulation by means of the interaction of the [Pt(en)(phen)]²⁺ complex and several derivatives by methylation of phen in different number and position and the d(GTCGAC)₂ DNA hexamer via groove binding using PM6-DH2 and DFT-D methods. The analysis of the geometries, electronic structure and energetics of the studied systems was compared to experimental works to gain insight into the relation structure-interaction for the studied systems with cytotoxicity. The trends are explained by means of the Non-Covalent Interaction (NCI) index, the Energy Decomposition Analysis (EDA) and solvation contributions. Our results are in agreement with the experiments, in which the methylation of position 4 of phen seems to favour the interaction via groove binding thus making the transition to the intercalation cytotoxic mode difficult. Looking at the NCI results, these interactions come not only from the CH/π and CH/n interactions of the methyl group in position 4 but also from the ethylenediamine (en) ligand, whose orientation in the Pt complex was found in such a way that it produces a high number of weak interactions with DNA, especially with the sugar and phosphate backbone.

New Insights on the Interaction of Phenanthroline Based Ligands and Metal Complexes and Polyoxometalates with Duplex DNA and G-Quadruplexes

This work provides new insights from our team regarding advances in targeting canonical and non-canonical nucleic acid structures. This modality of medical treatment is used as a form of molecular medicine specifically against the growth of cancer cells. Nevertheless, because of increasing concerns about bacterial antibiotic resistance, this medical strategy is also being explored in this field. Up to three strategies for the use of DNA as target have been studied in our research lines during the last few years: (1) the intercalation of phenanthroline derivatives with duplex DNA; (2) the interaction of metal complexes containing phenanthroline with G-quadruplexes; and (3) the activity of Mo polyoxometalates and other Mo-oxo species as artificial phosphoesterases to catalyze the hydrolysis of phosphoester bonds in DNA. We demonstrate some promising computational results concerning the favorable interaction of these small molecules with DNA that could correspond to cytotoxic effects against tumoral cells and microorganisms. Therefore, our results open the door for the pharmaceutical and medical applications of the compounds we propose.

The role of CH/π interactions in the high affinity binding of streptavidin and biotin

The streptavidin-biotin complex has an extraordinarily high affinity (Ka: 10¹⁵mol^-1) and contains one of the strongest non-covalent interactions known. This strong interaction is widely used in biological tools, including for affinity tags, detection, and immobilization of proteins. Although hydrogen bond networks and hydrophobic interactions have been proposed to explain this high affinity, the reasons for it remain poorly understood. Inspired by the deceased affinity of biotin observed for point mutations of streptavidin at tryptophan residues, we hypothesized that a CH/π interaction may also contribute to the strong interaction between streptavidin and biotin. CH/π interactions were explored and analyzed at the biotin-binding site and at the interface of the subunits by the fragment molecular orbital method (FMO) and extended applications: PIEDA and FMO4. The results show that CH/π interactions are involved in the high affinity for biotin at the binding site of streptavidin. We further suggest that the involvement of CH/π interactions at the subunit interfaces and an extended CH/π network play more critical roles in determining the high affinity, rather than involvement at the binding site.

Torsional and Electronic Factors Control the C-H⋅⋅⋅O Interaction

The precise role of non‐conventional hydrogen bonds such as the C−H⋅⋅⋅O interaction in influencing the conformation of small molecules remains unresolved. Here we survey a series of β‐turn mimetics using X‐ray crystallography and NMR spectroscopy in conjunction with quantum calculation, and conclude that favourable torsional and electronic effects are important for the population of states with conformationally influential C−H⋅⋅⋅O interactions. Our results also highlight the challenge in attempting to deconvolute a myriad of interdependent noncovalent interactions in order to focus on the contribution of a single one. Within a small molecule that is designed to resemble the complexity of the environment within peptides and proteins, the interplay of different steric burdens, hydrogen‐acceptor/‐donor properties and rotational profiles illustrate why unambiguous conclusions based solely on NMR chemical shift data are extremely challenging to rationalize.

A theoretical study of methylation and CH/π interactions in DNA intercalation: methylated 1,10-phenanthroline in adenine–thymine base pairs

This work shows that quality is better that quantity to estabilize the intercalation of methylated phen.

CH-π hydrogen bonds in biological macromolecules

This is a sequel to the previous Perspective "The CH-π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates", which featured in a PCCP themed issue on "Weak Hydrogen Bonds - Strong Effects?": Phys. Chem. Chem. Phys., 2011, 13, 13873-13900. Evidence that weak hydrogen bonds play an enormously important role in chemistry and biochemistry has now accumulated to an extent that the rigid classical concept of hydrogen bonds formulated by Pauling needs to be seriously revised and extended. The concept of a more generalized hydrogen bond definition is indispensable for understanding the folding mechanisms of proteins. The CH-π hydrogen bond, a weak molecular force occurring between a soft acid CH and a soft base π-electron system, among all is one of the most important and plays a functional role in defining the conformation and stability of 3D structures as well as in many molecular recognition events. This concept is also valuable in structure-based drug design efforts. Despite their frequent occurrence in organic molecules and bio-molecules, the importance of CH-π hydrogen bonds is still largely unknown to many chemists and biochemists. Here we present a review that deals with the evidence, nature, characteristics and consequences of the CH-π hydrogen bond in biological macromolecules (proteins, nucleic acids, lipids and polysaccharides). It is hoped that the present Perspective will show the importance of CH-π hydrogen bonds and stimulate interest in the interactions of biological macromolecules, one of the most fascinating fields in bioorganic chemistry. Implication of this concept is enormous and valuable in the scientific community.

B-DNA characteristics are preserved in double stranded d(A)3·d(T)3 and d(G)3·d(C)3 mini-helixes: conclusions from DFT/M06-2X study

We report the results of the first comprehensive DFT study on the d(A)3·d(T)3 and d(G)3·d(C)3 nucleic acid duplexes. The ability of mini-helixes to preserve the conformation of B-DNA in the gas phase and under the influence of such factors as: solvent, uncompensated charge, and counter-ions was evaluated using M06-2X functional with 6-31G(d,p) basis set. The accuracy of the models was ascertained based on their ability to reproduce key structural features of natural B-DNA. Analysis of the helicity suggests that the helical conformations adopt geometrical parameters which are close to those of the B-DNA form. The torsion angles fall somewhere between the values observed for BI/BII conformational classes. The comparative analysis of parameters of isolated Watson-Crick base pairs versus B-DNA-like conformations indicates the same tendency of base-pair polarization and hydration. Specifically, effects of polarization of nucleobases in continuum type dielectric medium mimicking water are stronger than those caused by the presence of backbone. Polar environment as well as the presence of counterions stabilizes duplexes, facilitating helix formation. Substantial conformational changes of nucleotides upon duplex formation decrease the binding energy. In spite of structural and energetic changes, the placement of a mini-helix into the gas phase does not lead to significant disruption of the structure. On the contrary, the duplex preserves its helicity and the strands remain bound.

Insights into a highly conserved network of hydrogen bonds in the agonist binding site of nicotinic acetylcholine receptors: a structural and theoretical study

Structural and theoretical studies on the geometrical features of a hydrogen-bond network occurring in the binding site of nicotinic acetylcholine receptors (nAChRs) and composed of interconnected WxPD (Trp-x-Pro-Asp) and SWyz (Ser-Trp-yz) sequences from loops A and B, respectively, have been carried out. Multiple sequence alignments using as template the sequence of the apoform of Aplysia californica acetylcholine binding protein (Ac-AChBP) show the strict conservation of serine and tryptophan residues of the loop B SWyz sequence. Considering a sample of 19 high resolution AChBP structures, the strong conformational preferences of the key tryptophan residue has been pointing out, whatever the form, free or bounded, of AChBP. The geometry of the motif hydrogen-bond network has been characterized through the analyses of seven distances. The robustness of the various hydrogen-bond interactions is pointed out, the one involving the aspartate carboxylate group and the serine residue being the shortest of the network. The role of a cooperative effect involving a NH(His145)…OH (Ser142) hydrogen bond is highlighted. Density functional theory calculations on several simplified models based on the motif hydrogen-bond network allow probing the importance of the various hydrogen-bond interactions. The removal of the Ser142 hydroxyl group induces strong structural rearrangements, in agreement with the structural observations. Molecular electrostatic potential calculations on model systems highlight the importance of a cooperative effect in the whole hydrogen-bond network. More precisely, the key role of the Ser142 hydroxyl group, involved in several hydrogen bonds, is underlined.

Mapping of the interaction sites of galanthamine: a quantitative analysis through pairwise potentials and quantum chemistry

A quantitative analysis of the interaction sites of the anti-Alzheimer drug galanthamine with molecular probes (water and benzene molecules) representative of its surroundings in the binding site of acetylcholinesterase (AChE) has been realized through pairwise potentials calculations and quantum chemistry. This strategy allows a full and accurate exploration of the galanthamine potential energy surface of interaction. Significantly different results are obtained according to the distances of approaches between the various molecular fragments and the conformation of the galanthamine N-methyl substituent. The geometry of the most relevant complexes has then been fully optimized through MPWB1K/6-31 + G(d,p) calculations, final energies being recomputed at the LMP2/aug-cc-pVTZ(-f) level of theory. Unexpectedly, galanthamine is found to interact mainly from its hydrogen-bond donor groups. Among those, CH groups in the vicinity of the ammonium group are prominent. The trends obtained provide rationales to the predilection of the equatorial orientation of the galanthamine N-methyl substituent for binding to AChE. The analysis of the interaction energies pointed out the independence between the various interaction sites and the rigid character of galanthamine. The comparison between the cluster calculations and the crystallographic observations in galanthamine-AChE co-crystals allows the validation of the theoretical methodology. In particular, the positions of several water molecules appearing as strongly conserved in galanthamine-AChE co-crystals are predicted by the calculations. Moreover, the experimental position and orientation of lateral chains of functionally important aminoacid residues are in close agreement with the ones predicted theoretically. Our study provides relevant information for a rational drug design of galanthamine based AChE inhibitors.

Preferential base pairing modes of T·T mismatches

It has long been recognized that T·T mismatches can adopt two different modes of exchangeable wobble base pairs in which no preferential pairing mode has been observed. In this study, we have performed a systematic nuclear magnetic resonance (NMR) investigation to study the sequence context effect on the pairing modes of T·T mismatches. Our results reveal for the first time that preferential pairing mode does exist in T·T mismatches with specific type of flanking base pairs.

Developing a computational model that accurately reproduces the structural features of a dinucleoside monophosphate unit within B-DNA

The ability of a dinucleoside monophosphate to mimic the conformation of B-DNA was tested using a combination of different phosphate models (anionic, neutral, counterion), environments (gas, water), and density functionals (B3LYP, MPWB1K, M06-2X) with the 6-31G(d,p) basis set. Three sequences (5'-GX(Py)-3', where X(Py) = T, U or (Br)U) were considered, which vary in the (natural or modified) 3' pyrimidine nucleobase (X(Py)). These bases were selected due to their presence in natural DNA, structural similarity to T and/or applications in radical-initiated anti-tumour therapies. The accuracy of each of the 54 model, method and sequence combinations was judged based on the ability to reproduce key structural features of natural B-DNA, including the stacked base-base orientation and important backbone torsion angles. B3LYP yields distorted or tilted relative base-base orientations, while many computational challenges were encountered for MPWB1K. Despite wide use in computational studies of DNA, the neutral (protonated) phosphate model could not consistently predict a stacked arrangement for all sequences. In contrast, stacked base-base arrangements were obtained for all sequences with M06-2X in conjunction with both the anionic and (sodium) counterion phosphate models. However, comparison of calculated and experimental backbone conformations reveals the charge-neutralized counterion phosphate model best mimics B-DNA. Structures optimized with implicit solvent (water) are comparable to gas-phase structures, which suggests similar results should be obtained in an environment of intermediate polarity. We recommend the use of either gas-phase or water M06-2X optimizations with the counterion phosphate model to study the structure and/or reactivity of natural or modified DNA with a dinucleoside monophosphate.

Importance of CH/π hydrogen bonds in recognition of the core motif in proline-recognition domains: an ab initio fragment molecular orbital study

We examined CH/π hydrogen bonds in protein/ligand complexes involving at least one proline residue using the ab initio fragment molecular orbital (FMO) method and the program CHPI. FMO calculations were carried out at the Hartree-Fock (HF)/6-31G*, HF/6-31G**, second-order Møller-Plesset perturbation (MP2)/6-31G*, and MP2/6-31G** levels for three Src homology 3 (SH3) domains and five proline-recognition domains (PRDs) complexed with their corresponding ligand peptides. PRDs use a conserved set of aromatic residues to recognize proline-rich sequences of specific ligands. Many CH/π hydrogen bonds were identified in these complexes. CH/π hydrogen bonds occurred, in particular, in the central part of the proline-rich motifs. Our results suggest that CH/π hydrogen bonds are important in the recognition of SH3 and PRDs by their ligand peptides and play a vital role in the signal transduction system. Combined use of the FMO method and CHPI analysis is a valuable tool for the study of protein/protein and protein/ligand interactions and may be useful in rational drug design.

The CH/π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates

The CH/π hydrogen bond is an attractive molecular force occurring between a soft acid and a soft base. Contribution from the dispersion energy is important in typical cases where aliphatic or aromatic CH groups are involved. Coulombic energy is of minor importance as compared to the other weak hydrogen bonds. The hydrogen bond nature of this force, however, has been confirmed by AIM analyses. The dual characteristic of the CH/π hydrogen bond is the basis for ubiquitous existence of this force in various fields of chemistry. A salient feature is that the CH/π hydrogen bond works cooperatively. Another significant point is that it works in nonpolar as well as polar, protic solvents such as water. The interaction energy depends on the nature of the molecular fragments, CH as well as π-groups: the stronger the proton donating ability of the CH group, the larger the stabilizing effect. This Perspective focuses on the consequence of this molecular force in the conformation of organic compounds and supramolecular chemistry. Implication of the CH/π hydrogen bond extends to the specificity of molecular recognition or selectivity in organic reactions, polymer science, surface phenomena and interactions involving proteins. Many problems, unsettled to date, will become clearer in the light of the CH/π paradigm.

Mutual relationship between stacking and hydrogen bonding in DNA. Theoretical study of guanine-cytosine, guanine-5-methylcytosine, and their dimers

The mutual relationship between stacking and hydrogen-bonding and the possible influence of stacking in the different behavior of cytosine (C) and 5-methylcytosine (C') in DNA have been studied through complete DFT optimization of different structures of G-C and G-C' dimers (i.e., G-C/C-G and G-C'/C'-G), using four different functionals. Our results show that stacking leads to an increase of the O(6)...H-N(4) hydrogen bond length and to a simultaneous decrease of the N(2)-H...O(2) one, in such a way that both lengths approach each other and, in some cases, an inversion occurs. These results suggest that stacking can be a factor to explain the disparity between theory and experiment on the relative strength of the two lateral hydrogen bonds. Regarding the difference between cytosine and 5-methylcytosine, we have shown that methylation enhances the stacking interactions, mainly due to the increase of polarizability. Methylation also favors the existence of slid structures which can produce local distortions of DNA.

Effects of the biological backbone on DNA-protein stacking interactions

The pi-pi stacking (face-to-face) interactions between the five natural DNA or RNA nucleobases and the four aromatic amino acids were compared using three different types of dimers: (1) a truncated nucleoside (nucleobase) stacked with a truncated amino acid; (2) a truncated nucleoside (nucleobase) stacked with an extended amino acid; and (3) a nucleoside (extended nucleobase) stacked with a truncated amino acid. Systematic (MP2/6-31G*(0.25)) potential energy surface scans reveal important information about the effects of the deoxyribose sugar and protein backbone on the structure and binding energy between truncated nucleobase and amino acid models that are typically implemented in the literature. Most notably, electrostatic and steric interactions arising from the bulkiness of the biological backbones can change the preferred relative orientations of DNA and protein pi-systems. More importantly, the protein backbone can strengthen the stacking energy (by up to 10 kJ mol(-1)), while the deoxyribose moiety can strengthen or weaken the stacking interaction depending on the positioning of the amino acid relative to the sugar residue. These effects are likely due to additional interactions between the amino acid or nucleobase ring and the backbone in the extended monomer rather than significant changes in the properties of the biological pi-systems upon model extension. Since the present work reveals that all calculated DNA-protein stacking interactions are significant and approach the strength of other noncovalent interactions between biomolecules, both pi-pi and backbone-pi interactions must be considered when attempting to gain a complete picture of DNA-protein binding.

An analysis of the different behavior of DNA and RNA through the study of the mutual relationship between stacking and hydrogen bonding

The mutual relationship between stacking and hydrogen bonding and the possible influence of stacking in the different behavior of DNA and RNA base pairs have been studied through complete DFT optimization of different structures of A-U and A-T dimers (i.e., A-U/U-A and A-T/T-A), using some functionals developed by the group of Truhlar. The results obtained in this work clearly show that stacking and hydrogen bonding are deeply connected. The different behavior of DNA and RNA when replacing uracil by thymine can be interpreted through the formation of a stabilizing CH/pi interaction between the methyl group of thymine and the five-member ring of adenine.

Towards the development of universal, fast and highly accurate docking/scoring methods: a long way to go

Accelerating the drug discovery process requires predictive computational protocols capable of reducing or simplifying the synthetic and/or combinatorial challenge. Docking-based virtual screening methods have been developed and successfully applied to a number of pharmaceutical targets. In this review, we first present the current status of docking and scoring methods, with exhaustive lists of these. We next discuss reported comparative studies, outlining criteria for their interpretation. In the final section, we describe some of the remaining developments that would potentially lead to a universally applicable docking/scoring method.