Aromatic-Aromatic Interactions Database, A(2)ID: an analysis of aromatic π-networks in proteins

The geometrical arrangement of the aromatic rings of phenylalanine, tyrosine, tryptophan and histidine has been analyzed at a database level using the X-ray crystal structure of proteins from PDB in order to find out the aromatic-aromatic (π-π) networks in proteins and to understand how these aromatic rings are connected with each-other in a specific π-π network. A stringent examination of the 7848 proteins indicates that close to 89% of the proteins have occurrence of at least a network of 2π or a higher π-π network. The occurrence of π-π networks in various protein superfamilies based on SCOP, CATH and EC classifiers has also been probed in the present work. In general, we find that multidomain and membrane proteins as well as lyases show a more number of these networks. Analysis of the distribution of angle between planes of two proximal aromatic rings (?) distribution indicates that at a larger cutoff distance (between centroid of two aromatic rings), above 5?, C-H?π interactions (T-shaped orientation) are more prevalent, while π-π interactions (stacked orientation) are more prevalent at a smaller cutoff distance. The connectivity patterns of π-π networks propose strong propensity of finding arrangement of aromatic residues as clusters rather than linear arrangement. We have also made a public domain database "Aromatic-Aromatic Interactions Database" (A(2)ID) comprising of all types of π-π networks and their connectivity pattern present in proteins. It can be accessed by url http://203.199.182.73/gnsmmg/databases/aidb/aidb.html.     

Platinum(II) imidazo[4,5-f]-1,10-phenanthroline chloride and thiolate complexes: Synthesis and crystal structures

The planar aromatic imidazo[4,5-f]-1,10-phenanthroline ligands have been used to prepare platinum(II) chloride and thiolate complexes. The X-ray structures of two thiolate compounds are reported, which show column-like packing in the solid state due to intermolecular aromatic π-π interactions. The compounds absorb moderately in the visible region, owing to {charge-transfer-to-diimine} electronic transition.     

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.     

Investigation of aggregation and assembly of alkali lignin using iodine as a probe

Molecular iodine has been introduced into the alkali lignin (AL) solutions to adjust the π-π aggregation, and the effect of lignin-iodine complexes on the aggregation and assembly characteristics of AL have been investigated by using fluorescence, UV-vis spectroscopy, light scattering, and viscometric techniques. Results show that AL form π-π aggregates (i.e., J-aggregates) in THF driven by the π-π interaction of the aromatic groups in AL, and the π-π aggregates undergo disaggregation in THF-I(2) media because of the formation of lignin-iodine charge-transfer complexes. By using iodine as a probe to investigate the aggregation behaviors and assembly characteristics, it is estimated that about 18 mol % aromatic groups of AL form π-π aggregates in AL molecular aggregates. When molecular iodine is introduced into the AL solutions, lignin-iodine complexes occur with charge-transfer transition from HOMO of the aromatic groups of AL to the LUMO of iodine. The formation of lignin-iodine complexes reduces the affinity of the aromatic groups approaching each other due to the electrostatic repulsion and then eliminates the π-π interaction of the aromatic groups. The disaggregation of the π-π aggregates brings a dissociation behavior of AL chains and a pronounced molecular expansion. This dissociation behavior and molecular expansion of AL in the dipping solutions induce a decrease in the adsorbed amount and an increase in the adsorption rate, when AL is transferred from the dipping solution to the self-assembled adsorbed films. Consequently, the adsorption behavior of AL can be controlled by adjusting the π-π aggregation. Above observations give insight into the occurrence of J-aggregation of the aromatic groups in the AL molecular aggregates and the disaggregation mechanism of AL aggregates induced by the lignin-iodine complexes for the first time. The understanding can provide an academic instruction in the efficient utilization of the alkali lignin from the waste liquor and also leads to further development in expanding functionalities of the aromatic compounds through manipulation of the π-π aggregation.     

Role of amino acid hydrophobicity, aromaticity, and molecular volume on IAPP(20-29) amyloid self-assembly

Aromatic amino acids strongly promote cross-β amyloid formation; whether the amyloidogenicity of aromatic residues is due to high hydrophobicity and β-sheet propensity or formation of stabilizing π-π interactions has been debated. To clarify the role of aromatic residues on amyloid formation, the islet amyloid polypeptide 20-29 fragment [IAPP(20-29)], which contains a single aromatic residue (Phe 23), was adopted as a model. The side chain of residue 23 does not self-associate in cross-β fibrils of IAPP(20-29) (Nielsen et al., Angew Chem Int Ed 2009;48:2118-2121), allowing investigation of the amyloidogenicity of aromatic amino acids in a context where direct π-π interactions do not occur. We prepared variants of IAPP(20-29) in which Tyr, Leu, Phe, pentafluorophenylalanine (F5-Phe), Trp, cyclohexylalanine (Cha), α-naphthylalanine (1-Nap), or β-naphthylalanine (2-Nap) (in order of increasing peptide hydrophobicity) were incorporated at position 23 (SNNXGAILSS-NH2), and the kinetic and thermodynamic effects of these mutations on cross-β self-assembly were assessed. The Tyr, Leu, and Trp 23 variants failed to readily self-assemble at concentrations up to 1.5 mM, while the Cha 23 mutant fibrillized with attenuated kinetics and similar thermodynamic stability relative to the wild-type Phe 23 peptide. Conversely, the F5-Phe, 1-Nap, and 2-Nap 23 variants self-assembled at enhanced rates, forming fibrils with greater thermodynamic stability than the wild-type peptide. These results indicate that the high amyloidogenicity of aromatic amino acids is a function of hydrophobicity, β-sheet propensity, and planar geometry and not the ability to form stabilizing or directing π-π bonds.     

Smoothing molecular interactions: the "kinetic buffer" effect of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) widely participate in molecular recognition and signaling processes in cells by interacting with other molecules. Compared with ordered proteins, IDPs usually possess stronger intermolecular interactions in binding. As a result, the interface structure of IDPs in complexes is distinct from that of ordered-protein complexes, and this difference may have essential effect on the response to various perturbations in a cell. In this study, we examined the perturbations of intermolecular interactions and temperature on the coupled folding and binding processes of pKID to KIX domains by performing molecular dynamics simulations. By comparing a series of virtual pKID systems with various degree of disorder, we found that the complex stability and the binding kinetics of the disordered systems were less sensitive to the perturbations than the ordered systems. The origin of the lower response sensitivity of IDPs was attributed to their higher flexibility in the complex interface, which was further supported by an analysis on protein complex structures. On the basis of our simulations and results from the literature, we speculate IDPs may not only interact with their biological partners with high specificity and low affinity but also may be resistant to the perturbations in the environment and transmit signals fast and smooth. We proposed to name it the "kinetic buffer" effect.     

Effect of water coordination on competition between π and non-π cation binding sites in aromatic amino acids: l-phenylalanine, l-tyrosine, and l-tryptophan Li+, Na+, and K+ complexes

Quantum chemistry methods have been applied to charged complexes of the alkali metals Li(+), Na(+), and K(+) with the aromatic amino acids (AAAs) phenylalanine (Phe), tyrosine (Tyr), and tryptophan (Trp). The geometries of 72 different complexes (Phe·M, Tyr·M, Trp·M, M is Li(+), Na(+), or K(+)) were completely optimized at the B3LYP/6-311+G(d,p) level of density functional theory. The solvent effect on the geometry and stability of individual complexes was studied by making use of a microsolvation model. The interaction enthalpies, entropies, and Gibbs energies of nine different complexes of the systems Phe·M, Tyr·M, and Trp·M (M is Li(+), Na(+), or K(+)) were also determined at the B3LYP density functional level of theory. The calculated Gibbs binding energies of the M(+)-AAA complexes follow the order Phe < Tyr < Trp for all three metal cations studied. Among the three AAAs studied, the indole ring of Trp is the best π donor for alkali metal cations. Our calculations demonstrated the existence of strong cation-π interactions between the alkali metals and the aromatic side chains of the three AAAs. These AAAs comprise about 8% of all known protein sequences. Thus, besides the potential for hydrogen-bond interaction, aromatic residues of Phe, Tyr, and Trp show great potential for π-donor interactions. The existence of cation-π interaction in proteins has also been demonstrated experimentally. However, more complex experimental studies of metal cation-π interaction in diverse biological systems will no doubt lead to more exact validation of these investigations.     

Very strong metal ligand aromatic cation-π interactions in transition metal complexes: intermolecular interaction in tetraphenylborate salts

Strong intermolecular interactions between ligands in cationic metal complexes and aromatic rings of tetraphenylborate anion, so-called metal ligand aromatic cation-π (MLACπ) interactions, were found by screening Cambridge Structural Database. Distances between phenyl ring and ligand are shorter in these structures than in previously reported MLACπ interactions by 0.2 Å.     

Probing the “fingers” domain binding pocket of Hepatitis C virus NS5B RdRp and D559G resistance mutation via molecular docking,molecular dynamics simulation and binding free energy calculations

The NS5B RdRp polymerase is a prominent enzyme for the replication of Hepatitis C virus (HCV). During the HCV replication, the template RNA binding takes place in the “fingers” sub-domain of NS5B. The “fingers” domain is a new emerging allosteric site for the HCV drug development. The inhibitors of the “fingers” sub-domain adopt a new antiviral mechanism called RNA intervention. The details of essential amino acid residues, binding mode of the ligand, and the active site intermolecular interactions of RNA intervention reflect that this mechanism is ambiguous in the experimental study. To elucidate these details, we performed molecular docking analysis of the fingers domain inhibitor quercetagetin (QGN) with NS5B polymerase. The detailed analysis of QGN-NS5B intermolecular interactions was carried out and found that QGN interacts with the binding pocket amino acid residues Ala97, Ala140, Ile160, Phe162, Gly283, Gly557, and Asp559; and also forms π?π stacking interaction with Phe162 and hydrogen bonding interaction with Gly283. These are found to be the essential interactions for the RNA intervention mechanism. Among the strong hydrogen bonding interactions, the QGN?Ala140 is a newly identified important hydrogen bonding interaction by the present work and this interaction was not resolved by the previously reported crystal structure. Since D559G mutation at the fingers domain was reported for reducing the inhibition percentage of QGN to sevenfold, we carried out molecular dynamics (MD) simulation for wild and D559G mutated complexes to study the stability of protein conformation and intermolecular interactions. At the end of 50?ns MD simulation, the π?π stacking interaction of Phe162 with QGN found in the wild-type complex is altered into T-shaped π stacking interaction, which reduces the inhibition strength. The origin of the D559G resistance mutation was studied using combined MD simulation, binding free energy calculations and principal component analysis. The results were compared with the wild-type complex. The mutation D559G reduces the binding affinity of the QGN molecule to the fingers domain. The free energy decomposition analysis of each residue of wild-type and mutated complexes revealed that the loss of non-polar energy contribution is the origin of the resistance.

Communicated by Ramaswamy H. Sarma     

The mighty arginine, the stable quaternary amines, the powerful aromatics, and the aggressive phosphate: their role in the noncovalent minuet

In the age of proteomics, the role of certain amino acid residues and some post-translational modifications in noncovalent complex formation are gaining in importance, as the understanding of interactions between biological molecules, is at the heart of the structure function relationship puzzle. In this work, mass spectrometry is used to highlight ammonium- or guanidinium-aromatic interactions through Cation-pi bonds and ammonium- or guanidinium-phosphate interactions through salt bridge formation. Such interactions are crucial factors in certain ligand-receptor interactions and receptor-receptor interactions. In addition, the ability of phosphorylated residues and phosphorylated lipids to form noncovalent complexes with guanidinium and quaternary ammonium (mostly through Coulombic interactions) is demonstrated, and could explain the stability of certain membrane embedded protein, or a possible role for phosphorylation in protein-protein interactions. Dougherty's work demonstrates cation-pi interactions in intra-protein interactions and folding, the present work explores inter-peptide interactions, i.e., the formation of noncovalent complexes between peptides' epitopes containing adjacent aromatic residues and ones containing adjacent Arg as a model to better understand the role of cation-pi complexes in protein-protein interaction. Complexes of peptides containing aromatic residues with quaternary amines as well as the interaction of aromatic compounds, with the guanidinium group of Arg are also investigated. Considering that an inordinate number of therapeutic compounds contain aromatic rings and quaternary amines, the above-described interactions could possibly be of great importance in better understanding their mechanism of action.     

Tuning β-sheet peptide self-assembly and hydrogelation behavior by modification of sequence hydrophobicity and aromaticity

Peptide self-assembly leading to cross-β amyloid structures is a widely studied phenomenon because of its role in amyloid pathology and the exploitation of amyloid as a functional biomaterial. The self-assembly process is governed by hydrogen bonding, hydrophobic, aromatic π-π, and electrostatic Coulombic interactions. A role for aromatic π-π interactions in peptide self-assembly leading to amyloid has been proposed, but the relative contributions of π-π versus general hydrophobic interactions in these processes are poorly understood. The Ac-(XKXK)(2)-NH(2) peptide was used to study the contributions of aromatic and hydrophobic interactions to peptide self-assembly. Position X was globally replaced by valine (Val), isoleucine (Ile), phenylalanine (Phe), pentafluorophenylalanine (F(5)-Phe), and cyclohexylalanine (Cha). At low pH, these peptides remain monomeric because of repulsion of charged lysine (Lys) residues. Increasing the solvent ionic strength to shield repulsive charge-charge interactions between protonated Lys residues facilitated cross-β fibril formation. It was generally found that as peptide hydrophobicity increased, the required ionic strength to induce self-assembly decreased. At [NaCl] ranging from 0 to 1000 mM, the Val sequence failed to assemble. Assembly of the Phe sequence commenced at 700 mM NaCl and at 300 mM NaCl for the less hydrophobic Ile variant, even though it displayed a mixture of random coil and β-sheet secondary structures over all NaCl concentrations. β-Sheet formation for F(5)-Phe and Cha sequences was observed at only 20 and 60 mM NaCl, respectively. Whereas self-assembly propensity generally correlated to peptide hydrophobicity and not aromatic character the presence of aromatic amino acids imparted unique properties to fibrils derived from these peptides. Nonaromatic peptides formed fibrils of 3-15 nm in diameter, whereas aromatic peptides formed nanotape or nanoribbon architectures of 3-7 nm widths. In addition, all peptides formed fibrillar hydrogels at sufficient peptide concentrations, but nonaromatic peptides formed weak gels, whereas aromatic peptides formed rigid gels. These findings clarify the influence of aromatic amino acids on peptide self-assembly processes and illuminate design principles for the inclusion of aromatic amino acids in amyloid-derived biomaterials.     

Binding of cationic conjugated polymers to DNA: atomistic simulations of adducts involving the Dickerson's dodecamer

We here describe the investigation at the atomistic level of the structure, stability, and dynamics of several complexes resulting from the interaction of oxidized poly(3,4-ethylenedioxythiophene) with the well-known Dickerson's dodecamer sequence. Four specific arrangements have been selected as referential structures for molecular dynamics simulations, and the resulting independent trajectories tend to converge in two distinguishable models with the strongest interactions. The first one presents a coiled DNA strand enveloping the oligomer chain, whereas in the second model, the conducting polymer chain and the disorganized DNA strand are facing side-by-side. Analysis of the intermolecular interactions indicates that the electrostatic interactions involving the negatively charged DNA phosphates and the positively charged units of the oligomer are much more frequent in the first model. In addition, aside from these electrostatic interactions, specific O · · · H and S · · · H hydrogen bonds, π-π stacking, and N-H · · · π interactions have been detected. Among all of these four specific interactions, we show that the π-π stacking is the most abundant and shows the best stability, whereas O · · · H hydrogen bonds are also frequent with long lifetimes. At the end, we have to underline that these specific interactions are predominant for the thymine and the guanine, which is in perfect agreement with previous experimental observations.     

Investigation of cation–π interactions in sugar-binding proteins

Cation–π interaction is a non-covalent binding force that plays a significant role in protein stability and drug–receptor interactions. In this work, we have investigated the structural role of cation–π interactions in sugar-binding proteins (SBPs). We observed 212 cation–π interactions in 53 proteins out of 59 SBPs in dataset. There is an average one energetically significant cation–π interaction for every 66 residues in SBPs. In addition, Arg is highly preferred to form cation–π interactions, and the average energy of Arg-Trp is high among six pairs. Long-range interactions are predominant in the analyzed cation–π interactions. Comparatively, all interaction pairs favor to accommodate in strand conformations. The analysis of solvent accessible area indicates that most of the aromatic residues are found on buried or partially buried whereas cationic residues were found mostly on the exposed regions of protein. The cation–π interactions forming residues were found that around 43% of cation–π residues had highly conserved with the conservation score ≥6. Almost cationic and π-residues equally share in the stabilization center. Sugar-binding site analysis in available complexes showed that the frequency of Trp and Arg is high, suggesting the potential role of these two residues in the interactions between proteins and sugar molecules. Our observations in this study could help to further understand the structural stability of SBPs.     

Cation-pi interactions in protein-protein interfaces

Arginine is an abundant residue in protein-protein interfaces. The importance of this residue relates to the versatility of its side chain in intermolecular interactions. Different classes of protein-protein interfaces were surveyed for cation-pi interactions. Approximately half of the protein complexes and one-third of the homodimers analyzed were found to contain at least one intermolecular cation-pi pair. Interactions between arginine and tyrosine were found to be the most abundant. The electrostatic interaction energy was calculated to be approximately 3 kcal/mol, on average. A distance-based search of guanidinium:aromatic interactions was also performed using the Macromolecular Structure Database (MSD). This search revealed that half of the guanidinium:aromatic pairs pack in a coplanar manner. Furthermore, it was found that the cationic group of the cation-pi pair is frequently involved in intermolecular hydrogen bonds. In this manner the arginine side chain can participate in multiple interactions, providing a mechanism for inter-protein specificity. Thus, the cation-pi interaction is established as an important contributor to protein-protein interfaces.     

In silico binding free energy predictability with π-π interaction energy-augmented scoring function: benzimidazole Raf inhibitors as a case study

The ability to estimate binding affinities of ligands precisely is of paramount importance in designing drugs. Docking programs are used primarily to predict the binding mode of ligands to receptors. However, current scoring functions as used in docking programs are not reliable enough to predict binding affinities of ligands without any further calculations. In the present study, we investigate the usefulness of adding π-π interaction energies between ring groups of residues and ligands to the scoring function for docking. It is found that such addition helps ranking ligand activities more correctly. LMP2 calculation is used to measure π-π interaction energies between ring groups. The result of this simple addition shows possibility of π-π interaction generalization in scoring functions.     

Co-operativity in non-covalent interactions in ternary complexes: a comprehensive electronic structure theory based investigation

The structure and stability of various ternary complexes in which an extended aromatic system such as coronene interacts with ions/atoms/molecules on opposite faces of the π-electron cloud were investigated using ab initio calculations. By characterizing the nature of the intermolecular interactions using an energy decomposition analysis, it was shown that there is an interplay between various types of interactions and that there are co-operativity effects, particularly when different types of interactions coexist in the same system.

Open image in new window

Graphical abstract Weak OH-π, π-π and van der Waals-π ternary systems are stabilized through dispersion interactions. Cation-π ternary systems are stabilized by through-space electrostatic interactions.

     

Molecular principles of the interactions of disordered proteins

Thorough knowledge of the molecular principles of protein-protein recognition is essential to our understanding of protein function at the cellular level. Whereas interactions of ordered proteins have been analyzed in great detail, complexes of intrinsically unstructured/disordered proteins (IUPs) have hardly been addressed so far. Here, we have collected a database of 39 complexes of experimentally verified IUPs, and compared their interfaces with those of 72 complexes of ordered, globular proteins. The characteristic differences found between the two types of complexes suggest that IUPs represent a distinct molecular implementation of the principles of protein-protein recognition. The interfaces do not differ in size, but those of IUPs cover a much larger part of the surface of the protein than for their ordered counterparts. Moreover, IUP interfaces are significantly more hydrophobic relative to their overall amino acid composition, but also in absolute terms. They rely more on hydrophobic-hydrophobic than on polar-polar interactions. Their amino acids in the interface realize more intermolecular contacts, which suggests a better fit with the partner due to induced folding upon binding that results in a better adaptation to the partner. The two modes of interaction also differ in that IUPs usually use only a single continuous segment for partner binding, whereas the binding sites of ordered proteins are more segmented. Probably, all these features contribute to the increased evolutionary conservation of IUP interface residues. These noted molecular differences are also manifested in the interaction energies of IUPs. Our approximation of these by low-resolution force-fields shows that IUPs gain much more stabilization energy from intermolecular contacts, than from folding, i.e. they use their binding energy for folding. Overall, our findings provide a structural rationale to the prior suggestions that many IUPs are specialized for functions realized by protein-protein interactions.     

Molecular determinants for ATP-binding in proteins: a data mining and quantum chemical analysis

Adenosine 5'-triphosphate (ATP) plays an essential role in all forms of life. Molecular recognition of ATP in proteins is a subject of great importance for understanding enzymatic mechanism and for drug design. We have carried out a large-scale data mining of the Protein Data Bank (PDB) to analyze molecular determinants for recognition of the adenine moiety of ATP by proteins. Non-bonded intermolecular interactions (hydrogen bonding, pi-pi stacking interactions, and cation-pi interactions) between adenine base and surrounding residues in its binding pockets are systematically analyzed for 68 non-redundant, high-resolution crystal structures of adenylate-binding proteins. In addition to confirming the importance of the widely known hydrogen bonding, we found out that cation-pi interactions between adenine base and positively charged residues (Lys and Arg) and pi-pi stacking interactions between adenine base and surrounding aromatic residues (Phe, Tyr, Trp) are also crucial for adenine binding in proteins. On average, there exist 2.7 hydrogen bonding interactions, 1.0 pi-pi stacking interactions, and 0.8 cation-pi interactions in each adenylate-binding protein complex. Furthermore, a high-level quantum chemical analysis was performed to analyze contributions of each of the three forms of intermolecular interactions (i.e. hydrogen bonding, pi-pi stacking interactions, and cation-pi interactions) to the overall binding force of the adenine moiety of ATP in proteins. Intermolecular interaction energies for representative configurations of intermolecular complexes were analyzed using the supermolecular approach at the MP2/6-311 + G* level, which resulted in substantial interaction strengths for all the three forms of intermolecular interactions. This work represents a timely undertaking at a historical moment when a large number of X-ray crystallographic structures of proteins with bound ATP ligands have become available, and when high-level quantum chemical analysis of intermolecular interactions of large biomolecular systems becomes computationally feasible. The establishment of the molecular basis for recognition of the adenine moiety of ATP in proteins will directly impact molecular design of ATP-binding site targeted enzyme inhibitors such as kinase inhibitors.     

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.