Effects of hydrogen-bonding and stacking interactions with amino acids on the acidity of uracil

The effects of hydrogen-bonding interactions with amino acids on the (N1) acidity of uracil are evaluated using (B3LYP) density functional theory. Many different binding arrangements of each amino acid to three uracil binding sites are considered. The effects on the uracil acidity are found to significantly depend upon the nature of the amino acid and the binding orientation, but weakly depend on the binding site. Our results reveal that in some instances small models for the amino acids can be used, while for other amino acids larger models are required to properly describe the binding to uracil. The gas-phase acidity of uracil is found to increase by up to approximately 60 kJ mol(-1) due to discrete hydrogen-bonding interactions. Although (MP2) stacking interactions with aromatic amino acids decrease the acidity of uracil, unexpected increases in the acidity are found when any of the aromatic amino acids, or the backbone, hydrogen bond to uracil. Consideration of enzymatic and aqueous environments leads to decreases in the effects of the amino acids on the acidity of uracil. However, we find that the magnitude of the decrease varies with the nature of the molecule bound, as well as the (gas-phase) binding orientations and strengths, and therefore solvation effects should be considered on a case-by-case basis in future work. Nevertheless, the effects of amino acid interactions within enzymatic environments are as much as approximately 35 kJ mol(-1). The present study has general implications for understanding the nature of active site amino acids in enzymes, such as DNA repair enzymes, that catalyze reactions involving anionic nucleobase intermediates.     

Computational comparison of the stacking interactions between the aromatic amino acids and the natural or (cationic) methylated nucleobases

The strongest gas-phase MP2/6-31G*(0.25) stacking energies between the aromatic amino acids and the natural or methylated nucleobases were considered. The potential energy surfaces of dimers were searched as a function of the vertical separation, angle of rotation and horizontal displacement between monomers stacked according to their centers of mass. Our calculations reveal that the stacking interactions of adducts for a given nucleobase are dependent on the methylation site (by up to 20 kJ mol(-1)), where the relative magnitudes of the interactions are determined by the dipole moments of the adducts and the proton affinities of nucleobase methylation sites. Nevertheless, the differences in the (gas-phase) stacking of methylated adducts are small compared with the differences between the stacking of the corresponding natural and methylated nucleobases. Indeed, methylation increases the stacking energy by up to 40 kJ mol(-1) (or 135%). Although immersing the dimers in different solvents decreases the gas-phase stacking energies with an increase in the polarity of the environment, base methylation still has a significant effect on the nucleobase stacking ability in solvents with large dipole moments, and, perhaps more importantly, environments that mimic enzyme active sites. Our results shed light on the workings of DNA repairs enzymes that selectively remove a wide variety of alkylated nucleobases over the natural bases.     

The importance of aromatic-type interactions in serotonin synthesis: protein-ligand interactions in tryptophan hydroxylase and aromatic amino acid decarboxylase

The biosynthesis of serotonin requires aromatic substrates to be bound in the active sites of the enzymes tryptophan hydroxylase and aromatic amino acid decarboxylase. These aromatic substrates are held in place partially by dispersion and induction interactions with the enzymes' aromatic amino acid residues. Mutations that decrease substrate binding can result in a decrease in serotonin production and thus can lead to depression and related disorders. We use optimized crystal structures of these two enzymes to examine pair-wise electronic interaction energies between aromatic residues in the active sites and the aromatic ligands. We also perform in silico mutations on the aromatic residues to determine the change in interaction energies as mutations occur. Our second-order Moller-Plessett perturbation theory calculations show that drastic changes in interaction energy can occur and, in light of our previous work, we are able to use these data to offer predictions on the loss of protein function and on the possibility of disease upon mutation. We also examine local and gradient corrected density functional theory methods to evaluate their ability to predict these induction/dispersion-dominated interaction energies. We find that the hybrid B3LYP cannot model these interactions well, whereas the GGA HCTH407 offers largely qualitatively correct results, and the local functional SVWN quantitatively mimics the MP2 results rather well.     

Twist‐dependent stacking energy of base‐pair steps in B‐DNA geometry: A density functional theory approach

Stacking energy of all the 10 unique DNA base‐pair steps (bp step) are calculated using density functional theory within the ultrasoft pseudopotential plane wave method and local density approximation for the exchange‐correlation functional. We have studied the dependence of stacking energy on twist angle, an aspect found difficult to explain using classical theory. We have found that the twist angle for different bp steps at stacking energy minimum matches extremely well with the values of average twist obtained from B‐DNA crystal structure data. This indicates that the use of a proper quantum chemical method to calculate the π‐π electronic interactions may explain stacking energy without incorporating hydrophobic interaction through solvent or effect of backbone through pseudobond. From the twist angle‐dependent stacking energy profile, we have also generated the probability distributions of twist for all the bp steps and calculated the variance of the distribution. Our calculated variances show similar trend to that of the experimental data for which sufficient numbers of data are available. The TA, AT, and CG doublets show large variances among the 10 possible bp steps, indicating their maximum flexibility. This might be the case of unusual deformation observed at the TATA‐box while binding to TBP protein. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Synthesis and X-ray single crystal structure of a bivalent glycocluster

The crystal structure of a bivalent glycocluster containing aromatic amides reveals that alkylation of secondary amides alters amide configuration and thus carbohydrate presentation. This also facilitates non covalent interactions (azide-azide, carbonyl-pyranose and aromatic-pyranose) and thus carbohydrate-carbohydrate stacking.     

Analysis of the pi-pi stacking interactions between the aminoglycoside antibiotic kinase APH(3')-IIIa and its nucleotide ligands

A key contact in the active site of an aminoglycoside phosphotransferase enzyme (APH(3')-IIIa) is a pi-pi stacking interaction between Tyr42 and the adenine ring of bound nucleotides. We investigated the prevalence of similar Tyr-adenine contacts and found that many different protein systems employ Tyr residues in the recognition of the adenine ring. The geometry of these stacking interactions suggests that electrostatics play a role in the attraction between these aromatic systems. Kinetic and calorimetric experiments on wild-type and mutant forms of APH(3')-IIIa yielded further experimental evidence of the importance of electrostatics in the adenine binding region and suggested that the stacking interaction contributes approximately 2 kcal/mol of binding energy. This type of information concerning the forces that govern nucleotide binding in APH(3')-IIIa will facilitate inhibitor design strategies that target the nucleotide binding site of APH-type enzymes.     

The effect of methylation on the hydrogen-bonding and stacking interaction of nucleic acid bases

Methylated nucleosides play an important role in DNA/RNA function, and may affect republication by interrupting the base-pairing and base-stacking. In order to investigate the effect of methylation on the interaction between nucleic acid bases, this work presents the hydrogen-bonding and stacking interactions between 5-methylcytosine and guanine (G), cytosine (C) and G, 1-methyladenine and thymine (T), as well as adenine and T. Geometry optimization and potential energy surface scan have been performed for the involved complexes by MP2 calculations. The interaction energies, which were corrected for the basis-set superposition error by the full Boys–Bernardi counterpoise correction scheme, were used to evaluate the interaction intensity of these nucleic acid bases. The atoms in molecules theory and natural bond orbital analysis have been performed to study the hydrogen bonds in these complexes. The result shows that the methyl substitute contributes the stability to these complexes because it enhances either the hydrogen bonding or the staking interaction between nucleic acid bases studied.     

How well can new-generation density functional methods describe stacking interactions in biological systems?

We compare the performance of four recently developed DFT methods (MPW1B95, MPWB1K, PW6B95, and PWB6K) and two previous, generally successful DFT methods (B3LYP and B97-1) for the calculation of stacking interactions in six nucleic acid bases complexes and five amino acid pairs and for the calculation of hydrogen bonding interactions in two Watson-Crick type base pairs. We found that the four newly developed DFT methods give reasonable results for the stacking interactions in the DNA base pairs and amino acid pairs, whereas the previous DFT methods fail to describe interactions in these stacked complexes. We conclude that the new generation of DFT methods have greatly improved performance for stacking interaction as compared to previously available methods. We recommend the PWB6K method for investigating large DNA or protein systems where stacking plays an important role.     

Conformationally Constrained Peptidomimetics as Inhibitors of the Protein Arginine Methyl Transferases

Protein arginine N‐methyl transferases (PRMTs) belong to a family of enzymes that modulate the epigenetic code through modifications of histones. In the present study, peptides emerging from a phage display screening were modified in the search for PRMT inhibitors through substitution with non‐proteinogenic amino acids, N‐alkylation of the peptide backbone, and incorporation of constrained dipeptide mimics. One of the modified peptides ( 23 ) showed an increased inhibitory activity towards several PRMTs in the low μm range and the conformational preference of this peptide was investigated and compared with the original hit using circular dichroism and NMR spectroscopy. Introducing two constrained tryptophan residue mimics (l ‐Aia) spaced by a single amino acid was found to induce a unique turn structure stabilized by a hydrogen bond and aromatic π‐stacking interaction between the two l ‐Aia residues.     

Substrate specificity of fluoroacetate dehalogenase: an insight from crystallographic analysis, fluorescence spectroscopy, and theoretical computations

The high substrate specificity of fluoroacetate dehalogenase was explored by using crystallographic analysis, fluorescence spectroscopy, and theoretical computations. A crystal structure for the Asp104Ala mutant of the enzyme from Burkholderia sp. FA1 complexed with fluoroacetate was determined at 1.2 ? resolution. The orientation and conformation of bound fluoroacetate is different from those in the crystal structure of the corresponding Asp110Asn mutant of the enzyme from Rhodopseudomonas palustris CGA009 reported recently (J. Am. Chem. Soc. 2011, 133, 7461). The fluorescence of the tryptophan residues of the wild-type and Trp150Phe mutant enzymes from Burkholderia sp. FA1 incubated with fluoroacetate and chloroacetate was measured to gain information on the environment of the tryptophan residues. The environments of the tryptophan residues were found to be different between the fluoroacetate- and chloroacetate-bound enzymes; this would come from different binding modes of these two substrates in the active site. Docking simulations and QM/MM optimizations were performed to predict favorable conformations and orientations of the substrates. The F atom of the substrate is oriented toward Arg108 in the most stable enzyme-fluoroacetate complex. This is a stable but unreactive conformation, in which the small O-C-F angle is not suitable for the S(N)2 displacement of the F(-) ion. The cleavage of the C-F bond is initiated by the conformational change of the substrate to a near attack conformation (NAC) in the active site. The second lowest energy conformation is appropriate for NAC; the C-O distance and the O-C-F angle are reasonable for the S(N) 2 reaction. The activation energy is greatly reduced in this conformation because of three hydrogen bonds between the leaving F atom and surrounding amino acid residues. Chloroacetate cannot reach the reactive conformation, due to the longer C-Cl bond; this results in an increase of the activation energy despite the weaker C-Cl bond.     

Structure and electronic properties of amino acid ionic liquids

The interactions between eight amino acid based anions and four imidazolium-based cations have been investigated by density functional theory. The electronic and structural properties of the resulting amino acid ionic liquids (AAILs) have been unveiled by means of the atoms in molecules framework. The calculated interaction energy was found to increase in magnitude with decreasing alkyl chain length at imidazolium ring. Moreover, AAILs composed of an amino acid with some functional group such as aromatic ring had decreased interaction energy. Finally, several correlative relationships between glass transition temperature and interaction energy as well as density at bond critical point have been checked for 1-ethyl-3-methylimidazolium based ILs. Although the obtained correlations do not show excellent fits, a preliminary estimation of the glass transition temperature of different AAILs can be achieved by use of their electronic properties.     

Proton nuclear magnetic resonance study on the aromatic amino acid-guanine nucleotide system: effect of base methylation on the stacking interaction with tyrosine and phenylalanine

The stacking interactions of tyrosine methylester (TyrOMe)-guanosine-5'-monophosphate (GMP), TyrOMe-7-methylguanosine-5'-monophosphate (m7GMP), phenylalanine methylester (PheOMe)-GMP and PheOMe-m7GMP pairs in neutral buffer solution have been studied by proton nuclear magnetic resonance (1H-NMR). The H8 proton signal of GMP showed no noticeable temperature dependence, while the signals of other protons showed usual dependences arising from the ring stacking interaction with aromatic amino acids. The results can be interpreted in terms of the intramolecular C-H ... O hydrogen bonding and ring stacking. Complex formations in 1:1 molar ratio were deduced for all pairs from their Job plots. The association constant for each pair was obtained by analysis of the Scatchard plot. Further, the van't Hoff plot provided thermodynamic parameters of the complex structure. The analyses of these data suggested that albeit the N-quaternization of GMP strengthens the stacking interaction with aromatic amino acid, the bulky methyl group in m7GMP facilitates the dissociation from the amino acid with small environmental change. The possible conformations of GMP and m7GMP in the interaction states are discussed on the basis of the coupling constants.     

Aromatic interaction profile to understand the molecular basis of raltegravir resistance

Integrase (IN) is the enzyme of human immunodeficiency virus (HIV) which inserts the viral DNA (vDNA) into the host genome for successful viral replication leading to the infection. However, the chemical basis of HIV IN catalysis is speculative due to lack of complete co-crystal structure. Using the recently published prototype foamy virus IN crystal structure, we developed a model structure of HIV IN showing interaction of vDNA, the metal (Mg²⁺) cofactor, and raltegravir (RLT) in the active site. Molecular docking and dynamics simulations studies showed that RLT uses it core central ring with diketo motif for Mg²⁺ chelation and bridge interaction with DDE motif. The triple arene interactions mediated by RLT with neighboring molecular motifs (Y143, cytosine, and adenine) is maintained during long simulation in wild type (WT). The fluorobenzyl and oxadiazole moieties of RLT forms aromatic stacking with cytosine base (head stacking) aromatic side chain of Y143 (tail stacking), respectively, while central ring further establishes aromatic stacking with distorted adenine base of vDNA (central stacking). The novel triple stacking systems were further explored to understand the molecular basis of drug resistance by molecular simulation. The in silico mutation (N155H, Q148H, and Q148H + G140S) and simulation studies elucidated the structural mechanism of resistance to RLT. The simulation studies provided the molecular basis for interdependency observed for the primary and secondary (Q148H and G140S) mutations and also explained the mechanism of viral fitness regain. Our study reveals that triple stacking and its consequence in terms of VdW energetic profile acts as a critical point to understand the drug-resistance. Here, we demonstrate that the root mean square deviation of centroid system (aromatic stacking) can be used as a major determinant of RLT binding toward the fold resistance. This is first kind of report, which discloses a strategy to explore the molecular level of drug resistance profile using aromatic interactions.     

Crystal engineering for topochemical polymerization of muconic esters using halogen-halogen and CH/pi interactions as weak intermolecular interactions

We now report the molecular and crystal structure design of muconic ester derivatives on the basis of crystal engineering using halogen-halogen contacts and CH/pi interactions. The solid-state photoreaction pathway of the dibenzyl (Z,Z)-muconates as the 1,3-diene dicarboxylic acid monomers depends on the structure of the ester groups. The substitution of a halogen atom for the aromatic hydrogen of a benzyl group induces topochemical polymerization to produce stereoregular polymers in a crystalline form, whereas the unsubstituted benzyl derivative isomerizes to yield the corresponding E,E isomer under similar conditions. The topochemical polymerization process is directly confirmed by the fact that the single-crystal structures before and after the polymerization are very similar to each other. From the crystal structure analysis for a series of substituted benzyl (Z,Z)- and (E,E)-muconates, it has been revealed that the planar diene moieties are closely packed to form a columnar structure in the crystals. The stacking of the polymerizable monomers is characterized by a stacking distance of 4.9-5.2 A along the columns. This structure is supported by a halogen-halogen interaction between the chlorine or bromine atoms introduced at the p position of the benzyl groups in addition to an aromatic stacking due to the CH/pi interaction between the benzylic methylene hydrogens and aromatic rings. The design of a monomer packing corresponds to the type and position of the introduced halogen atom and also the polymorphs. To make a stacking distance of 5 A using both halogen-halogen and CH/pi interactions as supramolecular synthons is important for the molecular design of muconic ester derivatives appropriate for topochemical polymerization.     

Metal ligand aromatic cation-pi interactions in metalloproteins: ligands coordinated to metal interact with aromatic residues

Cation-pi interactions between aromatic residues and cationic amino groups in side chains and have been recognized as noncovalent bonding interactions relevant for molecular recognition and for stabilization and definition of the native structure of proteins. We propose a novel type of cation-pi interaction in metalloproteins; namely interaction between ligands coordinated to a metal cation--which gain positive charge from the metal--and aromatic groups in amino acid side chains. Investigation of crystal structures of metalloproteins in the Protein Data Bank (PDB) has revealed that there exist quite a number of metalloproteins in which aromatic rings of phenylalanine, tyrosine, and tryptophan are situated close to a metal center interacting with coordinated ligands. Among these ligands are amino acids such as asparagine, aspartate, glutamate, histidine, and threonine, but also water and substrates like ethanol. These interactions play a role in the stability and conformation of metalloproteins, and in some cases may also be directly involved in the mechanism of enzymatic reactions, which occur at the metal center. For the enzyme superoxide dismutase, we used quantum chemical computation to calculate that Trp163 has an interaction energy of 10.09 kcal mol(-1) with the ligands coordinated to iron.     

Impact of the Chirality and Curvature of Carbon Nanostructures on Their Interaction with Aromatics and Amino Acids

Understanding noncovalent interactions on the surfaces of carbon nanostructures (CNSs) is of fundamental importance and also has implications in nano‐ and biotechnology. The interactions of aromatic compounds such as benzene, naphthalene, and aromatic amino acids with CNSs of varying diameter, chirality, and curvature were systematically explored by using density functional theory. Planar graphene exhibits stronger binding affinity than curved carbon nanotubes (CNTs), whereas zigzag CNTs appear to show stronger binding affinity than armchair CNTs. For hydrocarbons, there exist two competing modes, namely, π–π stacking interactions and CH ??? π interactions, which bring the aromatic motifs into parallel and perpendicular dispositions with respect to the CNSs, respectively. Our results reveal that π–π stacking interactions override CH ??? π interactions in such cases. However, in the case of aromatic amino acids, π–π interactions can exist simultaneously along with a range of other interactions, including CH ??? π. The polarizability and HOMO energy of the CNSs were found to be the key factors that determine the binding energies. The HOMO–LUMO energy gaps of the CNSs were found to be undisturbed by the noncovalent functionalization of the aromatic molecules.     

Effects of the biological backbone on stacking interactions at DNA-protein interfaces: the interplay between the backbone···π and π···π components

The (gas-phase) MP2/6-31G*(0.25) π···π stacking interactions between the five natural bases and the aromatic amino acids calculated using (truncated) monomers composed of conjugated rings and/or (extended) monomers containing the biological backbone (either the protein backbone or deoxyribose sugar) were previously compared. Although preliminary energetic results indicated that the protein backbone strengthens, while the deoxyribose sugar either strengthens or weakens, the interaction calculated using truncated models, the reasons for these effects were unknown. The present work explains these observations by dissecting the interaction energy of the extended complexes into individual backbone···π and π···π components. Our calculations reveal that the total interaction energy of the extended complex can be predicted as a sum of the backbone···π and π···π components, which indicates that the biological backbone does not significantly affect the ring system through π-polarization. Instead, we find that the backbone can indirectly affect the magnitude of the π···π contribution by changing the relative ring orientations in extended dimers compared with truncated dimers. Furthermore, the strengths of the individual backbone···π contributions are determined to be significant (up to 18 kJ mol(-1)). Therefore, the origin of the energetic change upon model extension is found to result from a balance between an additional (attractive) backbone···π component and differences in the strength of the π···π interaction. In addition, to understand the effects of the biological backbone on the stacking interactions at DNA-protein interfaces in nature, we analyzed the stacking interactions found in select DNA-protein crystal structures, and verified that an additive approach can be used to examine the strength of these interactions in biological complexes. Interestingly, although the presence of attractive backbone···π contacts is qualitatively confirmed using the quantum theory of atoms in molecules (QTAIM), QTAIM electron density analysis is unable to quantitatively predict the additive relationship of these interactions. Most importantly, this work reveals that both the backbone···π and π···π components must be carefully considered to accurately determine the overall stability of DNA-protein assemblies.     

Intra- and intermolecular interactions between cyclic-AMP receptor protein and DNA: ab initio fragment molecular orbital study

The ab initio fragment molecular orbital (FMO) calculations were performed for the cAMP receptor protein (CRP) complexed with a cAMP and DNA duplex to elucidate their sequence-specific binding and the stability of the DNA duplex, as determined by analysis of their inter- and intramolecular interactions. Calculations were performed with the AMBER94 force field and at the HF and MP2 levels with several basis sets. The interfragment interaction energies (IFIEs) were analyzed for interactions of CRP-cAMP with each base pair, DNA duplex with each amino acid residue, and each base pair with each residue. In addition, base-base interactions were analyzed including hydrogen bonding and stacking of DNA. In the interaction between DNA and CRP-cAMP, there was a significant charge transfer (CT) from the DNA to CRP, and this CT interaction played an important role as well as the electrostatic interactions. It is necessary to apply a quantum mechanical approach beyond the "classical" force-field approach to describe the sequence specificity. In the DNA intramolecular interaction, the dispersion interactions dominated the stabilization of the base-pair stacking interactions. Strong, attractive 1,2-stacking interactions and weak, repulsive 1,3-stacking interactions were observed. Comparison of the intramolecular interactions of free and complexed DNA revealed that the base-pairing interactions were stronger, and the stacking interactions were weaker, in the complexed structure. Therefore, the DNA duplex stability appears to change due to both the electrostatic and the CT interactions that take place under conditions of DNA-CRP binding.     

A quantum chemistry study of binding carotenoids in the bacterial light-harvesting complexes

Carotenoids play the dual function of light harvesting and photoprotection in photosynthetic organisms. Despite their functional importance, the molecular basis for binding of carotenoids in the photosynthetic proteins is poorly understood. We have discovered that all carotenoids are surrounded either by aromatic residues or by chlorophylls in all known crystal structures of the photosynthetic pigment-protein complexes. The intermolecular pi-pi stacking interactions between carotenoids and the surrounding aromatic residues in the light-harvesting complex II (LH-II) of Rhodospirillum molischianum were analyzed by high level ab initio electronic structure calculations. Intermolecular interaction energies were calculated with the second-order M?ller-Plesset perturbation method (MP2) using the modified 6-31G*(0.25) basis set with diffuse d-polarization by Hobza and co-workers. The MP2/6-31G*(0.25) calculations yield a total stabilization energy of -15.66 kcal/mol between the carotenoid molecule and the four surrounding aromatic residues (alpha-Trp-23, beta-Phe-20, beta-Phe-24, beta-Phe-27). It is thus concluded that pi-pi stacking interactions between carotenoids and the aromatic residues play an essential role in binding carotenoids in the LH-II complex of Rhodospirillum molischianum. The physical nature of the pi-pi stacking interactions was further analyzed, and the dispersion interactions were found to be the dominant intermolecular attraction force. There is also a substantial electrostatic contribution to the overall intermolecular stabilization energy.     

X-ray diffraction and theoretical studies for the quantitative assessment of intermolecular arene-perfluoroarene stacking interactions

The arene-perfluoroarene stacking interaction was studied by experimental and theoretical methods. A series of compounds with different possibilities for formation of this recognition motif in the solid state were synthesized, and their crystal structures determined by single-crystal X-ray diffraction. The crystal packing of these compounds, as well as the packing of related compounds retrieved from crystallographic databases, were analyzed with quantitative crystal potentials: total lattice energies and the cohesive energies of closest molecular pairs in the crystals were calculated. The arene-perfluoroarene recognition motif emerges as a dominant interaction in the non-hydrogen-bonding compounds studied here, to the point that asymmetric dimers formed over the stacking motif carry over to asymmetric units made of two molecules in the crystal both for pure compounds and for molecular complexes; however, inter-ring distances and angles range from 3.70 to 4.85 A and from 5 to 21 degrees , respectively. Pixel energy partitioning reveals that whenever aromatic rings stack, the largest cohesive energy contribution comes from dispersion, which roughly amounts to 20 kJ mol(-1) per phenyl ring, while the coulombic term is minor but significant enough to make a difference between the arene-arene or perfluoroarene-perfluoroarene interactions on the one hand, and arene-perfluoroarene interactions on the other, whereby the latter are favored by about 10 kJ mol(-1) per phenyl ring. No evidence of special interaction which can be attributed to HF confrontation was recognizable.