skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Ensemble-based enzyme design can recapitulate the effects of laboratory directed evolution in silico
Abstract The creation of artificial enzymes is a key objective of computational protein design. Although de novo enzymes have been successfully designed, these exhibit low catalytic efficiencies, requiring directed evolution to improve activity. Here, we use room-temperature X-ray crystallography to study changes in the conformational ensemble during evolution of the designed Kemp eliminase HG3 (kcat/KM146 M−1s−1). We observe that catalytic residues are increasingly rigidified, the active site becomes better pre-organized, and its entrance is widened. Based on these observations, we engineer HG4, an efficient biocatalyst (kcat/KM103,000 M−1s−1) containing key first and second-shell mutations found during evolution. HG4 structures reveal that its active site is pre-organized and rigidified for efficient catalysis. Our results show how directed evolution circumvents challenges inherent to enzyme design by shifting conformational ensembles to favor catalytically-productive sub-states, and suggest improvements to the design methodology that incorporate ensemble modeling of crystallographic data.  more » « less
Award ID(s):
1231306
PAR ID:
10587742
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Nature Communications
Date Published:
Journal Name:
Nature Communications
Volume:
11
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract De novo enzyme design has sought to introduce active sites and substrate-binding pockets that are predicted to catalyse a reaction of interest into geometrically compatible native scaffolds1,2, but has been limited by a lack of suitable protein structures and the complexity of native protein sequence–structure relationships. Here we describe a deep-learning-based ‘family-wide hallucination’ approach that generates large numbers of idealized protein structures containing diverse pocket shapes and designed sequences that encode them. We use these scaffolds to design artificial luciferases that selectively catalyse the oxidative chemiluminescence of the synthetic luciferin substrates diphenylterazine3and 2-deoxycoelenterazine. The designed active sites position an arginine guanidinium group adjacent to an anion that develops during the reaction in a binding pocket with high shape complementarity. For both luciferin substrates, we obtain designed luciferases with high selectivity; the most active of these is a small (13.9 kDa) and thermostable (with a melting temperature higher than 95 °C) enzyme that has a catalytic efficiency on diphenylterazine (kcat/Km = 106 M−1 s−1) comparable to that of native luciferases, but a much higher substrate specificity. The creation of highly active and specific biocatalysts from scratch with broad applications in biomedicine is a key milestone for computational enzyme design, and our approach should enable generation of a wide range of luciferases and other enzymes. 
    more » « less
  2. Abstract The roles of local interactions in the laboratory evolution of a highly active, computationally designed retroaldolase (RA) are examined. Partial Order Optimum Likelihood (POOL) is used to identify catalytically important amino acid interactions in several RA95 enzyme variants. The series RA95.5, RA95.5–5, RA95.5–8, and RA95.5–8F, representing progress along an evolutionary trajectory with increasing activity, is examined. Computed measures of coupling between charged states of residues show that, as evolution proceeds and higher activities are achieved, electrostatic coupling between the biochemically active amino acids and other residues is increased. In silico residue scanning suggests multiple coupling partners for the catalytic lysine K83. The effects of two predicted partners, Y51 and E85, are tested using site‐directed mutagenesis and kinetic analysis of the variants Y51F and E85Q. The Y51F variants show decreases inkcatrelative to wild type, with the greatest losses observed for the more evolved constructs; they also exhibit significant decreases inkcat/KMacross the series. Only modest decreases inkcat/KMare observed for the E85Q variants with little effect onkcat. Computed metrics of the degree of coupling between protonation states rise significantly as evolution proceeds and catalytic turnover rate increases. Specifically, the charge state of the catalytic lysine K83 becomes more strongly coupled to those of other amino acids as the enzyme evolves to a better catalyst. 
    more » « less
  3. AbstractRecently, we reported the discovery of a novel endoglucanase of the glycoside hydrolase family 12 (GH12), designated IfCelS12A, from the haloalkaliphilic anaerobic bacteriumIocasia fonsfrigidaestrain SP3-1, which was isolated from a hypersaline pond in the Samut Sakhon province of Thailand (ca. 2017). IfCelS12A exhibits high substrate specificity on carboxymethyl cellulose and amorphous cellulose but low substrate specificity on b-1,3;1,4-glucan. Unlike some endoglucanases of the GH12 family, IfCelS12A does not exhibit hydrolytic activity on crystalline cellulose (i.e., Avicel™). High-Pressure Liquid Chromatography (HPLC) and Thin Layer Chromatography (TLC) analyses of products resulting from IfCelS12-mediated hydrolysis indicate mode of action for this enzyme. Notably, IfCelS12A preferentially hydrolyzes cellotetraoses, cellopentaoses, and cellohexaoses with negligible activity on cellobiose or cellotriose. Kinetic analysis with cellopentaose and barely b-d-glucan as cellulosic substrates were conducted. On cellopentaose, IfCelS12A demonstrates a 16-fold increase in activity (KM = 0.27 mM;kcat = 0.36 s−1;kcat/KM = 1.34 mM−1s−1) compared to the enzymatic hydrolysis of barley b-d-glucan (KM: 0.04 mM,kcat: 0.51 s−1,kcat/KM = 0.08 mM−1s−1). Moreover, IfCelS12A enzymatic efficacy is stable in hypersaline sodium chlorids (NaCl) solutions (up to 10% NaCl). Specifically, IfCel12A retains notable activity after 24 h at 2M NaCl (10% saline solution). IfCelS12A used as a cocktail component with other cellulolytic enzymes and in conjunction with mobile sequestration platform technology offers additional options for deconstruction of ionic liquid–pretreated cellulosic feedstock. Key points•IfCelS12A from an anaerobic alkaliphile Iocasia fronsfrigidae shows salt tolerance•IfCelS12A in cocktails with other enzymes efficiently degrades cellulosic biomass•IfCelS12A used with mobile enzyme sequestration platforms enhances hydrolysis 
    more » « less
  4. Rosenzweig, Amy (Ed.)
    Metformin is the first-line treatment for type II diabetes patients and a pervasive pollutant with more than 180 million kg ingested globally and entering wastewater. The drug’s direct mode of action is currently unknown but is linked to effects on gut microbiomes and may involve specific gut microbial reactions to the drug. In wastewater treatment plants, metformin is known to be transformed by microbes to guanylurea, although genes encoding this metabolism had not been elucidated. In the present study, we revealed the function of two genes responsible for metformin decomposition (mfmAandmfmB) found in isolated bacteria from activated sludge. MfmA and MfmB form an active heterocomplex (MfmAB) and are members of the ureohydrolase protein superfamily with binuclear metal-dependent activity. MfmAB is nickel-dependent and catalyzes the hydrolysis of metformin to dimethylamine and guanylurea with a catalytic efficiency (kcat/KM) of 9.6 × 103M−1s−1and KMfor metformin of 0.82 mM. MfmAB shows preferential activity for metformin, being able to discriminate other close substrates by several orders of magnitude. Crystal structures of MfmAB show coordination of binuclear nickel bound in the active site of the MfmA subunit but not MfmB subunits, indicating that MfmA is the active site for the MfmAB complex. Mutagenesis of residues conserved in the MfmA active site revealed those critical to metformin hydrolase activity and its small substrate binding pocket allowed for modeling of bound metformin. This study characterizes the products of themfmABgenes identified in wastewater treatment plants on three continents, suggesting that metformin hydrolase is widespread globally in wastewater. 
    more » « less
  5. The aerobic hyperthermophile“Fervidibacter sacchari”catabolizes diverse polysaccharides and is the only cultivated member of the class“Fervidibacteria”within the phylumArmatimonadota. It encodes 117 putative glycoside hydrolases (GHs), including two from GH family 50 (GH50). In this study, we expressed, purified, and functionally characterized one of these GH50 enzymes, Fsa16295Glu. We show that Fsa16295Glu is a β-1,3-endoglucanase with optimal activity on carboxymethyl curdlan (CM-curdlan) and only weak agarase activity, despite most GH50 enzymes being described as β-agarases. The purified enzyme has a wide temperature range of 4–95°C (optimal 80°C), making it the first characterized hyperthermophilic representative of GH50. The enzyme is also active at a broad pH range of at least 5.5–11 (optimal 6.5–10). Fsa16295Glu possesses a relatively highkcat/KMof 1.82 × 107 s−1M−1with CM-curdlan and degrades CM-curdlan nearly completely to sugar monomers, indicating preferential hydrolysis of glucans containing β-1,3 linkages. Finally, a phylogenetic analysis of Fsa16295Glu and all other GH50 enzymes revealed that Fsa16295Glu is distant from other characterized enzymes but phylogenetically related to enzymes from thermophilic archaea that were likely acquired horizontally from“Fervidibacteria.”Given its functional and phylogenetic novelty, we propose that Fsa16295Glu represents a new enzyme subfamily, GH50_3. 
    more » « less