Low barrier hydrogen bonds in protein structure and function

Dynamic interchange between two protonation states is characteristic of active sites of cholinesterases

Cholinesterases are well-known and widely studied enzymes crucial to human health and involved in neurology, Alzheimer's, and lipid metabolism. The protonation pattern of active sites of cholinesterases influences all the chemical processes within, including reaction, covalent inhibition by nerve agents, and reactivation. Despite its significance, our comprehension of the fine structure of cholinesterases remains limited. In this study, we employed enhanced-sampling quantum-mechanical/molecular-mechanical calculations to show that cholinesterases predominantly operate as dynamic mixtures of two protonation states. The proton transfer between two non-catalytic glutamate residues follows the Grotthuss mechanism facilitated by a mediator water molecule. We show that this uncovered complexity of active sites presents a challenge for classical molecular dynamics simulations and calls for special treatment. The calculated proton transfer barrier of 1.65 kcal/mol initiates a discussion on the potential existence of two coupled low-barrier hydrogen bonds in the inhibited form of butyrylcholinesterase. These findings expand our understanding of structural features expressed by highly evolved enzymes and guide future advances in cholinesterase-related protein and drug design studies.

Insight into the autoproteolysis mechanism of the RsgI9 anti-σ factor from Clostridium thermocellum

Clostridium thermocellum is a potential microbial platform to convert abundant plant biomass to biofuels and other renewable chemicals. It efficiently degrades lignocellulosic biomass using a surface displayed cellulosome, a megadalton sized multienzyme containing complex. The enzymatic composition and architecture of the cellulosome is controlled by several transmembrane biomass-sensing RsgI-type anti-σ factors. Recent studies suggest that these factors transduce signals from the cell surface via a conserved RsgI extracellular (CRE) domain (also called a periplasmic domain) that undergoes autoproteolysis through an incompletely understood mechanism. Here we report the structure of the autoproteolyzed CRE domain from the C. thermocellum RsgI9 anti-σ factor, revealing that the cleaved fragments forming this domain associate to form a stable α/β/α sandwich fold. Based on AlphaFold2 modeling, molecular dynamics simulations, and tandem mass spectrometry, we propose that a conserved Asn-Pro bond in RsgI9 autoproteolyzes via a succinimide intermediate whose formation is promoted by a conserved hydrogen bond network holding the scissile peptide bond in a strained conformation. As other RsgI anti-σ factors share sequence homology to RsgI9, they likely autoproteolyze through a similar mechanism.

Revisiting the reaction pathways for phospholipid hydrolysis catalyzed by phospholipase A2 with QM/MM methods

Secreted phospholipase A2 (sPLA2) is a Ca2+-dependent, widely distributed enzyme superfamily in almost all mammalian tissues and bacteria. It is also a critical component of the venom of nearly all snakes, as well as many invertebrate species. In non-venomous contexts, sPLA2 assumes significance in cellular signaling pathways by binding cell membranes and catalyzing ester bond hydrolysis at the sn-2 position of phospholipids. Elevated levels of GIIA sPLA2 have been detected in the synovial fluid of arthritis patients, where it exhibits a pro-inflammatory function. Consequently, identifying sPLA2 inhibitors holds promise for creating highly effective pharmaceutical treatments. Beyond arthritis, the similarities among GIIA sPLA2s offer an opportunity for developing treatments against snakebite envenoming, the deadliest neglected tropical disease. Despite decades of study, the details of PLA2 membrane-binding, substrate-binding, and reaction mechanism remain elusive, demanding a comprehensive understanding of the sPLA2 catalytic mechanism. This study explores two reaction mechanism hypotheses, involving one or two water molecules, and distinct roles for the Ca2+ cofactor. Our research focuses on the human synovial sPLA2 enzyme bound to lipid bilayers of varying phospholipid compositions, and employing adiabatic QM/MM and QM/MM MD umbrella sampling methods to energetically and geometrically characterize the structures found along both reaction pathways. Our studies demonstrate the higher frequency of productive conformations within the single-water pathway, also revealing a lower free energy barrier for hydrolyzing POPC. Furthermore, we observe that the TS of this concerted one-step reaction closely resembles transition state geometries observed in X-ray crystallography complexes featuring high-affinity transition state analogue inhibitors.

Identification of natural product-based effective inhibitors of spleen tyrosine kinase (SYK) through virtual screening and molecular dynamics simulation approaches

Spleen tyrosine kinase (SYK) is a non-receptor tyrosine kinase that plays an essential role in signal transduction across different cell types. In the context of allergy and autoimmune disorders, it is a crucial regulator of immune receptor signaling in inflammatory cells such as B cells, mast cells, macrophages, and neutrophils. Developing SYK kinase inhibitors has gained significant interest for potential therapeutic applications in neurological and cancer-related conditions. The clinical use of the most advanced SYK inhibitor, Fostamatinib, has been limited due to its unwanted side effects. Thus, a more targeted approach to SYK inhibition would provide a more comprehensive treatment window. In this study, we used a virtual screening approach to identify potential SYK inhibitors from natural compounds from the IMPPAT database. We identified two compounds, Isolysergic acid and Michelanugine, which showed strong affinity and specificity for the SYK binding pocket. All-atom molecular dynamics (MD) simulations were also performed to explore the stability, conformational changes, and interaction mechanism of SYK in complexes with the identified compounds. The identified compounds might have the potential to be developed into promising SYK inhibitors for the treatment of various diseases, including autoimmune disorders, cancer, and inflammatory diseases. This work aims to identify potential phytochemicals to develop a new protein kinase inhibitor for treating advanced malignancies by providing an updated understanding of the role of SYK.Communicated by Ramaswamy H. Sarma.

Theoretical Improvements in Enzyme Efficiency Associated with Noisy Rate Constants and Increased Dissipation

Previous studies have revealed the extraordinarily large catalytic efficiency of some enzymes. High catalytic proficiency is an essential accomplishment of biological evolution. Natural selection led to the increased turnover number, kcat, and enzyme efficiency, kcat/KM, of uni-uni enzymes, which convert a single substrate into a single product. We added or multiplied random noise with chosen rate constants to explore the correlation between dissipation and catalytic efficiency for ten enzymes: beta-galactosidase, glucose isomerase, β-lactamases from three bacterial strains, ketosteroid isomerase, triosephosphate isomerase, and carbonic anhydrase I, II, and T200H. Our results highlight the role of biological evolution in accelerating thermodynamic evolution. The catalytic performance of these enzymes is proportional to overall entropy production-the main parameter from irreversible thermodynamics. That parameter is also proportional to the evolutionary distance of β-lactamases PC1, RTEM, and Lac-1 when natural or artificial evolution produces the optimal or maximal possible catalytic efficiency. De novo enzyme design and attempts to speed up the rate-limiting catalytic steps may profit from the described connection between kinetics and thermodynamics.

Low-Barrier Hydrogen Bond Determines Target-Binding Affinity and Specificity of the Antitubercular Drug Bedaquiline

The role of short strong hydrogen bonds (SSHBs) in ligand-target binding remains largely unexplored, thereby hindering a potentially important avenue in rational drug design. Here we investigate the interaction between the antituberculosis drug bedaquiline (Bq) and the mycobacterial ATP synthase to unravel the role of a specific hydrogen bond to a conserved acidic residue in the target affinity and specificity. Our ab initio molecular dynamics simulations reveal that this bond belongs to the SSHB category and accounts for a substantial fraction of the target binding free energy. We also demonstrate that the presence of an extra acidic residue, i.e., aspartic acid at position 32 (D32), found exclusively in mycobacteria, cooperatively enhances the HB strength, ensuring specificity for the mycobacterial target. Consistently, we show that the removal of D32 markedly weakens the affinity, leading to Bq resistance associated with mutations of D32 to nonacidic residues. By designing simple Bq analogs, we then explore the possibility to overcome the resistance and potentially broaden the Bq antimicrobial spectrum by making the SSHB independent of the presence of the extra acidic residue.

Structure of a ribonucleotide reductase R2 protein radical

Aerobic ribonucleotide reductases (RNRs) initiate synthesis of DNA building blocks by generating a free radical within the R2 subunit; the radical is subsequently shuttled to the catalytic R1 subunit through proton-coupled electron transfer (PCET). We present a high-resolution room temperature structure of the class Ie R2 protein radical captured by x-ray free electron laser serial femtosecond crystallography. The structure reveals conformational reorganization to shield the radical and connect it to the translocation path, with structural changes propagating to the surface where the protein interacts with the catalytic R1 subunit. Restructuring of the hydrogen bond network, including a notably short O-O interaction of 2.41 angstroms, likely tunes and gates the radical during PCET. These structural results help explain radical handling and mobilization in RNR and have general implications for radical transfer in proteins.

Tubulin degradation: Principles, agents, and applications

The microtubule system plays an important role in the mitosis and growth of eukaryotic cells, and it is considered as an appealing and highly successful molecular target for cancer treatment. In fact, microtubule targeting agents, such as paclitaxel and vinblastine, have been approved by FDA for tumor therapy, which have achieved significant therapeutic effects and sales performance. At present, microtubule targeting agents mainly include microtubule-destabilizing agents, microtubule-stabilizing agents, and a few tubulin degradation agents. Although there are few reports about tubulin degradation agents at present, tubulin degradation agents show great potential in overcoming multidrug resistance and reducing neurotoxicity. In addition, some natural drugs could specifically degrade tubulin in tumor cells, but have no effect in normal cells, thus showing a good biosafety profile. Therefore, tubulin degradation agents might exhibit a better application. Currently, some small molecules have been designed to promote tubulin degradation with potent antiproliferative activities, showing the potential for cancer treatment. In this work, we reviewed the reports on tubulin degradation, and focused on the degradation mechanism and important functional groups of chemically synthesized compounds, hoping to provide help for the degradation design of tubulin.

Chemical features and machine learning assisted predictions of protein-ligand short hydrogen bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 Å closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Chemical Features and Machine Learning Assisted Predictions of Protein-Ligand Short Hydrogen Bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 A closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the "gold standard" coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller-Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

QM/MM Study of the Reaction Mechanism of Thermophilic Glucuronoyl Esterase for Biomass Treatment

Hydrolysis of lignocellulosic biomass, composed of a lignin-carbohydrate-complex (LCC) matrix, is critical for producing bioethanol from glucose. However, current methods for LCC processing require costly and polluting processes. The fungal Thermothelomyces thermophila glucuronoyl esterase (TtGE) is a promising thermophilic enzyme that hydrolyses LCC ester bonds. This study describes the TtGE catalytic mechanism using QM/MM methods. Two nearly-degenerate rate-determining transition states were found, with barriers of 16 and 17 kcal ⋅ mol^-1 , both with a zwitterionic nature that results from a proton interplay from His346 to either the Ser213-hydroxyl or the lignin leaving group and the rehybridisation of the ester moiety of the substrate to an alkoxide. An oxyanion hole, characteristic of esterases, was provided by the conserved Arg214 through its backbone and sidechain. Our work further suggests that a mutation of Glu267 to a non-negative residue will decrease the energetic barrier in ca. -5 kcal ⋅ mol^-1 , improving the catalytic rate of TtGE.

Advances in Protein-Based Nanocarriers of Bioactive Compounds: From Microscopic Molecular Principles to Macroscopical Structural and Functional Attributes

Many proteins can be used to fabricate nanocarriers for encapsulation, protection, and controlled release of nutraceuticals. This review examined the protein-based nanocarriers from microscopic molecular characteristics to the macroscopical structural and functional attributes. Structural, physical, and chemical properties of protein-based nanocarriers were introduced in detail. The spatial size, shape, water dispersibility, colloidal stability, etc. of protein-based nanocarriers were largely determined by the molecular physicochemical principles of protein. Different preparative techniques, including antisolvent precipitation, pH-driven, electrospray, and gelation methods, among others, can be used to fabricate different protein-based nanocarriers. Various modifications based on physical, chemical, and enzymatic approaches can be used to improve the functional performance of these nanocarriers. Protein is a natural resource with a wide range of sources, including plant, animal, and microbial, which are usually used to fabricate the nanocarriers. Protein-based nanocarriers have many advantages in aid of the application of bioactive ingredients to the medical, food, and cosmetic industries.

Advances in Rational Protein Engineering toward Functional Architectures and Their Applications in Food Science

Protein biomolecules including enzymes, cagelike proteins, and specific peptides have been continuously exploited as functional biomaterials applied in catalysis, nutrient delivery, and food preservation in food-related areas. However, natural proteins usually function well in physiological conditions, not industrial conditions, or may possess undesirable physical and chemical properties. Currently, rational protein design as a valuable technology has attracted extensive attention for the rational engineering or fabrication of ideal protein biomaterials with novel properties and functionality. This article starts with the underlying knowledge of protein folding and assembly and is followed by the introduction of the principles and strategies for rational protein design. Basic strategies for rational protein engineering involving experienced protein tailoring, computational prediction, computation redesign, and de novo protein design are summarized. Then, we focus on the recent progress of rational protein engineering or design in the application of food science, and a comprehensive summary ranging from enzyme manufacturing to cagelike protein nanocarriers engineering and antimicrobial peptides preparation is given. Overall, this review highlights the importance of rational protein engineering in food biomaterial preparation which could be beneficial for food science.

Structural insights into the catalytic and inhibitory mechanisms of the flavin transferase FmnB in Listeria monocytogenes

Listeria monocytogenes, a food-borne Gram-positive pathogen, often causes diseases such as gastroenteritis, bacterial sepsis, and meningitis. Newly discovered extracellular electron transfer (EET) from L. monocytogenes plays critical roles in the generation of redox molecules as electron carriers in bacteria. A Mg2+-dependent protein flavin mononucleotide (FMN) transferase (FmnB; UniProt: LMRG_02181) in EET is responsible for the transfer of electrons from intracellular to extracellular by hydrolyzing cofactor flavin adenine dinucleotide (FAD) and transferring FMN. FmnB homologs have been investigated in Gram-negative bacteria but have been less well studied in Gram-positive bacteria. In particular, the catalytic and inhibitory mechanisms of FmnB homologs remain elusive. Here, we report a series of crystal structures of apo-FmnB and FmnB complexed with substrate FAD, three inhibitors AMP, ADP, and ATP, revealing the unusual catalytic triad center (Asp301-Ser257-His273) of FmnB. The three inhibitors indeed inhibited the activity of FmnB in varying degrees by occupying the binding site of the FAD substrate. The key residue Arg262 of FmnB was profoundly affected by ADP but not AMP or ATP. Overall, our studies not only provide insights into the promiscuous ligand recognition behavior of FmnB but also shed light on its catalytic and inhibitory mechanisms.

Neutron crystallography reveals mechanisms used by Pseudomonas aeruginosa for host-cell binding

The opportunistic pathogen Pseudomonas aeruginosa, a major cause of nosocomial infections, uses carbohydrate-binding proteins (lectins) as part of its binding to host cells. The fucose-binding lectin, LecB, displays a unique carbohydrate-binding site that incorporates two closely located calcium ions bridging between the ligand and protein, providing specificity and unusually high affinity. Here, we investigate the mechanisms involved in binding based on neutron crystallography studies of a fully deuterated LecB/fucose/calcium complex. The neutron structure, which includes the positions of all the hydrogen atoms, reveals that the high affinity of binding may be related to the occurrence of a low-barrier hydrogen bond induced by the proximity of the two calcium ions, the presence of coordination rings between the sugar, calcium and LecB, and the dynamic behaviour of bridging water molecules at room temperature. These key structural details may assist in the design of anti-adhesive compounds to combat multi-resistance bacterial infections.

Pseudomonas aeruginosa employs lectins to bind to its host cells, and is known to be the major cause of lung infections. Lectin B (LecB) from Pseudomonas aeruginosa binds specifically to galactose and fucose and is important for pathogenicity, adhesion and biofilm formation. In this work, the neutron crystal structure (1.9 Å) of the deuterated LecB/Ca/fucose complex is reported. The structure, in combination with perdeuteration of the ligand and the receptor allowed the observation of hydrogen atoms, protonation states and hydrogen bonds involved in the interaction between pathogenic bacteria and host cells. Thus the study provides structural insights into the mechanism of high affinity binding of LecB to its targets.

Effective prediction of short hydrogen bonds in proteins via machine learning method

Short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms lie within 2.7 Å, exhibit prominent quantum mechanical characters and are connected to a wide range of essential biomolecular processes. However, exact determination of the geometry and functional roles of SHBs requires a protein to be at atomic resolution. In this work, we analyze 1260 high-resolution peptide and protein structures from the Protein Data Bank and develop a boosting based machine learning model to predict the formation of SHBs between amino acids. This model, which we name as machine learning assisted prediction of short hydrogen bonds (MAPSHB), takes into account 21 structural, chemical and sequence features and their interaction effects and effectively categorizes each hydrogen bond in a protein to a short or normal hydrogen bond. The MAPSHB model reveals that the type of the donor amino acid plays a major role in determining the class of a hydrogen bond and that the side chain Tyr-Asp pair demonstrates a significant probability of forming a SHB. Combining electronic structure calculations and energy decomposition analysis, we elucidate how the interplay of competing intermolecular interactions stabilizes the Tyr-Asp SHBs more than other commonly observed combinations of amino acid side chains. The MAPSHB model, which is freely available on our web server, allows one to accurately and efficiently predict the presence of SHBs given a protein structure with moderate or low resolution and will facilitate the experimental and computational refinement of protein structures.

Helium Nanodroplet Infrared Action Spectroscopy of the Proton-Bound Dimer of Hydrogen Sulfate and Formate: Examining Nuclear Quantum Effects

The proton-bound dimer of hydrogen sulfate and formate is an archetypal structure for ionic hydrogen-bonding complexes that contribute to biogenic aerosol nucleation. Of central importance for the structure and properties of this complex is the location of the bridging proton connecting the two conjugate base moieties. The potential energy surface for bridging proton translocation features two local minima, with the proton localized at either the formate or hydrogen sulfate moiety. However, electronic structure methods reveal a shallow potential energy surface governing proton translocation, with a barrier on the order of the zero-point energy. This shallow potential complicates structural assignment and necessitates a consideration of nuclear quantum effects. In this work, we probe the structure of this complex and its isotopologues, utilizing infrared (IR) action spectroscopy of ions captured in helium nanodroplets. The IR spectra indicate a structure in which a proton is shared between the hydrogen sulfate and formate moieties, HSO₄^-···H⁺···^-OOCH. However, because of the nuclear quantum effects and vibrational anharmonicities associated with the shallow potential for proton translocation, the extent of proton displacement from the formate moiety remains unclear, requiring further experiments or more advanced theoretical treatments for additional insight.

Protonation Equilibrium in the Active Site of the Photoactive Yellow Protein

The role and existence of low-barrier hydrogen bonds (LBHBs) in enzymatic and protein activity has been largely debated. An interesting case is that of the photoactive yellow protein (PYP). In this protein, two short HBs adjacent to the chromophore, p-coumaric acid (pCA), have been identified by X-ray and neutron diffraction experiments. However, there is a lack of agreement on the chemical nature of these H-bond interactions. Additionally, no consensus has been reached on the presence of LBHBs in the active site of the protein, despite various experimental and theoretical studies having been carried out to investigate this issue. In this work, we perform a computational study that combines classical and density functional theory (DFT)-based quantum mechanical/molecular mechanical (QM/MM) simulations to shed light onto this controversy. Furthermore, we aim to deepen our understanding of the chemical nature and dynamics of the protons involved in the two short hydrogen bonds that, in the dark state of PYP, connect pCA with the two binding pocket residues (E46 and Y42). Our results support the existence of a strong LBHB between pCA and E46, with the H fully delocalized and shared between both the carboxylic oxygen of E46 and the phenolic oxygen of pCA. Additionally, our findings suggest that the pCA interaction with Y42 can be suitably described as a typical short ionic H-bond of moderate strength that is fully localized on the phenolic oxygen of Y42.