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: Analysis of electrostatic coupling throughout the laboratory evolution of a designed retroaldolase
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
Award ID(s):
1905214
PAR ID:
10449690
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Protein Science
Volume:
30
Issue:
8
ISSN:
0961-8368
Page Range / eLocation ID:
p. 1617-1627
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Interactions in enzymes between catalytic and neighboring amino acids and how these interactions facilitate catalysis are examined. In examples from both natural and designed enzymes, it is shown that increases in catalytic rates may be achieved through elongation of the buffer range of the catalytic residues; such perturbations in the protonation equilibria are, in turn, achieved through enhanced coupling of the protonation equilibria of the active ionizable residues with those of other ionizable residues. The strongest coupling between protonation states for a pair of residues that deprotonate to form an anion (or a pair that accept a proton to form a cation) is achieved when the difference in the intrinsic pKas of the two residues is approximately within 1 pH unit. Thus, catalytic aspartates and glutamates are often coupled to nearby acidic residues. For an anion-forming residue coupled to a cation-forming residue, the elongated buffer range is achieved when the intrinsic pKa of the anion-forming residue is higher than the intrinsic pKa of the (conjugate acid of the) cation-forming residue. Therefore, the high pKa, anion-forming residues tyrosine and cysteine make good coupling partners for catalytic lysine residues. For the anion–cation pairs, the optimum difference in intrinsic pKas is a function of the energy of interaction between the residues. For the energy of interaction ε expressed in units of (ln 10)RT, the optimum difference in intrinsic pKas is within ∼1 pH unit of ε. 
    more » « less
  3. Two amino acid variants in soybean serine hydroxymethyltransferase 8 (SHMT8) are associated with resistance to the soybean cyst nematode (SCN), a devastating agricultural pathogen with worldwide economic impacts on soybean production. SHMT8 is a cytoplasmic enzyme that catalyzes the pyridoxal 5‐phosphate‐dependent conversion of serine and tetrahydrofolate (THF) to glycine and 5,10‐methylenetetrahydrofolate. A previous study of the P130R/N358Y double variant of SHMT8, identified in the SCN‐resistant soybean cultivar (cv.) Forrest, showed profound impairment of folate binding affinity and reduced THF‐dependent enzyme activity, relative to the highly active SHMT8 in cv. Essex, which is susceptible to SCN. Given the importance of SCN‐resistance in soybean agriculture, we report here the biochemical and structural characterization of the P130R and N358Y single variants to elucidate their individual effects on soybean SHMT8. We find that both single variants have reduced THF‐dependent catalytic activity relative to Essex SHMT8 (10‐ to 50‐fold decrease inkcat/Km) but are significantly more active than the P130R/N368Y double variant. The kinetic data also show that the single variants lack THF‐substrate inhibition as found in Essex SHMT8, an observation with implications for regulation of the folate cycle. Five crystal structures of the P130R and N358Y variants in complex with various ligands (resolutions from 1.49 to 2.30 Å) reveal distinct structural impacts of the mutations and provide new insights into allosterism. Our results support the notion that the P130R/N358Y double variant in Forrest SHMT8 produces unique and unexpected effects on the enzyme, which cannot be easily predicted from the behavior of the individual variants. 
    more » « less
  4. Brine shrimp (Artemia) are the only animals to thrive at sodium concentrations above 4 M. Salt excretion is powered by the Na+,K+-ATPase (NKA), a heterodimeric (αβ) pump that usually exports 3Na+in exchange for 2 K+per hydrolyzed ATP.Artemiaexpress several NKA catalytic α-subunit subtypes. High-salinity adaptation increases abundance of α2KK, an isoform that contains two lysines (Lys308 and Lys758 in transmembrane segments TM4 and TM5, respectively) at positions where canonical NKAs have asparagines (Xenopusα1’s Asn333 and Asn785). Using de novo transcriptome assembly and qPCR, we found thatArtemiaexpress two salinity-independent canonical α subunits (α1NNand α3NN), as well as two β variants, in addition to the salinity-controlled α2KK. These β subunits permitted heterologous expression of the α2KKpump and determination of its CryoEM structure in a closed, ion-free conformation, showing Lys758 residing within the ion-binding cavity. We used electrophysiology to characterize the function of α2KKpumps and compared it to that ofXenopusα1 (and its α2KK-mimicking single- and double-lysine substitutions). The double substitution N333K/N785K confers α2KK-like characteristics toXenopusα1, and mutant cycle analysis reveals energetic coupling between these two residues, illustrating how α2KK’s Lys308 helps to maintain high affinity for external K+when Lys758 occupies an ion-binding site. By measuring uptake under voltage clamp of the K+-congener86Rb+, we prove that double-lysine-substituted pumps transport 2Na+and 1 K+per catalytic cycle. Our results show how the two lysines contribute to generate a pump with reduced stoichiometry allowingArtemiato maintain steeper Na+gradients in hypersaline environments. 
    more » « less
  5. Carboxylases are biocatalysts that capture and convert carbon dioxide (CO2) under mild conditions and atmospheric concentrations at a scale of more than 400 Gt annually. However, how these enzymes bind and control the gaseous CO2molecule during catalysis is only poorly understood. One of the most efficient classes of carboxylating enzymes are enoyl-CoA carboxylases/reductases (Ecrs), which outcompete the plant enzyme RuBisCO in catalytic efficiency and fidelity by more than an order of magnitude. Here we investigated the interactions of CO2within the active site of Ecr fromKitasatospora setae. Combining experimental biochemistry, protein crystallography, and advanced computer simulations we show that 4 amino acids, N81, F170, E171, and H365, are required to create a highly efficient CO2-fixing enzyme. Together, these 4 residues anchor and position the CO2molecule for the attack by a reactive enolate created during the catalytic cycle. Notably, a highly ordered water molecule plays an important role in an active site that is otherwise carefully shielded from water, which is detrimental to CO2fixation. Altogether, our study reveals unprecedented molecular details of selective CO2binding and C–C-bond formation during the catalytic cycle of nature’s most efficient CO2-fixing enzyme. This knowledge provides the basis for the future development of catalytic frameworks for the capture and conversion of CO2in biology and chemistry. 
    more » « less