Synergistic Regulation and Ligand-Induced Conformational Changes of Tryptophan Synthase

Allosteric regulation of β-reaction stage I in tryptophan synthase upon the α-ligand binding

Tryptophan synthase (TRPS) is a bifunctional enzyme consisting of α- and β-subunits that catalyzes the last two steps of L-tryptophan (L-Trp) biosynthesis. The first stage of the reaction at the β-subunit is called β-reaction stage I, which converts the β-ligand from an internal aldimine [E(Ain)] to an α-aminoacrylate [E(A-A)] intermediate. The activity is known to increase 3-10-fold upon the binding of 3-indole-D-glycerol-3'-phosphate (IGP) at the α-subunit. The effect of α-ligand binding on β-reaction stage I at the distal β-active site is not well understood despite the abundant structural information available for TRPS. Here, we investigate the β-reaction stage I by carrying out minimum-energy pathway searches based on a hybrid quantum mechanics/molecular mechanics (QM/MM) model. The free-energy differences along the pathway are also examined using QM/MM umbrella sampling simulations with QM calculations at the B3LYP-D3/aug-cc-pVDZ level of theory. Our simulations suggest that the sidechain orientation of βD305 near the β-ligand likely plays an essential role in the allosteric regulation: a hydrogen bond is formed between βD305 and the β-ligand in the absence of the α-ligand, prohibiting a smooth rotation of the hydroxyl group in the quinonoid intermediate, whereas the dihedral angle rotates smoothly after the hydrogen bond is switched from βD305-β-ligand to βD305-βR141. This switch could occur upon the IGP-binding at the α-subunit, as evidenced by the existing TRPS crystal structures.

In silico evaluation of geroprotective phytochemicals as potential sirtuin 1 interactors

Sirtuin 1 (SIRT1) belongs to the histone deacetylase enzyme family and its activity regulates various signaling networks associated with aging. SIRT1 is widely involved in a large number of biological processes, including senescence, autophagy, inflammation, and oxidative stress. In addition, SIRT1 activation may improve lifespan and health in numerous experimental models. Therefore, SIRT1 targeting is a potential strategy to delay or reverse aging and age-related diseases. Although SIRT1 is activated by a wide array of small molecules, only a limited number of phytochemicals that directly interact with SIRT1 have been identified. Using the Geroprotectors.org database and a literature search, the aim of this study was to identify geroprotective phytochemicals that might interact with SIRT1. We performed molecular docking, density functional theory studies, molecular dynamic simulations (MDS), and absorption, distribution, metabolism, excretion, and toxicity (ADMET) prediction to screen potential candidates against SIRT1. After the initial screening of 70 phytochemicals, crocin, celastrol, hesperidin, taxifolin, vitexin, and quercetin had significant binding affinity scores. These six compounds established multiple hydrogen-bonding and hydrophobic interactions with SIRT1 and showed good drug-likeness and ADMET properties. In particular, crocin was further analyzed using MDS to study its complex with SIRT1 during simulation. Crocin has a high reactivity to SIRT1 and can form a stable complex with it, showing a good ability to fit into the binding pocket. Although further investigations are required, our results suggest that these geroprotective phytochemicals, especially crocin, are novel interacting partners of SIRT1.

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.

Fermentative Indole Production via Bacterial Tryptophan Synthase Alpha Subunit and Plant Indole-3-Glycerol Phosphate Lyase Enzymes

Indole is produced in nature by diverse organisms and exhibits a characteristic odor described as animal, fecal, and floral. In addition, it contributes to the flavor in foods, and it is applied in the fragrance and flavor industry. In nature, indole is synthesized either from tryptophan by bacterial tryptophanases (TNAs) or from indole-3-glycerol phosphate (IGP) by plant indole-3-glycerol phosphate lyases (IGLs). While it is widely accepted that the tryptophan synthase α-subunit (TSA) has intrinsically low IGL activity in the absence of the tryptophan synthase β-subunit, in this study, we show that Corynebacterium glutamicum TSA functions as a bona fide IGL and can support fermentative indole production in strains providing IGP. By bioprospecting additional bacterial TSAs and plant IGLs that function as bona fide IGLs were identified. Capturing indole in an overlay enabled indole production to titers of about 0.7 g L^-1 in fermentations using C. glutamicum strains expressing either the endogenous TSA gene or the IGL gene from wheat.

Computational Analysis on the Allostery of Tryptophan Synthase: Relationship between α/β-Ligand Binding and Distal Domain Closure

Tryptophan synthase (TRPS) is a bifunctional enzyme consisting of α and β-subunits and catalyzes the last two steps of l-tryptophan (L-Trp) biosynthesis, namely, cleavage of 3-indole-d-glycerol-3'-phosphate (IGP) into indole and glyceraldehyde-3-phosphate (G3P) in the α-subunit, and a pyridoxal phosphate (PLP)-dependent reaction of indole and l-serine (L-Ser) to produce L-Trp in the β-subunit. Importantly, the IGP binding at the α-subunit affects the β-subunit conformation and its ligand-binding affinity, which, in turn, enhances the enzymatic reaction at the α-subunit. The intersubunit communications in TRPS have been investigated extensively for decades because of the fundamental and pharmaceutical importance, while it is still difficult to answer how TRPS allostery is regulated at the atomic detail. Here, we investigate the allosteric regulation of TRPS by all-atom classical molecular dynamics (MD) simulations and analyze the potential of mean-force (PMF) along conformational changes of the α- and β-subunits. The present simulation has revealed a widely opened conformation of the β-subunit, which provides a pathway for L-Ser to enter into the β-active site. The IGP binding closes the α-subunit and induces a wide opening of the β-subunit, thereby enhancing the binding affinity of L-Ser to the β-subunit. Structural analyses have identified critical hydrogen bonds (HBs) at the interface of the two subunits (αG181-βS178, αP57-βR175, etc.) and HBs between the β-subunit (βT110 - βH115) and a complex of PLP and L-Ser (an α-aminoacrylate intermediate). The former HBs regulate the allosteric, β-subunit opening, whereas the latter HBs are essential for closing the β-subunit in a later step. The proposed mechanism for how the interdomain communication in TRPS is realized with ligand bindings is consistent with the previous experimental data, giving a general idea to interpret the allosteric regulations in multidomain proteins.

Discovery of antimicrobial agent targeting tryptophan synthase

Antibiotic resistance is a continually growing challenge in the treatment of various bacterial infections worldwide. New drugs and new drug targets are necessary to curb the threat of infectious diseases caused by multidrug-resistant pathogens. The tryptophan biosynthesis pathway is essential for bacterial growth but is absent in higher animals and humans. Drugs that can inhibit the bacterial biosynthesis of tryptophan offer a new class of antibiotics. In this work, we combined a structure-based strategy using in silico docking screening and molecular dynamics (MD) simulations to identify compounds targeting the α subunit of tryptophan synthase with experimental methods involving the whole-cell minimum inhibitory concentration (MIC) test, solution state NMR, and crystallography to confirm the inhibition of L-tryptophan biosynthesis. Screening 1,800 compounds from the National Cancer Institute Diversity Set I against α subunit revealed 28 compounds for experimental validation; four of the 28 hit compounds showed promising activity in MIC testing. We performed solution state NMR experiments to demonstrate that a one successful inhibitor, 3-amino-3-imino-2-phenyldiazenylpropanamide (Compound 1) binds to the α subunit. We also report a crystal structure of Salmonella enterica serotype Typhimurium tryptophan synthase in complex with Compound 1 which revealed a binding site at the αβ interface of the dimeric enzyme. MD simulations were carried out to examine two binding sites for the compound. Our results show that this small molecule inhibitor could be a promising lead for future drug development.

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.

Backbone assignments and conformational dynamics in the S. typhimurium tryptophan synthase α-subunit from solution-state NMR

Backbone assignments for the isolated α-subunit of Salmonella typhimurium tryptophan synthase (TS) are reported based on triple resonance solution-state NMR experiments on a uniformly 2H,13C,15N-labeled sample. From the backbone chemical shifts, secondary structure and random coil index order parameters (RCI-S2) are predicted. Titration with the 3-indole-D-glycerol 3'-phosphate analog, N-(4'-trifluoromethoxybenzenesulfonyl)-2-aminoethyl phosphate (F9), leads to chemical shift perturbations indicative of conformational changes from which an estimate of the dissociation constant is obtained. Comparisons of the backbone chemical-shifts, RCI-S2 values, and site-specific relaxation times with and without F9 reveal allosteric changes including modulation in secondary structures and loop rigidity induced upon ligand binding. A comparison is made to the X-ray crystal structure of the α-subunit in the full TS αββα bi-enzyme complex and to two new X-ray crystal structures of the isolated TS α-subunit reported in this work.

Pharmacophore model-based virtual screening, docking, biological evaluation and molecular dynamics simulations for inhibitors discovery against α-tryptophan synthase from Mycobacterium tuberculosis

Discovery of new anti-tuberculosis drugs with novel mode of action is urgently needed. The tryptophan synthase is a genetically validated enzyme that catalyzes last step of tryptophan biosynthetic pathway required for growth and survival of Mycobacterium tuberculosis. Here, a ligand-based pharmacophore model was built using molecular operating environment (MOE) software (version 2010.12) and validation of generated pharmacophoric features was done using active, inactive and decoy set of molecules. The generated pharmacophore model was used for screening of 7,523,972 drug-like molecules of ZINC database. The best matches (RMSD < 1) retrieved as a result of screening were subjected to molecular docking studies into active pocket of α-subunit of tryptophan synthase from M. tuberculosis. The five hits were selected and validated through anti-tuberculosis activity analysis. Finally, a new inhibitor ZINC09150898 has been identified with best binding score -32.07 kcal/mol, showing 100% growth inhibition of M. tuberculosis (H37Rv strain) at 50 µg/mL. This identified inhibitor-protein complex was further subjected to MD simulations studies (50 ns) involving root mean square deviation, root mean square fluctuation, secondary structure analysis and pocket interaction analysis to explore its binding mode stability inside active pocket. The binding free energies of inhibitor-protein complex through MM-PBSA analysis suggested that van der Waals interactions play a vital role for retention of identified inhibitor inside the protein pocket. All these analyses confirmed retention of ligand inside pocket and no unfolding in protein structure was observed over explored time scale.Communicated by Ramaswamy H. Sarma.

Evolutionary Morphing of Tryptophan Synthase: Functional Mechanisms for the Enzymatic Channeling of Indole

Tryptophan synthase (TrpS) is a heterotetrameric αββα enzyme that exhibits complex substrate channeling and allosteric mechanisms and is a model system in enzymology. In this work, we characterize proposed early and late evolutionary states of TrpS and show that they have distinct quaternary structures caused by insertions-deletions of sequence segments (indels) in the β-subunit. Remarkably, indole hydrophobic channels that connect α and β active sites have re-emerged in both TrpS types, yet they follow different paths through the β-subunit fold. Also, both TrpS geometries activate the α-subunit through the rearrangement of loops flanking the active site. Our results link evolutionary sequence changes in the enzyme subunits with channeling and allostery in the TrpS enzymes. The findings demonstrate that indels allow protein quaternary architectures to escape "minima" in the evolutionary landscape, thereby overcoming the conservational constraints imposed by existing functional interfaces and being free to morph into new mechanistic enzymes.

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.

Aqueous Molecular Dynamics Simulations of the M. tuberculosis Enoyl-ACP Reductase-NADH System and Its Complex with a Substrate Mimic or Diphenyl Ethers Inhibitors

Molecular dynamics (MD) simulations of 12 aqueous systems of the NADH-dependent enoyl-ACP reductase from Mycobacterium tuberculosis (InhA) were carried out for up to 20–40 ns using the GROMACS 4.5 package. Simulations of the holoenzyme, holoenzyme-substrate, and 10 holoenzyme-inhibitor complexes were conducted in order to gain more insight about the secondary structure motifs of the InhA substrate-binding pocket. We monitored the lifetime of the main intermolecular interactions: hydrogen bonds and hydrophobic contacts. Our MD simulations demonstrate the importance of evaluating the conformational changes that occur close to the active site of the enzyme-cofactor complex before and after binding of the ligand and the influence of the water molecules. Moreover, the protein-inhibitor total steric (E_LJ) and electrostatic (E_C) interaction energies, related to Gly96 and Tyr158, are able to explain 80% of the biological response variance according to the best linear equation, pK_i = 7.772 − 0.1885 × Gly96 + 0.0517 × Tyr158 (R² = 0.80; n = 10), where interactions with Gly96, mainly electrostatic, increase the biological response, while those with Tyr158 decrease. These results will help to understand the structure-activity relationships and to design new and more potent anti-TB drugs.

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

The importance of protonation states and proton transfer in pyridoxal 5'-phosphate (PLP)-chemistry can hardly be overstated. Although experimental approaches to investigate pKa values can provide general guidance for assigning proton locations, only static pictures of the chemical species are available. To obtain the overall protein dynamics for the interpretation of detailed enzyme catalysis in this study, guided by information from solid-state NMR, we performed molecular dynamics (MD) simulations for the PLP-dependent enzyme tryptophan synthase (TRPS), whose catalytic mechanism features multiple quasi-stable intermediates. The primary objective of this work is to elucidate how the position of a single proton on the reacting substrate affects local and global protein dynamics during the catalytic cycle. In general, proteins create a chemical environment and an ensemble of conformational motions to recognize different substrates with different protonations. The study of these interactions in TRPS shows that functional groups on the reacting substrate, such as the phosphoryl group, pyridine nitrogen, phenolic oxygen and carboxyl group, of each PLP-bound intermediate play a crucial role in constructing an appropriate molecular interface with TRPS. In particular, the protonation states of the ionizable groups on the PLP cofactor may enhance or weaken the attractions between the enzyme and substrate. In addition, remodulation of the charge distribution for the intermediates may help generate a suitable environment for chemical reactions. The results of our study enhance knowledge of protonation states for several PLP intermediates and help to elucidate their effects on protein dynamics in the function of TRPS and other PLP-dependent enzymes.

Switches of hydrogen bonds during ligand-protein association processes determine binding kinetics

Revealing the processes of ligand-protein associations deepens our understanding of molecular recognition and binding kinetics. Hydrogen bonds (H-bonds) play a crucial role in optimizing ligand-protein interactions and ligand specificity. In addition to the formation of stable H-bonds in the final bound state, the formation of transient H-bonds during binding processes contributes binding kinetics that define a ligand as a fast or slow binder, which also affects drug action. However, the effect of forming the transient H-bonds on the kinetic properties is little understood. Guided by results from coarse-grained Brownian dynamics simulations, we used classical molecular dynamics simulations in an implicit solvent model and accelerated molecular dynamics simulations in explicit waters to show that the position and distribution of the H-bond donor or acceptor of a drug result in switching intermolecular and intramolecular H-bond pairs during ligand recognition processes. We studied two major types of HIV-1 protease ligands: a fast binder, xk263, and a slow binder, ritonavir. The slow association rate in ritonavir can be attributed to increased flexibility of ritonavir, which yields multistep transitions and stepwise entering patterns and the formation and breaking of complex H-bond pairs during the binding process. This model suggests the importance of conversions of spatiotemporal H-bonds during the association of ligands and proteins, which helps in designing inhibitors with preferred binding kinetics.

Catalytic roles of βLys87 in tryptophan synthase: (15)N solid state NMR studies

The proposed mechanism for tryptophan synthase shows βLys87 playing multiple catalytic roles: it bonds to the PLP cofactor, activates C4' for nucleophilic attack via a protonated Schiff base nitrogen, and abstracts and returns protons to PLP-bound substrates (i.e. acid-base catalysis). ε-¹⁵N-lysine TS was prepared to access the protonation state of βLys87 using ¹⁵N solid-state nuclear magnetic resonance (SSNMR) spectroscopy for three quasi-stable intermediates along the reaction pathway. These experiments establish that the protonation state of the ε-amino group switches between protonated and neutral states as the β-site undergoes conversion from one intermediate to the next during catalysis, corresponding to mechanistic steps where this lysine residue has been anticipated to play alternating acid and base catalytic roles that help steer reaction specificity in tryptophan synthase catalysis. This article is part of a Special Issue entitled: Cofactor-dependent proteins: evolution, chemical diversity and bio-applications. Guest Editors: Andrea Mozzarelli and Loredano Pollegioni.

Allostery and substrate channeling in the tryptophan synthase bienzyme complex: evidence for two subunit conformations and four quaternary states

The allosteric regulation of substrate channeling in tryptophan synthase involves ligand-mediated allosteric signaling that switches the α- and β-subunits between open (low activity) and closed (high activity) conformations. This switching prevents the escape of the common intermediate, indole, and synchronizes the α- and β-catalytic cycles. (19)F NMR studies of bound α-site substrate analogues, N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F6) and N-(4'-trifluoromethoxybenzenesulfonyl)-2-aminoethyl phosphate (F9), were found to be sensitive NMR probes of β-subunit conformation. Both the internal and external aldimine F6 complexes gave a single bound peak at the same chemical shift, while α-aminoacrylate and quinonoid F6 complexes all gave a different bound peak shifted by +1.07 ppm. The F9 complexes exhibited similar behavior, but with a corresponding shift of -0.12 ppm. X-ray crystal structures show the F6 and F9 CF3 groups located at the α-β subunit interface and report changes in both the ligand conformation and the surrounding protein microenvironment. Ab initio computational modeling suggests that the change in (19)F chemical shift results primarily from changes in the α-site ligand conformation. Structures of α-aminoacrylate F6 and F9 complexes and quinonoid F6 and F9 complexes show the α- and β-subunits have closed conformations wherein access of ligands into the α- and β-sites from solution is blocked. Internal and external aldimine structures show the α- and β-subunits with closed and open global conformations, respectively. These results establish that β-subunits exist in two global conformational states, designated open, where the β-sites are freely accessible to substrates, and closed, where the β-site portal into solution is blocked. Switching between these conformations is critically important for the αβ-catalytic cycle.

Long-range interactions in the α subunit of tryptophan synthase help to coordinate ligand binding, catalysis, and substrate channeling

The α-subunit of tryptophan synthase (αTS) catalyzes the conversion of indole-3-glycerol phosphate to d-glyceraldehyde-3-phosphate and indole. We propose that allosteric networks intrinsic to αTS are modulated by the binding of the β-subunit to regulate αTS function. Understanding these long-range amino acid networks in αTS thus gives insight into the coordination of the two active sites within TS. In this study, we have used Ala residues as probes for structural and dynamic changes of αTS throughout its catalytic cycle, in the absence of the β-subunit. Projection analysis of the chemical shift changes by site-specific amino acid substitutions and ligand titrations indicates that αTS has three important conformational states: ligand-free, glyceraldehyde-3-phosphate-bound(like), and the active states. The amino acid networks within these conformations are different, as suggested by chemical shift correlation analysis. In particular, there are long-range connections, only in the active state, between Ala47, which reports on structural and dynamic changes associated with the general acid/base Glu49, and residues within the β2α2 loop, which contains the catalytically important Asp60 residue. These long-range interactions are likely important for coordinating chemical catalysis. In the free state, but not in the active state, there are connections between the β2α2 and β6α6 loops that likely help to coordinate substrate binding. Changes in the allosteric networks are also accompanied by protein dynamic changes. During catalytic turnover, the protein becomes more rigid on the millisecond timescale and the active-site dynamics are driven to a faster nanosecond timescale.

Allosteric regulation of substrate channeling and catalysis in the tryptophan synthase bienzyme complex

The tryptophan synthase α2β2 bi-enzyme complex catalyzes the last two steps in the synthesis of l-tryptophan (l-Trp). The α-subunit catalyzes cleavage of 3-indole-d-glycerol 3'-phosphate (IGP) to give indole and d-glyceraldehyde 3'-phosphate (G3P). Indole is then transferred (channeled) via an interconnecting 25Å-long tunnel, from the α-subunit to the β-subunit where it reacts with l-Ser in a pyridoxal 5'-phosphate-dependent reaction to give l-Trp and a water molecule. The efficient utilization of IGP and l-Ser by tryptophan synthase to synthesize l-Trp utilizes a system of allosteric interactions that (1) function to switch the α-site on and off at different stages of the β-subunit catalytic cycle, and (2) prevent the escape of the channeled intermediate, indole, from the confines of the α- and β-catalytic sites and the interconnecting tunnel. This review discusses in detail the chemical origins of the allosteric interactions responsible both for switching the α-site on and off, and for triggering the conformational changes between open and closed states which prevent the escape of indole from the bienzyme complex.

Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

The role of oligomerization and cooperative regulation in protein function: the case of tryptophan synthase

The oligomerization/co-localization of protein complexes and their cooperative regulation in protein function is a key feature in many biological systems. The synergistic regulation in different subunits often enhances the functional properties of the multi-enzyme complex. The present study used molecular dynamics and Brownian dynamics simulations to study the effects of allostery, oligomerization and intermediate channeling on enhancing the protein function of tryptophan synthase (TRPS). TRPS uses a set of α/β–dimeric units to catalyze the last two steps of L-tryptophan biosynthesis, and the rate is remarkably slower in the isolated monomers. Our work shows that without their binding partner, the isolated monomers are stable and more rigid. The substrates can form fairly stable interactions with the protein in both forms when the protein reaches the final ligand–bound conformations. Our simulations also revealed that the α/β–dimeric unit stabilizes the substrate–protein conformation in the ligand binding process, which lowers the conformation transition barrier and helps the protein conformations shift from an open/inactive form to a closed/active form. Brownian dynamics simulations with a coarse-grained model illustrate how protein conformations affect substrate channeling. The results highlight the complex roles of protein oligomerization and the fine balance between rigidity and dynamics in protein function.

Conformational changes of enzymes are often related to regulating and creating an optimal environment for efficient chemistry. An increasing number of evidences also indicate that oligomerization/co-localization of proteins contributes to the efficiency of metabolic pathways. Although static structures have been available for many multi-enzyme complexes, their efficiency is also governed by the synergistic regulation between the multi-units. Our study applies molecular dynamics and Brownian dynamics simulations to the model system, the tryptophan synthase complex. The multi-enzyme complex is a bienzyme nanomachine and its catalytic activity is intimately related to the allosteric signaling and the metabolite transfer between its α– and β–subunits connected by a 25-Å long channel. Our studies suggest that the binding partner is crucial for the ligand binding processes. Although the isolated monomers are stable in the ligand–free state and can form stable interaction if the substrate is in the final bound conformation, it has higher energy barrier when ligand binds to the active site. We also show that the channel does not always exist, but it may be blocked before the enzyme reaches its final bound conformation. The results highlight the importance of forming protein complexes and the cooperative changes during different states.

T-Analyst: a program for efficient analysis of protein conformational changes by torsion angles

T-Analyst is a user-friendly computer program for analyzing trajectories from molecular modeling. Instead of using Cartesian coordinates for protein conformational analysis, T-Analyst is based on internal bond-angle-torsion coordinates in which internal torsion angle movements, such as side-chain rotations, can be easily detected. The program computes entropy and automatically detects and corrects angle periodicity to produce accurate rotameric states of dihedrals. It also clusters multiple conformations and detects dihedral rotations that contribute hinge-like motions. Correlated motions between selected dihedrals can also be observed from the correlation map. T-Analyst focuses on showing changes in protein flexibility between different states and selecting representative protein conformations for molecular docking studies. The program is provided with instructions and full source code in Perl.