Strength of CH···π interactions in the C-terminal subdomain of villin headpiece

The relative free energies of the folded structures of the seven model peptides with PLX (X = W, Y, F, H, and A) and ALX (X = W and A) sequences to the corresponding extended structures are calculated using the density functional methods in water to evaluate the relative strengths of CH···π interactions, especially proline···aromatic interactions for the PLX motif of the C-terminal subdomain of villin headpiece. It has been found that the Pro···π contacts for the folded structures of the PLW, PLY, PLF, and PLH peptides have in common a geometric pattern having the edge of the Pro ring interacting with the face of the aromatic ring, as found for functionally important Pro residues in proteins. At the M06-2X/cc-pVTZ//SMD M06-2X/6-31+G(d) level of theory, the relative stabilities of the folded structures to the extended structures are obtained in the order PLW > ALW > PLA > PLH > PLY > ALA > PLF by the conformational Gibbs free energies in water, which is reasonably consistent with the observed results from the CD thermal analysis for wild-type and mutants of the C-terminal subdomains of villin headpieces. Although the interaction energies excluding the solvation free energies play a role in determining the relative stabilities of the PLX and ALX peptides, the solvation and entropic terms are found to be of consequence, too. In particular, it has been known that ～40% of the total interaction energy of the PLW peptide is ascribed to the CH···π interactions of the contacting side chains for Pro and Trp residues, in which the dispersion terms play a role.     

The square conformation of phenylglycine-incorporated ascidiacyclamide is stabilized by CH/π interactions between amino acid side chains

We designed a phenylglycine (Phg)-incorporated ascidiacyclamide (ASC) analogue, cyclo(-Phg-oxazoline-d-Val-thiazole-Ile-oxazoline-d-Val-thiazole- ([Phg]ASC), with the aim of stabilizing the square conformation of ASC through interactions between amino acid side chains. X-ray diffraction analysis showed that [Phg]ASC has a square structure, similar to ASC, in which the sec-butyl group of Ile and the benzene ring of Phg are in close proximity. Consistent with that finding, 1H NMR experiments revealed significant high-field shifts in the sec-butyl group of Ile, which suggests a potential for CH/π interactions between the sec-butyl group of Ile and the benzene ring of Phg. The CD spectra of [Phg]ASC were less affected by TFE titration or increasing temperature than those of ASC. In addition, [Phg]ASC showed approximately three times greater toxicity toward HL-60 cells than ASC. Thus the potently cytotoxic conformation of [Phg]ASC may be stabilized by CH/π interactions between the side chains of the Ile and Phg residues.     

The origin of the generalized anomeric effect: possibility of CH/n and CH/π hydrogen bonds

Ab initio MO calculations were carried out at the MP4/6-311++G(3df,3pd)//MP2/6-311++G(3df,3pd) level to investigate the conformational Gibbs energy of a series of methyl ethers CH₃O-CH₂-X (X = OH, OCH₃, F, Cl, Br, CN, CCH, C₆H₅, CHO). It was found that the Gibbs energy of the gauche conformers is lower in every case than that of the corresponding anti conformers. In the more stable gauche conformers, the interatomic distance between X and the hydrogen atom was shorter than the sum of the van der Waals radii. The natural bonding orbital (NBO) charges of group X were more negative in the gauche conformers than in the anti conformers. We suggest that the CH/n and CH/π hydrogen bonds play an important role in stabilizing the gauche conformation of these compounds.     

Contributions of cation-π interactions to the collagen triple helix stability

Cation–π interactions are found to be an important noncovalent force in proteins. Collagen is a right-handed triple helix composed of three left-handed PPII helices, in which (X–Y-Gly) repeats dominate in the sequence. Molecular modeling indicates that cation–π interactions could be formed between the X and Y positions in adjacent collagen strands. Here, we used a host–guest peptide system: (Pro-Hyp-Gly)₃-(Pro-Y-Gly-X-Hyp-Gly)-(Pro-Hyp-Gly)₃, where X is an aromatic residue and Y is a cationic residue, to study the cation–π interaction in the collagen triple helix. Circular dichroism (CD) measurements and T_m data analysis show that the cation–π interactions involving Arg have a larger contribution to the conformational stability than do those involving Lys, and Trp forms a weaker cation–π interaction with cationic residues than expected as a result of steric effects. The results also show that the formation of cation–π interactions between Arg and Phe depends on their relative positions in the strand. Moreover, the fluorinated and methylated Phe substitutions show that an electron-withdrawing or electron-donating substituent on the aromatic ring can modulate its π–electron density and the cation–π interaction in collagen. Our data demonstrate that the cation–π interaction could play an important role in stabilizing the collagen triple helix.     

Self-assembly of short collagen-related peptides into fibrils via cation-π interactions

Introduction of a cationic residue at the N-terminus and an aromatic residue at the C-terminus of a collagen-related peptide can generate favorable cation-π interactions between the termini of collagen triple helices. The experimental results indicate that such cation-π interactions can promote the self-assembly of collagen triple helices into a higher-order structure in a head-to-tail manner. Our current work shows that cation-π interactions can serve as an effective force in preparing collagen-related biomaterials.     

Nature of cation-π interactions and their role in structural stability of immunoglobulin proteins

Cation-π interactions are known to be important contributors to protein stability and ligand-protein interactions. In this study, we have analyzed the influence of cation-π interactions in single chain immunoglobulin proteins. We observed 87 cation-π interactions in a data set of 33 proteins. These interactions are mainly formed by long-range contacts, and there is preference of Arg over Lys in these interactions. Arg-Tyr interactions are predominant among the various pairs analyzed. Despite the scarcity of interactions involving Trp, the average energy for Trp-cation interactions is quite high. This information suggests that the cation-π interactions involving Trp might be of high relevance to the proteins. Secondary structure analysis reveals that cation-π interactions are formed preferably between residues in which at least one is in β-strand. Proteins having β-strand regions have the highest number of cation-π interaction-forming residues.     

Cation-π and π-π stacking interactions allow selective inhibition of butyrylcholinesterase by modified quinine and cinchonidine alkaloids

Scaffold varied quaternized quinine and cinchonidine alkaloid derivatives were evaluated for their selective butyrylcholinesterase (BChE) inhibitory potential. K_i values were between 0.4–260.5 μM (non-competitive inhibition) while corresponding K_ivalues to acetylcholinesterase (AChE) ranged from 7.0–400 μM exhibiting a 250-fold selectivity for BChE.Docking arrangements (GOLD, PLANT) revealed that the extended aromatic moieties and the quaternized nitrogen of the inhibitors were responsible for specific π–π stacking and π–cation interactions with the choline binding site and the peripheral anionic site of BChE’s active site.     

Silver(I) complexes of dibenzo-18-crown-6-ether having cation-π and π-π interactions

Three silver(I) complexes of dibenzo-18-crown-6-ether (DB[18]C6), [Ag(DB[18]C6)(ClO₄)](THF) (1), [Ag(DB[18]6)(CF₃SO₃)]₂(acetone)₂ (2) and [Ag(DB[18]C6)(CF₃COO)]₂(AgCF₃COO)₂ (3) have been synthesized in different solvents and characterized structurally. In each complex, silver ions prefer an octahedral coordination geometry and form close dinuclear complex with DB[18]C6 based on cation-π interaction in η²-fashion. In particular, the coordination unit involving σ bonding at an oxygen group and π-π bonding between two benzene rings is quite unique.     

Cation-π interactions induce kinking of a molecular hinge in the RNA polymerase bridge-helix domain

RNAPs (RNA polymerases) are complex molecular machines that contain a highly conserved catalytic site surrounded by conformationally flexible domains. High-throughput mutagenesis in the archaeal model system Methanocaldococcus jannaschii has demonstrated that the nanomechanical properties of one of these domains, the bridge-helix, exert a key regulatory role on the rate of the NAC (nucleotide-addition cycle). Mutations that increase the probability and/or half-life of kink formation in the BH-HC (bridge-helix C-terminal hinge) cause a substantial increase in specific activity ('superactivity'). Fully atomistic molecular dynamics simulations show that kinking of the BH-HC appears to be driven by cation-π interactions and involve amino acid side chains that are exceptionally highly conserved in all prokaryotic and eukaryotic species.     

Cation-π interactions in lipocalins: structural and functional implications

The cation-π interaction impacts protein folding, structural stability, specificity, and molecular recognition. Cation-π interactions have been overlooked in the lipocalin family. To fill this gap, these interactions were analyzed in the 113 crystal and solution structures from the lipocalin family. The cation-π interactions link previously identified structurally conserved regions and reveal new motifs, which are beyond the reach of a sequence alignment algorithm. Functional and structural significance of the interactions were tested experimentally in human tear lipocalin (TL). TL, a prominent and promiscuous lipocalin, has a key role in lipid binding at the ocular surface. Ligand binding modulation through the loop AB at the "open" end of the barrel has been erroneously attributed solely to electrostatic interactions. Data revealed that the interloop cation-π interaction in the pair Phe28-Lys108 contributes significantly to stabilize the holo-conformation of the loop AB. Numerous energetically significant and conserved cation-π interactions were uncovered in TL and throughout the lipocalin family. Cation-π interactions, such as the highly conserved Trp17-Arg118 pair in TL, were educed in low temperature experiments of mutants with Trp to Tyr substitutions.     

Non-canonical H-bonds in β-lactamases: importance of C–H···π interactions

β-Lactamase production is the common mechanism of resistance of β-lactam antibiotics. Knowledge of inter-residue interactions in protein structures increases our understanding of protein structure and stability. We have systematically analysed the contribution of C–H···π interactions to the stability of β-lactamases. Most of the interactions are long range and most of the interacting residues are evolutionarily conserved. The occurrence of C–H···π interactions in active sites and metal binding sites is very low in β-lactamases. Hence, C–H···π interactions are important determinants of stability in β-lactamases and they may not play a significant role in specificity. The results from this study provide valuable insights for understanding the stability patterns of β-lactamases and their relation to various other environmental preferences.     

Roles of electrostatic interaction and dispersion in CH···CH,CH···π, and π···π ethylene dimers

The structural, mechanical, electronic, and optical properties of orthorhombic Bi₂S₃ and Bi₂Se₃ compounds have been investigated by means of first principles calculations. The calculated lattice parameters and internal coordinates are in very good agreement with the experimental findings. The elastic constants are obtained, then the secondary results such as bulk modulus, shear modulus, Young’s modulus, Poisson’s ratio, anisotropy factor, and Debye temperature of polycrystalline aggregates are derived, and the relevant mechanical properties are also discussed. Furthermore, the band structures and optical properties such as real and imaginary parts of dielectric functions, energy-loss function, the effective number of valance electrons, and the effective optical dielectric constant have been computed. We also calculated some nonlinearities for Bi₂S₃ and Bi₂Se₃ (tensors of elasto-optical coefficients) under pressure.

Computational study on C−H…π interactions of acetylene with benzene, 1,3,5-trifluorobenzene and coronene

Meta-hybrid density functional theory calculations using M06-2X/6-31+G(d,p) and M06-2X/6-311+G(d,p) levels of theory have been performed to understand the strength of C?H^…π interactions of two possible types for benzene-acetylene, 1,3,5-trifluorobenzene-acetylene and coronene-acetylene complexes. Our study reveals that the C?H^...π interaction complex where acetylene located above to the center of benzene ring (classical T-shaped) is the lowest energy structure. This structure is twice more stable than the configuration characterized by H atom of benzene interacting with the π-cloud of acetylene. The binding energy of 2.91 kcal/mol calculated at the M06-2X/6-311+G(d,p) level for the lowest energy configuration (1A) is in very good agreement with the experimental binding energy of 2.7?±?0.2 kcal/mol for benzene-acetylene complex. Interestingly, the C?H^...π interaction of acetylene above to the center of the aromatic ring is not the lowest energy configuration for 1,3,5-trifluorobenzene-acetylene and coronene-acetylene complexes. The lowest energy configuration (2A) for the former complex possesses both C?H^...π interaction and C?H^...F hydrogen bond, while the lowest energy structure for the coronene-acetylene complex involves both π-π and C?H^...π interactions. C?H stretching vibrational frequencies and the frequency shifts are reported and analyzed for all of the configurations. We observed red-shift of the vibrational frequency for the stretching mode of the C-H bond that interacts with the π-cloud. Acetylene in the lowest-energy structures of the complexes exhibits significant red-shift of the C?H stretching frequency and change in intensity of the corresponding vibrational frequency, compared to bare acetylene. We have examined the molecular electrostatic potential on the surfaces of benzene, 1,3,5-trifluorobenzene, coronene and acetylene to explain the binding strengths of various complexes studied here.     

A survey of aspartate-phenylalanine and glutamate-phenylalanine interactions in the protein data bank: searching for anion-π pairs

Protein structures are stabilized using noncovalent interactions. In addition to the traditional noncovalent interactions, newer types of interactions are thought to be present in proteins. One such interaction, an anion-π pair, in which the positively charged edge of an aromatic ring interacts with an anion, forming a favorable anion-quadrupole interaction, has been previously proposed [Jackson, M. R., et al. (2007) J. Phys. Chem. B111, 8242-8249]. To study the role of anion-π interactions in stabilizing protein structure, we analyzed pairwise interactions between phenylalanine (Phe) and the anionic amino acids, aspartate (Asp) and glutamate (Glu). Particular emphasis was focused on identification of Phe-Asp or -Glu pairs separated by less than 7 ? in the high-resolution, nonredundant Protein Data Bank. Simplifying Phe to benzene and Asp or Glu to formate molecules facilitated in silico analysis of the pairs. Kitaura-Morokuma energy calculations were performed on roughly 19000 benzene-formate pairs and the resulting energies analyzed as a function of distance and angle. Edgewise interactions typically produced strongly stabilizing interaction energies (-2 to -7.3 kcal/mol), while interactions involving the ring face resulted in weakly stabilizing to repulsive interaction energies. The strongest, most stabilizing interactions were identified as preferentially occurring in buried residues. Anion-π pairs are found throughout protein structures, in helices as well as β strands. Numerous pairs also had nearby cation-π interactions as well as potential π-π stacking. While more than 1000 structures did not contain an anion-π pair, the 3134 remaining structures contained approximately 2.6 anion-π pairs per protein, suggesting it is a reasonably common motif that could contribute to the overall structural stability of a protein.     

Investigations on unconventional hydrogen bonds in RNA binding proteins: The role of CH⋯OC interactions

We have investigated the roles played by CH⋯OC interactions in RNA binding proteins. There was an average of 78 CH⋯OC interactions per protein and also there was an average of one significant CH⋯OC interaction for every 6 residues in the 59 RNA binding proteins studied. Main chain–Main chain (MM) CH⋯OC interactions are the predominant type of interactions in RNA binding proteins. The donor atom contribution to CH⋯OC interactions was mainly from aliphatic residues. The acceptor atom contribution for MM CH⋯OC interactions was mainly from Val, Phe, Leu, Ile, Arg and Ala. The secondary structure preference analysis of CH⋯OC interacting residues showed that, Arg, Gln, Glu and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Ile, Leu, Lys, Met, Phe, Trp and Val preferred to be in strand conformation. Most of the CH⋯OC interacting polar amino acid residues were solvent exposed while, majority of the CH⋯OC interacting non polar residues were excluded from the solvent. Long and medium-range CH⋯OC interactions are the predominant type of interactions in RNA binding proteins. More than 50% of CH⋯OC interacting residues had a higher conservation score. Significant percentage of CH⋯OC interacting residues had one or more stabilization centers. Sixty-six percent of the theoretically predicted stabilizing residues were also involved in CH⋯OC interactions and hence these residues may also contribute additional stability to RNA binding proteins.     

Preparation and biodistribution of 99mTc chelate of strophanthidin—a cardiac aglycone

Cardiac glycosides are well known to exert a specific and powerful effect on myocardial tissue, and there is a possibility that this class of compound with a ^99mTc radiolabel may behave as a superior myocardial imaging agent in comparison to ²⁰¹T1 which is at present used clinically. Because of the extreme chemical complexity of cardiac glycosides a simpler aglycone, strophanthidin was selected as the pilot compound for preliminary labelling and in vivo distribution studies.Strophanthidin was converted to 19-thiosemicarbazone which nicely accomodated ^99mTc to produce a pure radiopharmaceutical of high specific activity. The distribution pattern in animal models was studied which is in accordance with the metabolic studies performed earlier with the ligand itself.     

Discovery of novel inhibitors of LEDGF/p75-IN protein–protein interactions

Human lens epithelium-derived growth factor (LEDGF)/p75 plays an important role in the HIV life cycle by stimulating integrase (IN)-led viral DNA integration into cellular chromosomes. Mechanistic studies show the majority of IN inhibitors chelate magnesium ions in the catalytic active site, a region topologically distant from the LEDGF/p75 binding site. Compounds disrupting the formation of LEDGF/p75 and IN complexes serve as a novel mechanistic approach different from current antiretroviral therapies. We previously built pharmacophore models mimicking LEDGF/p75 residues and identified four classes of LEDGF/p75-IN inhibitors. Substructure and similarity searches yielded additional LEDGF/p75-IN inhibitors containing an acylhydrazone moiety. The most potent of the acylhydrazones inhibited LEDGF/p75-IN interaction with an IC₅₀ value of 400 nM. We explored structure–activity relationships (SAR) and identified new acylhydrazones, hydrazines, and diazenes as lead molecules for further optimization. Two lead LEDGF/p75-IN inhibitors showed antiviral activity.     

The role of OH…O and CH…O hydrogen bonds and H…H interactions in ethanol/methanol–water heterohexamers

Bioethanol is one of the world’s most extensively produced biofuels. However, it is difficult to purify due to the formation of the ethanol–water azeotrope. Knowledge of the azeotrope structure at the molecular level can help to improve existing purification methods. In order to achieve a better understanding of this azeotrope structure, the characterization of (ethanol)₅–water heterohexamers was carried out by analyzing the results of electronic structure calculations performed at the B3LYP/6-31+G(d) level. Hexamerization energies were found to range between ?36.8 and ?25.8 kcal/mol. Topological analysis of the electron density confirmed the existence of primary (OH…O) hydrogen bonds (HBs), secondary (CH…O) HBs, and H…H interactions in these clusters. Comparison with three different solvated alcohol systems featuring the same types of atom–atom interactions permitted the following order of stability to be determined: (methanol)₅–water > (methanol)₆ > (ethanol)₅–water > (ethanol)₆. These findings, together with accompanying geometric and spectroscopic analyses, show that similar cooperative effects exist among the primary HBs for structures with the same arrangement of primary HBs, regardless of the nature of the molecules involved. This result provides an indication that the molecular ratio can be considered to determine the unusual behavior of the ethanol–water system. The investigation also highlights the presence of several types of weak interaction in addition to primary HBs.

Open image in new window

Graphical Abstract Water-ethanol clusters exhibit a variety of interaction types between their atoms, such as primary OH...O (blue), secondary CH...O (green) and H...H (yellow) interactions as revealed by Quantum Chemical Topology