Origins of the Stability of Imidazole–Imidazole, Benzene–Imidazole, and Benzene–Indole Dimers: CCSD(T)/CBS and SAPT Calculations

Histidine in Proteins: pH-Dependent Interplay between π-π, Cation-π, and CH-π Interactions

Histidine (His) stands out as the most versatile natural amino acid due to its side chain's facile propensity to protonate at physiological pH, leading to a transition from aromatic to cationic characteristics and thereby enabling diverse biomolecular interactions. In this study, our objective was to quantify the energetics and geometries of pairwise interactions involving His at varying pH levels. Through quantum chemical calculations, we discovered that His exhibits robust participation in both π-π and cation-π interactions, underscoring its ability to adopt a π or cationic nature, akin to other common residues. Of particular note, we found that the affinity of protonated His for aromatic residues (via cation-π interactions) is greater than the affinity of neutral His for either cationic residues (also via cation-π interactions) or aromatic residues (via π-π interactions). Furthermore, His frequently engages in CH-π interactions, and notably, depending on its protonation state, we found that some instances of hydrogen bonding by His exhibit greater stability than is typical for interamino acid hydrogen bonds. The strength of the pH-dependent pairwise energies of His with aromatic residues is supported by the abundance of pairwise interactions with His of low and high predicted pKa values. Overall, our findings illustrate the contribution of His interactions to protein stability and its potential involvement in conformational changes despite its relatively low abundance in proteins.

Cucurbit[7]uril-mediated Histidine Dimerization: Exploring the Structure and Binding Mechanism

Cucurbit[7,8]urils are known to form inclusion complexes with hydrophobic amino acids such as Trp, Tyr, Phe, and Met, as well as peptides containing these residues at the N-terminus. Despite their widespread use in protein purification, the affinity of histidine (His) for cucurbit[7,8]urils has not been extensively explored. In this study, X-ray diffraction experiments were conducted to investigate the binding of two histidine moieties to the cucurbit[7]uril (CB7) cavity, resulting in a network of π-π and hydrogen bonds. This assembly was found to induce a His pKa shift of ΔpKa=-4. Histidine weakly bound to CB7 or CB8; however, isothermal titration calorimetry revealed micromolar equilibrium dissociation constant values for CB7 and CB8 when bound to dipeptides containing His at the C-terminus. Conversely, dipeptides with His at the N-terminus exhibited millimolar values. Additionally, the His-Gly-Gly tripeptide formed a 2 : 1 complex with CB7. These findings suggest the potential use of histidine and histidine-containing tags in conjunction with CB7 for various biological applications.

Influence of Anion Species on Liquid-Liquid Phase Separation in [EMIm+][X-]/Benzene Mixtures

When a molar excess of benzene is mixed with an ionic liquid (IL), liquid-liquid phase separation may appear with a pure liquid phase almost composed of only benzene molecules separated from a liquid clathrate phase with benzene molecules dissolved in the IL. Our previous study (J. Phys. Chem. B, 124, 7929, 2020) on long-chain IL/benzene systems has concluded that benzene molecules, as planar nonpolar ones, majorly dissolve in the IL nonpolar domains consisting of cationic alkyl side chains. Nevertheless, the above mechanism is inadequate for explaining the experimental observations that benzene can also dissolve in IL systems with very short alkyl side chains. In this study, by molecular dynamics simulation of the [EMIm+][X-]/benzene mixtures with X- being Cl-, NO3-, PF6-, or Tf2N-, we still observe liquid-liquid phase separation of the pure benzene phase from a liquid clathrate (mixed IL/benzene) phase where benzene molecules are almost equally distributed near imidazolium rings through π-stacking or near alkyl side chains. The anion species strongly influences the solubility of benzene and the ratio of the two liquid phases via the alteration of anionic charge density, which tunes the strength of the electrostatic interaction among ions and thus the probability of benzene molecules interacting with both imidazolium rings and alkyl side chains: a larger anionic charge density corresponds to a lower solubility of benzene.

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.

Unraveling the Role of π-Stacking Interactions in Ligand Binding to the Thiamine Pyrophosphate Riboswitch with High-level Quantum Chemical Calculations and Docking Study

The thiamine pyrophosphate (TPP) riboswitch has emerged as the new target for designing new ligands for antibiotic purpose. Binding of the natural ligand TPP to the TPP riboswitch causes downregulation of the genes responsible for its biosynthesis. We have reported the role of π-stacking energy contributions to ligand binding with a TPP riboswitch. In conjunction with the docking study, the higher-level quantum chemical calculations performed with the wB97XD and Def2TZVPP basis set in the aqueous phase revealed that the optimum ring size is crucial to attain the effective binding efficiency of ligands with a TPP riboswitch. The π-stacking energy contributions observed for the ligands studied are largely similar; however, the cases studied with higher π-stacking energies with larger rings have a weaker ability to displace the radiolabeled thiamine from the riboswitch. The EDA and NCI analyses suggest the role of larger dispersive interactions in stabilizing the π-stacking rings. The contribution from hydrogen-bonding interactions of the hydrogen-bond donor groups on the A ring augments the binding affinity of the ligand. This study sheds light on various factors that contribute to the design of new ligands for efficient binding with a TPP riboswitch and inhibition of gene expression.

On the propagation of the OH radical produced by Cu-amyloid beta peptide model complexes. Insight from molecular modelling

Oxidative stress and metal dyshomeostasis are considered as crucial factors in the pathogenesis of Alzheimer's disease (AD). Indeed, transition metal ions such as Cu(ii) can generate Reactive Oxygen Species (ROS) via O₂ Fenton-like reduction, catalyzed by Cu(ii) coordinated to the Amyloid beta (Aβ) peptide. Despite intensive effort, the mechanisms of ROS-induced molecular damage remain poorly understood. In the present paper, we investigate on the basis of molecular modelling computations the mechanism of OH radical propagation toward the Aβ peptide, starting from the end-product of OH radical generation by Cu(ii)·Aβ. We evaluate (i) the OH oxidative capacity, as well as the energetics of the possible Aβ oxidation target residues, by quantum chemistry Density Functional Theory (DFT) on coordination models of Cu(ii)/OH/Aβ and (ii) the motion of the OH˙ approaching the Aβ target residues by classical Molecular Dynamics (MD) on the full peptide Cu(ii)/OH/Aβ(1-16). The results show that the oxidative capacity of OH coordinated Cu(ii)Aβ is significantly lower than that of the free OH radical and that propagation toward Aβ Asp and His residues is favoured over Tyr residues. These results are discussed on the basis of the recent literature on in vitro Aβ metal-catalyzed oxidation and on the possible implications for the AD oxidative stress mechanism.

Probing the Nature of Noncovalent Interactions in Dimers of Linear Tyrosine-Based Dipeptides

Tyrosine-based dipeptides self-assemble to form higher order structures. To gain insights into the nature of intermolecular interactions contributing to the early stages of the self-assembly of aromatic dipeptides, we study the dimers of linear dityrosine (YY) and tryptophan-tyrosine (WY) using quantum-chemical methods with dispersion corrections and universal solvation model based on density in combination with energy decomposition and natural bond orbital (NBO) analyses. We find that hydrogen bonding is a dominant stabilizing force. The lowest energy structure for the linear YY dimer is characterized by Ocarboxyl···H(O)tyr. In contrast, the lowest energy dimer of linear WY is stabilized by Ocarboxyl···H(N)trp and πtyr···πtyr. The solvent plays a critical role as it can change the strength and nature of interactions. The lowest energy for linear WY dimer in acetone is stabilized by Ocarboxyl···H(O)tyr, πtrp···H(C), and πtrp···H(N). The ΔG of dimerization and stabilization energies of solvated dipeptides reveal that the dipeptide systems are more stable in the solvent phase than in gas phase. NBO confirms increased magnitudes for donor-acceptor interaction for the solvated dipeptides.

Local structure and hydrogen bond characteristics of imidazole molecules for proton conduction in acid and base proton-conducting composite materials

Protons in composite materials of acidic polymers and imidazole molecules transport with rotational motion of imidazole in hydrogen bonds.

Structure, Spectroscopy, and Bonding within the Zn(q+)-Imidazole(n) (q = 0, 1, 2; n = 1-4) Clusters and Implications for Zeolitic Imidazolate Frameworks and Zn-Enzymes

Using density functional theory (DFT) with dispersion correction and ab initio post Hartree-Fock methods, we treat the bonding, the structure, the stability, and the spectroscopy of the complexes between Zn(q+) and imidazole (Im), Zn(q+)Imn (where q = 0, 1 and 2; n = 1-4). These entities are subunits of zeolitic imidazolate frameworks (ZIFs) and Zn-enzymes, which possess relevant roles in industrial and biological domains, respectively. We also investigate the Imn (n = 2-4) clusters for comparison. For each species, we determine several new structures that were not found previously. Our calculations show a competition between atomic metal solvation, by either σ-type interactions or π-stacking type interaction, and proton transfer through hydrogen bonding (H-bonding) in charged species. This results in several geometrical environments around the metal. These are connected with structural properties and the functional role of Zn cation within ZIFs and Zn-enzymes. Moreover, we show that the Zn(2+)Imn subunits do not absorb in the visible domain, which may be related to the photostability of ZIFs. Our findings are important for the development of new applications of ZIFs and metalloenzymes.

Hydrophobic noncovalent interactions of inosine-phenylalanine: a theoretical model for investigating the molecular recognition of nucleobases

Understanding the molecular recognition process of nucleobases is one of the greatest challenges for both computational chemistry and biophysics fields. In fact, our results point out that it is a hard task to take into account the hydrophobic interactions, such as π-π and T-stacking interactions, by theoretical calculations using conventional force fields due to quantum effects of hyperconjugation and electronic correlation. In this line, our findings put in evidence that simple modifications in the Lennard-Jones potential can improve theoretical predictions in scenarios where hydrophobic interactions can drive the molecular recognition.

CH-π "T-shape" interaction with histidine explains binding of aromatic galactosides to Pseudomonas aeruginosa lectin LecA

The galactose specific lectin LecA mediates biofilm formation in the opportunistic pathogen P. aeruginosa . The interaction between LecA and aromatic β-galactoside biofilm inhibitors involves an intermolecular CH-π T-shape interaction between C(ε1)-H of residue His50 in LecA and the aromatic ring of the galactoside aglycone. The generality of this interaction was tested in a diverse family of β-galactosides. LecA binding to aromatic β-galactosides (KD ∼ 8 μM) was consistently stronger than to aliphatic β-galactosides (KD ∼ 36 μM). The CH-π interaction was observed in the X-ray crystal structures of six different LecA complexes, with shorter than the van der Waals distances indicating productive binding. Related XH/cation/π-π interactions involving other residues were identified in complexes of aromatic glycosides with a variety of carbohydrate binding proteins such as concanavalin A. Exploiting such interactions might be generally useful in drug design against these targets.

Selecting DFT methods for use in optimizations of enzyme active sites: applications to ONIOM treatments of DNA glycosylases

When using a hybrid methodology to treat an enzymatic reaction, many factors contribute to selecting the method for the high-level region, which can be complicated by the presence of dispersion-driven interactions such as π–π stacking. In addition, the proper treatment of the reaction center often requires a large number of heavy atoms to be included in the high-level region, precluding the use of ab initio methods such as MP2 as well as large basis sets, in the optimization step. In the present work, popular DFT methods were tested to identify an appropriate functional for treating the high-level region in ONIOM optimizations of reactions catalyzed by nonmetalloenzymes. Eight different DFT methods (B3LYP, B97-2, MPW1K, MPWB1K, BB1K, B1B95, M06-2X, and ωB97X-D) in combination with four double-ζ quality Pople basis sets were tested for their ability to optimize noncovalent interactions (hydrogen bonding and π–π) and characterize reactions (proton transfer, SN2 hydrolysis, and unimolecular cleavage). Although the primary focus of this study is accurate structure determination, energetics were also examined at both the optimization level of theory, and with triple-ζ quality basis set and select (M06-2X or ωB97X-D) methods. If dispersion-driven interactions exist within the active site, then MPWB1K/6-31G(d,p) or M06-2X/6-31+G(d,p) are recommended for the optimization step with subsequent triple-ζ quality single-point energies. However, since dispersion-corrected functionals (M06-2X and ωB97X-D) generally require diffuse functions to yield appropriate geometries, the possible size of the high-level region is greatly limited with these methods. In contrast, if the model is large enough to recover steric constraints on π–π interactions, then B3LYP with a small basis set performs comparatively well for the optimization step and is significantly less computationally expensive. Interestingly, the functionals that afford the best geometries often do not yield the best energetics, which emphasizes the importance of structural benchmark studies.

Modulating the Strength of Hydrogen Bonds through Beryllium Bonds

The mutual influence between beryllium bonds and inter- or intramolecular hydrogen bonds (HBs) has been investigated at the B3LYP/6-311++G(3df,2p) level of theory, using the complexes between imidazole dimer and malonaldehyde with BeH2 and BeF2 as suitable model systems. Imidazole and its dimer form very strong beryllium bonds with both BeH2 and BeF2, accompanied by a significant geometry distortion of the Lewis acid. More importantly, we have found a clear cooperativity between these two noncovalent interactions, since the intermolecular HB between the two imidazole molecules in the dimer-BeX2 complex becomes much stronger than in the isolated dimer, whereas the beryllium bond becomes also stronger in the dimer-BeX2 complex, with respect to that found in the imidazole-BeX2 complex. The effects of beryllium bonds are also dramatic on the strength of intramolecular HBs. Depending on to which center the BeX2 is attached, the intramolecular HB becomes much stronger or much weaker. The first situation is found when the beryllium derivative is attached to the HB donor, whereas the second occurs if it is attached to the HB acceptor. The first effect can be so strong as to produce a spontaneous proton transfer, as it is actually the case of the malonaldehyde-BeF2 complex.