Conservation of Dynamics Associated with Biological Function in an Enzyme Superfamily

Ancestral ribonucleases back in motion for evolutionary-dynamics guided protein design

The dynamics behavior of a protein is essential for its functionality. Here, Doucet et al. demonstrate how the evolutionary analysis of conformational pathways within a protein family serves to identify common core scaffolds that accommodate branch-specific functional regions controlled by flexibility switches, offering a model for evolutionary-dynamics based protein design.

Ancestral sequence reconstruction dissects structural and functional differences among eosinophil ribonucleases

Evolutionarily conserved structural folds can give rise to diverse biological functions, yet predicting atomic-scale interactions that contribute to the emergence of novel activities within such folds remains challenging. Pancreatic-type ribonucleases illustrate this complexity, sharing a core structure that has evolved to accommodate varied functions. In this study, we used ancestral sequence reconstruction to probe evolutionary and molecular determinants that distinguish biological activities within eosinophil members of the RNase 2/3 subfamily. Our investigation unveils functional, structural, and dynamical behaviors that differentiate the evolved ancestral ribonuclease (AncRNase) from its contemporary eosinophil RNase orthologs. Leveraging the potential of ancestral reconstruction for protein engineering, we used AncRNase predictions to design a minimal 4-residue variant that transforms human RNase 2 into a chimeric enzyme endowed with the antimicrobial and cytotoxic activities of RNase 3 members. This work provides unique insights into mutational and evolutionary pathways governing structure, function, and conformational states within the eosinophil RNase subfamily, offering potential for targeted modulation of RNase-associated functions.

Sequence - dynamics - function relationships in protein tyrosine phosphatases

Protein tyrosine phosphatases (PTPs) are crucial regulators of cellular signaling. Their activity is regulated by the motion of a conserved loop, the WPD-loop, from a catalytically inactive open to a catalytically active closed conformation. WPD-loop motion optimally positions a catalytically critical residue into the active site, and is directly linked to the turnover number of these enzymes. Crystal structures of chimeric PTPs constructed by grafting parts of the WPD-loop sequence of PTP1B onto the scaffold of YopH showed WPD-loops in a wide-open conformation never previously observed in either parent enzyme. This wide-open conformation has, however, been observed upon binding of small molecule inhibitors to other PTPs, suggesting the potential of targeting it for drug discovery efforts. Here, we have performed simulations of both enzymes and show that there are negligible energetic differences in the chemical step of catalysis, but significant differences in the dynamical properties of the WPD-loop. Detailed interaction network analysis provides insight into the molecular basis for this population shift to a wide-open conformation. Taken together, our study provides insight into the links between loop dynamics and chemistry in these YopH variants specifically, and how WPD-loop dynamic can be engineered through modification of the internal protein interaction network.

A multi-layered computational structural genomics approach enhances domain-specific interpretation of Kleefstra syndrome variants in EHMT1

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretations of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within the domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism for variant dysfunction.

A Multi-Layered Computational Structural Genomics Approach Enhances Domain-Specific Interpretation of Kleefstra Syndrome Variants in EHMT1

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretation of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within this domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism of variant dysfunction.

Identifying structural and dynamic changes during the Biliverdin Reductase B catalytic cycle

Biliverdin Reductase B (BLVRB) is an NADPH-dependent reductase that catalyzes the reduction of multiple substrates and is therefore considered a critical cellular redox regulator. In this study, we sought to address whether both structural and dynamics changes occur between different intermediates of the catalytic cycle and whether these were relegated to just the active site or the entirety of the enzyme. Through X-ray crystallography, we determined the apo BLVRB structure for the first time, revealing subtle global changes compared to the holo structure and identifying the loss of a critical hydrogen bond that "clamps" the R78-loop over the coenzyme. Amide and Cα chemical shift perturbations were used to identify environmental and secondary structural changes between intermediates, with more distant global changes observed upon coenzyme binding compared to substrate interactions. NMR relaxation rate measurements provided insights into the dynamic behavior of BLVRB during the catalytic cycle. Specifically, the inherently dynamic R78-loop that becomes ordered upon coenzyme binding persists through the catalytic cycle while similar regions experience dynamic exchange. However, the dynamic exchange processes were found to differ through the catalytic cycle with several groups of residues exhibiting similar dynamic responses. Finally, both local and distal structural and dynamic changes occur within BLVRB that are dependent solely on the oxidative state of the coenzyme. Thus, through a comprehensive analysis here, this study revealed structural and dynamic alterations in BLVRB through its catalytic cycle that are not simply relegated to the active site, but instead, are allosterically coupled throughout the enzyme.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Conformational exchange divergence along the evolutionary pathway of eosinophil-associated ribonucleases

The evolutionary role of conformational exchange in the emergence and preservation of function within structural homologs remains elusive. While protein engineering has revealed the importance of flexibility in function, productive modulation of atomic-scale dynamics has only been achieved on a finite number of distinct folds. Allosteric control of unique members within dynamically diverse structural families requires a better appreciation of exchange phenomena. Here, we examined the functional and structural role of conformational exchange within eosinophil-associated ribonucleases. Biological and catalytic activity of various EARs was performed in parallel to mapping their conformational behavior on multiple timescales using NMR and computational analyses. Despite functional conservation and conformational seclusion to a specific domain, we show that EARs can display similar or distinct motional profiles, implying divergence rather than conservation of flexibility. Comparing progressively more distant enzymes should unravel how this subfamily has evolved new functions and/or altered their behavior at the molecular level.

Insights into the importance of WPD-loop sequence for activity and structure in protein tyrosine phosphatases

Protein tyrosine phosphatases (PTPs) possess a conserved mobile catalytic loop, the WPD-loop, which brings an aspartic acid into the active site where it acts as an acid/base catalyst. Prior experimental and computational studies, focused on the human enzyme PTP1B and the PTP from Yersinia pestis, YopH, suggested that loop conformational dynamics are important in regulating both catalysis and evolvability. We have generated a chimeric protein in which the WPD-loop of YopH is transposed into PTP1B, and eight chimeras that systematically restored the loop sequence back to native PTP1B. Of these, four chimeras were soluble and were subjected to detailed biochemical and structural characterization, and a computational analysis of their WPD-loop dynamics. The chimeras maintain backbone structural integrity, with somewhat slower rates than either wild-type parent, and show differences in the pH dependency of catalysis, and changes in the effect of Mg²⁺. The chimeric proteins' WPD-loops differ significantly in their relative stability and rigidity. The time required for interconversion, coupled with electrostatic effects revealed by simulations, likely accounts for the activity differences between chimeras, and relative to the native enzymes. Our results further the understanding of connections between enzyme activity and the dynamics of catalytically important groups, particularly the effects of non-catalytic residues on key conformational equilibria.

Capture of activated dioxygen intermediates at the copper-active site of a lytic polysaccharide monooxygenase

Metalloproteins perform a diverse array of redox-related reactions facilitated by the increased chemical functionality afforded by their metallocofactors. Lytic polysaccharide monooxygenases (LPMOs) are a class of copper-dependent enzymes that are responsible for the breakdown of recalcitrant polysaccharides via oxidative cleavage at the glycosidic bond. The activated copper-oxygen intermediates and their mechanism of formation remains to be established. Neutron protein crystallography which permits direct visualization of protonation states was used to investigate the initial steps of oxygen activation directly following active site copper reduction in Neurospora crassa LPMO9D. Herein, we cryo-trap an activated dioxygen intermediate in a mixture of superoxo and hydroperoxo states, and we identify the conserved second coordination shell residue His157 as the proton donor. Density functional theory calculations indicate that both superoxo and hydroperoxo active site states are stable. The hydroperoxo formed is potentially an early LPMO catalytic reaction intermediate or the first step in the mechanism of hydrogen peroxide formation in the absence of substrate. We observe that the N-terminal amino group of the copper coordinating His1 remains doubly protonated directly following molecular oxygen reduction by copper. Aided by molecular dynamics and mining minima free energy calculations we establish that the conserved second-shell His161 in MtPMO3* displays conformational flexibility in solution and that this flexibility is also observed, though to a lesser extent, in His157 of NcLPMO9D. The imidazolate form of His157 observed in our structure following oxygen intermediate protonation can be attributed to abolished His157 flexibility due steric hindrance in the crystal as well as the solvent-occluded active site environment due to crystal packing. A neutron crystal structure of NcLPMO9D at low pH further supports occlusion of the active site since His157 remains singly protonated even at acidic conditions.

Conserved conformational dynamics determine enzyme activity

Homologous enzymes often exhibit different catalytic rates despite a fully conserved active site. The canonical view is that an enzyme sequence defines its structure and function and, more recently, that intrinsic protein dynamics at different time scales enable and/or promote catalytic activity. Here, we show that, using the protein tyrosine phosphatase PTP1B, residues surrounding the PTP1B active site promote dynamically coordinated chemistry necessary for PTP1B function. However, residues distant to the active site also undergo distinct intermediate time scale dynamics and these dynamics are correlated with its catalytic activity and thus allow for different catalytic rates in this enzyme family. We identify these previously undetected motions using coevolutionary coupling analysis and nuclear magnetic resonance spectroscopy. Our findings strongly indicate that conserved dynamics drives the enzymatic activity of the PTP family. Characterization of these conserved dynamics allows for the identification of novel regulatory elements (therapeutic binding pockets) that can be leveraged for the control of enzymes.

In Search of a Dynamical Vocabulary: A Pipeline to Construct a Basis of Shared Traits in Large-Scale Motions of Proteins

The paradigmatic sequence–structure–dynamics–function relation in proteins is currently well established in the scientific community; in particular, a large effort has been made to probe the first connection, indeed providing convincing evidence of its strength and rationalizing it in a quantitative and general framework. In contrast, however, the role of dynamics as a link between structure and function has eluded a similarly clear-cut verification and description. In this work, we propose a pipeline aimed at building a basis for the quantitative characterization of the large-scale dynamics of a set of proteins, starting from the sole knowledge of their native structures. The method hinges on a dynamics-based clusterization, which allows a straightforward comparison with structural and functional protein classifications. The resulting basis set, obtained through the application to a group of related proteins, is shown to reproduce the salient large-scale dynamical features of the dataset. Most interestingly, the basis set is shown to encode the fluctuation patterns of homologous proteins not belonging to the initial dataset, thus highlighting the general applicability of the pipeline used to build it.

Mechanism of nucleotide discrimination by the translesion synthesis polymerase Rev1

Rev1 is a translesion DNA synthesis (TLS) polymerase involved in the bypass of adducted-guanine bases and abasic sites during DNA replication. During damage bypass, Rev1 utilizes a protein-template mechanism of DNA synthesis, where the templating DNA base is evicted from the Rev1 active site and replaced by an arginine side chain that preferentially binds incoming dCTP. Here, we utilize X-ray crystallography and molecular dynamics simulations to obtain structural insight into the dCTP specificity of Rev1. We show the Rev1 R324 protein-template forms sub-optimal hydrogen bonds with incoming dTTP, dGTP, and dATP that prevents Rev1 from adopting a catalytically competent conformation. Additionally, we show the Rev1 R324 protein-template forms optimal hydrogen bonds with incoming rCTP. However, the incoming rCTP adopts an altered sugar pucker, which prevents the formation of a catalytically competent Rev1 active site. This work provides novel insight into the mechanisms for nucleotide discrimination by the TLS polymerase Rev1.

Integrating dynamics into enzyme engineering

Enzyme engineering has become a widely adopted practice in research labs and industry. In parallel, the past decades have seen tremendous strides in characterizing the dynamics of proteins, using a growing array of methodologies. Importantly, links have been established between the dynamics of proteins and their function. Characterizing the dynamics of an enzyme prior to, and following, its engineering is beginning to inform on the potential of 'dynamic engineering', i.e. the rational modification of protein dynamics to alter enzyme function. Here we examine the state of knowledge at the intersection of enzyme engineering and protein dynamics, describe current challenges and highlight pioneering work in the nascent area of dynamic engineering.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Adaptability and specificity: how do proteins balance opposing needs to achieve function?

Many proteins select from a small repertoire of 3-dimensional folds retained over evolutional timescales and recruited for different functions, with changes in local structure and sequence to enable specificity. Recent studies have revealed the evolutionary constraints on protein dynamics to achieve function. The significance of protein dynamics in simultaneously satisfying conformational flexibility/malleability and stability/precision requirements becomes clear upon dissecting the spectrum of equilibrium motions accessible to fold families. Accessibility to highly conserved global modes of motions shared by family members, to low-to-intermediate-frequency modes that distinguish subfamilies and confer specificity, and to conserved high-frequency modes ensuring chemical precision and core stability underlies functional specialization while exploiting highly versatile folds. These design principles are illustrated for the family of PDZ domains.

Enzyme dynamics: Looking beyond a single structure

Conventional understanding of how enzymes function strongly emphasizes the role of structure. However, increasing evidence clearly indicates that enzymes do not remain fixed or operate exclusively in or close to their native structure. Different parts of the enzyme (from individual residues to full domains) undergo concerted motions on a wide range of time-scales, including that of the catalyzed reaction. Information obtained on these internal motions and conformational fluctuations has so far uncovered and explained many aspects of enzyme mechanisms, which could not have been understood from a single structure alone. Although there is wide interest in understanding enzyme dynamics and its role in catalysis, several challenges remain. In addition to technical difficulties, the vast majority of investigations are performed in dilute aqueous solutions, where conditions are significantly different than the cellular milieu where a large number of enzymes operate. In this review, we discuss recent developments, several challenges as well as opportunities related to this topic. The benefits of considering dynamics as an integral part of the enzyme function can also enable new means of biocatalysis, engineering enzymes for industrial and medicinal applications.

Allosteric pluripotency as revealed by protein kinase A

The functional response of a signaling system to an allosteric stimulus often depends on subcellular conditions, a phenomenon known as pluripotent allostery. For example, a single allosteric modulator, Rp-cAMPS, of the prototypical protein kinase A (PKA) switches from antagonist to agonist depending on MgATP levels. However, the mechanism underlying such pluripotent allostery has remained elusive for decades. Using nuclear magnetic resonance spectroscopy, ensemble models, kinase assays, and molecular dynamics simulations, we show that allosteric pluripotency arises from surprisingly divergent responses of highly homologous tandem domains. The differential responses perturb domain-domain interactions and remodel the free-energy landscape of inhibitory excited states sampled by the regulatory subunit of PKA. The resulting activation threshold values are comparable to the effective free energy of regulatory and catalytic subunit binding, which depends on metabolites, substrates, and mutations, explaining pluripotent allostery and warranting a general redefinition of allosteric targets to include specific subcellular environments.

Insights into Structural and Dynamical Changes Experienced by Human RNase 6 upon Ligand Binding

Ribonuclease 6 (RNase 6) is one of eight catalytically active human pancreatic-type RNases that belong to a superfamily of rapidly evolving enzymes. Like some of its human homologues, RNase 6 exhibits host defense properties such as antiviral and antibacterial activities. Recently solved crystal structures of this enzyme in its nucleotide-free form show the conservation of the prototypical kidney-shaped fold preserved among vertebrate RNases, in addition to revealing the presence of a unique secondary active site. In this study, we determine the structural and conformational properties experienced by RNase 6 upon binding to substrate and product analogues. We present the first crystal structures of RNase 6 bound to a nucleotide ligand (adenosine 5'-monophosphate), in addition to RNase 6 bound to phosphate ions. While the enzyme preserves B2 subsite ligand preferences, our results show a lack of typical B2 subsite interactions normally observed in homologous ligand-bound RNases. A comparison of the dynamical properties of RNase 6 in its apo-, substrate-, and product-bound states highlight the unique dynamical properties experienced on time scales ranging from nano- to milliseconds. Overall, our results confirm the specific evolutionary adaptation of RNase 6 relative to its unique catalytic and biological activities.

How good are comparative models in the understanding of protein dynamics?

The 3D structure of a protein is essential to understand protein dynamics. If experimentally determined structure is unavailable, comparative models could be used to infer dynamics. However, the effectiveness of comparative models, compared to experimental structures, in inferring dynamics is not clear. To address this, we compared dynamics features of ~800 comparative models with their crystal structures using normal mode analysis. Average similarity in magnitude, direction, and correlation of residue motions is >0.8 (where value 1 is identical) indicating that the dynamics of models and crystal structures are highly similar. Accuracy of 3D structure and dynamics is significantly higher for models built on multiple and/or high sequence identity templates (>40%). Three-dimensional (3D) structure and residue fluctuations of models are closer to that of crystal structures than to templates (TM score 0.9 vs 0.7 and square inner product 0.92 vs 0.88). Furthermore, long-range molecular dynamics simulations on comparative models of RNase 1 and Angiogenin showed significant differences in the conformational sampling of conserved active-site residues that characterize differences in their activity levels. Similar analyses on two EGFR kinase variant models highlight the effect of mutations on the functional state-specific αC helix motions and these results corroborate with the previous experimental observations. Thus, our study adds confidence to the use of comparative models in understanding protein dynamics.

The functional diversity of structural disorder in plant proteins

Structural disorder in proteins is a widespread feature distributed in all domains of life, particularly abundant in eukaryotes, including plants. In these organisms, intrinsically disordered proteins (IDPs) perform a diversity of functions, participating as integrators of signaling networks, in transcriptional and post-transcriptional regulation, in metabolic control, in stress responses and in the formation of biomolecular condensates by liquid-liquid phase separation. Their roles impact the perception, propagation and control of various developmental and environmental cues, as well as the plant defense against abiotic and biotic adverse conditions. In this review, we focus on primary processes to exhibit a broad perspective of the relevance of IDPs in plant cell functions. The information here might help to incorporate this knowledge into a more dynamic view of plant cells, as well as open more questions and promote new ideas for a better understanding of plant life.

Patterns of Dynamics Comprise a Conserved Evolutionary Trait

The importance of protein dynamics in function may suggest an evolutionary selection on large-scale protein motions. Here we systematically studied the dynamic characteristics in 2221 protein domains (58477 sequences) of the Pfam database. We defined the patterns of dynamics (PODs) based on the estimated NMR order parameters and the predicted degree of disorder, and found a significant correlation between them in families of both structured and disordered protein domains. We demonstrate that conservation of dynamic patterns frequently exceeds conservation of sequence and is comparable to the patterns of hydropathy and nonspecific interaction potential. Similarity of dynamic patterns is weakly correlated to structure similarity and to the degree of disorder. We illustrate that POD alignments could be applied to sequentially divergent or intrinsically disordered regions. We propose that patterns of dynamics comprise a conserved evolutionary trait, which could be used to infer evolutionary relationships as an alternative to sequence and structure.

Evolutionary Trends in RNA Base Selectivity Within the RNase A Superfamily

There is a growing interest in the pharmaceutical industry to design novel tailored drugs for RNA targeting. The vertebrate-specific RNase A superfamily is nowadays one of the best characterized family of enzymes and comprises proteins involved in host defense with specific cytotoxic and immune-modulatory properties. We observe within the family a structural variability at the substrate-binding site associated to a diversification of biological properties. In this work, we have analyzed the enzyme specificity at the secondary base binding site. Towards this end, we have performed a kinetic characterization of the canonical RNase types together with a molecular dynamic simulation of selected representative family members. The RNases' catalytic activity and binding interactions have been compared using UpA, UpG and UpI dinucleotides. Our results highlight an evolutionary trend from lower to higher order vertebrates towards an enhanced discrimination power of selectivity for adenine respect to guanine at the secondary base binding site (B2). Interestingly, the shift from guanine to adenine preference is achieved in all the studied family members by equivalent residues through distinct interaction modes. We can identify specific polar and charged side chains that selectively interact with donor or acceptor purine groups. Overall, we observe selective bidentate polar and electrostatic interactions: Asn to N1/N6 and N6/N7 adenine groups in mammals versus Glu/Asp and Arg to N1/N2, N1/O6 and O6/N7 guanine groups in non-mammals. In addition, kinetic and molecular dynamics comparative results on UpG versus UpI emphasize the main contribution of Glu/Asp interactions to N1/N2 group for guanine selectivity in lower order vertebrates. A close inspection at the B2 binding pocket also highlights the principal contribution of the protein ß6 and L4 loop regions. Significant differences in the orientation and extension of the L4 loop could explain how the same residues can participate in alternative binding modes. The analysis suggests that within the RNase A superfamily an evolution pressure has taken place at the B2 secondary binding site to provide novel substrate-recognition patterns. We are confident that a better knowledge of the enzymes' nucleotide recognition pattern would contribute to identify their physiological substrate and eventually design applied therapies to modulate their biological functions.

Engineering Dynamic Surface Peptide Networks on ButyrylcholinesteraseG117H for Enhanced Organophosphosphorus Anticholinesterase Catalysis

The single residue mutation of butyrylcholinesterase (BChE_G117H) hydrolyzes a number of organophosphosphorus (OP) anticholinesterases. Whereas other BChE active site/proximal mutations have been investigated, none are sufficiently active to be prophylactically useful. In a fundamentally different computer simulations driven strategy, we identified a surface peptide loop (residues 278-285) exhibiting dynamic motions during catalysis and modified it via residue insertions. We evaluated these loop mutants using computer simulations, substrate kinetics, resistance to inhibition, and enzyme reactivation assays using both the choline ester and OP substrates. A slight but significant increase in reactivation was noted with paraoxon with one of the mutants, and changes in K_M and catalytic efficiency were noted in others. Simulations suggested weaker interactions between OP versus choline substrates and the active site of all engineered versions of the enzyme. The results indicate that an improvement of OP anticholinesterase hydrolysis through surface loop engineering may be a more effective strategy in an enzyme with higher intrinsic OP compound hydrolase activity.

Nucleotide substrate binding characterization in human pancreatic-type ribonucleases

Human genome contains a group of more than a dozen similar genes with diverse biological functions including antiviral, antibacterial and angiogenesis activities. The characterized gene products of this group show significant sequence similarity and a common structural fold associated with binding and cleavage of ribonucleic acid (RNA) substrates. Therefore, these proteins have been categorized as members of human pancreatic-type ribonucleases (hRNases). hRNases differ in cell/tissue localization and display distinct substrate binding preferences and a wide range of ribonucleolytic catalytic efficiencies. Limited information is available about structural and dynamical properties that influence this diversity among these homologous RNases. Here, we use computer simulations to characterize substrate interactions, electrostatics and dynamical properties of hRNases 1-7 associated with binding to two nucleotide substrates (ACAC and AUAU). Results indicate that even with complete conservation of active-site catalytic triad associated with ribonucleolytic activity, these enzymes show significant differences in substrate interactions. Detailed characterization suggests that in addition to binding site electrostatic and van der Waals interactions, dynamics of distal regions may also play a role in binding. Another key insight is that a small difference in temperature of 300 K (used in experimental studies) and 310 K (physiological temperature) shows significant changes in enzyme-substrate interactions.

gRINN: a tool for calculation of residue interaction energies and protein energy network analysis of molecular dynamics simulations

Atomistic molecular dynamics (MD) simulations generate a wealth of information related to the dynamics of proteins. If properly analyzed, this information can lead to new insights regarding protein function and assist wet-lab experiments. Aiming to identify interactions between individual amino acid residues and the role played by each in the context of MD simulations, we present a stand-alone software called gRINN (get Residue Interaction eNergies and Networks). gRINN features graphical user interfaces (GUIs) and a command-line interface for generating and analyzing pairwise residue interaction energies and energy correlations from protein MD simulation trajectories. gRINN utilizes the features of NAMD or GROMACS MD simulation packages and automatizes the steps necessary to extract residue-residue interaction energies from user-supplied simulation trajectories, greatly simplifying the analysis for the end-user. A GUI, including an embedded molecular viewer, is provided for visualization of interaction energy time-series, distributions, an interaction energy matrix, interaction energy correlations and a residue correlation matrix. gRINN additionally offers construction and analysis of Protein Energy Networks, providing residue-based metrics such as degrees, betweenness-centralities, closeness centralities as well as shortest path analysis. gRINN is free and open to all users without login requirement at http://grinn.readthedocs.io.

The Structural Dynamics of Engineered β-Lactamases Vary Broadly on Three Timescales yet Sustain Native Function

Understanding the principles of protein dynamics will help guide engineering of protein function: altering protein motions may be a barrier to success or may be an enabling tool for protein engineering. The impact of dynamics on protein function is typically reported over a fraction of the full scope of motional timescales. If motional patterns vary significantly at different timescales, then only by monitoring motions broadly will we understand the impact of protein dynamics on engineering functional proteins. Using an integrative approach combining experimental and in silico methodologies, we elucidate protein dynamics over the entire span of fast to slow timescales (ps to ms) for a laboratory-engineered system composed of five interrelated β-lactamases: two natural homologs and three laboratory-recombined variants. Fast (ps-ns) and intermediate (ns-µs) dynamics were mostly conserved. However, slow motions (µs-ms) were few and conserved in the natural homologs yet were numerous and widely dispersed in their recombinants. Nonetheless, modified slow dynamics were functionally tolerated. Crystallographic B-factors from high-resolution X-ray structures were partly predictive of the conserved motions but not of the new slow motions captured in our solution studies. Our inspection of protein dynamics over a continuous range of timescales vividly illustrates the complexity of dynamic impacts of protein engineering as well as the functional tolerance of an engineered enzyme system to new slow motions.

Conserved Functional Dynamics: I Like to Move It, Move It!

The catalytic activities of enzymes are typically associated with their unique structures. However, structural dynamics may also play a functional role in enzyme catalysis. In support of this idea, Narayanan et al. (2018) show in this issue of Structure that enzymes are characterized by the presence of evolutionarily conserved dynamic traits.

A Biophysical Perspective on Enzyme Catalysis

Even after a century of investigation, our understanding of how enzymes work remains far from complete. In particular, several factors that enable enzymes to achieve high catalytic efficiencies remain only poorly understood. A number of theories have been developed, which propose or reaffirm that enzymes work as structural scaffolds, serving to bring together and properly orient the participants so that the reaction can proceed; therefore, leading to enzymes being viewed as only passive participants in the catalyzed reaction. A growing body of evidence shows that enzymes are not rigid structures but are constantly undergoing a wide range of internal motions and conformational fluctuations. In this Perspective, on the basis of studies from our group, we discuss the emerging biophysical model of enzyme catalysis that provides a detailed understanding of the interconnection among internal protein motions, conformational substates, enzyme mechanisms, and the catalytic efficiency of enzymes. For a number of enzymes, networks of conserved residues that extend from the surface of the enzyme all the way to the active site have been discovered. These networks are hypothesized to serve as pathways of energy transfer that enables thermodynamical coupling of the surrounding solvent with enzyme catalysis and play a role in promoting enzyme function. Additionally, the role of enzyme structure and electrostatic effects has been well acknowledged for quite some time. Collectively, the recent knowledge gained about enzyme mechanisms suggests that the conventional paradigm of enzyme structure encoding function is incomplete and needs to be extended to structure encodes dynamics, and together these enzyme features encode function including catalytic rate acceleration.

Differential Substrate Recognition by Maltose Binding Proteins Influenced by Structure and Dynamics

The genome of the hyperthermophile Thermotoga maritima contains three isoforms of maltose binding protein (MBP) that are high-affinity receptors for di-, tri-, and tetrasaccharides. Two of these proteins (tmMBP1 and tmMBP2) share significant sequence identity, approximately 90%, while the third (tmMBP3) shares less than 40% identity. MBP from Escherichia coli (ecMBP) shares 35% sequence identity with the tmMBPs. This subset of MBP isoforms offers an interesting opportunity to investigate the mechanisms underlying the evolution of substrate specificity and affinity profiles in a genome where redundant MBP genes are present. In this study, the X-ray crystal structures of tmMBP1, tmMBP2, and tmMBP3 are reported in the absence and presence of oligosaccharides. tmMBP1 and tmMBP2 have binding pockets that are larger than that of tmMBP3, enabling them to bind to larger substrates, while tmMBP1 and tmMBP2 also undergo substrate-induced hinge bending motions (∼52°) that are larger than that of tmMBP3 (∼35°). Small-angle X-ray scattering was used to compare protein behavior in solution, and computer simulations provided insights into dynamics of these proteins. Comparing quantitative protein-substrate interactions and dynamical properties of tmMBPs with those of the promiscuous ecMBP and disaccharide selective Thermococcus litoralis MBP provides insights into the features that enable selective binding. Collectively, the results provide insights into how the structure and dynamics of tmMBP homologues enable them to differentiate between a myriad of chemical entities while maintaining their common fold.

Mechanism of Selective Enzyme Inhibition through Uncompetitive Regulation of an Allosteric Agonist

Classical uncompetitive inhibitors are potent pharmacological modulators of enzyme function. Since they selectively target enzyme-substrate complexes (E:S), their inhibitory potency is amplified by increasing substrate concentrations. Recently, an unconventional uncompetitive inhibitor, called CE3F4R, was discovered for the exchange protein activated by cAMP isoform 1 (EPAC1). Unlike conventional uncompetitive inhibitors, CE3F4R is uncompetitive with respect to an allosteric effector, cAMP, as opposed to the substrate (i.e., CE3F4R targets the E:cAMP rather than the E:S complex). However, the mechanism of CE3F4R as an uncompetitive inhibitor is currently unknown. Here, we elucidate the mechanism of CE3F4R's action using NMR spectroscopy. Due to limited solubility and line broadening, which pose major challenges for traditional structural determination approaches, we resorted to a combination of protein- and ligand-based NMR experiments to comparatively analyze EPAC mutations, inhibitor analogs, and cyclic nucleotide derivatives that trap EPAC at different stages of activation. We discovered that CE3F4R binds within the EPAC cAMP-binding domain (CBD) at a subdomain interface distinct from the cAMP binding site, acting as a wedge that stabilizes a cAMP-bound mixed-intermediate. The mixed-intermediate includes attributes of both the apo/inactive and cAMP-bound/active states. In particular, the intermediate targeted by CE3F4R traps a CBD's hinge helix in its inactive conformation, locking EPAC into a closed domain topology that restricts substrate access to the catalytic domain. The proposed mechanism of action also explains the isoform selectivity of CE3F4R in terms of a single EPAC1 versus EPAC2 amino acid difference that destabilizes the active conformation of the hinge helix.

Ligand-Induced Variations in Structural and Dynamical Properties Within an Enzyme Superfamily

Enzyme catalysis is a complex process involving several steps along the reaction coordinates, including substrate recognition and binding, chemical transformation, and product release. Evidence continues to emerge linking the functional and evolutionary role of conformational exchange processes in optimal catalytic activity. Ligand binding changes the conformational landscape of enzymes, inducing long-range conformational rearrangements. Using functionally distinct members of the pancreatic ribonuclease superfamily as a model system, we characterized the structural and conformational changes associated with the binding of two mononucleotide ligands. By combining NMR chemical shift titration experiments with the chemical shift projection analysis (CHESPA) and relaxation dispersion experiments, we show that biologically distinct members of the RNase superfamily display discrete chemical shift perturbations upon ligand binding that are not conserved even in structurally related members. Amino acid networks exhibiting coordinated chemical shift displacements upon binding of the two ligands are unique to each of the RNases analyzed. Our results reveal the contribution of conformational rearrangements to the observed chemical shift perturbations. These observations provide important insights into the contribution of the different ligand binding specificities and effects of conformational exchange on the observed perturbations associated with ligand binding for functionally diverse members of the pancreatic RNase superfamily.