Unravelling the Graded Millisecond Allosteric Activation Mechanism of Imidazole Glycerol Phosphate Synthase

Deciphering the molecular mechanisms of enzymatic allosteric regulation requires the structural characterization of key functional states and also their time evolution toward the formation of the allosterically activated ternary complex. The transient nature and usually slow millisecond timescale interconversion between these functional states hamper their detailed experimental and computational characterization. Here, we design a computational strategy tailored to reconstruct millisecond timescale events to describe the graded allosteric activation of imidazole glycerol phosphate synthase (IGPS) in the ternary complex. IGPS is a heterodimeric bienzyme complex responsible for the hydrolysis of glutamine to glutamate in the HisH subunit and delivering ammonia for the cyclase activity in HisF. Despite significant advances in understanding the underlying allosteric mechanism, essential molecular details of the long-range millisecond allosteric activation pathway of wild-type IGPS remain hidden. Without using a priori information of the active state, our simulations uncover how IGPS, with the allosteric effector bound in HisF, spontaneously captures glutamine in a catalytically inactive HisH conformation, subsequently attains a closed HisF:HisH interface, and finally forms the oxyanion hole in HisH for efficient glutamine hydrolysis. We show that effector binding in HisF dramatically decreases the conformational barrier associated with the oxyanion hole formation in HisH, in line with the experimentally observed 4500-fold activity increase in glutamine production. The formation of the allosterically active state is controlled by time-evolving dynamic communication networks connecting the effector and substrate binding sites. This computational strategy can be generalized to study other unrelated enzymes undergoing millisecond timescale allosteric transitions.

Ensemble–function relationships to evaluate catalysis in the ketosteroid isomerase oxyanion hole

Following decades of insights from structure–function studies, there is now a need to progress from a static to dynamic view of enzymes. Comparison of prior cryo X-ray structures suggested that deleterious effects from ketosteroid isomerase (KSI) mutants arise from misalignment of the oxyanion hole catalytic residue, Y16. However, multi-conformer models from room temperature X-ray diffraction revealed an ensemble of Y16 conformers indistinguishable from WT for Y32F/Y57F KSI and a distinct, non-native ensemble for Y16 in Y57F KSI. Functional analyses suggested rate effects arise from weakened hydrogen bonding, due to disruption of the Y16/Y57/Y32 hydrogen bond network, and repositioning of the general base. In general, catalytic changes can be deconvoluted into effects on the probability of occupying a state (P-effects) and the reactivity of each state (k-effects). Our results underscore the need for ensemble–function analysis to decipher enzyme function and ultimately manipulate their extraordinary capabilities.

Structure of a reaction intermediate mimic in t6A biosynthesis bound in the active site of the TsaBD heterodimer from Escherichia coli

The tRNA modification N6-threonylcarbamoyladenosine (t6A) is universally conserved in all organisms. In bacteria, the biosynthesis of t6A requires four proteins (TsaBCDE) that catalyze the formation of t6A via the unstable intermediate l-threonylcarbamoyl-adenylate (TC-AMP). While the formation and stability of this intermediate has been studied in detail, the mechanism of its transfer to A37 in tRNA is poorly understood. To investigate this step, the structure of the TsaBD heterodimer from Escherichia coli has been solved bound to a stable phosphonate isosteric mimic of TC-AMP. The phosphonate inhibits t6A synthesis in vitro with an IC50 value of 1.3 μM in the presence of millimolar ATP and L-threonine. The inhibitor binds to TsaBD by coordination to the active site Zn atom via an oxygen atom from both the phosphonate and the carboxylate moieties. The bound conformation of the inhibitor suggests that the catalysis exploits a putative oxyanion hole created by a conserved active site loop of TsaD and that the metal essentially serves as a binding scaffold for the intermediate. The phosphonate bound crystal structure should be useful for the rational design of potent, drug-like small molecule inhibitors as mechanistic probes or potentially novel antibiotics.

Application of the quantum theory of atoms in molecules (QTAIM) to the study of the enzymatic kinetic resolution of propranolol, an amino alcohol with pharmaceutical applications

Propranolol, ((R,S)-1-iso-propylamino-3-(1-naphthoxy)-2-propanol), is a β-adrenergic antagonist and is commercially available as a racemic mixture. Only the S-enantiomer has the desired therapeutic effect. Therefore, many researchers have been working on strategies to obtain S-propranolol with high enantiomeric purity. One approach to carry out the acetylation of (R,S)-Propranolol using Candida antarctica lipase B, CalB. This reaction leads to an enantiomeric purity of 96% at a relatively low conversion rate of 30 %. In our research group, we have been studying this reaction. The CalB active site is composed by the triad catalytic (ASP 187, HIS 224 and SER 105) and oxyanion hole (GLN 106 and THR 40). In a previous work, a QM/MM (Quantum Mechanics / Molecular Mechanics) study was carried out, using a QM region consisting only of the catalytic triad of CalB and (R,S)-propranolol [1]. In the present study, we investigate the effect of expanding the quantum region to include the oxyanion hole and to comprehend the effect of intermolecular hydrogen bonds present between the (R,S)-propranolol and the CalB active site. The electronic structure was analyzed using the Quantum Theory of Atoms In Molecules, QTAIM. Our results show that: 1. the studied reactions are more exothermic with the inclusion of the oxyanion hole than with only the catalytic triad. 2. the intermolecular interactions between (R,S)-propranolol and the CalB active site are dominated by hydrogen bonds (HB). Among those HBs, only one between propranolol and HIS 224, and another one between THR 40 and the carbonyl oxygen of acetylated SER 105 play an important role.

Probing the dynamic structure-function and structure-free energy relationships of the corona virus main protease with Biodynamics theory

The SARS-CoV-2 Main protease (Mpro) is of major interest as an anti-viral drug target. Structure-based virtual screening efforts, fueled by a growing list of apo and inhibitor-bound SARS-CoV/CoV-2 Mpro crystal structures, are underway in many labs. However, little is known about the dynamic enzyme mechanism, which is needed to inform both structure-based design and assay development. Here, we apply Biodynamics theory to characterize the structural dynamics of substrate-induced Mpro activation, and explore the implications thereof for efficacious inhibition under non-equilibrium conditions. The catalytic cycle (including tetrahedral intermediate formation and hydrolysis) is governed by concerted dynamic structural rearrangements of domain 3 and the m-shaped loop (residues 132-147) on which Cys145 (comprising the thiolate nucleophile and one-half of the oxyanion hole) and Gly143 reside (comprising the other half of the oxyanion hole). In particular: Domain 3 undergoes dynamic rigid-body rotations about the domain 2-3 linker, alternately visiting two conformational states (denoted as ).The Gly143-containing crest of the m-shaped loop (denoted as crest B) undergoes up and down translations in concert with the domain 3 rotations (denoted as , whereas the Cys145-containing crest (denoted as crest A) remains statically in the up position. The crest B translations are driven by conformational transitions within the rising leg of the loop (Lys137-Asn142).We propose that substrates associate to the state, which promotes the state, dimerization (denoted as -substrate), and catalysis. The structure resets to the dynamic monomeric form upon dissociation of the N-terminal product. We describe the energetics of the aforementioned state transitions, and address the implications of our proposed mechanism for efficacious Mpro inhibition under native-like conditions.

A multifunctional surfactant catalyst inspired by hydrolases

The remarkable power of enzymes to undertake catalysis frequently stems from their grouping of multiple, complementary chemical units within close proximity around the enzyme active site. Motivated by this, we report here a bioinspired surfactant catalyst that incorporates a variety of chemical functionalities common to hydrolytic enzymes. The textbook hydrolase active site, the catalytic triad, is modeled by positioning the three groups of the triad (-OH, -imidazole, and -CO2H) on a single, trifunctional surfactant molecule. To support this, we recreate the hydrogen bond donating arrangement of the oxyanion hole by imparting surfactant functionality to a guanidinium headgroup. Self-assembly of these amphiphiles in solution drives the collection of functional headgroups into close proximity around a hydrophobic nano-environment, affording hydrolysis of a model ester at rates that challenge α-chymotrypsin. Structural assessment via NMR and XRD, paired with MD simulation and QM calculation, reveals marked similarities of the co-micelle catalyst to native enzymes.

Role of the I16-D194 ionic interaction in the trypsin fold

Activity in trypsin-like proteases is the result of proteolytic cleavage at R15 followed by an ionic interaction that ensues between the new N terminus of I16 and the side chain of the highly conserved D194. This mechanism of activation, first proposed by Huber and Bode, organizes the oxyanion hole and primary specificity pocket for substrate binding and catalysis. Using the clotting protease thrombin as a relevant model, we unravel contributions of the I16-D194 ionic interaction to Na+ binding, stability of the transition state and the allosteric E*-E equilibrium of the trypsin fold. The I16T mutation abolishes the I16-D194 interaction and compromises the architecture of the oxyanion hole. The D194A mutation also abrogates the I16-D194 interaction but, surprisingly, has no effect on the architecture of the oxyanion hole that remains intact through a new H-bond established between G43 and G193. In both mutants, loss of the I16-D194 ionic interaction compromises Na+ binding, reduces stability of the transition state, collapses the 215–217 segment into the primary specific pocket and abrogates the allosteric E*-E equilibrium in favor of a rigid conformation that binds ligand at the active site according to a simple lock-and-key mechanism. These findings refine the structural role of the I16-D194 ionic interaction in the Huber-Bode mechanism of activation and reveal a functional linkage with the allosteric properties of the trypsin fold like Na+ binding and the E*-E equilibrium.