Analysis of allosteric communication in a multienzyme complex by ancestral sequence reconstruction

NRIMD, a Web Server for Analyzing Protein Allosteric Interactions Based on Molecular Dynamics Simulation

Long-range allosteric communication between distant sites and active sites in proteins is central to biological regulation but still poorly characterized, limiting the development of protein engineering and drug design. Addressing this gap, NRIMD is an open-access web server for analyzing long-range interactions in proteins from molecular dynamics (MD) simulations, such as the effect of mutations at distal sites or allosteric ligand binding at allosteric sites on the active center. Based on our recent works on neural relational inference using graph neural networks, this cloud-based web server accepts MD simulation data on any length of residues in the alpha-carbon skeleton format from mainstream MD software. The input trajectory data are validated at the frontend deployed on the cloud and then processed on the backend deployed on a high-performance computer system with a collection of complementary tools. The web server provides a one-stop-shop MD analysis platform to predict long-range interactions and their paths between distant sites and active sites. It provides a user-friendly interface for detailed analysis and visualization. To the best of our knowledge, NRIMD is the first-of-its-kind online service to provide comprehensive long-range interaction analysis on MD simulations, which significantly lowers the barrier of predictions on protein long-range interactions using deep learning. The NRIMD web server is publicly available at https://nrimd.luddy.indianapolis.iu.edu/.

Ancestral Reconstruction and the Evolution of Protein Energy Landscapes

A protein's sequence determines its conformational energy landscape. This, in turn, determines the protein's function. Understanding the evolution of new protein functions therefore requires understanding how mutations alter the protein energy landscape. Ancestral sequence reconstruction (ASR) has proven a valuable tool for tackling this problem. In ASR, one phylogenetically infers the sequences of ancient proteins, allowing characterization of their properties. When coupled to biophysical, biochemical, and functional characterization, ASR can reveal how historical mutations altered the energy landscape of ancient proteins, allowing the evolution of enzyme activity, altered conformations, binding specificity, oligomerization, and many other protein features. In this article, we review how ASR studies have been used to dissect the evolution of energy landscapes. We also discuss ASR studies that reveal how energy landscapes have shaped protein evolution. Finally, we propose that thinking about evolution from the perspective of an energy landscape can improve how we approach and interpret ASR studies.

Harnessing conformational dynamics in enzyme catalysis to achieve nature-like catalytic efficiencies: the shortest path map tool for computational enzyme redesign

Enzymes exhibit diverse conformations, as represented in the free energy landscape (FEL). Such conformational diversity provides enzymes with the ability to evolve towards novel functions. The challenge lies in identifying mutations that enhance specific conformational changes, especially if located in distal sites from the active site cavity. The shortest path map (SPM) method, which we developed to address this challenge, constructs a graph based on the distances and correlated motions of residues observed in nanosecond timescale molecular dynamics (MD) simulations. We recently introduced a template based AlphaFold2 (tAF2) approach coupled with 10 nanosecond MD simulations to quickly estimate the conformational landscape of enzymes and assess how the FEL is shifted after mutation. In this study, we evaluate the potential of SPM when coupled with tAF2-MD in estimating conformational heterogeneity and identifying key conformationally-relevant positions. The selected model system is the beta subunit of tryptophan synthase (TrpB). We compare how the SPM pathways differ when integrating tAF2 with different MD simulation lengths from as short as 10 ns until 50 ns and considering two distinct Amber forcefield and water models (ff14SB/TIP3P versus ff19SB/OPC). The new methodology can more effectively capture the distal mutations found in laboratory evolution, thus showcasing the efficacy of tAF2-MD-SPM in rapidly estimating enzyme dynamics and identifying the key conformationally relevant hotspots for computational enzyme engineering.

The shortest path method (SPM) webserver for computational enzyme design

SPMweb is the online webserver of the Shortest Path Map (SPM) tool for identifying the key conformationally-relevant positions of a given enzyme structure and dynamics. The server is built on top of the DynaComm.py code and enables the calculation and visualization of the SPM pathways. SPMweb is easy-to-use as it only requires three input files: the three-dimensional structure of the protein of interest, and the two matrices (distance and correlation) previously computed from a Molecular Dynamics simulation. We provide in this publication information on how to generate the files for SPM construction even for non-expert users and discuss the most relevant parameters that can be modified. The tool is extremely fast (it takes less than one minute per job), thus allowing the rapid identification of distal positions connected to the active site pocket of the enzyme. SPM applications expand from computational enzyme design, especially if combined with other tools to identify the preferred substitution at the identified position, but also to rationalizing allosteric regulation, and even cryptic pocket identification for drug discovery. The simple user interface and setup make the SPM tool accessible to the whole scientific community. SPMweb is freely available for academia at http://spmosuna.com/.

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.

Impact of ancestral sequence reconstruction on mechanistic and structural enzymology

Ancestral sequence reconstruction (ASR) provides insight into the changes within a protein sequence across evolution. More specifically, it can illustrate how specific amino acid changes give rise to different phenotypes within a protein family. Over the last few decades it has established itself as a powerful technique for revealing molecular common denominators that govern enzyme function. Here, we describe the strength of ASR in unveiling catalytic mechanisms and emerging phenotypes for a range of different proteins, also highlighting biotechnological applications the methodology can provide.

Allostery, engineering and inhibition of tryptophan synthase

The final two steps of tryptophan biosynthesis are catalyzed by the enzyme tryptophan synthase (TS), composed of alpha (αTS) and beta (βTS) subunits. Recently, experimental and computational methods have mapped "allosteric networks" that connect the αTS and βTS active sites. In αTS, allosteric networks change across the catalytic cycle, which might help drive the conformational changes associated with its function. Directed evolution studies to increase catalytic function and expand the substrate profile of stand-alone βTS have also revealed the importance of αTS in modulating the conformational changes in βTS. These studies also serve as a foundation for the development of TS inhibitors, which can find utility against Mycobacterium tuberculosis and other bacterial pathogens.

Heme Pocket Structure and Its Functional Implications in an Ancestral Globin Protein

Proteins have undergone evolutionary processes to achieve optimal stability, increased functionality, and novel functions. Comparative analysis of existent and ancestral proteins provides insights into the factors that influence protein stability and function. Ancestral sequence reconstruction allows us to deduce the amino acid sequences of ancestral proteins. Here, we present the structural and functional characteristics of an ancestral protein, AncMH, reconstructed to be the last common ancestor of hemoglobins and myoglobins. Our findings reveal that AncMH harbors heme and that the heme binds oxygen. Furthermore, we demonstrate that the ferrous heme in AncMH is pentacoordinated, similar to that of human adult hemoglobin and horse myoglobin. A detailed comparison of the heme pocket structure indicates that the heme pocket in AncMH is more similar to that of hemoglobin than that of myoglobin. However, the autoxidation of AncMH is faster than that of both hemoglobin and myoglobin. Collectively, our results suggest that ancestral proteins of hemoglobins and myoglobins evolved in steps, including the hexa- to pentacoordination transition, followed by stabilization of the oxygen-bound form.

Simple mechanisms for the evolution of protein complexity

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Estimating conformational heterogeneity of tryptophan synthase with a template‐based Alphafold2 approach

The three‐dimensional structure of the enzymes provides very relevant information on the arrangement of the catalytic machinery and structural elements gating the active site pocket. The recent success of the neural network Alphafold2 in predicting the folded structure of proteins from the primary sequence with high levels of accuracy has revolutionized the protein design field. However, the application of Alphafold2 for understanding and engineering function directly from the obtained single static picture is not straightforward. Indeed, understanding enzymatic function requires the exploration of the ensemble of thermally accessible conformations that enzymes adopt in solution. In the present study, we evaluate the potential of Alphafold2 in assessing the effect of the mutations on the conformational landscape of the beta subunit of tryptophan synthase (TrpB). Specifically, we develop a template‐based Alphafold2 approach for estimating the conformational heterogeneity of several TrpB enzymes, which is needed for enhanced stand‐alone activity. Our results show the potential of Alphafold2, especially if combined with molecular dynamics simulations, for elucidating the changes induced by mutation in the conformational landscapes at a rather reduced computational cost, thus revealing its plausible application in computational enzyme design.

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.

Protein Conformational Space at the Edge of Allostery: Turning a Non-allosteric Malate Dehydrogenase into an "Allosterized" Enzyme using Evolution Guided Punctual Mutations

We unveil the intimate relationship between protein dynamics and allostery by following the trajectories of model proteins in their conformational and sequence spaces. Starting from a nonallosteric hyperthermophilic malate dehydrogenase, we have tracked the role of protein dynamics in the evolution of the allosteric capacity. Based on a large phylogenetic analysis of the malate (MalDH) and lactate dehydrogenase (LDH) superfamily, we identified two amino acid positions that could have had a major role for the emergence of allostery in LDHs, which we targeted for investigation by site-directed mutagenesis. Wild-type MalDH and the single and double mutants were tested with respect to their substrate recognition profiles. The double mutant displayed a sigmoid-shaped profile typical of homotropic activation in LDH. By using molecular dynamics simulations, we showed that the mutations induce a drastic change in the protein sampling of its conformational landscape, making transiently T-like (inactive) conformers, typical of allosteric LDHs, accessible. Our data fit well with the seminal key concept linking protein dynamics and evolvability. We showed that the selection of a new phenotype can be achieved by a few key dynamics-enhancing mutations causing the enrichment of low-populated conformational substates.

Design and engineering of allosteric communications in proteins

Allostery in proteins plays an important role in regulating protein activities and influencing many biological processes such as gene expression, enzyme catalysis, and cell signaling. The process of allostery takes place when a signal detected at a site on a protein is transmitted via a mechanical pathway to a functional site and, thus, influences its activity. The pathway of allosteric communication consists of amino acids that form a network with covalent and non-covalent bonds. By mutating residues in this allosteric network, protein engineers have successfully established novel allosteric pathways to achieve desired properties in the target protein. In this review, we highlight the most recent and state-of-the-art techniques for allosteric communication engineering. We also discuss the challenges that need to be overcome and future directions for engineering protein allostery.

A Model for Allosteric Communication in Drug Transport by the AcrAB-TolC Tripartite Efflux Pump

RND family efflux pumps are complex macromolecular machines involved in multidrug resistance by extruding antibiotics from the cell. While structural studies and molecular dynamics simulations have provided insights into the architecture and conformational states of the pumps, the path followed by conformational changes from the inner membrane protein (IMP) to the periplasmic membrane fusion protein (MFP) and to the outer membrane protein (OMP) in tripartite efflux assemblies is not fully understood. Here, we investigated AcrAB-TolC efflux pump's allostery by comparing resting and transport states using difference distance matrices supplemented with evolutionary couplings data and buried surface area measurements. Our analysis indicated that substrate binding by the IMP triggers quaternary level conformational changes in the MFP, which induce OMP to switch from the closed state to the open state, accompanied by a considerable increase in the interface area between the MFP subunits and between the OMPs and MFPs. This suggests that the pump's transport-ready state is at a more favourable energy level than the resting state, but raises the puzzle of how the pump does not become stably trapped in a transport-intermediate state. We propose a model for pump allostery that includes a downhill energetic transition process from a proposed 'activated' transport state back to the resting pump.

In Silico Identification and Experimental Validation of Distal Activity-Enhancing Mutations in Tryptophan Synthase

Allostery is a central mechanism for the regulation of multi-enzyme complexes. The mechanistic basis that drives allosteric regulation is poorly understood but harbors key information for enzyme engineering. In the present study, we focus on the tryptophan synthase complex that is composed of TrpA and TrpB subunits, which allosterically activate each other. Specifically, we develop a rational approach for identifying key amino acid residues of TrpB distal from the active site. Those residues are predicted to be crucial for shifting the inefficient conformational ensemble of the isolated TrpB to a productive ensemble through intra-subunit allosteric effects. The experimental validation of the conformationally driven TrpB design demonstrates its superior stand-alone activity in the absence of TrpA, comparable to those enhancements obtained after multiple rounds of experimental laboratory evolution. Our work evidences that the current challenge of distal active site prediction for enhanced function in computational enzyme design has become within reach.

Exploiting enzyme evolution for computational protein design

Recent years have seen an explosion of interest in understanding the physicochemical parameters that shape enzyme evolution, as well as substantial advances in computational enzyme design. This review discusses three areas where evolutionary information can be used as part of the design process: (i) using ancestral sequence reconstruction (ASR) to generate new starting points for enzyme design efforts; (ii) learning from how nature uses conformational dynamics in enzyme evolution to mimic this process in silico; and (iii) modular design of enzymes from smaller fragments, again mimicking the process by which nature appears to create new protein folds. Using showcase examples, we highlight the importance of incorporating evolutionary information to continue to push forward the boundaries of enzyme design studies.

Substitution of a Surface-Exposed Residue Involved in an Allosteric Network Enhances Tryptophan Synthase Function in Cells

Networks of noncovalent amino acid interactions propagate allosteric signals throughout proteins. Tryptophan synthase (TS) is an allosterically controlled bienzyme in which the indole product of the alpha subunit (αTS) is transferred through a 25 Å hydrophobic tunnel to the active site of the beta subunit (βTS). Previous nuclear magnetic resonance and molecular dynamics simulations identified allosteric networks in αTS important for its function. We show here that substitution of a distant, surface-exposed network residue in αTS enhances tryptophan production, not by activating αTS function, but through dynamically controlling the opening of the indole channel and stimulating βTS activity. While stimulation is modest, the substitution also enhances cell growth in a tryptophan-auxotrophic strain of Escherichia coli compared to complementation with wild-type αTS, emphasizing the biological importance of the network. Surface-exposed networks provide new opportunities in allosteric drug design and protein engineering, and hint at potential information conduits through which the functions of a metabolon or even larger proteome might be coordinated and regulated.

Exploring the Allosteric Territory of Protein Function

While the pervasiveness of allostery in proteins is commonly accepted, we further show the generic nature of allosteric mechanisms by analyzing here transmembrane ion-channel viroporin 3a and RNA-dependent RNA polymerase (RdRp) from SARS-CoV-2 along with metabolic enzymes isocitrate dehydrogenase 1 (IDH1) and fumarate hydratase (FH) implicated in cancers. Using the previously developed structure-based statistical mechanical model of allostery (SBSMMA), we share our experience in analyzing the allosteric signaling, predicting latent allosteric sites, inducing and tuning targeted allosteric response, and exploring the allosteric effects of mutations. This, yet incomplete list of phenomenology, forms a complex and unique allosteric territory of protein function, which should be thoroughly explored. We propose a generic computational framework, which not only allows one to obtain a comprehensive allosteric control over proteins but also provides an opportunity to approach the fragment-based design of allosteric effectors and drug candidates. The advantages of allosteric drugs over traditional orthosteric compounds, complemented by the emerging role of the allosteric effects of mutations in the expansion of the cancer mutational landscape and in the increased mutability of viral proteins, leave no choice besides further extensive studies of allosteric mechanisms and their biomedical implications.

Ancestral lysosomal enzymes with increased activity harbor therapeutic potential for treatment of Hunter syndrome

We show the successful application of ancestral sequence reconstruction to enhance the activity of iduronate-2-sulfatase (IDS), thereby increasing its therapeutic potential for the treatment of Hunter syndrome—a lysosomal storage disease caused by impaired function of IDS. Current treatment, enzyme replacement therapy with recombinant human IDS, does not alleviate all symptoms, and an unmet medical need remains. We reconstructed putative ancestral sequences of mammalian IDS and compared them with extant IDS. Some ancestral variants displayed up to 2-fold higher activity than human IDS in in vitro assays and cleared more substrate in ex vivo experiments in patient fibroblasts. This could potentially allow for lower dosage or enhanced therapeutic effect in enzyme replacement therapy, thereby improving treatment outcomes and cost efficiency, as well as reducing treatment burden. In summary, we showed that ancestral sequence reconstruction can be applied to lysosomal enzymes that function in concert with modern enzymes and receptors in cells.

•

Reconstruction of ancestral lysosomal enzymes that function in complex cellular context

•

Ancestral iduronate-2-sulfatases with increased activity compared with the human enzyme

•

Increased clearance of substrate in patient fibroblasts indicates therapeutic potential

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.

Tryptophan Synthase: Biocatalyst Extraordinaire

Tryptophan synthase (TrpS) has emerged as a paragon of noncanonical amino acid (ncAA) synthesis and is an ideal biocatalyst for synthetic and biological applications. TrpS catalyzes an irreversible, C-C bond-forming reaction between indole and serine to make l-tryptophan; native TrpS complexes possess fairly broad specificity for indole analogues, but are difficult to engineer to extend substrate scope or to confer other useful properties due to allosteric constraints and their heterodimeric structure. Directed evolution freed the catalytically relevant TrpS β-subunit (TrpB) from allosteric regulation by its TrpA partner and has enabled dramatic expansion of the enzyme's substrate scope. This review examines the long and storied career of TrpS from the perspective of its application in ncAA synthesis and biocatalytic cascades.

Substrate Channeling via a Transient Protein-Protein Complex: The case of D-Glyceraldehyde-3-Phosphate Dehydrogenase and L-Lactate Dehydrogenase

Substrate channeling studies have frequently failed to provide conclusive results due to poor understanding of this subtle phenomenon. We analyzed the mechanism of NADH-channeling from D-glyceraldehyde-3-phosphate dehydrogenase (GAPDH) to L-lactate Dehydrogenase (LDH) using enzymes from different cells. Enzyme kinetics studies showed that LDH activity with free NADH and GAPDH-NADH complex always take place in parallel. The channeling is observed only in assays that mimic cytosolic conditions where free NADH concentration is negligible and the GAPDH-NADH complex is dominant. Molecular dynamics and protein-protein interaction studies showed that LDH and GAPDH can form a leaky channeling complex only at the limiting NADH concentrations. Surface calculations showed that positive electric field between the NAD(H) binding sites on LDH and GAPDH tetramers can merge in the LDH-GAPDH complex. NAD(H)-channeling within the LDH-GAPDH complex can be an extension of NAD(H)-channeling within each tetramer. In the case of a transient LDH-(GAPDH-NADH) complex, the relative contribution from the channeled and the diffusive paths depends on the overlap between the off-rates for the LDH-(GAPDH-NADH) complex and the GAPDH-NADH complex. Molecular evolution or metabolic engineering protocols can exploit substrate channeling for metabolic flux control by fine-tuning substrate-binding affinity for the key enzymes in the competing reaction paths.

Small design from big alignment: engineering proteins with multiple sequence alignment as the starting point

Multiple sequence alignment (MSA) is a fundamental way to gain information that cannot be obtained from the analysis of any individual sequence included in the alignment. It provides ways to investigate the relationship between sequence and function from a perspective of evolution. Thus, the MSA of proteins can be employed as a reference for protein engineering. In this paper, we reviewed the recent advances to highlight how protein engineering was benefited from the MSA of proteins. These methods include (1) engineering the thermostability or solubility of proteins by making it closer to the consensus sequence of the alignment through introducing site mutations; (2) structure-based engineering proteins with comparative modeling; (3) creating paleoenzymes featured with high thermostability and promiscuity by constructing the ancestral sequences derived from multiple sequence alignment; and (4) incorporating site-mutations targeting the evolutionarily coupled sites identified from multiple sequence alignment.