Cooperativity and flexibility in enzyme evolution

In Silico Structural Insights into a Glucanase from Clostridium perfringens and Prediction of Structural Stability Improvement Through Hydrophobic Interaction Network and Aromatic Interaction

Glucanases are widely applied in industrial applications such as brewing, biomass conversion, food, and animal feed. Glucanases catalyze the hydrolysis of glucan to produce the sugar hemiacetal through hydrolytic cleavage of glycosidic bonds. Current study aimed to investigate structural insights of a glucanase from Clostridium perfringens through blind molecular docking, site-specific molecular docking, molecular dynamics (MD) simulation, and binding energy calculation. Furthermore, we aimed to enhance structural stabilization through formation of hydrophobic interaction network. The molecular docking results illustrated that residues Glu222 and Asp187 may act as nucleophile acid/base catalyst. Moreover, the MM/PBSA results illustrated a high binding affinity of 108.71 ± 8.5 kJ/mol between glucanase and barely glucan during 100 ns simulation. The RMSF analysis illustrated a high flexible surface loop with the highest mobility at position D130. Therefore, the structural engineering was carried out through introducing a double-mutant S125Y/D130P, and the structural stability was improved by forming the hydrophobic interaction network and one π-π aromatic interaction. The spatial distance between the mutation sites and the catalytic pocket attenuates their direct impact on binding interactions within the catalytic pocket.

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.

A Novel A105Y Mutant of CYP17A1 Exhibits Almost Perfect Regioselectivity in the Production of 17α-Hydroxyprogesterone

17α-Hydroxyprogesterone (17α-OHP) is a steroid hormone with significant biological activity that can be obtained by catalyzing progesterone (PROG), the main product of sitosterol, through CYP17A1. However, increasing the catalytic specificity of HCYP17A1 for C17 hydroxylation of progesterone (PROG) poses a formidable challenge due to the close proximity of the C16 and C17 positions. In this study, a rational design was utilized to alter the spatial configuration of the substrate channel, leading to the complete abolition of its 16-hydroxylation activity. Subsequent molecular dynamics simulations revealed that the A105Y mutation heightened the rigidity of the G95-I112 region of CYP17A1, consequently regulating the direction of the entry of PROG into the catalytic pocket. Moreover, the establishment of hydrogen bonding between Y105 and R239, coupled with Pi-stacking of A105Y with F114, effectively immobilizes the substrate PROG in a fixed position, explaining the practically perfect regioselectivity observed in A105Y. Finally, a multifaceted enzymatic cascade system, incorporating A105Y, cytochrome P450 reductase (CPR), and glucose-6-phosphate dehydrogenase (ZWF) for NADPH cofactor regeneration, was constructed in Pichia pastoris GS115. The resulting biocatalyst produced 106 ± 3.2 mg L-1 17α-OHP, a 4.6-fold increase compared with A105Y alone. Thus, this study provides valuable insights for improving the regioselectivity and activity of P450 enzymes.

Changes in Active Site Loop Conformation Relate to the Transition toward a Novel Enzymatic Activity

Enzymatic promiscuity, the ability of enzymes to catalyze multiple, distinct chemical reactions, has been well documented and is hypothesized to be a major driver of the emergence of new enzymatic functions. Yet, the molecular mechanisms involved in the transition from one activity to another remain debated and elusive. Here, we evaluated the redesign of the active site binding cleft of lactonase SsoPox using structure-based design and combinatorial libraries. We created variants with largely improved catalytic abilities against phosphotriesters, the best ones being >1000-fold better compared to the wild-type enzyme. The observed shifts in activity specificity are large, and some variants completely lost their initial activity. The selected combinations of mutations have considerably reshaped the active site cavity via side chain changes but mostly through large rearrangements of the active site loops and changes to their conformations, as revealed by a suite of crystal structures. This suggests that a specific active site loop configuration is critical to the lactonase activity. Interestingly, analysis of high-resolution structures hints at the potential role of conformational sampling and its directionality in defining the enzyme activity profile.

Flexible active-site loops fine-tune substrate specificity of hyperthermophilic metallo-oxidases

Hyperthermophilic ('superheat-loving') archaea found in high-temperature environments such as Pyrobaculum aerophilum contain multicopper oxidases (MCOs) with remarkable efficiency for oxidizing cuprous and ferrous ions. In this work, directed evolution was used to expand the substrate specificity of P. aerophilum McoP for organic substrates. Six rounds of error-prone PCR and DNA shuffling followed by high-throughput screening lead to the identification of a hit variant with a 220-fold increased efficiency (kcat/Km) than the wild-type for 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid) (ABTS) without compromising its intrinsic activity for metal ions. The analysis of the X-ray crystal structure reveals four proximal mutations close to the T1Cu active site. One of these mutations is within the 23-residues loop that occludes this site, a distinctive feature of prokaryotic MCOs. The increased flexibility of this loop results in an enlarged tunnel and one additional pocket that facilitates bulky substrate-enzyme interactions. These findings underscore the synergy between mutations that modulate the dynamics of the active-site loop enabling enhanced catalytic function. This study highlights the potential of targeting loops close to the T1Cu for engineering improvements suitable for biotechnological applications.

Mathematical models describing oxygen binding by hemoglobin

Despite the fact that the investigation of the structural and functional properties of hemoglobin dates back more than 150 years, the topic has not lost its relevance today. The most important component of these studies is the development of mathematical models that formalize and generalize the mechanisms determining the cooperative binding of ligands based on data on the structural and functional state of the protein. In this work, we review the mathematical relationships describing oxygen binding by hemoglobin, ranging from the classical Hüfner, Hill, and Adair equations to the Szabo-Karplus and tertiary two-state mathematical models based on the Monod-Wyman-Changeux and Koshland-Némethy-Filmer concepts. The generality of the considered equations as mathematical functions, bearing in their basis a power dependence, is demonstrated. The problems and possible solutions related to approximation of experimental data by the oxygenation equations with correlated fitting parameters are noted. Attention is paid to empirical equations, extended versions of the Hill equation, where the coefficient of cooperation is modulated by Gauss and Lorentz distributions as functions of partial oxygen pressure.

Development of a novel single-chain l-glutamate oxidase from Streptomyces sp. X-119-6 by inserting flexible linkers

L-glutamate oxidase (LGOX, EC: 1.4.3.11) is an oxidoreductase that catalyzes L-glutamate deamination. LGOX from Streptomyces sp. X-119-6 is used widely for L-glutamate quantification in research and industrial applications. This enzyme encoded as a single precursor chain that undergoes post-translational cleavage to four fragments by an endogenous protease to become highly active. Efficient preparation of active LGOX by heterologous expression without proteolysis process should be indispensable for wide application of this enzyme. Thus, developing an LGOX that requires no protease treatment should expand the potential applications of recombinant LGOX. In this report, we succeeded in obtaining an active single-chain LGOX by connecting the four fragments of the mature form with insertion of flexible linkers. The most active single-chain mutant showed the similar activity to that of the mature form from Streptomyces sp. X-119-6. The structure of this mutant was determined at 2.9 Å resolution by X-ray crystallography. It was revealed that this single-stranded mutant had the similar conformation to that of mature form. This single-chain LGOX can be produced efficiently and should expand LGOX applications.

Substrate-dependent inactivation of recombinant paraoxonase 1 during catalytic dihydrocoumarin turnover and the protective properties of surfactants

Human paraoxonase-1 (PON1) is the most studied member of the paraoxonases (PONs) family and catalyzes the hydrolysis of various substrates (lactones, aryl esters, and paraoxon). Numerous studies link PON1 to oxidative stress-related diseases such as cardiovascular disease, diabetes, HIV infection, autism, Parkinson's, and Alzheimer's, where the kinetic behavior of an enzyme is characterized by initial rates or by modern methods that obtain enzyme kinetic parameters by fitting the computed curves over the entire time-courses of product formation (progress curves). In the analysis of progress curves, the behavior of PON1 during hydrolytically catalyzed turnover cycles is unknown. Hence, progress curves for enzyme-catalyzed hydrolysis of the lactone substrate dihydrocoumarin (DHC) by recombinant PON1 (rePON1) were analyzed to investigate the effect of catalytic DHC turnover on the stability of rePON1. Although rePON1 was significantly inactivated during the catalytic DHC turnover, its activity was not lost due to the product inhibition or spontaneous inactivation of rePON1 in the sample buffers. Examination of the progress curves of DHC hydrolysis by rePON1 led to the conclusion that rePON1 inactivates itself during catalytic DHC turnover hydrolysis. Moreover, human serum albumin or surfactants protected rePON1 from inactivation during this catalytic process, which is significant because the activity of PON1 in clinical samples is measured in the presence of albumin.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

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.

Metallo-Beta-Lactamase-like Encoding Genes in Candidate Phyla Radiation: Widespread and Highly Divergent Proteins with Potential Multifunctionality

The Candidate Phyla Radiation (CPR) was found to harbor a vast repertoire of genes encoding for enzymes with potential antibiotic resistance activity. Among these, as many as 3349 genes were predicted in silico to contain a metallo-beta-lactamase-like (MBL-like) fold. These proteins were subject to an in silico functional characterization by comparing their protein profiles (presence/absence of conserved protein domains) to other MBLs, including 24 already expressed in vitro, along with those of the beta-lactamase database (BLDB) (n = 761). The sequence similarity network (SSN) was then used to predict the functional clusters of CPR MBL-like sequences. Our findings showed that CPR MBL-like sequences were longer and more diverse than bacterial MBL sequences, with a high content of functional domains. Most CPR MBL-like sequences did not show any SSN connectivity with expressed MBLs, indicating the presence of many potential, yet unidentified, functions in CPR. In conclusion, CPR was shown to have many protein functions and a large sequence variability of MBL-like folds, exceeding all known MBLs. Further experimental and evolutionary studies of this superfamily of hydrolyzing enzymes are necessary to illustrate their functional annotation, origin, and expansion for adaptation or specialization within a given niche or compared to a specific substrate.

The Conformational Transitions and Dynamics of Burkholderia cepacia Lipase Regulated by Water-Oil Interfaces

Structural dynamics and conformational transitions are crucial for the activities of enzymes. As one of the most widely used industrial biocatalysts, lipase could be activated by the water-oil interfaces. The interface activations were believed to be dominated by the close-to-open transitions of the lid subdomains. However, the detailed mechanism and the roles of structure transitions are still under debate. In this study, the dynamic structures and conformational transitions of Burkholderia cepacia lipase (LipA) were investigated by combining all-atom molecular dynamics simulations, enhanced sampling simulation, and spectrophotometric assay experiments. The conformational transitions between the lid-open and lid-closed states of LipA in aqueous solution are directly observed by the computational simulation methods. The interactions between the hydrophobic residues on the two lid-subdomains are the driven forces for the LipA closing. Meanwhile, the hydrophobic environment provided by the oil interfaces would separate the interactions between the lid-subdomains and promote the structure opening of LipA. Moreover, our studies demonstrate the opening of the lids structure is insufficient to initiate the interfacial activation, providing explanations for the inability of interfacial activation of many lipases with lid structures.

Changes in Active Site Loop Conformation Relate to the Transition toward a Novel Enzymatic Activity

Enzymatic promiscuity, the ability of enzymes to catalyze multiple, distinct chemical reactions, has been well documented and is hypothesized to be a major driver for the emergence of new enzymatic functions. Yet, the molecular mechanisms involved in the transition from one activity to another remain debated and elusive. Here, we evaluated the redesign of the active site binding cleft of the lactonase SsoPox using structure-based design and combinatorial libraries. We created variants with largely improved catalytic abilities against phosphotriesters, the best ones being > 1,000-fold better compared to the wild-type enzyme. The observed shifts in activity specificity are large, ~1,000,000-fold and beyond, since some variants completely lost their initial activity. The selected combinations of mutations have considerably reshaped the active site cavity via side chain changes but mostly through large rearrangements of the active site loops, as revealed by a suite of crystal structures. This suggests that specific active site loop configuration is critical to the lactonase activity. Interestingly, analysis of high-resolution structures hints at the potential role of conformational sampling and its directionality in defining an enzyme activity profile.

Bypassing evolutionary dead ends and switching the rate-limiting step of a human immunotherapeutic enzyme

The Trp metabolite kynurenine (KYN) accumulates in numerous solid tumours and mediates potent immunosuppression. Bacterial kynureninases (KYNases), which preferentially degrade kynurenine, can relieve immunosuppression in multiple cancer models, but immunogenicity concerns preclude their clinical use, while the human enzyme (HsKYNase) has very low activity for kynurenine and shows no therapeutic effect. Using fitness selections, we evolved a HsKYNase variant with 27-fold higher activity, beyond which exploration of >30 evolutionary trajectories involving the interrogation of >109 variants led to no further improvements. Introduction of two amino acid substitutions conserved in bacterial KYNases reduced enzyme fitness but potentiated rapid evolution of variants with ~500-fold improved activity and reversed substrate specificity, resulting in an enzyme capable of mediating strong anti-tumour effects in mice. Pre-steady-state kinetics revealed a switch in rate-determining step attributable to changes in both enzyme structure and conformational dynamics. Apart from its clinical significance, our work highlights how rationally designed substitutions can potentiate trajectories that overcome barriers in protein evolution.

A New Model of Hemoglobin Oxygenation

The study of hemoglobin oxygenation, starting from the classical works of Hill, has laid the foundation for molecular biophysics. The cooperative nature of oxygen binding to hemoglobin has been variously described in different models. In the Adair model, which better fits the experimental data, the constants of oxygen binding at various stages differ. However, the physical meaning of the parameters in this model remains unclear. In this work, we applied Hill's approach, extending its interpretation; we obtained a good agreement between the theory and the experiment. The equation in which the Hill coefficient is modulated by the Lorentz distribution for oxygen partial pressure approximates the experimental data better than not only the classical Hill equation, but also the Adair equation.

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.

Modulation of global stability, ligand binding and catalytic properties of trypsin by anions

Specific salts effect is well-known on stability and solubility of proteins, however, relatively limited knowledge is known regarding the effect on catalytic properties of enzymes. Here, we examined the effect of four sodium anions on thermal stability and catalytic properties of trypsin and binding of the fluorescent probe, p-aminobenzamidine (PAB), to the enzyme. We show that the specific anions effect on trypsin properties agrees with the localization of the anions in the Hofmeister series. Thermal stability of trypsin, T_m, the affinity of the fluorescent probe to the binding site, K_d, and the rate constant, k_cat, of trypsin-catalyzed hydrolysis of the substrate N-benzoyl-L-arginine ethyl ester (BAEE) increase with increasing kosmotropic character of anions in the order: perchlorate<bromide<chloride<sulfate, while the value of Michaelis constant, K_M, decreases. Correlations between the values of T_m, K_d for PAB, k_cat, and K_M for BAEE in the presence of 1 M studied salts suggest interrelation among these parameters of the enzyme. Global stabilization as well as increased rigidity of trypsin is accompanied by strengthening of interaction with fluorescent probe PAB and in accordance with decreasing values of K_M for the substrate BAEE. Strong correlations between parameters characterizing the trypsin properties with the charge densities of anions clearly indicate direct electrostatic interaction as a basis of the specific anion effect on the conformational and functional properties of the enzyme.

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.

The Use of the Statistical Entropy in Some New Approaches for the Description of Biosystems

Some problems of describing biological systems with the use of entropy as a measure of the complexity of these systems are considered. Entropy is studied both for the organism as a whole and for its parts down to the molecular level. Correlation of actions of various parts of the whole organism, intercellular interactions and control, as well as cooperativity on the microlevel lead to a more complex structure and lower statistical entropy. For a multicellular organism, entropy is much lower than entropy for the same mass of a colony of unicellular organisms. Cooperativity always reduces the entropy of the system; a simple example of ligand binding to a macromolecule carrying two reaction centers shows how entropy is consistent with the ambiguity of the result in the Bernoulli test scheme. Particular attention is paid to the qualitative and quantitative relationship between the entropy of the system and the cooperativity of ligand binding to macromolecules. A kinetic model of metabolism. corresponding to Schrödinger's concept of the maintenance biosystems by "negentropy feeding", is proposed. This model allows calculating the nonequilibrium local entropy and comparing it with the local equilibrium entropy inherent in non-living matter.

Chemical Mapping Exposes the Importance of Active Site Interactions in Governing the Temperature Dependence of Enzyme Turnover

Uncovering the role of global protein dynamics in enzyme turnover is needed to fully understand enzyme catalysis. Recently, we have demonstrated that the heat capacity of catalysis, ΔC _P ^‡, can reveal links between the protein free energy landscape, global protein dynamics, and enzyme turnover, suggesting that subtle changes in molecular interactions at the active site can affect long-range protein dynamics and link to enzyme temperature activity. Here, we use a model promiscuous enzyme (glucose dehydrogenase from Sulfolobus solfataricus) to chemically map how individual substrate interactions affect the temperature dependence of enzyme activity and the network of motions throughout the protein. Utilizing a combination of kinetics, red edge excitation shift (REES) spectroscopy, and computational simulation, we explore the complex relationship between enzyme-substrate interactions and the global dynamics of the protein. We find that changes in ΔC _P ^‡ and protein dynamics can be mapped to specific substrate-enzyme interactions. Our study reveals how subtle changes in substrate binding affect global changes in motion and flexibility extending throughout the protein.

Unlocking the Stereoselectivity and Substrate Acceptance of Enzymes: Proline-Induced Loop Engineering Test

Protein stability and evolvability influence each other. Although protein dynamics play essential roles in various catalytically important properties, their high flexibility and diversity makes it difficult to incorporate such properties into rational engineering. Therefore, how to unlock the potential evolvability in a user-friendly rational design process remains a challenge. In this endeavor, we describe a method for engineering an enantioselective alcohol dehydrogenase. It enables synthetically important substrate acceptance for 4-chlorophenyl pyridine-2-yl ketone, and perfect stereocontrol of both (S)- and (R)-configured products. Thermodynamic analysis unveiled the subtle interaction between enzyme stability and evolvability, while computational studies provided insights into the origin of selectivity and substrate recognition. Preparative-scale synthesis of the (S)-product (73 % yield; >99 % ee) was performed on a gram-scale. This proof-of-principle study demonstrates that interfaced proline residues can be rationally engineered to unlock evolvability and thus provide access to new biocatalysts with highly improved catalytic performance.

Loops around the Heme Pocket Have a Critical Role in the Function and Stability of BsDyP from Bacillus subtilis

Bacillus subtilis BsDyP belongs to class I of the dye-decolorizing peroxidase (DyP) family of enzymes and is an interesting biocatalyst due to its high redox potential, broad substrate spectrum and thermostability. This work reports the optimization of BsDyP using directed evolution for improved oxidation of 2,6-dimethoxyphenol, a model lignin-derived phenolic. After three rounds of evolution, one variant was identified displaying 7-fold higher catalytic rates and higher production yields as compared to the wild-type enzyme. The analysis of X-ray structures of the wild type and the evolved variant showed that the heme pocket is delimited by three long conserved loop regions and a small α helix where, incidentally, the mutations were inserted in the course of evolution. One loop in the proximal side of the heme pocket becomes more flexible in the evolved variant and the size of the active site cavity is increased, as well as the width of its mouth, resulting in an enhanced exposure of the heme to solvent. These conformational changes have a positive functional role in facilitating electron transfer from the substrate to the enzyme. However, they concomitantly resulted in decreasing the enzyme’s overall stability by 2 kcal mol⁻¹, indicating a trade-off between functionality and stability. Furthermore, the evolved variant exhibited slightly reduced thermal stability compared to the wild type. The obtained data indicate that understanding the role of loops close to the heme pocket in the catalysis and stability of DyPs is critical for the development of new and more powerful biocatalysts: loops can be modulated for tuning important DyP properties such as activity, specificity and stability.

Two β-Lactamase Variants with Reduced Clavulanic Acid Inhibition Display Different Millisecond Dynamics

The β-lactamase of Mycobacterium tuberculosis, BlaC, is susceptible to inhibition by clavulanic acid. The ability of this enzyme to escape inhibition through mutation was probed using error-prone PCR combined with functional screening in Escherichia coli. The variant that was found to confer the most inhibitor resistance, K234R, as well as variant G132N that was found previously were characterized using X-ray crystallography and nuclear magnetic resonance (NMR) relaxation experiments to probe structural and dynamic properties. The G132N mutant exists in solution in two almost equally populated conformations that exchange with a rate of ca. 88 s-1. The conformational change affects a broad region of the enzyme. The crystal structure reveals that the Asn132 side chain forces the peptide bond between Ser104 and Ile105 in a cis-conformation. The crystal structure suggests multiple conformations for several side chains (e.g., Ser104 and Ser130) and a short loop (positions 214 to 216). In the K234R mutant, the active-site dynamics are significantly diminished with respect to the wild-type enzyme. These results show that multiple evolutionary routes are available to increase inhibitor resistance in BlaC and that active-site dynamics on the millisecond time scale are not required for catalytic function.

Metallo-β-lactamases in the Age of Multidrug Resistance: From Structure and Mechanism to Evolution, Dissemination, and Inhibitor Design

Antimicrobial resistance is one of the major problems in current practical medicine. The spread of genes coding for resistance determinants among bacteria challenges the use of approved antibiotics, narrowing the options for treatment. Resistance to carbapenems, last resort antibiotics, is a major concern. Metallo-β-lactamases (MBLs) hydrolyze carbapenems, penicillins, and cephalosporins, becoming central to this problem. These enzymes diverge with respect to serine-β-lactamases by exhibiting a different fold, active site, and catalytic features. Elucidating their catalytic mechanism has been a big challenge in the field that has limited the development of useful inhibitors. This review covers exhaustively the details of the active-site chemistries, the diversity of MBL alleles, the catalytic mechanism against different substrates, and how this information has helped developing inhibitors. We also discuss here different aspects critical to understand the success of MBLs in conferring resistance: the molecular determinants of their dissemination, their cell physiology, from the biogenesis to the processing involved in the transit to the periplasm, and the uptake of the Zn(II) ions upon metal starvation conditions, such as those encountered during an infection. In this regard, the chemical, biochemical and microbiological aspects provide an integrative view of the current knowledge of MBLs.

A reverse catalytic triad Asp containing loop shaping a wide substrate binding pocket of a feruloyl esterase from Lactobacillus plantarum

Feruloyl esterase is an indispensable biocatalyst in food processing, pesticide and pharmaceutical industries, catalyzing the cleavage of the ester bond cross-linked between the polysaccharide side chain of hemicellulose and ferulic acid in plant cell walls. LP_0796 from Lactobacillus plantarum was identified as a feruloyl esterase that may have potential applications in the food industry, but the lack of the substrate recognition and catalytic mechanisms limits its application. Here, LP_0796 showed the highest activity towards methyl caffeate at pH 6.6 and 40 °C. The crystal structure of LP_0796 was determined at 2.5 Å resolution and featured a catalytic triad Asp195-containing loop facing the opposite direction, thus forming a wider substrate binding pocket. Molecular docking simulation and site-directed mutagenesis studies further demonstrated that in addition to the catalytic triad (Ser94, Asp195, His225), Arg125 and Val128 played essential roles in the function of the active site. Our data also showed that Asp mutation of Ala23 and Ile198 increased the catalytic efficiency to 4- and 5-fold, respectively. Collectively, this work provided a better understanding of the substrate recognition and catalytic mechanisms of LP_0796 and may facilitate the future protein design of this important feruloyl esterase.

Comparative QM/MM Molecular Dynamics and Umbrella Sampling Simulations: Interaction of the Zinc-Bound Intermediate Gem-Diolate Trapoxin A Inhibitor and Acetyl-l-lysine Substrate with Histone Deacetylase 8

Targeting the genetic material without destruction is a priority to develop safe anticancer drugs. Histone deacetylase 8 (HDAC8), which is proved to be involved in carcinogenesis, is an enzyme associated with the chromatin for post-translational deacetylation of acetylated lysine. In this study, HDAC8 co-crystallized with the intermediate state tetrapeptide Trapoxin A (TA) inhibitor and the holoenzyme are utilized to find their conformational ensembles. Furthermore, the co-crystallized intermediate gem-diolate TA was used to find optimum interactions with the active site residues by conventional molecular dynamics (MD) simulation and QM/MM umbrella sampling. Finally, the intermediate state of the acetyl-l-lysine substrate was explored by QM/MM steered MD and compared to the binding of the intermediate state of the inhibitor. This research showed that HDAC8 is flexible and exists in conformational ensembles in its holoenzyme state. Binding of the intermediate state TA stabilizes its conformation. The optimum binding to the active site of HDAC8 for structures of gem-diolate TA (intermediate state) and acetyl-l-lysine (intermediate state) was determined according to the corresponding energy profiles. The use of these models will aid in the design of potentially reversible, potent, and selective inhibitors of HDAC8 for cancer treatment.

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

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

Pervasive cooperative mutational effects on multiple catalytic enzyme traits emerge via long-range conformational dynamics

Multidimensional fitness landscapes provide insights into the molecular basis of laboratory and natural evolution. To date, such efforts usually focus on limited protein families and a single enzyme trait, with little concern about the relationship between protein epistasis and conformational dynamics. Here, we report a multiparametric fitness landscape for a cytochrome P450 monooxygenase that was engineered for the regio- and stereoselective hydroxylation of a steroid. We develop a computational program to automatically quantify non-additive effects among all possible mutational pathways, finding pervasive cooperative signs and magnitude epistasis on multiple catalytic traits. By using quantum mechanics and molecular dynamics simulations, we show that these effects are modulated by long-range interactions in loops, helices and β-strands that gate the substrate access channel allowing for optimal catalysis. Our work highlights the importance of conformational dynamics on epistasis in an enzyme involved in secondary metabolism and offers insights for engineering P450s.

Connecting conformational dynamics and epistasis has so far been limited to a few proteins and a single fitness trait. Here, the authors provide evidence of positive epistasis on multiple catalytic traits in the evolution and dynamics of engineered cytochrome P450 monooxygenase, offering insights for in silico protein design.

A Single-Point Mutation in d-Arginine Dehydrogenase Unlocks a Transient Conformational State Resulting in Altered Cofactor Reactivity

Proteins are inherently dynamic, and proper enzyme function relies on conformational flexibility. In this study, we demonstrated how an active site residue changes an enzyme's reactivity by modulating fluctuations between conformational states. Replacement of tyrosine 249 (Y249) with phenylalanine in the active site of the flavin-dependent d-arginine dehydrogenase yielded an enzyme with both an active yellow FAD (Y249F-y) and an inactive chemically modified green FAD, identified as 6-OH-FAD (Y249F-g) through various spectroscopic techniques. Structural investigation of Y249F-g and Y249F-y variants by comparison to the wild-type enzyme showed no differences in the overall protein structure and fold. A closer observation of the active site of the Y249F-y enzyme revealed an alternative conformation for some active site residues and the flavin cofactor. Molecular dynamics simulations probed the alternate conformations observed in the Y249F-y enzyme structure and showed that the enzyme variant with FAD samples a metastable conformational state, not available to the wild-type enzyme. Hybrid quantum/molecular mechanical calculations identified differences in flavin electronics between the wild type and the alternate conformation of the Y249F-y enzyme. The computational studies further indicated that the alternate conformation in the Y249F-y enzyme is responsible for the higher spin density at the C6 atom of flavin, which is consistent with the formation of 6-OH-FAD in the variant enzyme. The observations in this study are consistent with an alternate conformational space that results in fine-tuning the microenvironment around a versatile cofactor playing a critical role in enzyme function.

Linking molecular evolution to molecular grafting

Molecular grafting is a strategy for the engineering of molecular scaffolds into new functional agents, such as next-generation therapeutics. Despite its wide use, studies so far have focused almost exclusively on demonstrating its utility rather than understanding the factors that lead to either poor or successful grafting outcomes. Here, we examine protein evolution and identify parallels between the natural process of protein functional diversification and the artificial process of molecular grafting. We discuss features of natural proteins that are correlated to innovability-the capacity to acquire new functions-and describe their implications to molecular grafting scaffolds. Disulfide-rich peptides are used as exemplars because they are particularly promising scaffolds onto which new functions can be grafted. This article provides a perspective on why some scaffolds are more suitable for grafting than others, identifying opportunities on how molecular grafting might be improved.

Temperature-Dependent Kinetic Isotope Effects in R67 Dihydrofolate Reductase from Path-Integral Simulations

Calculation of temperature-dependent kinetic isotope effects (KIE) in enzymes presents a significant theoretical challenge. Additionally, it is not trivial to identify enzymes with available experimental accurate intrinsic KIEs in a range of temperatures. In the current work, we present a theoretical study of KIEs in the primitive R67 dihydrofolate reductase (DHFR) enzyme and compare with experimental work. The advantage of R67 DHFR is its significantly lower kinetic complexity compared to more evolved DHFR isoforms. We employ mass-perturbation-based path-integral simulations in conjunction with umbrella sampling and a hybrid quantum mechanics-molecular mechanics Hamiltonian. We obtain temperature-dependent KIEs in good agreement with experiments and ascribe the temperature-dependent KIEs primarily to zero-point energy effects. The active site in the primitive enzyme is found to be poorly preorganized, which allows excessive water access to the active site and results in loosely bound reacting ligands.

Heme-binding enables allosteric modulation in an ancient TIM-barrel glycosidase

Glycosidases are phylogenetically widely distributed enzymes that are crucial for the cleavage of glycosidic bonds. Here, we present the exceptional properties of a putative ancestor of bacterial and eukaryotic family-1 glycosidases. The ancestral protein shares the TIM-barrel fold with its modern descendants but displays large regions with greatly enhanced conformational flexibility. Yet, the barrel core remains comparatively rigid and the ancestral glycosidase activity is stable, with an optimum temperature within the experimental range for thermophilic family-1 glycosidases. None of the ∼5500 reported crystallographic structures of ∼1400 modern glycosidases show a bound porphyrin. Remarkably, the ancestral glycosidase binds heme tightly and stoichiometrically at a well-defined buried site. Heme binding rigidifies this TIM-barrel and allosterically enhances catalysis. Our work demonstrates the capability of ancestral protein reconstructions to reveal valuable but unexpected biomolecular features when sampling distant sequence space. The potential of the ancestral glycosidase as a scaffold for custom catalysis and biosensor engineering is discussed.

Family 1 glycosidases (GH1) are present in the three domains of life and share classical TIM-barrel fold. Structural and biochemical analyses of a resurrected ancestral GH1 enzyme reveal heme binding, not known in its modern descendants. Heme rigidifies the TIM-barrel and allosterically enhances catalysis.

Conformational Dynamics Contribute to Substrate Selectivity and Catalysis in Human Kynureninase

Kynureninases (KYNases) are enzymes that play a key role in tryptophan catabolism through the degradation of intermediate kynurenine and 3'-hydroxy-kynurenine metabolites (KYN and OH-KYN, respectively). Bacterial KYNases exhibit high catalytic efficiency toward KYN and moderate activity toward OH-KYN, whereas animal KYNases are highly selective for OH-KYN, exhibiting only minimal activity toward the smaller KYN substrate. These differences reflect divergent pathways for KYN and OH-KYN utilization in the respective kingdoms. We examined the Homo sapiens and Pseudomonas fluorescens KYNases (HsKYNase and PfKYNase respectively) using pre-steady-state and hydrogen-deuterium exchange mass spectrometry (HDX-MS) methodologies. We discovered that the activity of HsKYNase critically depends on formation of hydrogen bonds with the hydroxyl group of OH-KYN to stabilize the entire active site and allow productive substrate turnover. With the preferred OH-KYN substrate, stabilization is observed at the substrate-binding site and the region surrounding the PLP cofactor. With the nonpreferred KYN substrate, less stabilization occurs, revealing a direct correlation with activity. This correlation holds true for PfKYNases; however there is only a modest stabilization at the substrate-binding site, suggesting that substrate discrimination is simply achieved by steric hindrance. We speculate that eukaryotic KYNases use dynamic mobility as a mechanism of substrate specificity to commit OH-KYN to nicotinamide synthesis and avoid futile hydrolysis of KYN. These findings have important ramifications for the engineering of HsKynase with high KYN activity as required for clinical applications in cancer immunotherapy. Our study shows how homologous enzymes with conserved active sites can use dynamics to discriminate between two highly similar substrates.

Enzyme Promiscuous Activity: How to Define it and its Evolutionary Aspects

Enzymes are among the most studied biological molecules because better understanding enzymes structure and activity will shed more light on their biological processes and regulation; from a biotechnological point of view there are many examples of enzymes used with the aim to obtain new products and/or to make industrial processes less invasive towards the environment. Enzymes are known for their high specificity in the recognition of a substrate but considering the particular features of an increasing number of enzymes this is not completely true, in fact, many enzymes are active on different substrates: this ability is called enzyme promiscuity. Usually, promiscuous activities have significantly lower kinetic parameters than to that of primary activity, but they have a crucial role in gene evolution. It is accepted that gene duplication followed by sequence divergence is considered a key evolutionary mechanism to generate new enzyme functions. In this way, promiscuous activities are the starting point to increase a secondary activity in the main activity and then get a new enzyme. The primary activity can be lost or reduced to a promiscuous activity. In this review we describe the differences between substrate and enzyme promiscuity, and its rule in gene evolution. From a practical point of view the knowledge of promiscuity can facilitate the in vitro progress of proteins engineering, both for biomedical and industrial applications. In particular, we report cases regarding esterases, phosphotriesterases and cytochrome P450.

Asymmetric dynamic coupling promotes alternative evolutionary pathways in an enzyme dimer

The importance of dynamic factors in enzyme evolution is gaining recognition. Here we study how the evolution of a new enzymatic activity exploits conformational tinkering and demonstrate that conversion of a dimeric phosphotriesterase to an arylesterase in Pseudomonas diminuta is accompanied by structural divergence between the two subunits. Deviations in loop conformations increase with promiscuity, leading to functionally distinct states, while they decrease during specialisation for the new function. We show that opposite loop movements in the two subunits are due to a dynamic coupling with the dimer interface, the importance of which is also corroborated by the co-evolution of the loop and interface residues. These results illuminate how protein dynamics promotes conformational heterogeneity in a dimeric enzyme, leading to alternative evolutionary pathways for the emergence of a new function.

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.

Thermal adaptation in the honeybee (Apis mellifera) via changes to the structure of malate dehydrogenase

In honeybees there are three alleles of cytosolic malate dehydrogenase gene: F, M and S. Allele frequencies are correlated with environmental temperature, suggesting that the alleles have temperature-dependent fitness benefits. We determined the enzyme activity of each allele across a range of temperatures in vitro The F and S alleles have higher activity and are less sensitive to high temperatures than the M allele, which loses activity after incubation at temperatures found in the thorax of foraging bees in hot climates. Next, we predicted the protein structure of each allele and used molecular dynamics simulations to investigate their molecular flexibility. The M allozyme is more flexible than the S and F allozymes at 50°C, suggesting a plausible explanation for its loss of activity at high temperatures, and has the greatest structural flexibility at 15°C, suggesting that it can retain some enzyme activity at cooler temperatures. MM bees recovered from 2 h of cold narcosis significantly better than all other genotypes. Combined, these results explain clinal variation in malate dehydrogenase allele frequencies in the honeybee at the molecular level.

Elucidating the potential neurotoxicity of chiral phenthoate: Molecular insight from experimental and computational studies

Chiral organophosphorus pollutants are existed ubiquitously in the ecological environment, but the enantioselective toxicities of these nerve agents to humans and their molecular bases have not been fully elucidated. Using experimental and computational approaches, this story was to explore the neurotoxic response process of the target acetylcholinesterase (AChE) to chiral phenthoate and further decipher the microscopic mechanism of such toxicological effect at the enantiomeric level. The results showed that the toxic reaction of AChE with chiral phenthoate exhibited significant enantioselectivity, and (R)-phenthoate (K＝1.486 × 10⁵ M^-1) has a bioaffinity for the nerve enzyme nearly three times that of (S)-phenthoate (K＝4.503 × 10⁴ M^-1). Dynamic research outcomes interpreted the wet experiments, and the inherent conformational flexibility of the target enzyme has a great influence on the enantioselective neurotoxicological action processes, especially reflected in the conformational changes of the three key loop regions (i.e. residues His-447, Gly-448, and Tyr-449; residues Gly-122, Phe-123, and Tyr-124; and residues Thr-75, Leu-76, and Tyr-77) around the reaction patch. This was supported by the quantitative results of conformational studies derived from circular dichroism spectroscopy (α-helix: 34.7%→30.2%/31.6%; β-sheet: 23.6%→19.5%/20.7%; turn: 19.2%→22.4%/21.9%; and random coil: 22.5%→27.9%/25.8%). Meanwhile, via analyzing the modes of toxic action and free energies, we can find that (R)-phenthoate has a strong inhibitory effect on the enzymatic activity of AChE, as compared with (S)-phenthoate, and electrostatic energy (-23.79/－17.77 kJ mol^-1) played a critical role in toxicological reactions. These points were the underlying causes of chiral phenthoate displaying different degrees of enantioselective neurotoxicity.

Resurrection of ancestral effector caspases identifies novel networks for evolution of substrate specificity

Apoptotic caspases evolved with metazoans more than 950 million years ago (MYA), and a series of gene duplications resulted in two subfamilies consisting of initiator and effector caspases. The effector caspase genes (caspases-3, -6, and -7) were subsequently fixed into the Chordata phylum more than 650 MYA when the gene for a common ancestor (CA) duplicated, and the three effector caspases have persisted throughout mammalian evolution. All caspases prefer an aspartate residue at the P1 position of substrates, so each caspase evolved discrete cellular roles through changes in substrate recognition at the P4 position combined with allosteric regulation. We examined the evolution of substrate specificity in caspase-6, which prefers valine at the P4 residue, compared with caspases-3 and -7, which prefer aspartate, by reconstructing the CA of effector caspases (AncCP-Ef1) and the CA of caspase-6 (AncCP-6An). We show that AncCP-Ef1 is a promiscuous enzyme with little distinction between Asp, Val, or Leu at P4. The specificity of caspase-6 was defined early in its evolution, where AncCP-6An demonstrates a preference for Val over Asp at P4. Structures of AncCP-Ef1 and of AncCP-6An show a network of charged amino acids near the S4 pocket that, when combined with repositioning a flexible active site loop, resulted in a more hydrophobic binding pocket in AncCP-6An. The ancestral protein reconstructions show that the caspase-hemoglobinase fold has been conserved for over 650 million years and that only three substitutions in the scaffold are necessary to shift substrate selection toward Val over Asp.

Harnessing Conformational Plasticity to Generate Designer Enzymes

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

Microsecond dynamics in proteins by two-dimensional ESR: Predictions

Two-dimensional electron-electron double resonance (2D-ELDOR) provides extensive insight into molecular motions. Recent developments permitting experiments at higher frequencies (95 GHz) provide molecular orientational resolution, enabling a clearer description of the nature of the motions. In this work, simulations are provided for the example of domain motions within proteins that are themselves slowly tumbling in solution. These show the nature of the exchange cross-peaks that are predicted to develop in real time from such domain motions. However, we find that the existing theoretical methods for computing 2D-ELDOR experiments over a wide motional range begin to fail seriously when applied to very slow motions characteristic of proteins in solution. One reason is the failure to obtain accurate eigenvectors and eigenvalues of the complex symmetric stochastic Liouville matrices describing the experiment when computed by the efficient Lanczos algorithm in the range of very slow motion. Another, perhaps more serious, issue is that these matrices are "non-normal," such that for the very slow motional range even rigorous diagonalization algorithms do not yield the correct eigenvalues and eigenvectors. We have employed algorithms that overcome both these issues and lead to valid 2D-ELDOR predictions even for motions approaching the rigid limit. They are utilized to describe the development of cross-peaks in 2D-ELDOR at 95 GHz for a particular case of domain motion.

Temperature, Dynamics, and Enzyme-Catalyzed Reaction Rates

We review the adaptations of enzyme activity to different temperatures. Psychrophilic (cold-adapted) enzymes show significantly different activation parameters (lower activation enthalpies and entropies) from their mesophilic counterparts. Furthermore, there is increasing evidence that the temperature dependence of many enzyme-catalyzed reactions is more complex than is widely believed. Many enzymes show curvature in plots of activity versus temperature that is not accounted for by denaturation or unfolding. This is explained by macromolecular rate theory: A negative activation heat capacity for the rate-limiting chemical step leads directly to predictions of temperature optima; both entropy and enthalpy are temperature dependent. Fluctuations in the transition state ensemble are reduced compared to the ground state. We show how investigations combining experiment with molecular simulation are revealing fundamental details of enzyme thermoadaptation that are relevant for understanding aspects of enzyme evolution. Simulations can calculate relevant thermodynamic properties (such as activation enthalpies, entropies, and heat capacities) and reveal the molecular mechanisms underlying experimentally observed behavior.

Dynamics, a Powerful Component of Current and Future in Silico Approaches for Protein Design and Engineering

Computational prediction has become an indispensable aid in the processes of engineering and designing proteins for various biotechnological applications. With the tremendous progress in more powerful computer hardware and more efficient algorithms, some of in silico tools and methods have started to apply the more realistic description of proteins as their conformational ensembles, making protein dynamics an integral part of their prediction workflows. To help protein engineers to harness benefits of considering dynamics in their designs, we surveyed new tools developed for analyses of conformational ensembles in order to select engineering hotspots and design mutations. Next, we discussed the collective evolution towards more flexible protein design methods, including ensemble-based approaches, knowledge-assisted methods, and provable algorithms. Finally, we highlighted apparent challenges that current approaches are facing and provided our perspectives on their further development.

Ligand-Binding-Site Structure Shapes Allosteric Signal Transduction and the Evolution of Allostery in Protein Complexes

The structure of ligand-binding sites has been shown to profoundly influence the evolution of function in homomeric protein complexes. Complexes with multichain binding sites (MBSs) have more conserved quaternary structure, more similar binding sites and ligands between homologs, and evolve new functions slower than homomers with single-chain binding sites (SBSs). Here, using in silico analyses of protein dynamics, we investigate whether ligand-binding-site structure shapes allosteric signal transduction pathways, and whether the structural similarity of binding sites influences the evolution of allostery. Our analyses show that: 1) allostery is more frequent among MBS complexes than in SBS complexes, particularly in homomers; 2) in MBS homomers, semirigid communities and critical residues frequently connect interfaces and thus they are characterized by signal transduction pathways that cross protein-protein interfaces, whereas SBS homomers usually not; 3) ligand binding alters community structure differently in MBS and SBS homomers; and 4) except MBS homomers, allosteric proteins are more likely to have homologs with similar binding site than nonallosteric proteins, suggesting that binding site similarity is an important factor driving the evolution of allostery.