Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors.

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MProWT) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373), and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ∼400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decrease in free energy (ΔG) and thermal stability (Tm), respectively, relative to MProWT. The H41A mutation results in a 20-fold increase in the dissociation constant (Kd), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and Tm, respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MProH41A, whereas no significant change is observed in binding to MProWT. Structures of the four inhibitor complexes with MPro1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MProWT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MProWT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MProC145A also results in a significant, but smaller, increase in Kd and decrease in ΔG and Tm, relative to NMV.

AI-Accelerated Design of Targeted Covalent Inhibitors for SARS-CoV-2.

Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic's evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (Mpro), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. This work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic "warhead" to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of Mpro with micromolar affinities (KI of 5.27 µM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.

Covalent inhibition of hAChE by organophosphates causes homodimer dissociation through long-range allosteric effects.

Acetylcholinesterase (EC 3.1.1.7), a key acetylcholine-hydrolyzing enzyme in cholinergic neurotransmission, is present in a variety of states in situ, including monomers, C-terminally disulfide-linked homodimers, homotetramers, and up to three tetramers covalently attached to structural subunits. Could oligomerization that ensures high local concentrations of catalytic sites necessary for efficient neurotransmission be affected by environmental factors? Using small-angle X-ray scattering (SAXS) and cryo-EM, we demonstrate that homodimerization of recombinant monomeric human acetylcholinesterase (hAChE) in solution occurs through a C-terminal four-helix bundle at micromolar concentrations. We show that diethylphosphorylation of the active serine in the catalytic gorge or isopropylmethylphosphonylation by the RP enantiomer of sarin promotes a 10-fold increase in homodimer dissociation. We also demonstrate the dissociation of organophosphate (OP)-conjugated dimers is reversed by structurally diverse oximes 2PAM, HI6, or RS194B, as demonstrated by SAXS of diethylphosphoryl-hAChE. However, binding of oximes to the native ligand-free hAChE, binding of high-affinity reversible ligands, or formation of an SP-sarin-hAChE conjugate had no effect on homodimerization. Dissociation monitored by time-resolved SAXS occurs in milliseconds, consistent with rates of hAChE covalent inhibition. OP-induced dissociation was not observed in the SAXS profiles of the double-mutant Y337A/F338A, where the active center gorge volume is larger than in wildtype hAChE. These observations suggest a key role of the tightly packed acyl pocket in allosterically triggered OP-induced dimer dissociation, with the potential for local reduction of acetylcholine-hydrolytic power in situ. Computational models predict allosteric correlated motions extending from the acyl pocket toward the four-helix bundle dimerization interface 25 Å away.

High-Throughput Virtual Screening and Validation of a SARS-CoV-2 Main Protease Noncovalent Inhibitor.

Despite the recent availability of vaccines against the acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the search for inhibitory therapeutic agents has assumed importance especially in the context of emerging new viral variants. In this paper, we describe the discovery of a novel noncovalent small-molecule inhibitor, MCULE-5948770040, that binds to and inhibits the SARS-Cov-2 main protease (Mpro) by employing a scalable high-throughput virtual screening (HTVS) framework and a targeted compound library of over 6.5 million molecules that could be readily ordered and purchased. Our HTVS framework leverages the U.S. supercomputing infrastructure achieving nearly 91% resource utilization and nearly 126 million docking calculations per hour. Downstream biochemical assays validate this Mpro inhibitor with an inhibition constant (Ki) of 2.9 µM (95% CI 2.2, 4.0). Furthermore, using room-temperature X-ray crystallography, we show that MCULE-5948770040 binds to a cleft in the primary binding site of Mpro forming stable hydrogen bond and hydrophobic interactions. We then used multiple µs-time scale molecular dynamics (MD) simulations and machine learning (ML) techniques to elucidate how the bound ligand alters the conformational states accessed by Mpro, involving motions both proximal and distal to the binding site. Together, our results demonstrate how MCULE-5948770040 inhibits Mpro and offers a springboard for further therapeutic design.

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor.

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ∼5-50 cm-1. This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, ß-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effects of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal mol-1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

Unusual zwitterionic catalytic site of SARS-CoV-2 main protease revealed by neutron crystallography.

The main protease (3CL Mpro) from SARS-CoV-2, the etiological agent of COVID-19, is an essential enzyme for viral replication. 3CL Mpro possesses an unusual catalytic dyad composed of Cys145 and His41 residues. A critical question in the field has been what the protonation states of the ionizable residues in the substrate-binding active-site cavity are; resolving this point would help understand the catalytic details of the enzyme and inform rational drug development against this pernicious virus. Here, we present the room-temperature neutron structure of 3CL Mpro, which allowed direct determination of hydrogen atom positions and, hence, protonation states in the protease. We observe that the catalytic site natively adopts a zwitterionic reactive form in which Cys145 is in the negatively charged thiolate state and His41 is doubly protonated and positively charged, instead of the neutral unreactive state usually envisaged. The neutron structure also identified the protonation states, and thus electrical charges, of all other amino acid residues and revealed intricate hydrogen-bonding networks in the active-site cavity and at the dimer interface. The fine atomic details present in this structure were made possible by the unique scattering properties of the neutron, which is an ideal probe for locating hydrogen positions and experimentally determining protonation states at near-physiological temperature. Our observations provide critical information for structure-assisted and computational drug design, allowing precise tailoring of inhibitors to the enzyme's electrostatic environment.

Rational design, synthesis, and evaluation of uncharged, "smart" bis-oxime antidotes of organophosphate-inhibited human acetylcholinesterase.

Organophosphate (OP) intoxications from nerve agent and OP pesticide exposures are managed with pyridinium aldoxime-based therapies whose success rates are currently limited. The pyridinium cation hampers uptake of OPs into the central nervous system (CNS). Furthermore, it frequently binds to aromatic residues of OP-inhibited acetylcholinesterase (AChE) in orientations that are nonproductive for AChE reactivation, and the structural diversity of OPs impedes efficient reactivation. Improvements of OP antidotes need to include much better access of AChE reactivators to the CNS and optimized orientation of the antidotes' nucleophile within the AChE active-center gorge. On the basis of X-ray structures of a CNS-penetrating reactivator, monoxime RS194B, reversibly bound to native and venomous agent X (VX)-inhibited human AChE, here we created seven uncharged acetamido bis-oximes as candidate antidotes. Both oxime groups in these bis-oximes were attached to the same central, saturated heterocyclic core. Diverse protonation of the heterocyclic amines and oxime groups of the bis-oximes resulted in equilibration among up to 16 distinct ionization forms, including uncharged forms capable of diffusing into the CNS and multiple zwitterionic forms optimal for reactivation reactions. Conformationally diverse zwitterions that could act as structural antidote variants significantly improved in vitro reactivation of diverse OP-human AChE conjugates. Oxime group reorientation of one of the bis-oximes, forcing it to point into the active center for reactivation, was confirmed by X-ray structural analysis. Our findings provide detailed structure-activity properties of several CNS-directed, uncharged aliphatic bis-oximes holding promise for use as protonation-dependent, conformationally adaptive, "smart" accelerated antidotes against OP toxicity.

Characterization and structural analysis of a thermophilic GH11 xylanase from compost metatranscriptome.

Xylanase is efficient for xylan degradation and widely applied in industries. We found a GH11 family xylanase (Xyn11A) with high thermostability and catalytic activity from compost metatranscriptome. This xylanase has the optimal reaction temperature at 80 °C with the activity of 2907.3 U/mg. The X-ray crystallographic structure shows a typical "right hand" architecture, which is the characteristics of the GH11 family enzymes. Comparing it with the mesophilic XYN II, a well-studied GH11 xylanase from Trichoderma reesei, Xyn11A is more compact with more H-bonds. Our mutagenic results show that the electrostatic interactions in the thumb and palm region of Xyn11A could result in its high thermostability and activity. Introducing a disulfide bond at the N-terminus further increased its optimal reaction temperature to 90 °C with augmented activity. KEY POINTS: • A hyperthermophilic xylanase with high activity was discovered using the metatranscriptomic method. • The mechanisms of thermophilicity and high activity were revealed using X-ray crystallography, mutagenesis, and molecular dynamics simulations. • The thermostability and activity were further improved by introducing a disulfide bond.

Productive reorientation of a bound oxime reactivator revealed in room temperature X-ray structures of native and VX-inhibited human acetylcholinesterase.

Exposure to organophosphorus compounds (OPs) may be fatal if untreated, and a clear and present danger posed by nerve agent OPs has become palpable in recent years. OPs inactivate acetylcholinesterase (AChE) by covalently modifying its catalytic serine. Inhibited AChE cannot hydrolyze the neurotransmitter acetylcholine leading to its build-up at the cholinergic synapses and creating an acute cholinergic crisis. Current antidotes, including oxime reactivators that attack the OP-AChE conjugate to free the active enzyme, are inefficient. Better reactivators are sought, but their design is hampered by a conformationally rigid portrait of AChE extracted exclusively from 100K X-ray crystallography and scarcity of structural knowledge on human AChE (hAChE). Here, we present room temperature X-ray structures of native and VX-phosphonylated hAChE with an imidazole-based oxime reactivator, RS-170B. We discovered that inhibition with VX triggers substantial conformational changes in bound RS-170B from a "nonproductive" pose (the reactive aldoxime group points away from the VX-bound serine) in the reactivator-only complex to a "semi-productive" orientation in the VX-modified complex. This observation, supported by concurrent molecular simulations, suggested that the narrow active-site gorge of hAChE may be significantly more dynamic than previously thought, allowing RS-170B to reorient inside the gorge. Furthermore, we found that small molecules can bind in the choline-binding site hindering approach to the phosphorous of VX-bound serine. Our results provide structural and mechanistic perspectives on the reactivation of OP-inhibited hAChE and demonstrate that structural studies at physiologically relevant temperatures can deliver previously overlooked insights applicable for designing next-generation antidotes.

Neutron structures of the Helicobacter pylori 5'-methylthioadenosine nucleosidase highlight proton sharing and protonation states.

MTAN (5'-methylthioadenosine nucleosidase) catalyzes the hydrolysis of the N-ribosidic bond of a variety of adenosine-containing metabolites. The Helicobacter pylori MTAN (HpMTAN) hydrolyzes 6-amino-6-deoxyfutalosine in the second step of the alternative menaquinone biosynthetic pathway. Substrate binding of the adenine moiety is mediated almost exclusively by hydrogen bonds, and the proposed catalytic mechanism requires multiple proton-transfer events. Of particular interest is the protonation state of residue D198, which possesses a pKa above 8 and functions as a general acid to initiate the enzymatic reaction. In this study we present three corefined neutron/X-ray crystal structures of wild-type HpMTAN cocrystallized with S-adenosylhomocysteine (SAH), Formycin A (FMA), and (3R,4S)-4-(4-Chlorophenylthiomethyl)-1-[(9-deaza-adenin-9-yl)methyl]-3-hydroxypyrrolidine (p-ClPh-Thio-DADMe-ImmA) as well as one neutron/X-ray crystal structure of an inactive variant (HpMTAN-D198N) cocrystallized with SAH. These results support a mechanism of D198 pKa elevation through the unexpected sharing of a proton with atom N7 of the adenine moiety possessing unconventional hydrogen-bond geometry. Additionally, the neutron structures also highlight active site features that promote the stabilization of the transition state and slight variations in these interactions that result in 100-fold difference in binding affinities between the DADMe-ImmA and ImmA analogs.

Neutron Crystallography Detects Differences in Protein Dynamics: Structure of the PKG II Cyclic Nucleotide Binding Domain in Complex with an Activator.

As one of the main receptors of a second messenger, cGMP, cGMP-dependent protein kinase (PKG) isoforms I and II regulate distinct physiological processes. The design of isoform-specific activators is thus of great biomedical importance and requires detailed structural information about PKG isoforms bound with activators, including accurate positions of hydrogen atoms and a description of the hydrogen bonding and water architecture. Here, we determined a 2.2 Å room-temperature joint X-ray/neutron (XN) structure of the human PKG II carboxyl cyclic nucleotide binding (CNB-B) domain bound with a potent PKG II activator, 8-pCPT-cGMP. The XN structure directly visualizes intermolecular interactions and reveals changes in hydrogen bonding patterns upon comparison to the X-ray structure determined at cryo-temperatures. Comparative analysis of the backbone hydrogen/deuterium exchange patterns in PKG II:8-pCPT-cGMP and previously reported PKG Iß:cGMP XN structures suggests that the ability of these agonists to activate PKG is related to how effectively they quench dynamics of the cyclic nucleotide binding pocket and the surrounding regions.

Direct evidence that an extended hydrogen-bonding network influences activation of pyridoxal 5'-phosphate in aspartate aminotransferase.

Pyridoxal 5'-phosphate (PLP) is a fundamental, multifunctional enzyme cofactor used to catalyze a wide variety of chemical reactions involved in amino acid metabolism. PLP-dependent enzymes optimize specific chemical reactions by modulating the electronic states of PLP through distinct active site environments. In aspartate aminotransferase (AAT), an extended hydrogen bond network is coupled to the pyridinyl nitrogen of the PLP, influencing the electrophilicity of the cofactor. This network, which involves residues Asp-222, His-143, Thr-139, His-189, and structural waters, is located at the edge of PLP opposite the reactive Schiff base. We demonstrate that this hydrogen bond network directly influences the protonation state of the pyridine nitrogen of PLP, which affects the rates of catalysis. We analyzed perturbations caused by single- and double-mutant variants using steady-state kinetics, high resolution X-ray crystallography, and quantum chemical calculations. Protonation of the pyridinyl nitrogen to form a pyridinium cation induces electronic delocalization in the PLP, which correlates with the enhancement in catalytic rate in AAT. Thus, PLP activation is controlled by the proximity of the pyridinyl nitrogen to the hydrogen bond microenvironment. Quantum chemical calculations indicate that Asp-222, which is directly coupled to the pyridinyl nitrogen, increases the pKa of the pyridine nitrogen and stabilizes the pyridinium cation. His-143 and His-189 also increase the pKa of the pyridine nitrogen but, more significantly, influence the position of the proton that resides between Asp-222 and the pyridinyl nitrogen. These findings indicate that the second shell residues directly enhance the rate of catalysis in AAT.

Direct determination of protonation states and visualization of hydrogen bonding in a glycoside hydrolase with neutron crystallography.

Glycoside hydrolase (GH) enzymes apply acid/base chemistry to catalyze the decomposition of complex carbohydrates. These ubiquitous enzymes accept protons from solvent and donate them to substrates at close to neutral pH by modulating the pKa values of key side chains during catalysis. However, it is not known how the catalytic acid residue acquires a proton and transfers it efficiently to the substrate. To better understand GH chemistry, we used macromolecular neutron crystallography to directly determine protonation and ionization states of the active site residues of a family 11 GH at multiple pD (pD=pH+0.4) values. The general acid glutamate (Glu) cycles between two conformations, upward and downward, but is protonated only in the downward orientation. We performed continuum electrostatics calculations to estimate the pKa values of the catalytic Glu residues in both the apo- and substrate-bound states of the enzyme. The calculated pKa of the Glu increases substantially when the side chain moves down. The energy barrier required to rotate the catalytic Glu residue back to the upward conformation, where it can protonate the glycosidic oxygen of the substrate, is 4.3 kcal/mol according to free energy simulations. These findings shed light on the initial stage of the glycoside hydrolysis reaction in which molecular motion enables the general acid catalyst to obtain a proton from the bulk solvent and deliver it to the glycosidic oxygen.

Joint neutron crystallographic and NMR solution studies of Tyr residue ionization and hydrogen bonding: Implications for enzyme-mediated proton transfer.

Human carbonic anhydrase II (HCA II) uses a Zn-bound OH(-)/H2O mechanism to catalyze the reversible hydration of CO2. This catalysis also involves a separate proton transfer step, mediated by an ordered solvent network coordinated by hydrophilic residues. One of these residues, Tyr7, was previously shown to be deprotonated in the neutron crystal structure at pH 10. This observation indicated that Tyr7 has a perturbed pKa compared with free tyrosine. To further probe the pKa of this residue, NMR spectroscopic measurements of [(13)C]Tyr-labeled holo HCA II (with active-site Zn present) were preformed to titrate all Tyr residues between pH 5.4-11.0. In addition, neutron studies of apo HCA II (with Zn removed from the active site) at pH 7.5 and holo HCA II at pH 6 were conducted. This detailed interrogation of tyrosines in HCA II by NMR and neutron crystallography revealed a significantly lowered pKa of Tyr7 and how pH and Tyr proximity to Zn affect hydrogen-bonding interactions.

Mannobiose Binding Induces Changes in Hydrogen Bonding and Protonation States of Acidic Residues in Concanavalin A As Revealed by Neutron Crystallography.

Plant lectins are carbohydrate-binding proteins with various biomedical applications. Concanavalin A (Con A) holds promise in treating cancerous tumors. To better understand the Con A carbohydrate binding specificity, we obtained a room-temperature neutron structure of this legume lectin in complex with a disaccharide Manα1-2Man, mannobiose. The neutron structure afforded direct visualization of the hydrogen bonding between the protein and ligand, showing that the ligand is able to alter both protonation states and interactions for residues located close to and distant from the binding site. An unprecedented low-barrier hydrogen bond was observed forming between the carboxylic side chains of Asp28 and Glu8, with the D atom positioned equidistant from the oxygen atoms having an O···D···O angle of 101.5°.

Neutron and Atomic Resolution X-ray Structures of a Lytic Polysaccharide Monooxygenase Reveal Copper-Mediated Dioxygen Binding and Evidence for N-Terminal Deprotonation.

A 1.1 Å resolution, room-temperature X-ray structure and a 2.1 Å resolution neutron structure of a chitin-degrading lytic polysaccharide monooxygenase domain from the bacterium Jonesia denitrificans (JdLPMO10A) show a putative dioxygen species equatorially bound to the active site copper. Both structures show an elongated density for the dioxygen, most consistent with a Cu(II)-bound peroxide. The coordination environment is consistent with Cu(II). In the neutron and X-ray structures, difference maps reveal the N-terminal amino group, involved in copper coordination, is present as a mixed ND2 and ND-, suggesting a role for the copper ion in shifting the pKa of the amino terminus.

Toward resolving the catalytic mechanism of dihydrofolate reductase using neutron and ultrahigh-resolution X-ray crystallography.

Dihydrofolate reductase (DHFR) catalyzes the NADPH-dependent reduction of dihydrofolate (DHF) to tetrahydrofolate (THF). An important step in the mechanism involves proton donation to the N5 atom of DHF. The inability to determine the protonation states of active site residues and substrate has led to a lack of consensus regarding the catalytic mechanism involved. To resolve this ambiguity, we conducted neutron and ultrahigh-resolution X-ray crystallographic studies of the pseudo-Michaelis ternary complex of Escherichia coli DHFR with folate and NADP(+). The neutron data were collected to 2.0-Å resolution using a 3.6-mm(3) crystal with the quasi-Laue technique. The structure reveals that the N3 atom of folate is protonated, whereas Asp27 is negatively charged. Previous mechanisms have proposed a keto-to-enol tautomerization of the substrate to facilitate protonation of the N5 atom. The structure supports the existence of the keto tautomer owing to protonation of the N3 atom, suggesting that tautomerization is unnecessary for catalysis. In the 1.05-Å resolution X-ray structure of the ternary complex, conformational disorder of the Met20 side chain is coupled to electron density for a partially occupied water within hydrogen-bonding distance of the N5 atom of folate; this suggests direct protonation of substrate by solvent. We propose a catalytic mechanism for DHFR that involves stabilization of the keto tautomer of the substrate, elevation of the pKa value of the N5 atom of DHF by Asp27, and protonation of N5 by water that gains access to the active site through fluctuation of the Met20 side chain even though the Met20 loop is closed.

Phosphoryl Transfer Reaction Snapshots in Crystals: INSIGHTS INTO THE MECHANISM OF PROTEIN KINASE A CATALYTIC SUBUNIT.

To study the catalytic mechanism of phosphorylation catalyzed by cAMP-dependent protein kinase (PKA) a structure of the enzyme-substrate complex representing the Michaelis complex is of specific interest as it can shed light on the structure of the transition state. However, all previous crystal structures of the Michaelis complex mimics of the PKA catalytic subunit (PKAc) were obtained with either peptide inhibitors or ATP analogs. Here we utilized Ca(2+) ions and sulfur in place of the nucleophilic oxygen in a 20-residue pseudo-substrate peptide (CP20) and ATP to produce a close mimic of the Michaelis complex. In the ternary reactant complex, the thiol group of Cys-21 of the peptide is facing Asp-166 and the sulfur atom is positioned for an in-line phosphoryl transfer. Replacement of Ca(2+) cations with Mg(2+) ions resulted in a complex with trapped products of ATP hydrolysis: phosphate ion and ADP. The present structural results in combination with the previously reported structures of the transition state mimic and phosphorylated product complexes complete the snapshots of the phosphoryl transfer reaction by PKAc, providing us with the most thorough picture of the catalytic mechanism to date.

Long-Range Electrostatics-Induced Two-Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site.

Neutron crystallography was used to directly locate two protons before and after a pH-induced two-proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV-1 protease. The two-proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low-pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.