Three-Photon Infrared Stimulation of Endogenous Neuroreceptors in Vivo.

To interrogate neural circuits and crack their codes, in vivo brain activity imaging must be combined with spatiotemporally precise stimulation in three dimensions using genetic or pharmacological specificity. This challenge requires deep penetration and focusing as provided by infrared light and multiphoton excitation, and has promoted two-photon photopharmacology and optogenetics. However, three-photon brain stimulation in vivo remains to be demonstrated. We report the regulation of neuronal activity in zebrafish larvae by three-photon excitation of a photoswitchable muscarinic agonist at 50 pM, a billion-fold lower concentration than used for uncaging, and with mid-infrared light of 1560 nm, the longest reported photoswitch wavelength. Robust, physiologically relevant photoresponses allow modulating brain activity in wild-type animals with spatiotemporal and pharmacological precision. Computational calculations predict that azobenzene-based ligands have high three-photon absorption cross-section and can be used directly with pulsed infrared light. The expansion of three-photon pharmacology will deeply impact basic neurobiology and neuromodulation phototherapies.

Different Structures-Similar Effect: Do Substituted 5-(4-Methoxyphenyl)-1H-indoles and 5-(4-Methoxyphenyl)-1H-imidazoles Represent a Common Pharmacophore for Substrate Selective Inhibition of Linoleate Oxygenase Activity of ALOX15?

Mammalian 15-lipoxygenases (ALOX15) are lipid peroxidizing enzymes that exhibit variable functionality in different cancer and inflammation models. The pathophysiological role of linoleic acid- and arachidonic acid-derived ALOX15 metabolites rendered this enzyme a target for pharmacological research. Several indole and imidazole derivatives inhibit the catalytic activity of rabbit ALOX15 in a substrate-specific manner, but the molecular basis for this allosteric inhibition remains unclear. Here, we attempt to define a common pharmacophore, which is critical for this allosteric inhibition. We found that substituted imidazoles induce weaker inhibitory effects when compared with the indole derivatives. In silico docking studies and molecular dynamics simulations using a dimeric allosteric enzyme model, in which the inhibitor occupies the substrate-binding pocket of one monomer, whereas the substrate fatty acid is bound at the catalytic center of another monomer within the ALOX15 dimer, indicated that chemical modification of the core pharmacophore alters the enzyme-inhibitor interactions, inducing a reduced inhibitory potency. In our dimeric ALOX15 model, the structural differences induced by inhibitor binding are translated to the hydrophobic dimerization cluster and affect the structures of enzyme-substrate complexes. These data are of particular importance since substrate-specific inhibition may contribute to elucidation of the putative roles of ALOX15 metabolites derived from different polyunsaturated fatty acids in mammalian pathophysiology.

Functional Characterization of Mouse and Human Arachidonic Acid Lipoxygenase 15B (ALOX15B) Orthologs and of Their Mutants Exhibiting Humanized and Murinized Reaction Specificities.

The arachidonic acid lipoxygenase 15B (ALOX15B) orthologs of men and mice form different reaction products when arachidonic acid is used as the substrate. Tyr603Asp+His604Val double mutation in mouse arachidonic acid lipoxygenase 15b humanized the product pattern and an inverse mutagenesis strategy murinized the specificity of the human enzyme. As the mechanistic basis for these functional differences, an inverse substrate binding at the active site of the enzymes has been suggested, but experimental proof for this hypothesis is still pending. Here we expressed wildtype mouse and human arachidonic acid lipoxygenase 15B orthologs as well as their humanized and murinized double mutants as recombinant proteins and analyzed the product patterns of these enzymes with different polyenoic fatty acids. In addition, in silico substrate docking studies and molecular dynamics simulation were performed to explore the mechanistic basis for the distinct reaction specificities of the different enzyme variants. Wildtype human arachidonic acid lipoxygenase 15B converted arachidonic acid and eicosapentaenoic acid to their 15-hydroperoxy derivatives but the Asp602Tyr+Val603His exchange murinized the product pattern. The inverse mutagenesis strategy in mouse arachidonic acid lipoxygenase 15b (Tyr603Asp+His604Val exchange) humanized the product pattern with these substrates, but the situation was different with docosahexaenoic acid. Here, Tyr603Asp+His604Val substitution in mouse arachidonic acid lipoxygenase 15b also humanized the specificity but the inverse mutagenesis (Asp602Tyr+Val603His) did not murinize the human enzyme. With linoleic acid Tyr603Asp+His604Val substitution in mouse arachidonic acid lipoxygenase 15b humanized the product pattern but the inverse mutagenesis in human arachidonic acid lipoxygenase 15B induced racemic product formation. Amino acid exchanges at critical positions of human and mouse arachidonic acid lipoxygenase 15B orthologs humanized/murinized the product pattern with C20 fatty acids, but this was not the case with fatty acid substrates of different chain lengths. Asp602Tyr+Val603His exchange murinized the product pattern of human arachidonic acid lipoxygenase 15B with arachidonic acid, eicosapentaenoic acid, and docosahexaenoic acid. An inverse mutagenesis strategy on mouse arachidonic acid lipoxygenase 15b (Tyr603Asp+His604Val exchange) did humanize the reaction products with arachidonic acid and eicosapentaenoic acid, but not with docosahexaenoic acid.

Hydroperoxidation of Docosahexaenoic Acid by Human ALOX12 and pigALOX15-mini-LOX.

Human lipoxygenase 12 (hALOX12) catalyzes the conversion of docosahexaenoic acid (DHA) into mainly 14S-hydroperoxy-4Z,7Z,10Z,12E,16Z,19Z-docosahexaenoic acid (14S-H(p)DHA). This hydroperoxidation reaction is followed by an epoxidation and hydrolysis process that finally leads to maresin 1 (MaR1), a potent bioactive specialized pro-resolving mediator (SPM) in chronic inflammation resolution. By combining docking, molecular dynamics simulations, and quantum mechanics/molecular mechanics calculations, we have computed the potential energy profile of DHA hydroperoxidation in the active site of hALOX12. Our results describe the structural evolution of the molecular system at each step of this catalytic reaction pathway. Noteworthy, the required stereospecificity of the reaction leading to MaR1 is explained by the configurations adopted by DHA bound to hALOX12, along with the stereochemistry of the pentadienyl radical formed after the first step of the mechanism. In pig lipoxygenase 15 (pigALOX15-mini-LOX), our calculations suggest that 14S-H(p)DHA can be formed, but with a stereochemistry that is inadequate for MaR1 biosynthesis.

Molecular Insights into the Regulation of 3-Phosphoinositide-Dependent Protein Kinase 1: Modeling the Interaction between the Kinase and the Pleckstrin Homology Domains.

The 3-phosphoinositide-dependent protein kinase 1 (PDK1) K465E mutant kinase can still activate protein kinase B (PKB) at the membrane in a phosphatidylinositol-3,4,5-trisphosphate (PIP3, PtdIns(3,4,5)P3) independent manner. To understand this new PDK1 regulatory mechanism, docking and molecular dynamics calculations were performed for the first time to simulate the wild-type kinase domain-pleckstrin homology (PH) domain complex with PH-in and PH-out conformations. These simulations were then compared to the PH-in model of the KD-PH(mutant K465E) PDK1 complex. Additionally, three KD-PH complexes were simulated, including a substrate analogue bound to a hydrophobic pocket (denominated the PIF-pocket) substrate-docking site. We find that only the PH-out conformation, with the PH domain well-oriented to interact with the cellular membrane, is active for wild-type PDK1. In contrast, the active conformation of the PDK1 K465E mutant is PH-in, being ATP-stable at the active site while the PIF-pocket is more accessible to the peptide substrate. We corroborate that both the docking-site binding and the catalytic activity are in fact enhanced in knock-in mouse samples expressing the PDK1 K465E protein, enabling the phosphorylation of PKB in the absence of PIP3 binding.

On the Computational Design of Azobenzene-Based Multi-State Photoswitches.

In order to theoretically design multi-state photoswitches with specific properties, an exhaustive computational study is first carried out for an azobenzene dimer that has been recently synthesized and experimentally studied. This study allows for a full comprehension of the factors that govern the photoactivated isomerization processes of these molecules so to provide a conceptual/computational protocol that can be applied to generic multi-state photoswitches. From this knowledge a new dimer with a similar chemical design is designed and also fully characterized. Our theoretical calculations predict that the new dimer proposed is one step further in the quest for a double photoswitch, where the four metastable isomers could be selectively interconverted through the use of different irradiation sequences.

Theoretical Characterization of the Step-by-Step Mechanism of Conversion of Leukotriene A4 to Leukotriene B4 Catalysed by the Enzyme Leukotriene A4 Hydrolase.

LTA4H is a bifunctional zinc metalloenzyme that converts leukotriene A4 (LTA4) into leukotriene B4 (LTB4), one of the most potent chemotactic agents involved in acute and chronic inflammatory diseases. In this reaction, LTA4H acts as an epoxide hydrolase with a unique and fascinating mechanism, which includes the stereoselective attachment of one water molecule to the carbon backbone of LTA4 several methylene units away from the epoxide moiety. By combining Molecular Dynamics simulations and Quantum Mechanics/Molecular Mechanics calculations, we obtained a very detailed molecular picture of the different consecutive steps of that mechanism. By means of a rather unusual 1,7-nucleophilic substitution through a clear SN1 mechanism, the epoxide opens and the triene moiety of the substrate twists in such a way that the bond C6-C7 adopts its cis (Z) configuration, thus exposing the R face of C12 to the addition of a water molecule hydrogen-bonded to ASP375. Thus, the two stereochemical features that are required for the bioactivity of LTB4 appear to be closely related. The noncovalent π-π stacking interactions between the triene moiety and two tyrosines (TYR267 and, especially, TYR378) that wrap the triene system along the whole reaction explain the preference for the cis configuration inside LTA4H.

The role of acetylated cyclooxygenase-2 in the biosynthesis of resolvin precursors derived from eicosapentaenoic acid.

Specialized pro-resolving lipid mediators (SPMs) are natural bioactive agents actively involved in inflammation resolution. SPMs act when uncontrolled inflammatory processes are developed, for instance, in patients of COVID-19 or other diseases. The so-called resolution pharmacology aims at developing new treatments based on the use of SPMs as agonists, which promote inflammation resolution without unwanted side effects. It has been shown that the biosynthesis of SPMs called eicosapentaenoic acid (EPA)-derived E-series resolvins is initiated by aspirin-acetylated COX-2 from EPA, leading to 18-hydroperoxy-eicosapentaenoic acid (18-HpEPE). However, there are many open questions concerning the intriguing role of aspirin in the molecular mechanism of resolvin formation. Our MD simulations, combined with QM/MM calculations, show that the potential energy barriers for the H16-abstraction from EPA, required for forming 18-HpEPE, are higher than for the H13-abstraction, thus explaining why 18-HpEPE is a marginal product of COX-2 catalysis. By contrast, in the aspirin-acetylated COX-2/EPA complex, the H16proS-abstraction energy barriers are somewhat lower than the H13proS energy barriers and much smaller than the H16-transfer barriers in the wild type COX-2/EPA system. Those results agree with the experimental observation that aspirin favours the synthesis of several SPMs known as aspirin-triggered resolvins. In the following step of the catalytic mechanism, the calculated O2 addition to C18 is preferred versus the addition to C14 which also agrees with 18R-HEPE and 18S-HEPE being the main products from EPA in aspirin-acetylated COX-2.

N-Substituted 5-(1H-Indol-2-yl)-2-methoxyanilines Are Allosteric Inhibitors of the Linoleate Oxygenase Activity of Selected Mammalian ALOX15 Orthologs: Mechanism of Action.

Here, we describe the first systematic study on the mechanism of substrate-selective inhibition of mammalian ALOX15 orthologs. For this purpose, we prepared a series of N-substituted 5-(1H-indol-2-yl)anilines and found that (N-(5-(1H-indol-2-yl)-2-methoxyphenyl)sulfamoyl)carbamates and their monofluorinated analogues are potent and selective inhibitors of the linoleate oxygenase activity of rabbit and human ALOX15. Introduction of a 2-methoxyaniline moiety into the core pharmacophore plays a crucial role in substrate-selective inhibition of ALOX15-catalyzed oxygenation of linoleic acid at submicromolar concentrations without affecting arachidonic acid oxygenation. Steady-state kinetics, mutagenesis studies, and molecular dynamics (MD) simulations suggested an allosteric mechanism of action. Using a dimer model of ALOX15, our MD simulations suggest that the binding of the inhibitor at the active site of one monomer induces conformational alterations in the other monomer so that the formation of a productive enzyme-linoleic acid complex is energetically compromised.

Predicting the Electronic Absorption Band Shape of Azobenzene Photoswitches.

Simulations based on molecular dynamics coupled to excitation energy calculations were used to generate simulated absorption spectra for a family of halide derivatives of azobenzene, a family of photoswitch molecules with a weak absorption band around 400-600 nm and potential uses in living tissue. This is a case where using the conventional approach in theoretical spectroscopy (estimation of absorption maxima based on the vertical transition from the potential energy minimum on the ground electronic state) does not provide valid results that explain how the observed band shape extends towards the low energy region of the spectrum. The method affords a reasonable description of the main features of the low-energy UV-Vis spectra of these compounds. A bathochromic trend was detected linked to the size of the halide atom. Analysis of the excitation reveals a correlation between the energy of the molecular orbital where excitation starts and the energy of the highest occupied atomic orbital of the free halide atom. This was put to the test with a new brominated compound with good results. The energy level of the highest occupied orbital on the free halide was identified as a key factor that strongly affects the energy gap in the photoswitch. This opens the way for the design of bathochromically shifted variants of the photoswitch with possible applications.

Accounting for the instantaneous disorder in the enzyme-substrate Michaelis complex to calculate the Gibbs free energy barrier of an enzyme reaction.

Many enzyme reactions present instantaneous disorder. These dynamic fluctuations in the enzyme-substrate Michaelis complexes generate a wide range of energy barriers that cannot be experimentally observed, but that determine the measured kinetics of the reaction. These individual energy barriers can be calculated using QM/MM methods, but then the problem is how to deal with this dispersion of energy barriers to provide kinetic information. So far, the most usual procedure has implied the so-called exponential average of the energy barriers. In this paper, we discuss the foundations of this method, and we use the free energy perturbation theory to derive an alternative equation to get the Gibbs free energy barrier of the enzyme reaction. In addition, we propose a practical way to implement it. We have chosen four enzyme reactions as examples. In particular, we have studied the hydrolysis of a glycosidic bond catalyzed by the enzyme Thermus thermophilus ß-glycosidase, and the mutant Y284P Ttb-gly, and the hydrogen abstraction reactions from C13 and C7 of arachidonic acid catalyzed by the enzyme rabbit 15-lipoxygenase-1.

Conformational Heterogeneity and Cooperative Effects of Mammalian ALOX15.

Arachidonic acid lipoxygenases (ALOXs) have been suggested to function as monomeric enzymes, but more recent data on rabbit ALOX15 indicated that there is a dynamic monomer-dimer equilibrium in aqueous solution. In the presence of an active site ligand (the ALOX15 inhibitor RS7) rabbit ALOX15 was crystalized as heterodimer and the X-ray coordinates of the two monomers within the dimer exhibit subtle structural differences. Using native polyacrylamide electrophoresis, we here observed that highly purified and predominantly monomeric rabbit ALOX15 and human ALOX15B are present in two conformers with distinct electrophoretic mobilities. In silico docking studies, molecular dynamics simulations, site directed mutagenesis experiments and kinetic measurements suggested that in aqueous solutions the two enzymes exhibit motional flexibility, which may impact the enzymatic properties.

Deciphering the Molecular Details of the Lipoxin Formation Mechanism in the 5(S),15(S)-DiHpETE Biosynthetic Pathway Catalyzed by Reticulocyte 15-Lipoxygenase-1.

Chronic inflammation is now widely recognized to play important roles in many commonly occurring diseases, including COVID-19. The resolution response to this chronic inflammation is an active process governed by specialized pro-resolving mediators (SPMs) like the lipid mediators known as lipoxins. The biosynthesis of lipoxins is catalyzed by several lipoxygenases (LOXs) from arachidonic acid. However, the molecular details of the mechanisms involved are not well known yet. In this paper, we have combined molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations to analyze how reticulocyte 15-LOX-1 catalyzes the production of lipoxins from 5(S),15(S)-diHpETE. Our results indicate that the dehydration mechanism from 5(S),15(S)-diHpETE, via the formation of an epoxide, presents huge energy barriers even though it was one of the two a priori synthetic proposals. This result is compatible with the fact that no epoxide has been directly detected as an intermediate in the catalytic formation of lipoxins from 5(S),15(S)-diHpETE. Conversely, the oxygenation of 5(S),15(S)-diHpETE at C14 is feasible because there is an open channel connecting the protein surface with this carbon atom, and the energy barrier for oxygen addition through this channel is small. The analysis of the following steps of this mechanism, leading to the corresponding hydroperoxide at the 15-LOX-1 active site, indicates that the oxygenation mechanism will lead to the formation of lipoxinB4 after the final action of a reductase. In contrast, our calculations are in agreement with experiments that lipoxinA4 cannot derive from 5(S),15(S)-diHpETE by either of the two proposed mechanisms and that 5(S),15(S)-diHETE is not an intermediate of lipoxin biosynthesis catalyzed by 15-LOX-1.

A protocol to obtain multidimensional quantum tunneling corrections derived from QM(DFT)/MM calculations for an enzyme reaction.

The significance of tunneling contributions to the rate constants of enzymatic reactions has been described at length using experimental evidence as well as theoretical computations. Within the frame of variational transition state theory (VTST), tunneling corrections are included using the so-called ground-state tunneling transmission coefficient. For the calculation of those enzymatic rate constants using the ensemble-averaged extension of VTST on QM/MM potential energy surfaces, the transmission coefficient at a given temperature is averaged over a set of coefficient values, each one obtained from an individual minimum energy path (MEP). However, the calculation of accurate QM/MM MEPs for tunneling calculations, also using a reliable QM method like DFT, is highly costly in enzyme models. For this reason, more affordable methodologies have been used. In this paper, we validate a feasible computational strategy to compute multidimensional tunneling corrections that describes better than cheaper alternatives the physics of the hydrogen abstraction from linoleic acid catalyzed by the enzyme 15-rLOX-1. Our recommendations to obtain better values of kinetic isotope effects and, especially, of rate constants are based on multidimensional small-curvature tunneling (SCT) coefficients derived from electrostatic embedding QM(DFT)/MM MEPs. The MEPs used must be calculated with a small enough step-size. Also, the number of gradients and Hessians along the reaction path must be checked to cover the whole tunneling region and to obtain converged adiabatic potential energy profiles. Distinguished reaction coordinates (DCPs) that are commonly used to describe enzyme reaction mechanisms are not adequate for tunneling calculations in such biological systems.

A role of Gln596 in fine-tuning mammalian ALOX15 specificity, protein stability and allosteric properties.

His596 of human ALOX12 has been suggested to interact with the COO--group of arachidonic acid during ALOX catalysis. In mammalian ALOX15 orthologs Gln596 occupies this position and this amino acid exchange might contribute to the functional differences between the two ALOX-isoforms. To explore the role of Gln596 for ALOX15 functionality we mutated this amino acid to different residues in rabbit and human ALOX15 and investigated the impact of these mutations on structural, catalytic and allosteric enzyme properties. To shed light on the molecular basis of the observed functional alterations we performed in silico substrate docking studies and molecular dynamics simulations and also explored the impact of Gln596 exchange on the protein structure. The combined theoretical and experimental data suggest that Gln596 may not directly interact with the COO--group of arachidonic acid. In contrast, mutations at Gln596 destabilize the secondary and tertiary structure of ALOX15 orthologs, which may be related to a disturbance of the electrostatic interaction network with other amino acids in the immediate surrounding. Moreover, our MD-simulations suggest that the geometry of the dimer interface depends on the structure of substrate bound inside the substrate-binding pocket and that Gln596Ala exchange impairs the allosteric properties of the enzyme. Taken together, these data indicate the structural and functional importance of Gln596 for ALOX15 catalysis.

Unraveling how the Gly526Ser mutation arrests prostaglandin formation from arachidonic acid catalyzed by cyclooxygenase-2: a combined molecular dynamics and QM/MM study.

Cyclooxygenases (COXs) are the enzymes responsible for the biosynthesis of prostaglandins, eicosanoids that play a major role in many physiological processes. Particularly, prostaglandins are known to trigger inflammation, and COX-2, the enzyme isoform associated with this inflammatory response, catalyzes the cyclooxidation of arachidonic acid, leading to prostaglandin G2. For this reason, COX-2 has been a very important pharmacological target for several decades now. The catalytic mechanism of COX-2, a so-called all-radical mechanism, consists of six chemical steps. One of the most intriguing aspects of this mechanism is how COX-2 manages to control the regio- and stereospecificity of the products formed at each step. Mutagenesis experiments have previously been performed in an attempt to find those hot-spot residues that make such control possible. In this context, it is worth mentioning that in experiments with the Gly526Ser COX-2 mutant, prostaglandins were not detected. In this paper, we have combined molecular dynamics simulations and quantum mechanics/molecular mechanics calculations to analyze how the COX-2 catalytic mechanism is modified in the Gly526Ser mutant. Therefore, this study provides new insights into the COX-2 catalytic function.

Mutations of Triad Determinants Changes the Substrate Alignment at the Catalytic Center of Human ALOX5.

For the specificity of ALOX15 orthologs of different mammals, the geometry of the amino acids Phe353, Ile418, Met419, and Ile593 ("triad determinants") is important, and mutagenesis of these residues altered the reaction specificity of these enzymes. Here we expressed wild-type human ALOX5 and its F359W/A424I/N425M/A603I mutant in Sf9 insect cells and characterized the catalytic differences of the two enzyme variants. We found that wild-type ALOX5 converted arachidonic acid mainly to 5(S)-hydroperoxyeicosatetraenoic acid (HpETE). In contrast, 15(S)- and 8(S)-H(p)ETE were formed by the mutant enzyme. In addition to arachidonic acid, wild-type ALOX5 accepted eicosapentaenoic acid (EPA) as substrate, but C18 fatty acids were not oxygenated. The quadruple mutant also accepted linoleic acid and α- and γ-linolenic acid as substrate. Structural analysis of the oxygenation products and kinetic studies with stereospecifically labeled 11(S)- and 11(R)-deutero-linoleic acid suggested alternative ways of substrate orientation at the active site. In silico docking studies, molecular dynamics simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations confirmed this hypothesis. These data indicate that "triad determinant" mutagenesis alters the catalytic properties of ALOX5 abolishing its leukotriene synthase activity but improving its biosynthetic capacity for pro-resolving lipoxins.

Deciphering the grounds of the suitability of acylhydrazones as efficient photoswitches.

Many efforts are currently being devoted to designing molecular photoswitches with specific properties. In this sense, a recent publication [D. J. van Dijken et al., J. Am. Chem. Soc., 2015, 137, 14982-14991] has synthesized and analyzed the photochromic properties of a large set of acylhydrazones (ACHs), a relatively unexploited class of potential photoswitches with two stable E and Z isomers. This study has revealed a very diverse and complex pattern of the absorption/emission properties of ACHs depending on the substituents attached to the ACH motif. In this work, high level theoretical calculations are performed on a representative set of the experimentally studied ACHs in order to analyze, at the molecular level, the reasons behind the different photochemistries experimentally observed. This systematic study allows for the classification of the full set of ACHs into just four categories. The two more common groups display a small, either positive or negative, shift of the maximum wavelength of absorption between the E and Z isomers. Less common, but far more interesting from a practical point of view, are the compounds that show a large (>100 nm) Stokes shift. This behavior may arise from two different situations. The most common one implies the possibility of an intramolecular proton transfer in the excited electronic state of the less stable Z isomer. The less likely scenario would also involve a loss of the azidic proton through an intermolecular proton transfer that would take place with the aid of the solvent.

Synthetic Photoswitchable Neurotransmitters Based on Bridged Azobenzenes.

Photoswitchable neurotransmitters of ionotropic kainate receptors were synthesized by tethering a glutamate moiety to disubstituted C2-bridged azobenzenes, which were prepared through a novel methodology that allows access to diazocines with higher yields and versatility. Because of the singular properties of these photochromes, photoisomerizable compounds were obtained with larger thermal stability for their inert cis isomer than for their biologically activity trans state. This enabled selective neuronal firing upon irradiation without background activity in the dark.

Discovery of processive catalysis by an exo-hydrolase with a pocket-shaped active site.

Substrates associate and products dissociate from enzyme catalytic sites rapidly, which hampers investigations of their trajectories. The high-resolution structure of the native Hordeum exo-hydrolase HvExoI isolated from seedlings reveals that non-covalently trapped glucose forms a stable enzyme-product complex. Here, we report that the alkyl ß-D-glucoside and methyl 6-thio-ß-gentiobioside substrate analogues perfused in crystalline HvExoI bind across the catalytic site after they displace glucose, while methyl 2-thio-ß-sophoroside attaches nearby. Structural analyses and multi-scale molecular modelling of nanoscale reactant movements in HvExoI reveal that upon productive binding of incoming substrates, the glucose product modifies its binding patterns and evokes the formation of a transient lateral cavity, which serves as a conduit for glucose departure to allow for the next catalytic round. This path enables substrate-product assisted processive catalysis through multiple hydrolytic events without HvExoI losing contact with oligo- or polymeric substrates. We anticipate that such enzyme plasticity could be prevalent among exo-hydrolases.