Empirical formulation and parameterization of cation-π interactions for protein modeling

Aromatic Residues in Proteins: Re-Evaluating the Geometry and Energetics of π-π, Cation-π, and CH-π Interactions

Aromatic residues can participate in various biomolecular interactions, such as π-π, cation-π, and CH-π interactions, which are essential for protein structure and function. Here, we re-evaluate the geometry and energetics of these interactions using quantum mechanical (QM) calculations, focusing on pairwise interactions involving the aromatic amino acids Phe, Tyr, and Trp and the cationic amino acids Arg and Lys. Our findings reveal that π-π interactions, while energetically favorable, are less abundant in structured proteins than commonly assumed and are often overshadowed by previously underappreciated, yet prevalent, CH-π interactions. Cation-π interactions, particularly those involving Arg, show strong binding energies and a specific geometric preference toward stacked conformations, despite the global QM minimum, suggesting that a rather perpendicular T-shape conformation should be more favorable. Our results support a more nuanced understanding of protein stabilization via interactions involving aromatic residues. On the one hand, our results challenge the traditional emphasis on π-π interactions in structured proteins by showing that CH-π and cation-π interactions contribute significantly to their structure. On the other hand, π-π interactions appear to be key stabilizers in solvated regions and thus may be particularly important to the stabilization of intrinsically disordered proteins.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

The phenoxyl group-modulated interplay of cation-π and σ-type interactions in the alkali metal series

The interaction potential energy surfaces (IPESs) of four alkaline metal cations (Na+, K+, Rb+ and Cs+) complexed with phenol and catechol were explored by accurate ab initio calculations to investigate the interplay of different noncovalent interactions and their behavior along the alkali metal series and upon -OH substitution. Selected one-dimensional interaction energy curves revealed two different minimum energy configurations for all phenol- and catechol-metal complexes, characterized either by cation-π or σ-type interactions. For each investigated complex several two-dimensional IPES maps were also computed, exploiting the computational advantages of the MP2mod approach. The size of the alkali cation was found to play a similar role in modulating both kinds of complexes, as the interaction strength always decreases along the metal series, from Na+ to Cs+. Conversely, the number of hydroxyl substituents markedly affected cation-π complexes vs. σ-type ones. As a most relevant finding, in catechol-metal complexes the strength of cation-π interactions is around half that of the σ-type ones. It is argued that the combined effect of cation dimensions and hydroxyl substitution in catechol-Na+ complexes makes σ-type configurations remarkably more stable and easily accessible than cation-π ones. Besides shedding new light on the origin of biological phenomena connected with underwater adhesion, the quantum mechanical interaction energy database provided herein may offer a useful reference for tuning accurate force fields, suitable for molecular dynamics simulations, where environmental effects might be also taken into account.

Using the Supermolecule Approach To Predict the Nonlinear Optics Potential of a Novel Asymmetric Azine

Organic molecules with electron acceptors or withdrawal substituents terminal at π-conjugated system are promising candidates to be explored as materials with high linear and nonlinear optical properties. On the basis of these features, a novel asymmetric azine ( 7E, 8E)-2-(3-methoxy-4-hydroxy-benzylidene)-1-(4-nitrobenzylidene)hydrazineC₁₅H₁₃N₃O₄ (NMZ) was synthesized. The molecular structure of NMZ was elucidated by X-ray crystallography and the supramolecular arrangement was analyzed from Hirshfeld surface methodology. An iterative electrostatic scheme using a super molecule approach, where neighboring molecules are represented by charge points, was employed to investigate optical dipole moment (μ), the linear polarization (α) and the first (β) and second (γ) hyperpolarizabilities. The NMZ crystallized in the centrosymetric space group P2₁/n and packs via combined O-H···O, C-H···O, and N···π interactions. The macroscopic property of third order χ⁽³⁾ found for the NMZ is 298.62 times greater than values reported for chalcone derivative (2 E)-1-(3-bromophenyl)-3-[4 (methylsulfanyl)phenyl]prop-2-en-1-one. The results for NMZ indicate a good nonlinear optical effect.

Predictions of Ligand Selectivity from Absolute Binding Free Energy Calculations

Binding selectivity is a requirement for the development of a safe drug, and it is a critical property for chemical probes used in preclinical target validation. Engineering selectivity adds considerable complexity to the rational design of new drugs, as it involves the optimization of multiple binding affinities. Computationally, the prediction of binding selectivity is a challenge, and generally applicable methodologies are still not available to the computational and medicinal chemistry communities. Absolute binding free energy calculations based on alchemical pathways provide a rigorous framework for affinity predictions and could thus offer a general approach to the problem. We evaluated the performance of free energy calculations based on molecular dynamics for the prediction of selectivity by estimating the affinity profile of three bromodomain inhibitors across multiple bromodomain families, and by comparing the results to isothermal titration calorimetry data. Two case studies were considered. In the first one, the affinities of two similar ligands for seven bromodomains were calculated and returned excellent agreement with experiment (mean unsigned error of 0.81 kcal/mol and Pearson correlation of 0.75). In this test case, we also show how the preferred binding orientation of a ligand for different proteins can be estimated via free energy calculations. In the second case, the affinities of a broad-spectrum inhibitor for 22 bromodomains were calculated and returned a more modest accuracy (mean unsigned error of 1.76 kcal/mol and Pearson correlation of 0.48); however, the reparametrization of a sulfonamide moiety improved the agreement with experiment.

Exploring Strong Interactions in Proteins with Quantum Chemistry and Examples of Their Applications in Drug Design

OBJECTIVES

Three strong interactions between amino acid side chains (salt bridge, cation-π, and amide bridge) are studied that are stronger than (or comparable to) the common hydrogen bond interactions, and play important roles in protein-protein interactions.

METHODS

Quantum chemical methods MP2 and CCSD(T) are used in calculations of interaction energies and structural optimizations.

RESULTS

The energies of three types of amino acid side chain interactions in gaseous phase and in aqueous solutions are calculated using high level quantum chemical methods and basis sets. Typical examples of amino acid salt bridge, cation-π, and amide bridge interactions are analyzed, including the inhibitor design targeting neuraminidase (NA) enzyme of influenza A virus, and the ligand binding interactions in the HCV p7 ion channel. The inhibition mechanism of the M2 proton channel in the influenza A virus is analyzed based on strong amino acid interactions.

CONCLUSION

(1) The salt bridge interactions between acidic amino acids (Glu- and Asp-) and alkaline amino acids (Arg+, Lys+ and His+) are the strongest residue-residue interactions. However, this type of interaction may be weakened by solvation effects and broken by lower pH conditions. (2) The cation- interactions between protonated amino acids (Arg+, Lys+ and His+) and aromatic amino acids (Phe, Tyr, Trp and His) are 2.5 to 5-fold stronger than common hydrogen bond interactions and are less affected by the solvation environment. (3) The amide bridge interactions between the two amide-containing amino acids (Asn and Gln) are three times stronger than hydrogen bond interactions, which are less influenced by the pH of the solution. (4) Ten of the twenty natural amino acids are involved in salt bridge, or cation-, or amide bridge interactions that often play important roles in protein-protein, protein-peptide, protein-ligand, and protein-DNA interactions.

Insight into a molecular interaction force supporting peptide backbones and its implication to protein loops and folding

Although not being classified as the most fundamental protein structural elements like α-helices and β-strands, the loop segment may play considerable roles for protein stability, flexibility, and dynamic activity. Meanwhile, the protein loop is also quite elusive; i.e. its interactions with the other parts of protein as well as its own shape-maintaining forces have still remained as a puzzle or at least not quite clear yet. Here, we report a molecular force, the so-called polar hydrogen-π interaction (Hp-π), which may play an important role in supporting the backbones of protein loops. By conducting the potential energy surface scanning calculations on the quasi π-plane of peptide bond unit, we have observed the following intriguing phenomena: (1) when the polar hydrogen atom of a peptide unit is perpendicularly pointing to the π-plane of other peptide bond units, a remarkable Hp-π interaction occurs; (2) the interaction is distance and orientation dependent, acting in a broad space, and belonging to the 'point-to-plane' one. The molecular force reported here may provide useful interaction concepts and insights into better understanding the loop's unique stability and flexibility feature, as well as the driving force of the protein global folding.

Rational design of a fluorescent poly(N-aryleneindole ether sulfone) switch by cation–π interactions

A fluorescence emission on–off switch is achieved by adjusting the assembly of poly(N-aryleneindole ether sulfone) (PESIN) and pyridine hydrochloride via the cation–π interactions.

In depth analysis on the binding sites of adamantane derivatives in HCV (hepatitis C virus) p7 channel based on the NMR structure

BACKGROUND

The recently solved solution structure of HCV (hepatitis C virus) p7 ion channel provides a solid structure basis for drug design against HCV infection. In the p7 channel the ligand amantadine (or rimantadine) was determined in a hydrophobic pocket. However the pharmocophore (-NH2) of the ligand was not assigned a specific binding site.

RESULTS

The possible binding sites for amino group of adamantane derivatives is studied based on the NMR structure of p7 channel using QM calculation and molecular modeling. In the hydrophobic cavity and nearby three possible binding sites are proposed: His17, Phe20, and Trp21. The ligand binding energies at the three binding sites are studied using high level QM method CCSD(T)/6-311+G(d,p) and AutoDock calculations, and the interaction details are analyzed. The potential application of the binding sites for rational inhibitor design are discussed.

CONCLUSIONS

Some useful viewpoints are concluded as follows. (1) The amino group (-NH2) of adamantane derivatives is protonated (-NH3+), and the positively charged cation may form cation-π interactions with aromatic amino acids. (2) The aromatic amino acids (His17, Phe20, and Trp21) are the possible binding sites for the protonated amino group (-NH3+) of adamantane derivatives, and the cation-π bond energies are 3 to 5 times stronger than the energies of common hydrogen bonds. (3) The higher inhibition potent of rimantadine than amantadine probably because of its higher pKa value (pKa = 10.40) and the higher positive charge in the amino group. The potential application of p7 channel structure for inhibitor design is discussed.

Cation-π interactions of methylated ammonium ions: a quantum mechanical study

Cation-π interactions of methylated ammonium ions play a key role in a broad range of biochemical systems. These include methyl-lysine binding proteins which bind to methylated sites on histone proteins, lysine demethylase enzymes which demethylate these sites, and neurotransmitter receptor complexes which bind choline-derived ligands. Recognition in these systems is achieved through an 'aromatic cage' motif in the binding site. Here we use high-level quantum mechanical calculations to address how cation-π interactions of methylated ammonium ions are modulated by a change in methylation state and interaction geometry. We survey methyl-lysine and choline-derived complexes in the Protein Databank to validate our results against available structural data. A quantitative description of cation-π interactions of methylated ammonium systems is critical to structure-based efforts to target methyl-lysine binding proteins and demethylase enzymes in the treatment of unregulated transcriptional control, and neurotransmitter receptors in the treatment of neurological disease. It is our hope that our work will serve as a benchmark for the development of physical chemistry based force fields that can accurately model the contribution of cation-π interactions to binding and specificity in these systems.

Theoretical study on the polar hydrogen-π (Hp-π) interactions between protein side chains

BACKGROUND

In the study of biomolecular structures and interactions the polar hydrogen-π bonds (Hp-π) are an extensive molecular interaction type. In proteins 11 of 20 natural amino acids and in DNA (or RNA) all four nucleic acids are involved in this type interaction.

RESULTS

The Hp-π in proteins are studied using high level QM method CCSD/6-311 + G(d,p) + H-Bq (ghost hydrogen basis functions) in vacuum and in solutions (water, acetonitrile, and cyclohexane). Three quantum chemical methods (B3LYP, CCSD, and CCSD(T)) and three basis sets (6-311 + G(d,p), TZVP, and cc-pVTZ) are compared. The Hp-π donors include R2NH, RNH2, ROH, and C6H5OH; and the acceptors are aromatic amino acids, peptide bond unit, and small conjugate π-groups. The Hp-π interaction energies of four amino acid pairs (Ser-Phe, Lys-Phe, His-Phe, and Tyr-Phe) are quantitatively calculated.

CONCLUSIONS

Five conclusion points are abstracted from the calculation results. (1) The common DFT method B3LYP fails in describing the Hp-π interactions. On the other hand, CCSD/6-311 + G(d,p) plus ghost atom H-Bq can yield better results, very close to the state-of-the-art method CCSD(T)/cc-pVTZ. (2) The Hp-π interactions are point to π-plane interactions, possessing much more interaction conformations and broader energy range than other interaction types, such as common hydrogen bond and electrostatic interactions. (3) In proteins the Hp-π interaction energies are in the range 10 to 30 kJ/mol, comparable or even larger than common hydrogen bond interactions. (4) The bond length of Hp-π interactions are in the region from 2.30 to 3.00 Å at the perpendicular direction to the π-plane, much longer than the common hydrogen bonds (~1.9 Å). (5) Like common hydrogen bond interactions, the Hp-π interactions are less affected by solvation effects.

Effects of the aromatic substitution pattern in cation-π sandwich complexes

A computational study investigating the effects of the aromatic substitution pattern on the structure and binding energies of cation-π sandwich complexes is reported. The correlation between the binding energies (Ebind) and Hammett substituent constants is approximately the same as what is observed for cation-π half-sandwich complexes. For cation-π sandwich complexes where both aromatics contain substituents the issue of relative conformation is a possible factor in the strength of the binding; however, the work presented here shows the Ebind values are approximately the same regardless of the relative conformation of the two substituted aromatics. Finally, recent computational work has shown conflicting results on whether cation-π sandwich Ebind values (Ebind,S) are approximately equal to twice the respective half-sandwich Ebind values (Ebind,HS), or if cation-π sandwich Ebind,S values are less than double the respective half-sandwich Ebind,HS values. The work presented here shows that for cation-π sandwich complexes involving substituted aromatics the Ebind,S values are less than twice the respective half-sandwich Ebind,HS values, and this is termed nonadditive. The extent to which the cation-π sandwich complexes investigated here are nonadditive is greater for B3LYP calculated values than for MP2 calculated values and for sandwich complexes with electron-donating substituents than those with electron-withdrawing groups.

The multiple roles of histidine in protein interactions

Background

Among the 20 natural amino acids histidine is the most active and versatile member that plays the multiple roles in protein interactions, often the key residue in enzyme catalytic reactions. A theoretical and comprehensive study on the structural features and interaction properties of histidine is certainly helpful.

Results

Four interaction types of histidine are quantitatively calculated, including: (1) Cation-π interactions, in which the histidine acts as the aromatic π-motif in neutral form (His), or plays the cation role in protonated form (His⁺); (2) π-π stacking interactions between histidine and other aromatic amino acids; (3) Hydrogen-π interactions between histidine and other aromatic amino acids; (4) Coordinate interactions between histidine and metallic cations. The energies of π-π stacking interactions and hydrogen-π interactions are calculated using CCSD/6-31+G(d,p). The energies of cation-π interactions and coordinate interactions are calculated using B3LYP/6-31+G(d,p) method and adjusted by empirical method for dispersion energy.

Conclusions

The coordinate interactions between histidine and metallic cations are the strongest one acting in broad range, followed by the cation-π, hydrogen-π, and π-π stacking interactions. When the histidine is in neutral form, the cation-π interactions are attractive; when it is protonated (His⁺), the interactions turn to repulsive. The two protonation forms (and pK_a values) of histidine are reversibly switched by the attractive and repulsive cation-π interactions. In proteins the π-π stacking interaction between neutral histidine and aromatic amino acids (Phe, Tyr, Trp) are in the range from -3.0 to -4.0 kcal/mol, significantly larger than the van der Waals energies.