Protonation states and catalysis: Molecular dynamics studies of intermediates in tryptophan synthase

Neutron diffraction from a microgravity-grown crystal reveals the active site hydrogens of the internal aldimine form of tryptophan synthase

Pyridoxal 5'-phosphate (PLP), the biologically active form of vitamin B6, is an essential cofactor in many biosynthetic pathways. The emergence of PLP-dependent enzymes as drug targets and biocatalysts, such as tryptophan synthase (TS), has underlined the demand to understand PLP-dependent catalysis and reaction specificity. The ability of neutron diffraction to resolve the positions of hydrogen atoms makes it an ideal technique to understand how the electrostatic environment and selective protonation of PLP regulates PLP-dependent activities. Facilitated by microgravity crystallization of TS with the Toledo Crystallization Box, we report the 2.1 Å joint X-ray/neutron (XN) structure of TS with PLP in the internal aldimine form. Positions of hydrogens were directly determined in both the α- and β-active sites, including PLP cofactor. The joint XN structure thus provides insight into the selective protonation of the internal aldimine and the electrostatic environment of TS necessary to understand the overall catalytic mechanism.

Molecular Dynamic Simulations of Aqueous Micellar Organometallic Catalysis: Methane Functionalization as a Case Study

Molecular Dynamics (MD) simulations constitute a powerful tool that provides a 3D perspective of the dynamical behavior of chemical systems. Herein the first MD study of the dynamics of a catalytic organometallic system, in micellar media, is presented. The challenging methane catalytic functionalization into ethyl propionate through a silver-catalyzed process has been targeted as the case study. The intimate nature of the micelles formed with the surfactants sodium dodecylsulfate (SDS) and potassium perfluorooctane sulfonate (PFOS) has been ascertained, as well as the relative distribution of the main actors in this transformation, namely methane, the diazo reagent and the silver catalyst, the latter in two different forms: the initial compound and a silver-carbene intermediate. Catalyst deactivation occurs with halide containing surfactants dodecyltrimethylammonium chloride (DTAC) and Triton X-100. Computed simulations allow explaining the experimental results, indicating that micelles behave differently regarding the degree of accumulation and the local distribution of the reactants and their effect in the molecular collisions leading to net reaction.

A Molecular Dynamics Simulation Study of the Effects of βGln114 Mutation on the Dynamic Behavior of the Catalytic Site of the Tryptophan Synthase

L-tryptophan (l-Trp), a vital amino acid for the survival of various organisms, is synthesized by the enzyme tryptophan synthase (TS) in organisms such as eubacteria, archaebacteria, protista, fungi, and plantae. TS, a pyridoxal 5'-phosphate (PLP)-dependent enzyme, comprises α and β subunits that typically form an α2β2 tetramer. The enzyme's activity is regulated by the conformational switching of its α and β subunits between the open (T state) and closed (R state) conformations. Many microorganisms rely on TS for growth and replication, making the enzyme and the l-Trp biosynthetic pathway potential drug targets. For instance, Mycobacterium tuberculosis, Chlamydiae bacteria, Streptococcus pneumoniae, Francisella tularensis, Salmonella bacteria, and Cryptosporidium parasitic protozoa depend on l-Trp synthesis. Antibiotic-resistant salmonella strains have emerged, underscoring the need for novel drugs targeting the l-Trp biosynthetic pathway, especially for salmonella-related infections. A single amino acid mutation can significantly impact enzyme function, affecting stability, conformational dynamics, and active or allosteric sites. These changes influence interactions, catalytic activity, and protein-ligand/protein-protein interactions. This study focuses on the impact of mutating the βGln114 residue on the catalytic and allosteric sites of TS. Extensive molecular dynamics simulations were conducted on E(PLP), E(AEX1), E(A-A), and E(C3) forms of TS using the WT, βQ114A, and βQ114N versions. The results show that both the βQ114A and βQ114N mutations increase protein backbone root mean square deviation fluctuations, destabilizing all TS forms. Conformational and hydrogen bond analyses suggest the significance of βGln114 drifting away from cofactor/intermediates and forming hydrogen bonds with water molecules necessary for l-Trp biosynthesis. The βQ114A mutation creates a gap between βAla114 and cofactor/intermediates, hindering hydrogen bond formation due to short side chains and disrupting β-sites. Conversely, the βQ114N mutation positions βAsn114 closer to cofactor/intermediates, forming hydrogen bonds with O3 of cofactors/intermediates and nearby water molecules, potentially disrupting the l-Trp biosynthetic mechanism.

Improving the Activity of Tryptophan Synthetase via a Nucleic Acid Scaffold

Tryptophan synthetase (TSase), which functions as a tetramer, is a typical enzyme with a substrate channel effect, and shows excellent performance in the production of non-standard amino acids, histamine, and other biological derivatives. Based on previous work, we fused a mutant CE protein (colistin of E. coli, a polypeptide with antibacterial activity) sequence with the sequence of TSase to explore whether its catalytic activity could be enhanced, and we also analyzed whether the addition of a DNA scaffold was a feasible strategy. Here, dCE (CE protein without DNase activity) protein tags were constructed and fused to the TrapA and TrapB subunits of TSase, and the whole cell was used for the catalytic reaction. The results showed that after the dCE protein tag was fused to the TrapB subunit, its whole cell catalytic activity increased by 50%. Next, the two subunits were expressed separately, and the proteins were bound in vitro to ensure equimolar combination between the two subunits. After the dCE label was fused to TrapB, the activity of TSase assembled with TrapA also improved. A series of experiments revealed that the enzyme fused with dCE9 showed higher activity than the wild-type protein. In general, the activity of assembly TSase was optimal when the temperature was 50 °C and the pH was about 9.0. After a long temperature treatment, the enzyme maintained good activity. With the addition of exogenous nucleic acid, the activity of the enzyme increased. The maximum yield was 0.58 g/L, which was almost three times that of the wild-type TSase (0.21 g/L). The recombinant TSase constructed in this study with dCE fusion had the advantages of higher heat resistance and higher activity, and confirmed the feasibility of adding a nucleic acid scaffold, providing a new idea for the improvement of structurally similar enzymes.

Information-Rich, Dual-Function 13C/2H-Isotopic Crosstalk NMR Assay for Human Serine Racemase (hSR) Provides a PLP-Enzyme "Partitioning Fingerprint" and Reveals Disparate Chemotypes for hSR Inhibition

The first dual-function assay for human serine racemase (hSR), the only bona fide racemase in human biology, is reported. The hSR racemization function is essential for neuronal signaling, as the product, d-serine (d-Ser), is a potent N-methyl d-aspartate (NMDA) coagonist, important for learning and memory, with dysfunctional d-Ser-signaling being observed in some neuronal disorders. The second hSR function is β-elimination and gives pyruvate; this activity is elevated in colorectal cancer. This new NMR-based assay allows one to monitor both α-proton-exchange chemistry and β-elimination using only the native l-Ser substrate and hSR and is the most sensitive such assay. The assay judiciously employs segregated dual 13C-labeling and 13C/2H crosstalk, exploiting both the splitting and shielding effects of deuterium. The assay is deployed to screen a 1020-compound library and identifies an indolo-chroman-2,4-dione inhibitor family that displays allosteric site binding behavior (noncompetitive inhibition vs l-Ser substrate; competitive inhibition vs adenosine 5'-triphosphate (ATP)). This assay also reveals important mechanistic information for hSR; namely, that H/D exchange is ∼13-fold faster than racemization, implying that K56 protonates the carbanionic intermediate on the si-face much faster than does S84 on the re-face. Moreover, the 13C NMR peak pattern seen is suggestive of internal return, pointing to K56 as the likely enamine-protonating residue for β-elimination. The 13C/2H-isotopic crosstalk assay has also been applied to the enzyme tryptophan synthase and reveals a dramatically different partition ratio in this active site (β-replacement: si-face protonation ∼6:1 vs β-elimination: si-face protonation ∼1:3.6 for hSR), highlighting the value of this approach for fingerprinting the pyridoxal phosphate (PLP) enzyme mechanism.

Cysteine Oxidation in Proteins: Structure, Biophysics, and Simulation

Cysteine side chains can exist in distinct oxidation states depending on the pH and redox potential of the environment, and cysteine oxidation plays important yet complex regulatory roles. Compared with the effects of post-translational modifications such as phosphorylation, the effects of oxidation of cysteine to sulfenic, sulfinic, and sulfonic acid on protein structure and function remain relatively poorly characterized. We present an analysis of the role of cysteine reactivity as a regulatory factor in proteins, emphasizing the interplay between electrostatics and redox potential as key determinants of the resulting oxidation state. A review of current computational approaches suggests underdeveloped areas of research for studying cysteine reactivity through molecular simulations.

Allosteric regulation of substrate channeling: Salmonella typhimurium tryptophan synthase

The regulation of the synthesis of L-tryptophan (L-Trp) in enteric bacteria begins at the level of gene expression where the cellular concentration of L-Trp tightly controls expression of the five enzymes of the Trp operon responsible for the synthesis of L-Trp. Two of these enzymes, trpA and trpB, form an αββα bienzyme complex, designated as tryptophan synthase (TS). TS carries out the last two enzymatic processes comprising the synthesis of L-Trp. The TS α-subunits catalyze the cleavage of 3-indole D-glyceraldehyde 3′-phosphate to indole and D-glyceraldehyde 3-phosphate; the pyridoxal phosphate-requiring β-subunits catalyze a nine-step reaction sequence to replace the L-Ser hydroxyl by indole giving L-Trp and a water molecule. Within αβ dimeric units of the αββα bienzyme complex, the common intermediate indole is channeled from the α site to the β site via an interconnecting 25 Å-long tunnel. The TS system provides an unusual example of allosteric control wherein the structures of the nine different covalent intermediates along the β-reaction catalytic path and substrate binding to the α-site provide the allosteric triggers for switching the αββα system between the open (T) and closed (R) allosteric states. This triggering provides a linkage that couples the allosteric conformational coordinate to the covalent chemical reaction coordinates at the α- and β-sites. This coupling drives the α- and β-sites between T and R conformations to achieve regulation of substrate binding and/or product release, modulation of the α- and β-site catalytic activities, prevention of indole escape from the confines of the active sites and the interconnecting tunnel, and synchronization of the α- and β-site catalytic activities. Here we review recent advances in the understanding of the relationships between structure, function, and allosteric regulation of the complex found in Salmonella typhimurium.

Atomic-resolution chemical characterization of (2x)72-kDa tryptophan synthase via four- and five-dimensional 1H-detected solid-state NMR

NMR chemical shifts provide detailed information on the chemical properties of molecules, thereby complementing structural data from techniques like X-ray crystallography and electron microscopy. Detailed analysis of protein NMR data, however, often hinges on comprehensive, site-specific assignment of backbone resonances, which becomes a bottleneck for molecular weights beyond 40 to 45 kDa. Here, we show that assignments for the (2x)72-kDa protein tryptophan synthase (665 amino acids per asymmetric unit) can be achieved via higher-dimensional, proton-detected, solid-state NMR using a single, 1-mg, uniformly labeled, microcrystalline sample. This framework grants access to atom-specific characterization of chemical properties and relaxation for the backbone and side chains, including those residues important for the catalytic turnover. Combined with first-principles calculations, the chemical shifts in the β-subunit active site suggest a connection between active-site chemistry, the electrostatic environment, and catalytically important dynamics of the portal to the β-subunit from solution.

Mutation of βGln114 to Ala Alters the Stabilities of Allosteric States in Tryptophan Synthase Catalysis

The tryptophan synthase (TS) bienzyme complexes found in bacteria, yeasts, and molds are pyridoxal 5'-phosphate (PLP)-requiring enzymes that synthesize l-Trp. In the TS catalytic cycle, switching between the open and closed states of the α- and β-subunits via allosteric interactions is key to the efficient conversion of 3-indole-d-glycerol-3'-phosphate and l-Ser to l-Trp. In this process, the roles played by β-site residues proximal to the PLP cofactor have not yet been fully established. βGln114 is one such residue. To explore the roles played by βQ114, we conducted a detailed investigation of the βQ114A mutation on the structure and function of tryptophan synthase. Initial steady-state kinetic and static ultraviolet-visible spectroscopic analyses showed the Q to A mutation impairs catalytic activity and alters the stabilities of intermediates in the β-reaction. Therefore, we conducted X-ray structural and solid-state nuclear magnetic resonance spectroscopic studies to compare the wild-type and βQ114A mutant enzymes. These comparisons establish that the protein structural changes are limited to the Gln to Ala replacement, the loss of hydrogen bonds among the side chains of βGln114, βAsn145, and βArg148, and the inclusion of waters in the cavity created by substitution of the smaller Ala side chain. Because the conformations of the open and closed allosteric states are not changed by the mutation, we hypothesize that the altered properties arise from the lost hydrogen bonds that alter the relative stabilities of the open (β^T state) and closed (β^R state) conformations of the β-subunit and consequently alter the distribution of intermediates along the β-subunit catalytic path.

Theoretical Studies of the Acid-Base Equilibria in a Model Active Site of the Human 20S Proteasome

The 20S proteasome is a macromolecule responsible for the chemical step in the ubiquitin-proteasome system of degrading unnecessary and unused proteins of the cell. It plays a central role both in the rapid growth of cancer cells and in viral infection cycles. Herein, we present a computational study of the acid-base equilibria in an active site of the human proteasome (caspase-like), an aspect which is often neglected despite the crucial role protons play in the catalysis. As example substrates, we take the inhibition by epoxy- and boronic acid-containing warheads. We have combined cluster quantum mechanical calculations, replica exchange molecular dynamics, and Bayesian optimization of nonbonded potential terms in the inhibitors. In relation to the latter, we propose an easily scalable approach for the reevaluation of nonbonded potentials making use of the hybrid quantum mechanics molecular mechanics dynamics information. Our results show that coupled acid-base equilibria need to be considered when modeling the inhibition mechanism. The coupling between a neighboring lysine and the reacting threonine is not affected by the presence of the studied inhibitors.

Computational investigations of allostery in aromatic amino acid biosynthetic enzymes

Allostery, in which binding of ligands to remote sites causes a functional change in the active sites, is a fascinating phenomenon observed in enzymes. Allostery can occur either with or without significant conformational changes in the enzymes, and the molecular basis of its mechanism can be difficult to decipher using only experimental techniques. Computational tools for analyzing enzyme sequences, structures, and dynamics can provide insights into the allosteric mechanism at the atomic level. Combining computational and experimental methods offers a powerful strategy for the study of enzyme allostery. The aromatic amino acid biosynthesis pathway is essential in microorganisms and plants. Multiple enzymes involved in this pathway are sensitive to feedback regulation by pathway end products and are known to use allostery to control their activities. To date, four enzymes in the aromatic amino acid biosynthesis pathway have been computationally investigated for their allosteric mechanisms, including 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase, anthranilate synthase, chorismate mutase, and tryptophan synthase. Here we review the computational studies and findings on the allosteric mechanisms of these four enzymes. Results from these studies demonstrate the capability of computational tools and encourage future computational investigations of allostery in other enzymes of this pathway.

Structural Basis of the Stereochemistry of Inhibition of Tryptophan Synthase by Tryptophan and Derivatives

We have examined the reaction of Salmonella enterica serovar typhimurium tryptophan (Trp) synthase α₂β₂ complex with l-Trp, d-Trp, oxindolyl-l-alanine (OIA), and dioxindolyl-l-alanine (DOA) in the presence of disodium (dl)-α-glycerol phosphate (GP), using stopped-flow spectrophotometry and X-ray crystallography. All structures contained the d-isomer of GP bound at the α-active site. (3S)-OIA reacts with the pyridoxal-5'-phosphate (PLP) of Trp synthase to form a mixture of external aldimine and quinonoid complexes. The α-carboxylate of OIA rotates about 90° to become planar with the PLP when the quinonoid complex is formed, resulting in a conformational change in the loop of residues 110-115. The COMM domain of the Trp synthase-OIA complex is found as a mixture of two conformations. The (3R)-diastereomer of DOA binds about 5-fold more tightly than (3S)-OIA and also forms a mixture of aldimine and quinonoid complexes. DOA forms an additional H-bond between the 3-OH of DOA and βLys-87. l-Trp does not form a covalent complex with the PLP of Trp synthase. However, d-Trp forms a mixture of two external aldimine complexes which differ in the orientation of the α-carboxylate. In one conformation, the α-carboxylate is in the plane of the PLP, while in the other conformation, the α-carboxylate is perpendicular to the PLP plane. These results confirm that the stereochemistry of the transient indolenine quinonoid intermediate in the mechanism of Trp synthase is (3S) and demonstrate the linkage between aldimine and quinonoid reaction intermediates in the β-active site and allosteric communications with the α-active site.

Engineering Improves Enzymatic Synthesis of L-Tryptophan by Tryptophan Synthase from Escherichia coli

To improve the thermostability of tryptophan synthase, the molecular modification of tryptophan synthase was carried out by rational molecular engineering. First, B-FITTER software was used to analyze the temperature factor (B-factor) of each amino acid residue in the crystal structure of tryptophan synthase. A key amino acid residue, G395, which adversely affected the thermal stability of the enzyme, was identified, and then, a mutant library was constructed by site-specific saturation mutation. A mutant (G395S) enzyme with significantly improved thermal stability was screened from the saturated mutant library. Error-prone PCR was used to conduct a directed evolution of the mutant enzyme (G395S). Compared with the parent, the mutant enzyme (G395S /A191T) had a K_m of 0.21 mM and a catalytic efficiency k_cat/K_m of 5.38 mM⁻¹∙s⁻¹, which was 4.8 times higher than that of the wild-type strain. The conditions for L-tryptophan synthesis by the mutated enzyme were a L-serine concentration of 50 mmol/L, a reaction temperature of 40 °C, pH of 8, a reaction time of 12 h, and an L-tryptophan yield of 81%. The thermal stability of the enzyme can be improved by using an appropriate rational design strategy to modify the correct site. The catalytic activity of tryptophan synthase was increased by directed evolution.

Direct dynamic nuclear polarization of 15N and 13C spins at 14.1 T using a trityl radical and magic angle spinning

We investigate solid-state dynamic nuclear polarization of 13C and 15N nuclei using monoradical trityl OX063 as a polarizing agent in a magnetic field of 14.1 T with magic angle spinning at ∼100 K. We monitored the field dependence of direct 13C and 15N polarization for frozen [13C, 15N] urea and achieved maximum absolute enhancement factors of 240 and 470, respectively. The field profiles are consistent with polarization of 15N spins via either the solid effect or the cross effect, and polarization of 13C spins via a combination of cross effect and solid effect. For microcrystalline, 15N-enriched tryptophan synthase sample containing trityl radical, a 1500-fold increase in 15N signal was observed under microwave irradiation. These results show the promise of trityl radicals and their derivatives for direct polarization of low gamma, spin-½ nuclei at high magnetic fields and suggest a novel approach for selectively polarizing specific moieties or for polarizing systems which have low levels of protonation.

Human Serine Racemase: Key Residues/Active Site Motifs and Their Relation to Enzyme Function

Serine racemase (SR) is the first racemase enzyme to be identified in human biology and converts L-serine to D-serine, an important neuronal signaling molecule that serves as a co-agonist of the NMDA (N-methyl-D-aspartate) receptor. This overview describes key molecular features of the enzyme, focusing on the side chains and binding motifs that control PLP (pyridoxal phosphate) cofactor binding as well as activity modulation through the binding of both divalent cations and ATP, the latter showing allosteric modulation. Discussed are catalytically important residues in the active site including K56 and S84—the si- and re-face bases, respectively,—and R135, a residue that appears to play a critical role in the binding of both negatively charged alternative substrates and inhibitors. The interesting bifurcated mechanism followed by this enzyme whereby substrate L-serine can be channeled either into D-serine (racemization pathway) or into pyruvate (β-elimination pathway) is discussed extensively, as are studies that focus on a key loop region (the so-called “triple serine loop”), the modification of which can be used to invert the normal in vitro preference of this enzyme for the latter pathway over the former. The possible cross-talk between the PLP enzymes hSR and hCBS (human cystathionine β-synthase) is discussed, as the former produces D-serine and the latter produces H₂S, both of which stimulate the NMDAR and both of which have been implicated in neuronal infarction pursuant to ischemic stroke. Efforts to gain a more complete mechanistic understanding of these PLP enzymes are expected to provide valuable insights for the development of specific small molecule modulators of these enzymes as tools to study their roles in neuronal signaling and in modulation of NMDAR function.

Investigation of Structural Dynamics of Enzymes and Protonation States of Substrates Using Computational Tools

This review discusses the use of molecular modeling tools, together with existing experimental findings, to provide a complete atomic-level description of enzyme dynamics and function. We focus on functionally relevant conformational dynamics of enzymes and the protonation states of substrates. The conformational fluctuations of enzymes usually play a crucial role in substrate recognition and catalysis. Protein dynamics can be altered by a tiny change in a molecular system such as different protonation states of various intermediates or by a significant perturbation such as a ligand association. Here we review recent advances in applying atomistic molecular dynamics (MD) simulations to investigate allosteric and network regulation of tryptophan synthase (TRPS) and protonation states of its intermediates and catalysis. In addition, we review studies using quantum mechanics/molecular mechanics (QM/MM) methods to investigate the protonation states of catalytic residues of β-Ketoacyl ACP synthase I (KasA). We also discuss modeling of large-scale protein motions for HIV-1 protease with coarse-grained Brownian dynamics (BD) simulations.

Visualizing the tunnel in tryptophan synthase with crystallography: Insights into a selective filter for accommodating indole and rejecting water

Four new X-ray structures of tryptophan synthase (TS) crystallized with varying numbers of the amphipathic N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F6) molecule are presented. These structures show one of the F6 ligands threaded into the tunnel from the β-site and reveal a distinct hydrophobic region. Over this expanse, the interactions between F6 and the tunnel are primarily nonpolar, while the F6 phosphoryl group fits into a polar pocket of the β-subunit active site. Further examination of TS structures reveals that one portion of the tunnel (T1) binds clusters of water molecules, whereas waters are not observed in the nonpolar F6 binding region of the tunnel (T2). MD simulation of another TS structure with an unobstructed tunnel also indicates the T2 region of the tunnel excludes water, consistent with a dewetted state that presents a significant barrier to the transfer of water into the closed β-site. We conclude that hydrophobic molecules can freely diffuse between the α- and β-sites via the tunnel, while water does not. We propose that exclusion of water serves to inhibit reaction of water with the α-aminoacrylate intermediate to form ammonium ion and pyruvate, a deleterious side reaction in the αβ-catalytic cycle. Finally, while most TS structures show βPhe280 partially blocking the tunnel between the α- and β-sites, new structures show an open tunnel, suggesting the flexibility of the βPhe280 side chain. Flexible docking studies and MD simulations confirm that the dynamic behavior of βPhe280 allows unhindered transfer of indole through the tunnel, therefore excluding a gating role for this residue.