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

Dynamic Allostery: Evolution's Double-Edged Sword in Protein Function and Disease

Allostery is a core mechanism in biology that allows proteins to communicate and regulate activity over long structural distances. While classical models of allostery focus on conformational changes triggered by ligand binding, dynamic allostery-where protein function is modulated through alterations in thermal fluctuations without major conformational shifts-has emerged as a critical evolutionary mechanism. This review explores how evolution leverages dynamic allostery to fine-tune protein function through subtle mutations at distal sites, preserving core structural architecture while dramatically altering functional properties. Using a combination of computational approaches including Dynamic Flexibility Index (DFI), Dynamic Coupling Index (DCI), and vibrational density of states (VDOS) analysis, we demonstrate that functional adaptations in proteins often involve "hinge-shift" mechanisms, where redistribution of rigid and flexible regions modulates collective motions without changing the overall fold. This evolutionary principle is a double-edged sword: the same mechanisms that enable functional innovation also create vulnerabilities that can be exploited in disease states. Disease-associated variants frequently occur at positions highly coupled to functional sites despite being physically distant, forming Dynamic Allosteric Residue Couples (DARC sites). We demonstrate applications of these principles in understanding viral evolution, drug resistance, and capsid assembly dynamics. Understanding dynamic allostery provides critical insights into protein evolution and offers new avenues for therapeutic interventions targeting allosteric regulation.

Biosynthetic Pathways of Alaremycin and Its Derivative: Inhibitors of Porphobilinogen Synthase in Porphyrin Biosynthesis from Streptomyces sp. A012304

The antibiotic alaremycin (5-acetamido-4-oxo-5-hexenoic acid, 1), isolated from Streptomyces sp. A012304, structurally resembles 5-aminolevulinic acid (ALA), a precursor in porphyrin biosynthesis, and inhibits porphobilinogen synthase, the enzyme responsible for catalyzing the first common step of this pathway. In our previous study, the biosynthetic gene cluster responsible for alaremycin production-composed of almA (ALA synthase homologue), almB (N-acetyltransferase), almC (oxidoreductase), and almE (MFS-type transporter)-was identified, and a potential biosynthetic pathway was proposed. In this study, the biosynthetic pathway of 1 was confirmed by detecting intermediates using the liquid chromatography-mass spectrometry/MS (LC-MS/MS) analysis of extracts from Escherichia coli cells transformed with the biosynthetic genes, followed by in vitro reconstitution of the biosynthetic reactions using purified enzymes. AlmA catalyzed the condensation of l-serine and succinyl-CoA to produce 5-amino-6-hydroxy-4-oxohexanoic acid (2), AlmB catalyzed the N-acetylation of 2 to produce 5-acetamido-6-hydroxy-4-oxohexanoic acid (3), and AlmC catalyzed the dehydration of 3 to form 1. The AlmC-catalyzed reaction may involve a two-step mechanism including reduction by NADH and oxidation by Fe3+. Additionally, a novel derivative of 1 was identified in the culture broth of the producer strain, and its structure was determined as 5,6-dihydroalaremycin (5-acetamido-4-oxohexanoic acid, 4). It was revealed that 4 is synthesized via the same biosynthetic pathway but with AlmA and AlmB utilizing l-alanine as the amino acid precursor instead of l-serine.

How enzyme functions evolve: genetic, structural, and kinetic perspectives

Understanding the emergence or loss of enzyme functions comprises several approaches, such as genetic, structural, and kinetic studies. Promiscuous enzyme activities have been proposed as starting points for the emergence of novel enzyme functions, for example, through genetic models such as neofunctionalization and subfunctionalization. In both cases, neutral evolution would fix gene redundancy, critical in relaxing functional constraints and allowing specific mutations to drive innovation. The evolution of enzyme activities has a structural basis, with genetic mutations modifying the active site architecture, conformational dynamics, or interaction networks, which leads to the creation, enhancement, or restriction of enzyme functions where epistatic interactions are crucial. These structural changes impact the described kinetic mechanisms like ground-state stabilization (affinity), transition-state stabilization (catalysis), or a combination of both. Case studies across diverse enzyme families illustrate these principles, emphasizing the interplay between genetic, structural, and kinetic approaches. Finally, we discuss the importance of understanding evolutionary mechanisms and their impact on protein engineering and drug design for biomedical and industrial applications. However, these studies highlight that further experimental evolutionary data collection is necessary to enable the training of advanced machine learning models for use in biotechnological applications.

Tailoring industrial enzymes for thermostability and activity evolution by the machine learning-based iCASE strategy

The pursuit of obtaining enzymes with high activity and stability remains a grail in enzyme evolution due to the stability-activity trade-off. Here, we develop an isothermal compressibility-assisted dynamic squeezing index perturbation engineering (iCASE) strategy to construct hierarchical modular networks for enzymes of varying complexity. Molecular mechanism analysis elucidates that the peak of adaptive evolution is reached through a structural response mechanism among variants. Furthermore, this dynamic response predictive model using structure-based supervised machine learning is established to predict enzyme function and fitness, demonstrating robust performance across different datasets and reliable prediction for epistasis. The universality of the iCASE strategy is validated by four sorts of enzymes with different structures and catalytic types. This machine learning-based iCASE strategy provides guidance for future research on the fitness evolution of enzymes.

Global hinge sites of proteins as target sites for drug binding

Hinge sites of proteins play a key role in mediating conformational mechanics. Among them, those involved in the most collective modes of motion, also called global hinges, are of particular interest, as they support cooperative rearrangements that are often functional. Yet, the utility of targeting global hinges for modulating function remains to be established. We present here a systematic study of a series of proteins resolved in drug-bound forms to examine the probabilistic occurrence of spatial overlaps between hinge sites and drug-binding pockets. Our analysis reveals a high propensity of drug binding to hinge sites compared to random. Notably, one-third of currently approved drugs are colocalized with hinge sites. These mechanosensitive sites are predictable by simple models such as the Gaussian Network Model. Their targeting thus emerges as a viable strategy for developing a new class of drugs that would exploit and modulate the target proteins' intrinsic dynamics, and potentially alleviate drug-resistance when used in combination with orthosteric or allosteric drugs.

Engineering the native ensemble to tune protein function: Diverse mutational strategies and interlinked molecular mechanisms

Natural proteins are fragile entities, intrinsically sensitive to perturbations both at the level of sequence and their immediate environment. Here, we highlight the diverse strategies available for engineering function through mutations influencing backbone conformational entropy, charge-charge interactions, and in the loops and hinge regions, many of which are located far from the active site. It thus appears that there are potentially numerous ways to microscopically vary the identity of residues and the constituent interactions to tune function. Functional modulation could occur via changes in native-state stability, altered thermodynamic coupling extents within the folded structure, redistributed dynamics, or through modulation of the population of conformational substates. As these mechanisms are intrinsically linked and given the pervasive long-range effects of mutations, it is crucial to consider the interaction network as a whole and fully map the native conformational landscape to place mutational effects in the context of allostery and protein evolution.

Transition Path Sampling Based Free Energy Calculations of Evolution's Effect on Rates in β-Lactamase: The Contributions of Rapid Protein Dynamics to Rate

β-Lactamases are one of the primary enzymes responsible for antibiotic resistance and have existed for billions of years. The structural differences between a modern class A TEM-1 β-lactamase compared to a sequentially reconstructed Gram-negative bacteria β-lactamase are minor. Despite the similar structures and mechanisms, there are different functions between the two enzymes. We recently identified differences in dynamics effects that result from evolutionary changes that could potentially account for the increase in substrate specificity and catalytic rate. In this study, we used transition path sampling-based calculations of free energies to identify how evolutionary changes found between an ancestral β-lactamase, and its extant counterpart TEM-1 β-lactamase affect rate.

Dynamics-based protein network features accurately discriminate neutral and rheostat positions

In some proteins, a unique class of nonconserved positions is characterized by their ability to generate diverse functional outcomes through single amino acid substitutions. Due to their ability to tune protein function, accurately identifying such "rheostat" positions is crucial for protein design, for understanding the impact of mutations observed in humans, and for predicting the evolution of pathogen drug resistance. However, identifying rheostat positions has been challenging, due-in part-to the absence of a clear structural relationship with binding sites. In this study, experimental data from our previous study of the Escherichia coli lactose repressor protein (LacI) was used to identify rheostat positions for which mutations tune in vivo EC50 for the allosteric ligand "IPTG." We next used the rheostat assignments to test the hypothesis that rheostat positions have unique dynamic features that will enable their identification. To that end, we integrated all-atom molecular dynamics simulations with perturbation residue response analysis. Results first revealed distinct dynamic behavior in IPTG-bound LacI compared with apo LacI, which was consistent with IPTG's role as an allosteric inducer. Next, we used a variety of dynamic features to build a classification model that discriminates experimentally characterized rheostat positions in LacI from positions with other types of substitution outcomes. In parallel, we built a second classifier model based on the 3D structural "static" network features of LacI. In comparative studies, the dynamic model better identified rheostat positions that were >8 Å from the binding site. In summary, our study provides insights into the dynamic characteristics of rheostat positions and suggests that models built on dynamic features may be useful for predicting the locations of rheostat positions in a wide range of proteins.

Functional and Mechanistic Dissection of Protein Glutaminase PG3 and Its Rational Engineering for Enhanced Modification of Myofibrillar Proteins

Protein glutaminases (PG; EC = 3.5.1.44) are enzymes known for enhancing protein functionality. In this study, we cloned and expressed the gene chryb3 encoding protein glutaminase PG3, exhibiting 39.4 U/mg specific activity. Mature-PG3 featured a substrate channel surrounded by aromatic and hydrophobic amino acids at positions 38-45 and 78-84, with Val81 playing a pivotal role in substrate affinity. The dynamic opening and closing motions between Gly65, Thr66, and Cys164 at the catalytic cleft greatly influence substrate binding and product release. Redesigning catalytic pocket and cocatalytic region produced combinatorial mutant MT6 showing a 2.69-fold increase in specific activity and a 2.99-fold increase at t65 °C1/2. Furthermore, MT6 boosted fish myofibrillar protein (MP) solubility without NaCl. Key residues such as Thr3, Asn54, Val81, Tyr82, Asn107, and Ser108 were vital for PG3-myosin interaction, particularly Asn54 and Asn107. This study sheds light on the catalytic mechanism of PG3 and guided its rational engineering and utilization in low-salt fish MP product production.

The Evolution of the Acylation Mechanism in β-Lactamase and Rapid Protein Dynamics

β-Lactamases are a class of well-studied enzymes that are known to have existed since billions of years ago, starting as a defense mechanism to stave off competitors and are now enzymes responsible for antibiotic resistance. Using ancestral sequence reconstruction, it is possible to study the crystal structure of a laboratory resurrected 2-3 billion year-old β-lactamase. Comparing the ancestral enzyme to its modern counterpart, a TEM-1 β-lactamase, the structural changes are minor, and it is probable that dynamic effects play an important role in the evolution of function. We used molecular dynamics simulations and employed transition path sampling methods to identify the presence of rate-enhancing dynamics at the femtosecond level in both systems, found that these fast motions are more efficiently coordinated in the modern enzyme, and examined how specific dynamics can pinpoint evolutionary effects that are essential for improving enzymatic catalysis.

Friends and relatives: insight into conformational regulation from orthologues and evolutionary lineages using KIF and KIN

Noncovalent interaction networks provide a powerful means to represent and analyze protein structure. Such networks can represent both static structures and dynamic conformational ensembles. We have recently developed two tools for analyzing such interaction networks and generating hypotheses for protein engineering. Here, we apply these tools to the conformational regulation of substrate specificity in class A β-lactamases, particularly the evolutionary development from generalist to specialist catalytic function and how that can be recapitulated or reversed by protein engineering. These tools, KIF and KIN, generate a set of prioritized residues and interactions as targets for experimental protein engineering.

Evidence that the cold- and menthol-sensing functions of the human TRPM8 channel evolved separately

Transient receptor potential melastatin 8 (TRPM8) is a temperature- and menthol-sensitive ion channel that contributes to diverse physiological roles, including cold sensing and pain perception. Clinical trials targeting TRPM8 have faced repeated setbacks predominantly due to the knowledge gap in unraveling the molecular underpinnings governing polymodal activation. A better understanding of the molecular foundations between the TRPM8 activation modes may aid the development of mode-specific, thermal-neutral therapies. Ancestral sequence reconstruction was used to explore the origins of TRPM8 activation modes. By resurrecting key TRPM8 nodes along the human evolutionary trajectory, we gained valuable insights into the trafficking, stability, and function of these ancestral forms. Notably, this approach unveiled the differential emergence of cold and menthol sensitivity over evolutionary time, providing a fresh perspective on complex polymodal behavior. These studies provide a paradigm for understanding polymodal behavior in TRPM8 and other proteins with the potential to enhance our understanding of sensory receptor biology and pave the way for innovative therapeutic interventions.

Simultaneous enhancement of multiple functional properties using evolution-informed protein design

A major challenge in protein design is to augment existing functional proteins with multiple property enhancements. Altering several properties likely necessitates numerous primary sequence changes, and novel methods are needed to accurately predict combinations of mutations that maintain or enhance function. Models of sequence co-variation (e.g., EVcouplings), which leverage extensive information about various protein properties and activities from homologous protein sequences, have proven effective for many applications including structure determination and mutation effect prediction. We apply EVcouplings to computationally design variants of the model protein TEM-1 β-lactamase. Nearly all the 14 experimentally characterized designs were functional, including one with 84 mutations from the nearest natural homolog. The designs also had large increases in thermostability, increased activity on multiple substrates, and nearly identical structure to the wild type enzyme. This study highlights the efficacy of evolutionary models in guiding large sequence alterations to generate functional diversity for protein design applications.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 Spike (S) protein. With this approach, we first identified Candidate Adaptive Polymorphisms (CAPs) of the SARS-CoV-2 Spike protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

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.

Allosteric regulatory control in dihydrofolate reductase is revealed by dynamic asymmetry

We investigated the relationship between mutations and dynamics in Escherichia coli dihydrofolate reductase (DHFR) using computational methods. Our study focused on the M20 and FG loops, which are known to be functionally important and affected by mutations distal to the loops. We used molecular dynamics simulations and developed position-specific metrics, including the dynamic flexibility index (DFI) and dynamic coupling index (DCI), to analyze the dynamics of wild-type DHFR and compared our results with existing deep mutational scanning data. Our analysis showed a statistically significant association between DFI and mutational tolerance of the DHFR positions, indicating that DFI can predict functionally beneficial or detrimental substitutions. We also applied an asymmetric version of our DCI metric (DCIasym ) to DHFR and found that certain distal residues control the dynamics of the M20 and FG loops, whereas others are controlled by them. Residues that are suggested to control the M20 and FG loops by our DCIasym metric are evolutionarily nonconserved; mutations at these sites can enhance enzyme activity. On the other hand, residues controlled by the loops are mostly deleterious to function when mutated and are also evolutionary conserved. Our results suggest that dynamics-based metrics can identify residues that explain the relationship between mutation and protein function or can be targeted to rationally engineer enzymes with enhanced activity.

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.

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.

Second-Shell Residues Contribute to Catalysis by Predominately Preorganizing the Apo State in PafA

Residues beyond the first coordination shell are often observed to make considerable cumulative contributions in enzymes. Due to typically indirect perturbations of multiple physicochemical properties of the active site, however, their individual and specific roles in enzyme catalysis and disease-causing mutations remain difficult to predict and understand at the molecular level. Here we analyze the contributions of several second-shell residues in phosphate-irrepressible alkaline phosphatase of flavobacterium (PafA), a representative system as one of the most efficient enzymes. By adopting a multifaceted approach that integrates quantum-mechanical/molecular-mechanical free energy computations, molecular-mechanical molecular dynamics simulations, and density functional theory cluster model calculations, we probe the rate-limiting phosphoryl transfer step and structural properties of all relevant enzyme states. In combination with available experimental data, our computational results show that mutations of the studied second-shell residues impact catalytic efficiency mainly by perturbation of the apo state and therefore substrate binding, while they do not affect the ground state or alter the nature of phosphoryl transfer transition state significantly. Several second-shell mutations also modulate the active site hydration level, which in turn influences the energetics of phosphoryl transfer. These mechanistic insights also help inform strategies that may improve the efficiency of enzyme design and engineering by going beyond the current focus on the first coordination shell.

Simultaneous enhancement of multiple functional properties using evolution-informed protein design

Designing optimized proteins is important for a range of practical applications. Protein design is a rapidly developing field that would benefit from approaches that enable many changes in the amino acid primary sequence, rather than a small number of mutations, while maintaining structure and enhancing function. Homologous protein sequences contain extensive information about various protein properties and activities that have emerged over billions of years of evolution. Evolutionary models of sequence co-variation, derived from a set of homologous sequences, have proven effective in a range of applications including structure determination and mutation effect prediction. In this work we apply one of these models (EVcouplings) to computationally design highly divergent variants of the model protein TEM-1 β-lactamase, and characterize these designs experimentally using multiple biochemical and biophysical assays. Nearly all designed variants were functional, including one with 84 mutations from the nearest natural homolog. Surprisingly, all functional designs had large increases in thermostability and most had a broadening of available substrates. These property enhancements occurred while maintaining a nearly identical structure to the wild type enzyme. Collectively, this work demonstrates that evolutionary models of sequence co-variation (1) are able to capture complex epistatic interactions that successfully guide large sequence departures from natural contexts, and (2) can be applied to generate functional diversity useful for many applications in protein design.

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

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.

β-Lactams from the Ocean

The title of this essay is as much a question as it is a statement. The discovery of the β-lactam antibiotics-including penicillins, cephalosporins, and carbapenems-as largely (if not exclusively) secondary metabolites of terrestrial fungi and bacteria, transformed modern medicine. The antibiotic β-lactams inactivate essential enzymes of bacterial cell-wall biosynthesis. Moreover, the ability of the β-lactams to function as enzyme inhibitors is of such great medical value, that inhibitors of the enzymes which degrade hydrolytically the β-lactams, the β-lactamases, have equal value. Given this privileged status for the β-lactam ring, it is therefore a disappointment that the exemplification of this ring in marine secondary metabolites is sparse. It may be that biologically active marine β-lactams are there, and simply have yet to be encountered. In this report, we posit a second explanation: that the value of the β-lactam to secure an ecological advantage in the marine environment might be compromised by its close structural similarity to the β-lactones of quorum sensing. The steric and reactivity similarities between the β-lactams and the β-lactones represent an outside-of-the-box opportunity for correlating new structures and new enzyme targets for the discovery of compelling biological activities.

Achieving Ultrasensitive Chromogenic Probes for Rapid, Direct Detection of Carbapenemase-Producing Bacteria in Sputum

Carbapenemase-producing bacteria (CPB) stand as the most dangerous "superbugs" in the clinic. Rapid point-of-care (POC) detection of CPB in clinical samples is key to timely and effective infection management. We herein report the first ultrasensitive chromogenic probe that allows direct POC detection of CPB in clinical sputum samples at a sample-to-result time of less than 15 min. This chromogenic probe is modularly designed by conjugating the carbapenem core with a benzene derivative bearing an electronegativity-tunable substituent. Unexpectedly high sensitivity was achieved simply by choosing strong electron-withdrawing substituents, such as -N⁺(CH₃)₃, without resorting to complex molecular design. Through integrating the probes with a portable paper chip, 24 out of 80 clinical sputum samples from sepsis patients with lung infections were quickly diagnosed as CPB-positive, exhibiting 100% clinical sensitivity and specificity. This low-cost paper chip assay can be readily performed on-site, breaking through the dilemma of rapid CPB detection, especially in resource-limited settings.

DHFR Mutants Modulate Their Synchronized Dynamics with the Substrate by Shifting Hydrogen Bond Occupancies

Antibiotic resistance is a global health problem in which mutations occurring in functional proteins render drugs ineffective. The working mechanisms of the arising mutants are seldom apparent; a methodology to decipher these mechanisms systematically would render devising therapies to control the arising mutational pathways possible. Here we utilize C_α-C_β bond vector relaxations obtained from moderate length MD trajectories to determine conduits for functionality of the resistance conferring mutants of Escherichia coli dihydrofolate reductase. We find that the whole enzyme is synchronized to the motions of the substrate, irrespective of the mutation introducing gain-of-function or loss-of function. The total coordination of the motions suggests changes in the hydrogen bond dynamics with respect to the wild type as a possible route to determine and classify the mode-of-action of individual mutants. As a result, nine trimethoprim-resistant point mutations arising frequently in evolution experiments are categorized. One group of mutants that display the largest occurrence (L28R, W30G) work directly by modifying the dihydrofolate binding region. Conversely, W30R works indirectly by the formation of the E139-R30 salt bridge which releases energy resulting from tight binding by distorting the binding cavity. A third group (D27E, F153S, I94L) arising as single, resistance invoking mutants in evolution experiment trajectories allosterically and dynamically affects a hydrogen bonding motif formed at residues 59-69-71 which in turn modifies the binding site dynamics. The final group (I5F, A26T, R98P) consists of those mutants that have properties most similar to the wild type; these only appear after one of the other mutants is fixed on the protein structure and therefore display clear epistasis. Thus, we show that the binding event is governed by the entire enzyme dynamics while the binding site residues play gating roles. The adjustments made in the total enzyme in response to point mutations are what make quantifying and pinpointing their effect a hard problem. Here, we show that hydrogen bond dynamics recorded on sub-μs time scales provide the necessary fingerprints to decipher the various mechanisms at play.

Investigating the allosteric response of the PICK1 PDZ domain to different ligands with all-atom simulations

The PDZ family is comprised of small modular domains that play critical roles in the allosteric modulation of many cellular signaling processes by binding to the C-terminal tail of different proteins. As dominant modular proteins that interact with a diverse set of peptides, it is of particular interest to explore how different binding partners induce different allosteric effects on the same PDZ domain. Because the PICK1 PDZ domain can bind different types of ligands, it is an ideal test case to answer this question and explore the network of interactions that give rise to dynamic allostery. Here, we use all-atom molecular dynamics simulations to explore dynamic allostery in the PICK1 PDZ domain by modeling two PICK1 PDZ systems: PICK1 PDZ-DAT and PICK1 PDZ-GluR2. Our results suggest that ligand binding to the PICK1 PDZ domain induces dynamic allostery at the αA helix that is similar to what has been observed in other PDZ domains. We found that the PICK1 PDZ-ligand distance is directly correlated with both dynamic changes of the αA helix and the distance between the αA helix and βB strand. Furthermore, our work identifies a hydrophobic core between DAT/GluR2 and I35 as a key interaction in inducing such dynamic allostery. Finally, the unique interaction patterns between different binding partners and the PICK1 PDZ domain can induce unique dynamic changes to the PICK1 PDZ domain. We suspect that unique allosteric coupling patterns with different ligands may play a critical role in how PICK1 performs its biological functions in various signaling networks.

General Theory of Specific Binding: Insights from a Genetic-Mechano-Chemical Protein Model

Proteins need to selectively interact with specific targets among a multitude of similar molecules in the cell. However, despite a firm physical understanding of binding interactions, we lack a general theory of how proteins evolve high specificity. Here, we present such a model that combines chemistry, mechanics, and genetics and explains how their interplay governs the evolution of specific protein-ligand interactions. The model shows that there are many routes to achieving molecular discrimination-by varying degrees of flexibility and shape/chemistry complementarity-but the key ingredient is precision. Harder discrimination tasks require more collective and precise coaction of structure, forces, and movements. Proteins can achieve this through correlated mutations extending far from a binding site, which fine-tune the localized interaction with the ligand. Thus, the solution of more complicated tasks is enabled by increasing the protein size, and proteins become more evolvable and robust when they are larger than the bare minimum required for discrimination. The model makes testable, specific predictions about the role of flexibility and shape mismatch in discrimination, and how evolution can independently tune affinity and specificity. Thus, the proposed theory of specific binding addresses the natural question of "why are proteins so big?". A possible answer is that molecular discrimination is often a hard task best performed by adding more layers to the protein.

Directed reconstruction of a novel ancestral alcohol dehydrogenase featuring shifted pH-profile, enhanced thermostability and expanded substrate spectrum

Ancestral enzymes are promising for industrial biotechnology due to high stability and catalytic promiscuity. An effective protocol was developed for the directed resurrection of ancestral enzymes. Employing genome mining with diaryl alcohol dehydrogenase KpADH as the probe, descendant enzymes D10 and D11 were firstly identified. Then through ancestral sequence reconstruction, A64 was resurrected with a specific activity of 4.3 U·mg^-1. The optimum pH of A64 was 7.5, distinct from 5.5 of D10. The T15 50 and T_m values of A64 were 57.5 °C and 61.7 °C, significantly higher than those of the descendant counterpart. Substrate spectrum of A64 was quantitively characterized with a Shannon-Wiener index of 2.38, more expanded than D10, especially, towards bulky ketones in Group A and B. A64 also exhibited higher enantioselectivity. This study provides an effective protocol for constructing of ancestral enzymes and an efficient ancestral enzyme of industrial relevance for asymmetric synthesis of chiral alcohols.

The Statistical Trends of Protein Evolution: A Lesson from AlphaFold Database

The recent development of artificial intelligence provides us with new and powerful tools for studying the mysterious relationship between organism evolution and protein evolution. In this work, based on the AlphaFold Protein Structure Database (AlphaFold DB), we perform comparative analyses of the proteins of different organisms. The statistics of AlphaFold-predicted structures show that, for organisms with higher complexity, their constituent proteins will have larger radii of gyration, higher coil fractions, and slower vibrations, statistically. By conducting normal mode analysis and scaling analyses, we demonstrate that higher organismal complexity correlates with lower fractal dimensions in both the structure and dynamics of the constituent proteins, suggesting that higher functional specialization is associated with higher organismal complexity. We also uncover the topology and sequence bases of these correlations. As the organismal complexity increases, the residue contact networks of the constituent proteins will be more assortative, and these proteins will have a higher degree of hydrophilic-hydrophobic segregation in the sequences. Furthermore, by comparing the statistical structural proximity across the proteomes with the phylogenetic tree of homologous proteins, we show that, statistical structural proximity across the proteomes may indirectly reflect the phylogenetic proximity, indicating a statistical trend of protein evolution in parallel with organism evolution. This study provides new insights into how the diversity in the functionality of proteins increases and how the dimensionality of the manifold of protein dynamics reduces during evolution, contributing to the understanding of the origin and evolution of lives.

Slow Protein Dynamics Elicits New Enzymatic Functions by Means of Epistatic Interactions

Protein evolution depends on the adaptation of these molecules to different functional challenges. This occurs by tuning their biochemical, biophysical, and structural traits through the accumulation of mutations. While the role of protein dynamics in biochemistry is well recognized, there are limited examples providing experimental evidence of the optimization of protein dynamics during evolution. Here we report an NMR study of four variants of the CTX-M β-lactamases, in which the interplay of two mutations outside the active site enhances the activity against a cephalosporin substrate, ceftazidime. The crystal structures of these enzymes do not account for this activity enhancement. By using NMR, here we show that the combination of these two mutations increases the backbone dynamics in a slow timescale and the exposure to the solvent of an otherwise buried β-sheet. The two mutations located in this β-sheet trigger conformational changes in loops located at the opposite side of the active site. We postulate that the most active variant explores alternative conformations that enable binding of the more challenging substrate ceftazidime. The impact of the mutations in the dynamics is context-dependent, in line with the epistatic effect observed in the catalytic activity of the different variants. These results reveal the existence of a dynamic network in CTX-M β-lactamases that has been exploited in evolution to provide a net gain-of-function, highlighting the role of alternative conformations in protein evolution.

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.

Efficient Base-Catalyzed Kemp Elimination in an Engineered Ancestral Enzyme

The routine generation of enzymes with completely new active sites is a major unsolved problem in protein engineering. Advances in this field have thus far been modest, perhaps due, at least in part, to the widespread use of modern natural proteins as scaffolds for de novo engineering. Most modern proteins are highly evolved and specialized and, consequently, difficult to repurpose for completely new functionalities. Conceivably, resurrected ancestral proteins with the biophysical properties that promote evolvability, such as high stability and conformational diversity, could provide better scaffolds for de novo enzyme generation. Kemp elimination, a non-natural reaction that provides a simple model of proton abstraction from carbon, has been extensively used as a benchmark in de novo enzyme engineering. Here, we present an engineered ancestral β-lactamase with a new active site that is capable of efficiently catalyzing Kemp elimination. The engineering of our Kemp eliminase involved minimalist design based on a single function-generating mutation, inclusion of an extra polypeptide segment at a position close to the de novo active site, and sharply focused, low-throughput library screening. Nevertheless, its catalytic parameters (kcat/KM~2·105 M-1 s-1, kcat~635 s-1) compare favorably with the average modern natural enzyme and match the best proton-abstraction de novo Kemp eliminases that are reported in the literature. The general implications of our results for de novo enzyme engineering are discussed.

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.

A review on the immobilization of pepsin: A Lys-poor enzyme that is unstable at alkaline pH values

Pepsin is a protease used in many different applications, and in many instances, it is utilized in an immobilized form to prevent contamination of the reaction product. This enzyme has two peculiarities that make its immobilization complex. The first one is related to the poor presence of primary amino groups on its surface (just one Lys and the terminal amino group). The second one is its poor stability at alkaline pH values. Both features make the immobilization of this enzyme to be considered a complicated goal, as most of the immobilization protocols utilize primary amino groups for immobilization. This review presents some of the attempts to get immobilized pepsin biocatalyst and their applications. The high density of anionic groups (Asp and Glu) make the anion exchange of the enzyme simpler, but this makes many of the strategies utilized to immobilize the enzyme (e.g., amino-glutaraldehyde supports) more related to a mixed ion exchange/hydrophobic adsorption than to real covalent immobilization. Finally, we propose some possibilities that can permit not only the covalent immobilization of this enzyme, but also their stabilization via multipoint covalent attachment.

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.

The role of conformational dynamics on the activity of polymer-conjugated CalB in organic solvents

A perennial interest in enzyme catalysis has been expanding its applicability from aqueous phase where enzymes are naturally evolved to organic solvents in which the majority of industrial chemical synthesis...

The Interplay of Electrostatics and Chemical Positioning in the Evolution of Antibiotic Resistance in TEM β-Lactamases

The interplay of enzyme active site electrostatics and chemical positioning is important for understanding the origin(s) of enzyme catalysis and the design of novel catalysts. We reconstruct the evolutionary trajectory of TEM-1 β-lactamase to TEM-52 toward extended-spectrum activity to better understand the emergence of antibiotic resistance and to provide insights into the structure-function paradigm and noncovalent interactions involved in catalysis. Utilizing a detailed kinetic analysis and the vibrational Stark effect, we quantify the changes in rates and electric fields in the Michaelis and acyl-enzyme complexes for penicillin G and cefotaxime to ascertain the evolutionary role of electric fields to modulate function. These data are combined with MD simulations to interpret and quantify the substrate-dependent structural changes during evolution. We observe that this evolutionary trajectory utilizes a large preorganized electric field and substrate-dependent chemical positioning to facilitate catalysis. This governs the evolvability, substrate promiscuity, and protein fitness landscape in TEM β-lactamase antibiotic resistance.

Allosteric Control of Structural Mimicry and Mutational Escape in the SARS-CoV-2 Spike Protein Complexes with the ACE2 Decoys and Miniprotein Inhibitors: A Network-Based Approach for Mutational Profiling of Binding and Signaling

We developed a computational framework for comprehensive and rapid mutational scanning of binding energetics and residue interaction networks in the SARS-CoV-2 spike protein complexes. Using this approach, we integrated atomistic simulations and conformational landscaping of the SARS-CoV-2 spike protein complexes with ensemble-based mutational screening and network modeling to characterize mechanisms of structure-functional mimicry and resilience toward mutational escape by the ACE2 protein decoy and de novo designed miniprotein inhibitors. A detailed analysis of structural plasticity of the SARS-CoV-2 spike proteins obtained from atomistic simulations of conformational landscapes and sequence-based profiling of the disorder propensities revealed the intrinsically flexible regions that harbor key functional sites targeted by circulating variants. The conservation of collective dynamics in the SARS-CoV-2 spike protein complexes showed that mutational escape positions are important for modulation of functional motions and that mutational changes in these sites can alter allosteric interaction networks. Through mutational profiling of binding and allosteric propensities in the SARS-CoV-2 spike protein complexes, we identified the key binding and regulatory hotspots that collectively determine functional response and resilience of miniproteins to mutational variants. The results suggest that binding affinities and allosteric signatures of the SARS-CoV-2 complexes can be determined by dynamic crosstalk between structurally stable regulatory centers and conformationally adaptable allosteric hotspots that collectively control the resilience toward mutational escape. This may underlie a mechanism in which moderate perturbations in the mutational escape positions can induce global allosteric changes and alter functional protein response by modulating signaling in the residue interaction networks.

Computational strategies for protein conformational ensemble detection

Protein function is constrained by the three-dimensional structure but is delineated by its dynamics. This framework must satisfy specificity of function along with adaptability to changing environments and evolvability under external constraints. The accessibility of the available regions of the energy landscape for a set of conditions and shifts in the populations upon their modulation have effects propagating across scales, from biomolecular interactions, to organisms, to populations. Developing the ability to detect and juggle protein conformations supplemented by a physics-based understanding has implications for not only in vivo problems but also for resistance impeding drug discovery and bionano-sensor design.