Cooperativity and flexibility in enzyme evolution

Structural and biochemical consequences of pyridoxine-dependent epilepsy mutations that target the aldehyde binding site of aldehyde dehydrogenase ALDH7A1

In humans, certain mutations in the gene encoding aldehyde dehydrogenase 7A1 are associated with pyridoxine-dependent epilepsy (PDE). Understanding the impact of PDE-causing mutations on the structure and activity of ALDH7A1 could allow for the prediction of symptom-severity and aid the development of patient-specific medical treatments. Herein, we investigate the biochemical and structural consequences of PDE missense mutations targeting residues in the aldehyde substrate binding site: N167S, P169S, A171V, G174V, and W175G. All but G174V could be purified for biochemical and X-ray crystallographic analysis. W175G has a relatively mild kinetic defect, exhibiting a fivefold decrease in k_cat with no change in K_m . P169S and N167S have moderate defects, characterized by catalytic efficiencies of 20- and 100-times lower than wild-type, respectively. A171V has a profound functional defect, with catalytic efficiency 2000-times lower than wild-type. The crystal structures of the variants are the first for any PDE-associated mutant of ALDH7A1. The structures show that missense mutations that decrease the steric bulk of the side chain tend to create a cavity in the active site. The protein responds by relaxing into the vacant space, and this structural perturbation appears to cause misalignment of the aldehyde substrate in W175G and N167S. The P169S structure is nearly identical to that of the wild-type enzyme; however, analysis of B-factors suggests the catalytic defect may result from altered protein dynamics. The A171V structure suggests that the potential for steric clash with Val171 prevents Glu121 from ion pairing with the amino group of the aldehyde substrate. ENZYMES: Aldehyde dehydrogenase 7A1 (EC1.2.1.31). DATABASES: Coordinates have been deposited in the Protein Data Bank under the following accession codes: 6O4B, 6O4C, 6O4D, 6O4E, 6O4F, 6O4G, 6O4H.

NAD+ -Dependent Enzymatic Route for the Epimerization of Hydroxysteroids

Epimerization of cholic and chenodeoxycholic acid (CA and CDCA, respectively) is a notable conversion for the production of ursodeoxycholic acid (UDCA). Two enantiocomplementary hydroxysteroid dehydrogenases (7α- and 7β-HSDHs) can carry out this transformation fully selectively by specific oxidation of the 7α-OH group of the substrate and subsequent reduction of the keto intermediate to the final product (7β-OH). With a view to developing robust and active biocatalysts, novel NADH-active 7β-HSDH species are necessary to enable a solely NAD+ -dependent redox-neutral cascade for UDCA production. A wild-type NADH-dependent 7β-HSDH from Lactobacillus spicheri (Ls7β-HSDH) was identified, recombinantly expressed, purified, and biochemically characterized. Using this novel NAD+ -dependent 7β-HSDH enzyme in combination with 7α-HSDH from Stenotrophomonas maltophilia permitted the biotransformations of CA and CDCA in the presence of catalytic amounts of NAD+ , resulting in high yields (>90 %) of UCA and UDCA.

Protein Flexibility and Stiffness Enable Efficient Enzymatic Catalysis

The enormous rate accelerations observed for many enzyme catalysts are due to strong stabilizing interactions between the protein and reaction transition state. The defining property of these catalysts is their specificity for binding the transition state with a much higher affinity than substrate. Experimental results are presented which show that the phosphodianion-binding energy of phosphate monoester substrates is used to drive conversion of their protein catalysts from flexible and entropically rich ground states to stiff and catalytically active Michaelis complexes. These results are generalized to other enzyme-catalyzed reactions. The existence of many enzymes in flexible, entropically rich, and inactive ground states provides a mechanism for utilization of ligand-binding energy to mold these catalysts into stiff and active forms. This reduces the substrate-binding energy expressed at the Michaelis complex, while enabling the full and specific expression of large transition-state binding energies. Evidence is presented that the complexity of enzyme conformational changes increases with increases in the enzymatic rate acceleration. The requirement that a large fraction of the total substrate-binding energy be utilized to drive conformational changes of floppy enzymes is proposed to favor the selection and evolution of protein folds with multiple flexible unstructured loops, such as the TIM-barrel fold. The effect of protein motions on the kinetic parameters for enzymes that undergo ligand-driven conformational changes is considered. The results of computational studies to model the complex ligand-driven conformational change in catalysis by triosephosphate isomerase are presented.

Local frustration around enzyme active sites

Conflicting biological goals often meet in the specification of protein sequences for structure and function. Overall, strong energetic conflicts are minimized in folded native states according to the principle of minimal frustration, so that a sequence can spontaneously fold, but local violations of this principle open up the possibility to encode the complex energy landscapes that are required for active biological functions. We survey the local energetic frustration patterns of all protein enzymes with known structures and experimentally annotated catalytic residues. In agreement with previous hypotheses, the catalytic sites themselves are often highly frustrated regardless of the protein oligomeric state, overall topology, and enzymatic class. At the same time a secondary shell of more weakly frustrated interactions surrounds the catalytic site itself. We evaluate the conservation of these energetic signatures in various family members of major enzyme classes, showing that local frustration is evolutionarily more conserved than the primary structure itself.

Protein engineers turned evolutionists-the quest for the optimal starting point

The advent of laboratory directed evolution yielded a fruitful crosstalk between the disciplines of molecular evolution and bio-engineering. Here, we outline recent developments in both disciplines with respect to how one can identify the best starting points for directed evolution, such that highly efficient and robust tailor-made enzymes can be obtained with minimal optimization. Directed evolution studies have highlighted essential features of engineer-able enzymes: highly stable, mutationally robust enzymes with the capacity to accept a broad range of substrates. Robust, evolvable enzymes can be inferred from the natural sequence record. Broad substrate spectrum relates to conformational plasticity and can also be predicted by phylogenetic analyses and/or by computational design. Overall, an increasingly powerful toolkit is becoming available for identifying optimal starting points including network analyses of enzyme superfamilies and other bioinformatics methods.

Loop Motion in Triosephosphate Isomerase Is Not a Simple Open and Shut Case

Conformational changes are crucial for the catalytic action of many enzymes. A prototypical and well-studied example is loop opening and closure in triosephosphate isomerase (TIM), which is thought to determine the rate of catalytic turnover in many circumstances. Specifically, TIM loop 6 "grips" the phosphodianion of the substrate and, together with a change in loop 7, sets up the TIM active site for efficient catalysis. Crystal structures of TIM typically show an open or a closed conformation of loop 6, with the tip of the loop moving ∼7 Å between conformations. Many studies have interpreted this motion as a two-state, rigid-body transition. Here, we use extensive molecular dynamics simulations, with both conventional and enhanced sampling techniques, to analyze loop motion in apo and substrate-bound TIM in detail, using five crystal structures of the dimeric TIM from Saccharomyces cerevisiae. We find that loop 6 is highly flexible and samples multiple conformational states. Empirical valence bond simulations of the first reaction step show that slight displacements away from the fully closed-loop conformation can be sufficient to abolish most of the catalytic activity; full closure is required for efficient reaction. The conformational change of the loops in TIM is thus not a simple "open and shut" case and is crucial for its catalytic action. Our detailed analysis of loop motion in a highly efficient enzyme highlights the complexity of loop conformational changes and their role in biological catalysis.

Tuning site-specific dynamics to drive allosteric activation in a pneumococcal zinc uptake regulator

MarR (multiple antibiotic resistance repressor) family proteins are bacterial repressors that regulate transcription in response to a wide range of chemical signals. Although specific features of MarR family function have been described, the role of atomic motions in MarRs remains unexplored thus limiting insights into the evolution of allostery in this ubiquitous family of repressors. Here, we provide the first experimental evidence that internal dynamics play a crucial functional role in MarR proteins. Streptococcus pneumoniae AdcR (adhesin-competence repressor) regulates Zn^II homeostasis and Zn^II functions as an allosteric activator of DNA binding. Zn^II coordination triggers a transition from somewhat independent domains to a more compact structure. We identify residues that impact allosteric activation on the basis of Zn^II-induced perturbations of atomic motions over a wide range of timescales. These findings appear to reconcile the distinct allosteric mechanisms proposed for other MarRs and highlight the importance of conformational dynamics in biological regulation.

Molecular modeling of conformational dynamics and its role in enzyme evolution

With increasing computational power, biomolecular simulations have become an invaluable tool for understanding enzyme mechanisms and the origins of enzyme catalysis. More recently, computational studies have started to focus on understanding how enzyme activity itself evolves, both in terms of enhancing the native or new activities on existing enzyme scaffolds, or completely de novo on previously non-catalytic scaffolds. In this context, both experiment and molecular modeling provided strong evidence for an important role of conformational dynamics in the evolution of enzyme functions. This contribution will present a brief overview of the current state of the art for computationally exploring enzyme conformational dynamics in enzyme evolution, and, using several showcase studies, illustrate the ways molecular modeling can be used to shed light on how enzyme function evolves, at the most fundamental molecular level.

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

Conformational Sampling of the Intrinsically Disordered C-Terminal Tail of DERA Is Important for Enzyme Catalysis

2-Deoxyribose-5-phosphate aldolase (DERA) catalyzes the reversible conversion of acetaldehyde and glyceraldehyde-3-phosphate into deoxyribose-5-phosphate. DERA is used as a biocatalyst for the synthesis of drugs such as statins and is a promising pharmaceutical target due to its involvement in nucleotide catabolism. Despite previous biochemical studies suggesting the catalytic importance of the C-terminal tyrosine residue found in several bacterial DERAs, the structural and functional basis of its participation in catalysis remains elusive because the electron density for the last eight to nine residues (i.e., the C-terminal tail) is absent in all available crystal structures. Using a combination of NMR spectroscopy and molecular dynamics simulations, we conclusively show that the rarely studied C-terminal tail of E. coli DERA (ecDERA) is intrinsically disordered and exists in equilibrium between open and catalytically relevant closed states, where the C-terminal tyrosine (Y259) enters the active site. Nuclear Overhauser effect distance restraints, obtained due to the presence of a substantial closed state population, were used to derive the solution-state structure of the ecDERA closed state. Real-time NMR hydrogen/deuterium exchange experiments reveal that Y259 is required for efficiency of the proton abstraction step of the catalytic reaction. Phosphate titration experiments show that, in addition to the phosphate-binding residues located near the active site, as observed in the available crystal structures, ecDERA contains previously unknown auxiliary phosphate-binding residues on the C-terminal tail which could facilitate in orienting Y259 in an optimal position for catalysis. Thus, we present significant insights into the structural and mechanistic importance of the ecDERA C-terminal tail and illustrate the role of conformational sampling in enzyme catalysis.

Biotechnological and protein-engineering implications of ancestral protein resurrection

Approximations to the sequences of ancestral proteins can be derived from the sequences of their modern descendants. Proteins encoded by such reconstructed sequences can be prepared in the laboratory and subjected to experimental scrutiny. These 'resurrected' ancestral proteins often display remarkable properties, reflecting ancestral adaptations to intra-cellular and extra-cellular environments that differed from the environments hosting modern/extant proteins. Recent experimental and computational work has specifically discussed high stability, substrate and catalytic promiscuity, conformational flexibility/diversity and altered patterns of interaction with other sub-cellular components. In this review, we discuss these remarkable properties as well as recent attempts to explore their biotechnological and protein-engineering potential.