Ancient thioredoxins evolved to modern-day stability-function requirement by altering native state ensemble

Thr-to-Ala Mutation Leads to a Larger Aromatic Pair and Reduced Packing Density in α1,α3-Helices during Thioredoxin Cold Adaptation

This study investigates the impact of aromatic-aromatic interactions on the cold adaptation of thioredoxin (Trx), a small redox protein with a conserved Trx-fold structure. Two Trx orthologs, one from the psychrophilic Arctic bacterium Sphingomonas sp. (SpTrx) and the other from the mesophilic Escherichia coli (EcTrx), display distinct aromatic interactions in their α1,α3-helices. SpTrx features a larger Trp11-Phe69 pair, while EcTrx employs a smaller Phe12-Tyr70 pair along with an additional Asp9-Thr66 hydrogen bond. Smaller aromatic residues in SpTrx (Phe-Phe or Phe-Tyr pair) lead to decreased thermal and thermodynamic stabilities, increased conformational flexibility, and reduced enzyme activity. In contrast, EcTrx's thermal stability is primarily influenced by the larger Trp residue, especially in the more hydrophobic Trp-Phe pair compared to the Trp-Tyr pair. Both SpTrx and EcTrx exhibit a strengthening of the Asp-Thr hydrogen bond by a Phe-Tyr pair and a weakening by a Trp-Phe pair. Additionally, the Asp8-Thr65 hydrogen bond in SpTrx contributes to the destabilization of the Phe-Phe pair. Molecular dynamics simulations of SpTrx indicate that a smaller aromatic pair or the Asp-Thr hydrogen bond in the α1,α3-helices further destabilizes the α2-helix across the central β-sheet. Our results suggest that the Thr-to-Ala mutation destabilizes the α1,α3-helices, resulting in a larger aromatic pair and reduced packing density in psychrophilic Trxs during cold adaptation. These findings enhance our understanding of Trx's adaptation to colder temperatures.

Structural conservation in the glutathione binding in Sphingomonas sp. glutaredoxin Grx3 and variations for cold adaptation

Glutaredoxin 3 (Grx3), a redox protein with a thioredoxin-fold structure, maintains structural integrity and glutathione (GSH) binding capabilities across varying habitat temperatures. The cis-Pro loop, essential for GSH binding, relies on the Arg-Asp salt bridge (α2-α3) and Gln-His hydrogen bond (β3-β4) for its conformation. In some psychrophilic Grx3 variants, Arg in α2 is replaced with Tyr, and His in β4 is replaced with Phe. This study examines the roles of these bonds in Grx3's structure, function, and cold adaptation, using SpGrx3 from the Arctic bacterium Sphingomonas sp. Despite its cold habitat, SpGrx3 maintains the Arg51-Asp69 salt bridge and Gln56-His63 hydrogen bond. The R51Y substitution disrupts the α2-α3 salt bridge, while the H63F and H63Y substitutions hinder the salt bridge through cation-π interactions with Arg51, involving Phe63/Tyr63, thereby enhancing flexibility. Conversely, mutations that disrupt the hydrogen bond (Q56A, H63A, and H63F) reduce thermal stability. In the psychrophilic Grx3 configuration A48T/R51Y/H63F, a Thr48-Gln56 hydrogen bond stabilizes the cis-Pro loop, enhancing flexibility by disrupting both bonds. Furthermore, all mutants exhibit reduced α-helical content and catalytic efficiency. In summary, the highly conserved Arg51-Asp69 salt bridge and Gln56-His63 hydrogen bond are crucial for stabilizing the cis-Pro loop and catalytic activity in SpGrx3. His63 is favored as it avoids cation-π interactions with Arg51, unlike Phe63/Tyr63. Psychrophilic Grx3 variants have adapted to cold environments by reducing GSH binding and increasing structural flexibility. These findings deepen our understanding of the structural conservation in Grx3 for GSH binding and the critical alterations required for cold adaptation.

Role of thioredoxin in chronic obstructive pulmonary disease (COPD): a promising future target

INTRODUCTION

Thioredoxin (Trx) is a secretory protein that acts as an antioxidant, redox regulator, anti-allergic, and anti-inflammatory molecule. It has been used to treat dermatitis and inflammation of the digestive tract. In the lungs, Trx has a significant anti-inflammatory impact. On the other hand, Chronic Obstructive Pulmonary Disease (COPD) is one of the significant causes of death in the developed world, with a tremendous individual and socioeconomic impact. Despite new initiatives and endless treatment trials, COPD incidence and death will likely escalate in the coming decades.

AREAS COVERED

COPD is a chronic inflammatory disease impacting the airways, lung parenchyma, and pulmonary vasculature. Oxidative stress and protease-antiprotease imbalances are thought to be involved in the process. The most popular respiratory inflammatory and allergic disorders therapies are corticosteroids and β-receptor agonists. These medications are helpful but have some drawbacks, such as infection and immunosuppression; thus, addressing Trx signalling treatments may be a viable COPD treatment approach. This review shall cover the pathophysiology of COPD, the pharmacognosy of anti-COPD drugs, including the assets and liabilities of each, and the role and mechanism of Trx in COPD treatment.

EXPERT OPINION

Limited research has targeted the thioredoxin system as an anti-COPD drug. Spectating the increase in the mortality rates of COPD, this review article would be an interesting one to research.

Conformational Tuning Shapes the Balance between Functional Promiscuity and Specialization in Paralogous Plasmodium Acyl-CoA Binding Proteins

Paralogous proteins confer enhanced fitness to organisms via complex sequence-conformation codes that shape functional divergence, specialization, or promiscuity. Here, we dissect the underlying mechanism of promiscuous binding versus partial subfunctionalization in paralogues by studying structurally identical acyl-CoA binding proteins (ACBPs) from Plasmodium falciparum that serve as promising drug targets due to their high expression during the protozoan proliferative phase. Combining spectroscopic measurements, solution NMR, SPR, and simulations on two of the paralogues, A16 and A749, we show that minor sequence differences shape nearly every local and global conformational feature. A749 displays a broader and heterogeneous native ensemble, weaker thermodynamic coupling and cooperativity, enhanced fluctuations, and a larger binding pocket volume compared to A16. Site-specific tryptophan probes signal a graded reduction in the sampling of substates in the holo form, which is particularly apparent in A749. The paralogues exhibit a spectrum of binding affinities to different acyl-CoAs with A749, the more promiscuous and hence the likely ancestor, binding 1000-fold stronger to lauroyl-CoA under physiological conditions. We thus demonstrate how minor sequence changes modulate the extent of long-range interactions and dynamics, effectively contributing to the molecular evolution of contrasting functional repertoires in paralogues.

Engineering gain-of-function mutants of a WW domain by dynamics and structural analysis

Proteins gain optimal fitness such as foldability and function through evolutionary selection. However, classical studies have found that evolutionarily designed protein sequences alone cannot guarantee foldability, or at least not without considering local contacts associated with the initial folding steps. We previously showed that foldability and function can be restored by removing frustration in the folding energy landscape of a model WW domain protein, CC16, which was designed based on Statistical Coupling Analysis (SCA). Substitutions ensuring the formation of five local contacts identified as "on-path" were selected using the closest homolog native folded sequence, N21. Surprisingly, the resulting sequence, CC16-N21, bound to Group I peptides, while N21 did not. Here, we identified single-point mutations that enable N21 to bind a Group I peptide ligand through structure and dynamic-based computational design. Comparison of the docked position of the CC16-N21/ligand complex with the N21 structure showed that residues at positions 9 and 19 are important for peptide binding, whereas the dynamic profiles identified position 10 as allosterically coupled to the binding site and exhibiting different dynamics between N21 and CC16-N21. We found that swapping these positions in N21 with matched residues from CC16-N21 recovers nature-like binding affinity to N21. This study validates the use of dynamic profiles as guiding principles for affecting the binding affinity of small proteins.

Unraveling the roles of aromatic cluster side-chain interactions on the structural stability and functional significance of psychrophilic Sphingomonas sp. glutaredoxin 3

This study investigates the impact of aromatic cluster side-chain interactions in Grx3 (SpGrx3) from the psychrophilic Arctic bacterium Sphingomonas sp. Grx3 is a class I oxidoreductase with a unique parallel arrangement of aromatic residues in its aromatic cluster, unlike the tetrahedral geometry observed in Trxs. Hydrophilic-to-hydrophobic substitutions were made in the aromatic cluster, in β1 (E5V and Y7F), adjacent β2 (Y32F and Y32L), both β1 and β2 (E5V/Y32L), and short α2 (R47F). The hydrophobic substitutions, particularly those at or near Tyr7 (E5V, Y7F, Y32F, and R47F), increased melting temperatures and conformational stability, whereas disrupting β1-β2 interactions (Y32L and E5V/Y32L) led to structural instability of SpGrx3. However, excessive hydrophobic interactions (Y7F and E5V/Y32L) caused protein aggregation at elevated temperatures. All mutations resulted in a reduction in α-helical content and an increase in β-strand content. The R47F mutant, which formed dimers and exhibited the highest β-strand content, showed increased conformational flexibility and a significant decrease in catalytic rate due to the disturbance of β1-α2 interactions. In summary, the configuration of the aromatic cluster, especially Tyr7 in the buried β1 and Arg47 in the short α2, played crucial roles in maintaining the active conformation of SpGrx3 and preventing its protein aggregation. These modifications, reducing hydrophobicity in the central β-sheet, distinguish Grx3 from other Trx-fold proteins, highlighting evolutionary divergence within the Trx-fold superfamily and its functional versatility.

Protein dynamics provide mechanistic insights about epistasis among common missense polymorphisms

Sequencing of the protein coding genome has revealed many different missense mutations of human proteins and different population frequencies of corresponding haplotypes, which consist of different sets of those mutations. Here, we present evidence for pairwise intramolecular epistasis (i.e., nonadditive interactions) between many such mutations through an analysis of protein dynamics. We suggest that functional compensation for conserving protein dynamics is a likely evolutionary mechanism that maintains high-frequency mutations that are individually nonneutral but epistatically compensating within proteins. This analysis is the first of its type to look at human proteins with specific high population frequency mutations and examine the relationship between mutations that make up that observed high-frequency protein haplotype. Importantly, protein dynamics revealed a separation between high and low frequency haplotypes within a target protein cytochrome P450 2A7, with the high-frequency haplotypes showing behavior closer to the wild-type protein. Common protein haplotypes containing two mutations display dynamic compensation in which one mutation can correct for the dynamic effects of the other. We also utilize a dynamics-based metric, EpiScore, that evaluates the epistatic interactions and allows us to see dynamic compensation within many other proteins.

Correlated Evolution of Low-Frequency Vibrations and Function in Enzymes

Previous studies of the flexibility of ancestral proteins suggest that proteins evolve their function by altering their native state ensemble. Here, we propose a more direct method to analyze such changes during protein evolution by comparing thermally activated vibrations at frequencies below 6 THz, which report on the dynamics of collective protein modes. We analyzed the backbone vibrational density of states of ancestral and extant β-lactamases and thioredoxins and observed marked changes in the vibrational spectrum in response to evolution. Coupled with previously observed changes in protein flexibility, the observed shifts of vibrational mode densities suggest that protein dynamics and dynamical allostery are critical factors for the evolution of enzymes with specialized catalytic and biophysical properties.

Structure-guided isoform identification for the human transcriptome

Recently developed methods to predict three-dimensional protein structure with high accuracy have opened new avenues for genome and proteome research. We explore a new hypothesis in genome annotation, namely whether computationally predicted structures can help to identify which of multiple possible gene isoforms represents a functional protein product. Guided by protein structure predictions, we evaluated over 230,000 isoforms of human protein-coding genes assembled from over 10,000 RNA sequencing experiments across many human tissues. From this set of assembled transcripts, we identified hundreds of isoforms with more confidently predicted structure and potentially superior function in comparison to canonical isoforms in the latest human gene database. We illustrate our new method with examples where structure provides a guide to function in combination with expression and evolutionary evidence. Additionally, we provide the complete set of structures as a resource to better understand the function of human genes and their isoforms. These results demonstrate the promise of protein structure prediction as a genome annotation tool, allowing us to refine even the most highly curated catalog of human proteins. More generally we demonstrate a practical, structure-guided approach that can be used to enhance the annotation of any genome.

Structural and Dynamic Effects of PTEN C-Terminal Tail Phosphorylation

The phosphatase and tensin homologue deleted on chromosome 10 (PTEN) tumor suppressor gene encodes a tightly regulated dual-specificity phosphatase that serves as the master regulator of PI3K/AKT/mTOR signaling. The carboxy-terminal tail (CTT) is key to regulation and harbors multiple phosphorylation sites (Ser/Thr residues 380-385). CTT phosphorylation suppresses the phosphatase activity by inducing a stable, closed conformation. However, little is known about the mechanisms of phosphorylation-induced CTT-deactivation dynamics. Using explicit solvent microsecond molecular dynamics simulations, we show that CTT phosphorylation leads to a partially collapsed conformation, which alters the secondary structure of PTEN and induces long-range conformational rearrangements that encompass the active site. The active site rearrangements prevent localization of PTEN to the membrane, precluding lipid phosphatase activity. Notably, we have identified phosphorylation-induced allosteric coupling between the interdomain region and a hydrophobic site neighboring the active site in the phosphatase domain. Collectively, the results provide a mechanistic understanding of CTT phosphorylation dynamics and reveal potential druggable allosteric sites in a previously believed clinically undruggable protein.

Dynamic coupling of residues within proteins as a mechanistic foundation of many enigmatic pathogenic missense variants

Many pathogenic missense mutations are found in protein positions that are neither well-conserved nor fall in any known functional domains. Consequently, we lack any mechanistic underpinning of dysfunction caused by such mutations. We explored the disruption of allosteric dynamic coupling between these positions and the known functional sites as a possible mechanism for pathogenesis. In this study, we present an analysis of 591 pathogenic missense variants in 144 human enzymes that suggests that allosteric dynamic coupling of mutated positions with known active sites is a plausible biophysical mechanism and evidence of their functional importance. We illustrate this mechanism in a case study of β-Glucocerebrosidase (GCase) in which a vast majority of 94 sites harboring Gaucher disease-associated missense variants are located some distance away from the active site. An analysis of the conformational dynamics of GCase suggests that mutations on these distal sites cause changes in the flexibility of active site residues despite their distance, indicating a dynamic communication network throughout the protein. The disruption of the long-distance dynamic coupling caused by missense mutations may provide a plausible general mechanistic explanation for biological dysfunction and disease.

Dynamic Allostery Highlights the Evolutionary Differences between the CoV-1 and CoV-2 Main Proteases

The SARS-CoV-2 coronavirus has become one of the most immediate and widely-studied systems since its identification and subsequent global outbreak from 2019-2021. In an effort to understand the biophysical changes as a result of mutations, the mechanistic details of multiple different proteins within the SARS-CoV-2 virus have been studied and compared with SARS-CoV-1. Focusing on the main protease (mPro), we first explored the long-range dynamics using the Dynamic Coupling Index (DCI) to investigate the dynamic coupling between the catalytic site residues and the rest of the protein, both inter and intra chain, for the CoV-1 and CoV-2 mPro. We found that there is significant cross-chain coupling between these active sites and specific distal residues in the CoV-2 mPro not present in CoV-1. The enhanced long distance interactions, particularly between the two chains, suggest subsequently enhanced cooperativity for CoV-2. A further comparative analysis of the dynamic flexibility using the Dynamic Flexibility Index (DFI) between the CoV-1 and CoV-2 mPros shows that the inhibitor binding near active sites induces change in flexibility to a distal region of the protein, opposite in behavior between the two systems; this region becomes more flexible upon inhibitor binding in CoV-1 while it becomes less flexible in the CoV-2 mPro. Upon inspection, we show that, on average, the dynamic flexibility of the sites substituted from CoV-1 to CoV-2 changes significantly less than the average calculated across all residues within the structure, indicating that the differences in behaviors between the two systems is likely the result of allosteric influence, where the new substitutions in CoV-2 induce flexibility and dynamical changes elsewhere in the structure.

Design of novel cyanovirin-N variants by modulation of binding dynamics through distal mutations

We develop integrated co-evolution and dynamic coupling (ICDC) approach to identify, mutate, and assess distal sites to modulate function. We validate the approach first by analyzing the existing mutational fitness data of TEM-1 β-lactamase and show that allosteric positions co-evolved and dynamically coupled with the active site significantly modulate function. We further apply ICDC approach to identify positions and their mutations that can modulate binding affinity in a lectin, cyanovirin-N (CV-N), that selectively binds to dimannose, and predict binding energies of its variants through Adaptive BP-Dock. Computational and experimental analyses reveal that binding enhancing mutants identified by ICDC impact the dynamics of the binding pocket, and show that rigidification of the binding residues compensates for the entropic cost of binding. This work suggests a mechanism by which distal mutations modulate function through dynamic allostery and provides a blueprint to identify candidates for mutagenesis in order to optimize protein function.

Essential roles of buried phenylalanine in the structural stability of thioredoxin from a psychrophilic Arctic bacterium Sphingomonas sp

Thioredoxin (Trx), a small redox protein, exhibits thermal stability at high temperatures regardless of its origin, including psychrophiles. Trxs have a common structure consisting of the central β-sheet flanked by an aliphatic cluster on one side and an aromatic cluster on the other side. Although the roles of aromatic amino acids in the folding and stability of proteins have been studied extensively, the contributions of aromatic residues to the stability and function of Trx, particularly Trxs from cold-adapted organisms, have not been fully elucidated. This study examined the roles of aromatic amino acids in the aromatic cluster of a Trx from the psychrophilic Arctic bacterium Sphingomonas sp. PAMC 26621 (SpTrx). The aromatic cluster of SpTrx was comprised of W11, F26, F69, and F80, in which F26 at the β2 terminus was buried inside. The substitution of tyrosine for F26 changed the SpTrx conformation substantially compared to that of F69 and F80. Further biochemical and spectroscopic investigations on F26 showed that the F26Y, F26W, and F26A mutants resulted in structural instability of SpTrx in both urea- and temperature-induced unfolding and lower insulin reduction activities. The Trx reductase (SpTR) showed lower catalytic efficiencies against F26 mutants compared to the wild-type SpTrx. These results suggest that buried F26 is essential for maintaining the active-site conformation of SpTrx as an oxidoreductase and its structural stability for interactions with SpTR at colder temperatures.

Can thiol-based redox systems be utilized as parts for synthetic biology applications?

OBJECTIVES

Synthetic biology has emerged from molecular biology and engineering approaches and aims to develop novel, biologically-inspired systems for industrial and basic research applications ranging from biocomputing to drug production. Surprisingly, redoxin (thioredoxin, glutaredoxin, peroxiredoxin) and other thiol-based redox systems have not been widely utilized in many of these synthetic biology applications.

METHODS

We reviewed thiol-based redox systems and the development of synthetic biology applications that have used thiol-dependent parts.

RESULTS

The development of circuits to facilitate cytoplasmic disulfide bonding, biocomputing and the treatment of intestinal bowel disease are amongst the applications that have used thiol-based parts. We propose that genetically encoded redox sensors, thiol-based biomaterials and intracellular hydrogen peroxide generators may also be valuable components for synthetic biology applications.

DISCUSSION

Thiol-based systems play multiple roles in cellular redox metabolism, antioxidant defense and signaling and could therefore offer a vast and diverse portfolio of components, parts and devices for synthetic biology applications. However, factors limiting the adoption of redoxin systems for synthetic biology applications include the orthogonality of thiol-based components, limitations in the methods to characterize thiol-based systems and an incomplete understanding of the design principles of these systems.

Hydroperoxide-Reducing Enzymes in the Regulation of Free-Radical Processes

The review presents current concepts of the molecular mechanisms of oxidative stress development and describes main stages of the free-radical reactions in oxidative stress. Endogenous and exogenous factors of the oxidative stress development, including dysfunction of cell oxidoreductase systems, as well as the effects of various external physicochemical factors, are discussed. The review also describes the main components of the antioxidant defense system and stages of its evolution, with a special focus on peroxiredoxins, glutathione peroxidases, and glutathione S-transferases, which share some phylogenetic, structural, and catalytic properties. The substrate specificity, as well as the similarities and differences in the catalytic mechanisms of these enzymes, are discussed in detail. The role of peroxiredoxins, glutathione peroxidases, and glutathione S-transferases in the regulation of hydroperoxide-mediated intracellular and intercellular signaling and interactions of these enzymes with receptors and non-receptor proteins are described. An important contribution of hydroperoxide-reducing enzymes to the antioxidant protection and regulation of such cell processes as growth, differentiation, and apoptosis is demonstrated.

Substitutions at Nonconserved Rheostat Positions Modulate Function by Rewiring Long-Range, Dynamic Interactions

Amino acid substitutions at nonconserved protein positions can have noncanonical and "long-distance" outcomes on protein function. Such outcomes might arise from changes in the internal protein communication network, which is often accompanied by changes in structural flexibility. To test this, we calculated flexibilities and dynamic coupling for positions in the linker region of the lactose repressor protein. This region contains nonconserved positions for which substitutions alter DNA-binding affinity. We first chose to study 11 substitutions at position 52. In computations, substitutions showed long-range effects on flexibilities of DNA-binding positions, and the degree of flexibility change correlated with experimentally measured changes in DNA binding. Substitutions also altered dynamic coupling to DNA-binding positions in a manner that captured other experimentally determined functional changes. Next, we broadened calculations to consider the dynamic coupling between 17 linker positions and the DNA-binding domain. Experimentally, these linker positions exhibited a wide range of substitution outcomes: Four conserved positions tolerated hardly any substitutions ("toggle"), ten nonconserved positions showed progressive changes from a range of substitutions ("rheostat"), and three nonconserved positions tolerated almost all substitutions ("neutral"). In computations with wild-type lactose repressor protein, the dynamic couplings between the DNA-binding domain and these linker positions showed varied degrees of asymmetry that correlated with the observed toggle/rheostat/neutral substitution outcomes. Thus, we propose that long-range and noncanonical substitutions outcomes at nonconserved positions arise from rewiring long-range communication among functionally important positions. Such calculations might enable predictions for substitution outcomes at a range of nonconserved positions.

Diverse Native Ensembles Dictate the Differential Functional Responses of Nuclear Receptor Ligand-Binding Domains

Native states of folded proteins are characterized by a large ensemble of conformations whose relative populations and interconversion dynamics determine the functional output. This is more apparent in transcription factors that have evolved to be inherently sensitive to small perturbations, thus fine-tuning gene expression. To explore the extent to which such functional features are imprinted on the folding landscape of transcription factor ligand-binding domains (LBDs), we characterize paralogous LBDs of the nuclear receptor (NR) family employing an energetically detailed and ensemble-based Ising-like statistical mechanical model. We find that the native ensembles of the LBDs from glucocorticoid receptor, PPAγ, and thyroid hormone receptor display a remarkable diversity in the width of the native wells, the number and nature of partially structured states, and hence the degree of conformational order. Monte Carlo simulations employing the full state representation of the ensemble highlight that many of the functional conformations coexist in equilibrium, whose relative populations are sensitive to both temperature and the strength of ligand binding. Allosteric modulation of the degree of structure at a coregulator binding site on ligand binding is shown to arise via a redistribution of populations in the native ensembles of glucocorticoid and PPAγ LBDs. Our results illustrate how functional requirements can drive the evolution of conformationally diverse native ensembles in paralogs.

Hinge-shift mechanism as a protein design principle for the evolution of β-lactamases from substrate promiscuity to specificity

TEM-1 β-lactamase degrades β-lactam antibiotics with a strong preference for penicillins. Sequence reconstruction studies indicate that it evolved from ancestral enzymes that degraded a variety of β-lactam antibiotics with moderate efficiency. This generalist to specialist conversion involved more than 100 mutational changes, but conserved fold and catalytic residues, suggesting a role for dynamics in enzyme evolution. Here, we develop a conformational dynamics computational approach to rationally mold a protein flexibility profile on the basis of a hinge-shift mechanism. By deliberately weighting and altering the conformational dynamics of a putative Precambrian β-lactamase, we engineer enzyme specificity that mimics the modern TEM-1 β-lactamase with only 21 amino acid replacements. Our conformational dynamics design thus re-enacts the evolutionary process and provides a rational allosteric approach for manipulating function while conserving the enzyme active site.

The Role of Rigid Residues in Modulating TEM-1 β-Lactamase Function and Thermostability

The relationship between protein motions (i.e., dynamics) and enzymatic function has begun to be explored in β-lactamases as a way to advance our understanding of these proteins. In a recent study, we analyzed the dynamic profiles of TEM-1 (a ubiquitous class A β-lactamase) and several ancestrally reconstructed homologues. A chief finding of this work was that rigid residues that were allosterically coupled to the active site appeared to have profound effects on enzyme function, even when separated from the active site by many angstroms. In the present work, our aim was to further explore the implications of protein dynamics on β-lactamase function by altering the dynamic profile of TEM-1 using computational protein design methods. The Rosetta software suite was used to mutate amino acids surrounding either rigid residues that are highly coupled to the active site or to flexible residues with no apparent communication with the active site. Experimental characterization of ten designed proteins indicated that alteration of residues surrounding rigid, highly coupled residues, substantially affected both enzymatic activity and stability; in contrast, native-like activities and stabilities were maintained when flexible, uncoupled residues, were targeted. Our results provide additional insight into the structure-function relationship present in the TEM family of β-lactamases. Furthermore, the integration of computational protein design methods with analyses of protein dynamics represents a general approach that could be used to extend our understanding of the relationship between dynamics and function in other enzyme classes.

Protein folding stability and binding interactions through the lens of evolution: a dynamical perspective

While the function of a protein depends heavily on its ability to fold into a correct 3D structure, billions of years of evolution have tailored proteins from highly stable objects to flexible molecules as they adapted to environmental changes. Nature maintains the fine balance of protein folding and stability while still evolving towards new function through generations of fine-tuning necessary interactions with other proteins and small molecules. Here we focus on recent computational and experimental studies that shed light onto how evolution molds protein folding and the functional landscape from a conformational dynamics' perspective. Particularly, we explore the importance of dynamic allostery throughout protein evolution and discuss how the protein anisotropic network can give rise to allosteric and epistatic interactions.

Slow Folding of a Helical Protein: Large Barriers, Strong Internal Friction, or a Shallow, Bumpy Landscape?

The rate at which a protein molecule folds is determined by opposing energetic and entropic contributions to the free energy that shape the folding landscape. Delineating the extent to which they impact the diffusional barrier-crossing events, including the magnitude of internal friction and barrier height, has largely been a challenging task. In this work, we extract the underlying thermodynamic and dynamic contributions to the folding rate of an unusually slow-folding helical DNA-binding domain, PurR, which shares the characteristics of ultrafast downhill-folding proteins but nonetheless appears to exhibit an apparent two-state equilibrium. We combine equilibrium spectroscopy, temperature-viscosity-dependent kinetics, statistical mechanical modeling, and coarse-grained simulations to show that the conformational behavior of PurR is highly heterogeneous characterized by a large spread in melting temperatures, marginal thermodynamic barriers, and populated partially structured states. PurR appears to be at the threshold of disorder arising from frustrated electrostatics and weak packing that in turn slows down folding due to a shallow, bumpy landscape and not due to large thermodynamic barriers or strong internal friction. Our work highlights how a strong temperature dependence on the pre-exponential could signal a shallow landscape and not necessarily a slow-folding diffusion coefficient, thus determining the folding timescales of even millisecond folding proteins and hints at possible structural origins for the shallow landscape.

Subcellular compartmentalization of the plant antioxidant system: an integrated overview

The antioxidant system (AOS) maintains the optimal concentration of reactive oxygen species (ROS) in a cell and protects it against oxidative stress. In plants, the AOS consists of seven main classes of antioxidant enzymes, low-molecular antioxidants (e.g., ascorbate, glutathione, and their oxidized forms) and thioredoxin/glutaredoxin systems which can serve as reducing agents for antioxidant enzymes. The number of genes encoding AOS enzymes varies between classes, and same class enzymes encoded by different gene copies may have different subcellular localizations, functional loads and modes of evolution. These facts hereafter reinforce the complex nature of AOS regulation and functioning. Further studies can describe new trends in the behavior and functioning of systems components, and provide new fundamental knowledge about systems regulation. The system is revealed to have a lot of interactions and interplay pathways between its components at the subcellular level (antioxidants, enzymes, ROS level, and hormonal and transcriptional regulation). These facts should be taken into account in further studies during the AOS modeling by describing the main pathways of generating and utilizing ROS, as well as the associated signaling processes and regulation of the system on cellular and organelle levels, which is a complicated and ambitious task. Another objective for studying the phenomenon of the AOS is related to the influence of cell dynamics and circadian rhythms on it. Therefore, the AOS requires an integrated and multi-level approach to study. We focused this review on the existing scientific background and experimental data used for the systems biology research of the plant AOS.

Harnessing Conformational Plasticity to Generate Designer Enzymes

Recent years have witnessed an explosion of interest in understanding the role of conformational dynamics both in the evolution of new enzymatic activities from existing enzymes and in facilitating the emergence of enzymatic activity de novo on scaffolds that were previously non-catalytic. There are also an increasing number of examples in the literature of targeted engineering of conformational dynamics being successfully used to alter enzyme selectivity and activity. Despite the obvious importance of conformational dynamics to both enzyme function and evolvability, many (although not all) computational design approaches still focus either on pure sequence-based approaches or on using structures with limited flexibility to guide the design. However, there exist a wide variety of computational approaches that can be (re)purposed to introduce conformational dynamics as a key consideration in the design process. Coupled with laboratory evolution and more conventional existing sequence- and structure-based approaches, these techniques provide powerful tools for greatly expanding the protein engineering toolkit. This Perspective provides an overview of evolutionary studies that have dissected the role of conformational dynamics in facilitating the emergence of novel enzymes, as well as advances in computational approaches that allow one to target conformational dynamics as part of enzyme design. Harnessing conformational dynamics in engineering studies is a powerful paradigm with which to engineer the next generation of designer biocatalysts.

Allostery and Epistasis: Emergent Properties of Anisotropic Networks

Understanding the underlying mechanisms behind protein allostery and non-additivity of substitution outcomes (i.e., epistasis) is critical when attempting to predict the functional impact of mutations, particularly at non-conserved sites. In an effort to model these two biological properties, we extend the framework of our metric to calculate dynamic coupling between residues, the Dynamic Coupling Index (DCI) to two new metrics: (i) EpiScore, which quantifies the difference between the residue fluctuation response of a functional site when two other positions are perturbed with random Brownian kicks simultaneously versus individually to capture the degree of cooperativity of these two other positions in modulating the dynamics of the functional site and (ii) DCI_asym, which measures the degree of asymmetry between the residue fluctuation response of two sites when one or the other is perturbed with a random force. Applied to four independent systems, we successfully show that EpiScore and DCI_asym can capture important biophysical properties in dual mutant substitution outcomes. We propose that allosteric regulation and the mechanisms underlying non-additive amino acid substitution outcomes (i.e., epistasis) can be understood as emergent properties of an anisotropic network of interactions where the inclusion of the full network of interactions is critical for accurate modeling. Consequently, mutations which drive towards a new function may require a fine balance between functional site asymmetry and strength of dynamic coupling with the functional sites. These two tools will provide mechanistic insight into both understanding and predicting the outcome of dual mutations.

Molecular origins of folding rate differences in the thioredoxin family

Thioredoxins are a family of conserved oxidoreductases responsible for maintaining redox balance within cells. They have also served as excellent model systems for protein design and engineering studies particularly through ancestral sequence reconstruction methods. The recent work by Gamiz-Arco et al. [Biochem J (2019) 476, 3631–3647] answers fundamental questions on how specific sequence differences can contribute to differences in folding rates between modern and ancient thioredoxins but also among a selected subset of modern thioredoxins. They surprisingly find that rapid unassisted folding, a feature of ancient thioredoxins, is not conserved in the modern descendants suggestive of co-evolution of better folding machinery that likely enabled the accumulation of mutations that slow-down folding. The work thus provides an interesting take on the expected folding-stability-function constraint while arguing for additional factors that contribute to sequence evolution and hence impact folding efficiency.

The Role of Conformational Dynamics and Allostery in Modulating Protein Evolution

Advances in sequencing techniques and statistical methods have made it possible not only to predict sequences of ancestral proteins but also to identify thousands of mutations in the human exome, some of which are disease associated. These developments have motivated numerous theories and raised many questions regarding the fundamental principles behind protein evolution, which have been traditionally investigated horizontally using the tip of the phylogenetic tree through comparative studies of extant proteins within a family. In this article, we review a vertical comparison of the modern and resurrected ancestral proteins. We focus mainly on the dynamical properties responsible for a protein's ability to adapt new functions in response to environmental changes. Using the Dynamic Flexibility Index and the Dynamic Coupling Index to quantify the relative flexibility and dynamic coupling at a site-specific, single-amino-acid level, we provide evidence that the migration of hinges, which are often functionally critical rigid sites, is a mechanism through which proteins can rapidly evolve. Additionally, we show that disease-associated mutations in proteins often result in flexibility changes even at positions distal from mutational sites, particularly in the modulation of active site dynamics.

Intrinsic dynamics is evolutionarily optimized to enable allosteric behavior

Allosteric behavior is central to the function of many proteins, enabling molecular machinery, metabolism, signaling and regulation. Recent years have shown that the intrinsic dynamics of allosteric proteins defined by their 3-dimensional architecture or by the topology of inter-residue contacts favors cooperative motions that bear close similarity to structural changes they undergo during their allosteric actions. These conformational motions are usually driven by energetically favorable or soft modes at the low frequency end of the mode spectrum, and they are evolutionarily conserved among orthologs. These observations brought into light evolutionary adaptation mechanisms that help maintain, optimize or regulate allosteric behavior as the evolution from bacterial to higher organisms introduces sequential heterogeneities and structural complexities.

Functional analyses of ancestral thioredoxins provide insights into their evolutionary history

Thioredoxin (Trx) is a conserved, cytosolic reductase in all known organisms. The enzyme receives two electrons from NADPH via thioredoxin reductase (TrxR) and passes them on to multiple cellular reductases via disulfide exchange. Despite the ubiquity of thioredoxins in all taxa, little is known about the functions of resurrected ancestral thioredoxins in the context of a modern mesophilic organism. Here, we report on functional in vitro and in vivo analyses of seven resurrected Precambrian thioredoxins, dating back 1–4 billion years, in the Escherichia coli cytoplasm. Using synthetic gene constructs for recombinant expression of the ancestral enzymes, along with thermodynamic and kinetic assays, we show that all ancestral thioredoxins, as today's thioredoxins, exhibit strongly reducing redox potentials, suggesting that thioredoxins served as catalysts of cellular reduction reactions from the beginning of evolution, even before the oxygen catastrophe. A detailed, quantitative characterization of their interactions with the electron donor TrxR from Escherichia coli and the electron acceptor methionine sulfoxide reductase, also from E. coli, strongly hinted that thioredoxins and thioredoxin reductases co-evolved and that the promiscuity of thioredoxins toward downstream electron acceptors was maintained during evolution. In summary, our findings suggest that thioredoxins evolved high specificity for their sole electron donor TrxR while maintaining promiscuity to their multiple electron acceptors.

Protein ensembles link genotype to phenotype

Classically, phenotype is what is observed, and genotype is the genetic makeup. Statistical studies aim to project phenotypic likelihoods of genotypic patterns. The traditional genotype-to-phenotype theory embraces the view that the encoded protein shape together with gene expression level largely determines the resulting phenotypic trait. Here, we point out that the molecular biology revolution at the turn of the century explained that the gene encodes not one but ensembles of conformations, which in turn spell all possible gene-associated phenotypes. The significance of a dynamic ensemble view is in understanding the linkage between genetic change and the gained observable physical or biochemical characteristics. Thus, despite the transformative shift in our understanding of the basis of protein structure and function, the literature still commonly relates to the classical genotype-phenotype paradigm. This is important because an ensemble view clarifies how even seemingly small genetic alterations can lead to pleiotropic traits in adaptive evolution and in disease, why cellular pathways can be modified in monogenic and polygenic traits, and how the environment may tweak protein function.

Mutations Utilize Dynamic Allostery to Confer Resistance in TEM-1 β-lactamase

β-lactamases are enzymes produced by bacteria to hydrolyze β-lactam antibiotics as a common mechanism of resistance. Evolution in such enzymes has been rendering a wide variety of antibiotics impotent, therefore posing a major threat. Clinical and in vitro studies of evolution in TEM-1 β-lactamase have revealed a large number of single point mutations that are responsible for driving resistance to antibiotics and/or inhibitors. The distal locations of these mutations from the active sites suggest that these allosterically modulate the antibiotic resistance. We investigated the effects of resistance driver mutations on the conformational dynamics of the enzyme to provide insights about the mechanism of their long-distance interactions. Through all-atom molecular dynamics (MD) simulations, we obtained the dynamic flexibility profiles of the variants and compared those with that of the wild type TEM-1. While the mutational sites in the variants did not have any direct van der Waals interactions with the active site position S70 and E166, we observed a change in the flexibility of these sites, which play a very critical role in hydrolysis. Such long distance dynamic interactions were further confirmed by dynamic coupling index (DCI) analysis as the sites involved in resistance driving mutations exhibited high dynamic coupling with the active sites. A more exhaustive dynamic analysis, using a selection pressure for ampicillin and cefotaxime resistance on all possible types of substitutions in the amino acid sequence of TEM-1, further demonstrated the observed mechanism. Mutational positions that play a crucial role for the emergence of resistance to new antibiotics exhibited high dynamic coupling with the active site irrespective of their locations. These dynamically coupled positions were neither particularly rigid nor particularly flexible, making them more evolvable positions. Nature utilizes these sites to modulate the dynamics of the catalytic sites instead of mutating the highly rigid positions around the catalytic site.

Biotechnological and protein-engineering implications of ancestral protein resurrection

Approximations to the sequences of ancestral proteins can be derived from the sequences of their modern descendants. Proteins encoded by such reconstructed sequences can be prepared in the laboratory and subjected to experimental scrutiny. These 'resurrected' ancestral proteins often display remarkable properties, reflecting ancestral adaptations to intra-cellular and extra-cellular environments that differed from the environments hosting modern/extant proteins. Recent experimental and computational work has specifically discussed high stability, substrate and catalytic promiscuity, conformational flexibility/diversity and altered patterns of interaction with other sub-cellular components. In this review, we discuss these remarkable properties as well as recent attempts to explore their biotechnological and protein-engineering potential.