Protein design through systematic catalytic loop exchange in the (beta/alpha)8 fold

Phthalate hydrolase: distribution, diversity and molecular evolution

The alpha/beta-fold superfamily of hydrolases is rapidly becoming one of the largest groups of structurally related enzymes with diverse catalytic functions. In this superfamily of enzymes, esterase deserves special attention because of their wide distribution in biological systems and importance towards environmental and industrial applications. Among various esterases, phthalate hydrolases are the key alpha/beta enzymes involved in the metabolism of structurally diverse estrogenic phthalic acid esters, ubiquitously distributed synthetic chemicals, used as plasticizer in plastic manufacturing processes. Although they vary both at the sequence and functional levels, these hydrolases use a similar acid-base-nucleophile catalytic mechanism to catalyse reactions on structurally different substrates. The current review attempts to present insights on phthalate hydrolases, describing their sources, structural diversities, phylogenetic affiliations and catalytically different types or classes of enzymes, categorized as diesterase, monoesterase and diesterase-monoesterase, capable of hydrolysing phthalate diester, phthalate monoester and both respectively. Furthermore, available information on in silico analyses and site-directed mutagenesis studies revealing structure-function integrity and altered enzyme kinetics have been highlighted along with the possible scenario of their evolution at the molecular level.

A Novel Glutathione S-Transferase Gtt2 Class (VpGSTT2) Is Found in the Genome of the AHPND/EMS Vibrio parahaemolyticus Shrimp Pathogen

Glutathione S-transferases are a family of detoxifying enzymes that catalyze the conjugation of reduced glutathione (GSH) with different xenobiotic compounds using either Ser, Tyr, or Cys as a primary catalytic residue. We identified a novel GST in the genome of the shrimp pathogen V. parahaemolyticus FIM- S1708⁺, a bacterial strain associated with Acute Hepatopancreatic Necrosis Disease (AHPND)/Early Mortality Syndrome (EMS) in cultured shrimp. This new GST class was named Gtt2. It has an atypical catalytic mechanism in which a water molecule instead of Ser, Tyr, or Cys activates the sulfhydryl group of GSH. The biochemical properties of Gtt2 from Vibrio parahaemolyticus (VpGSTT2) were characterized using kinetic and crystallographic methods. Recombinant VpGSTT2 was enzymatically active using GSH and CDNB as substrates, with a specific activity of 5.7 units/mg. Low affinity for substrates was demonstrated using both Michaelis–Menten kinetics and isothermal titration calorimetry. The crystal structure showed a canonical two-domain structure comprising a glutathione binding G-domain and a hydrophobic ligand H domain. A water molecule was hydrogen-bonded to residues Thr9 and Ser 11, as reported for the yeast Gtt2, suggesting a primary role in the reaction. Molecular docking showed that GSH could bind at the G-site in the vicinity of Ser11. G-site mutationsT9A and S11A were analyzed. S11A retained 30% activity, while T9A/S11A showed no detectable activity. VpGSTT2 was the first bacterial Gtt2 characterized, in which residues Ser11 and Thr9 coordinated a water molecule as part of a catalytic mechanism that was characteristic of yeast GTT2. The GTT2 family has been shown to provide protection against metal toxicity; in some cases, excess heavy metals appear in shrimp ponds presenting AHPND/EMS. Further studies may address whether GTT2 in V. parahaemolyticus pathogenic strains may provide a competitive advantage as a novel detoxification mechanism.

Deletion and Randomization of Structurally Variable Regions in B. subtilis Lipase A (BSLA) Alter Its Stability and Hydrolytic Performance Against Long Chain Fatty Acid Esters

The continuous search for novel enzyme backbones and the engineering of already well studied enzymes for biotechnological applications has become an increasing challenge, especially by the increasing potential diversity space provided by directed enzyme evolution approaches and the demands of experimental data generated by rational design of enzymes. In this work, we propose a semi-rational mutational strategy focused on introducing diversity in structurally variable regions in enzymes. The identified sequences are subjected to a progressive deletion of two amino acids and the joining residues are subjected to saturation mutagenesis using NNK degenerate codons. This strategy offers a novel library diversity approach while simultaneously decreasing enzyme size in the variable regions. In this way, we intend to identify and reduce variable regions found in enzymes, probably resulting from neutral drift evolution, and simultaneously studying the functional effect of said regions. This strategy was applied to Bacillus. subtilis lipase A (BSLA), by selecting and deleting six variable enzyme regions (named regions 1 to 6) by the deletion of two amino acids and additionally randomizing the joining amino acid residues. After screening, no active variants were found in libraries 1% and 4%, 15% active variants were found in libraries 2% and 3%, and 25% for libraries 5 and 6 (n = 3000 per library, activity detected using tributyrin agar plates). Active variants were assessed for activity in microtiter plate assay (pNP-butyrate), thermal stability, substrate preference (pNP-butyrate, -palmitate), and compared to wildtype BSLA. From these analyses, variant P5F3 (F41L-ΔW42-ΔD43-K44P), from library 3 was identified, showing increased activity towards longer chain p-nitrophenyl fatty acid esters, when compared to BSLA. This study allowed to propose the targeted region 3 (positions 40-46) as a potential modulator for substrate specificity (fatty acid chain length) in BSLA, which can be further studied to increase its substrate spectrum and selectivity. Additionally, this variant showed a decreased thermal resistance but interestingly, higher isopropanol and Triton X-100 resistance. This deletion-randomization strategy could help to expand and explore sequence diversity, even in already well studied and characterized enzyme backbones such as BSLA. In addition, this strategy can contribute to investigate and identify important non-conserved regions in classic and novel enzymes, as well as generating novel biocatalysts with increased performance in specific processes, such as enzyme immobilization.

Synthesis of Human Milk Oligosaccharides: Protein Engineering Strategies for Improved Enzymatic Transglycosylation

Human milk oligosaccharides (HMOs) signify a unique group of oligosaccharides in breast milk, which is of major importance for infant health and development. The functional benefits of HMOs create an enormous impetus for biosynthetic production of HMOs for use as additives in infant formula and other products. HMO molecules can be synthesized chemically, via fermentation, and by enzymatic synthesis. This treatise discusses these different techniques, with particular focus on harnessing enzymes for controlled enzymatic synthesis of HMO molecules. In order to foster precise and high-yield enzymatic synthesis, several novel protein engineering approaches have been reported, mainly concerning changing glycoside hydrolases to catalyze relevant transglycosylations. The protein engineering strategies for these enzymes range from rationally modifying specific catalytic residues, over targeted subsite -1 mutations, to unique and novel transplantations of designed peptide sequences near the active site, so-called loop engineering. These strategies have proven useful to foster enhanced transglycosylation to promote different types of HMO synthesis reactions. The rationale of subsite -1 modification, acceptor binding site matching, and loop engineering, including changes that may alter the spatial arrangement of water in the enzyme active site region, may prove useful for novel enzyme-catalyzed carbohydrate design in general.

βαβ Super-Secondary Motifs: Sequence, Structural Overview, and Pursuit of Potential Autonomously Folding βαβ Sequences from (β/α)(8)/TIM Barrels

βαβ super-secondary structures constitute the basic building blocks of (β/α)₈ class of proteins. Despite the success in designing super-secondary structures, till date, there is not a single example of a natural βαβ sequence known to fold in isolation. In this chapter, to address the finding the "needles" in the haystack scenario, we have combined the sequence preferences and structural features of independent βαβ motifs, dictated by natural selection, with rationally derived parameters from a designed βαβ motif adopting stable fold in solution. Guided by this approach, a set of potential βαβ sequences from (β/α)₈/TIM barrels are proposed as likely candidates for autonomously folding based on the assessment of their foldability.

LoopX: A Graphical User Interface-Based Database for Comprehensive Analysis and Comparative Evaluation of Loops from Protein Structures

Due to their crucial role in function, folding, and stability, protein loops are being targeted for grafting/designing to create novel or alter existing functionality and improve stability and foldability. With a view to facilitate a thorough analysis and effectual search options for extracting and comparing loops for sequence and structural compatibility, we developed, LoopX a comprehensively compiled library of sequence and conformational features of ∼700,000 loops from protein structures. The database equipped with a graphical user interface is empowered with diverse query tools and search algorithms, with various rendering options to visualize the sequence- and structural-level information along with hydrogen bonding patterns, backbone φ, ψ dihedral angles of both the target and candidate loops. Two new features (i) conservation of the polar/nonpolar environment and (ii) conservation of sequence and conformation of specific residues within the loops have also been incorporated in the search and retrieval of compatible loops for a chosen target loop. Thus, the LoopX server not only serves as a database and visualization tool for sequence and structural analysis of protein loops but also aids in extracting and comparing candidate loops for a given target loop based on user-defined search options.

Correlation of fitness landscapes from three orthologous TIM barrels originates from sequence and structure constraints

Sequence divergence of orthologous proteins enables adaptation to environmental stresses and promotes evolution of novel functions. Limits on evolution imposed by constraints on sequence and structure were explored using a model TIM barrel protein, indole-3-glycerol phosphate synthase (IGPS). Fitness effects of point mutations in three phylogenetically divergent IGPS proteins during adaptation to temperature stress were probed by auxotrophic complementation of yeast with prokaryotic, thermophilic IGPS. Analysis of beneficial mutations pointed to an unexpected, long-range allosteric pathway towards the active site of the protein. Significant correlations between the fitness landscapes of distant orthologues implicate both sequence and structure as primary forces in defining the TIM barrel fitness landscape and suggest that fitness landscapes can be translocated in sequence space. Exploration of fitness landscapes in the context of a protein fold provides a strategy for elucidating the sequence-structure-fitness relationships in other common motifs.

The TIM barrel fold is an evolutionarily conserved motif found in proteins with a variety of enzymatic functions. Here the authors explore the fitness landscape of the TIM barrel protein IGPS and uncover evolutionary constraints on both sequence and structure, accompanied by long range allosteric interactions.

Active Hydrogen Bond Network (AHBN) and Applications for Improvement of Thermal Stability and pH-Sensitivity of Pullulanase from Bacillus naganoensis

A method, so called “active hydrogen bond network” (AHBN), is proposed for site-directed mutations of hydrolytic enzymes. In an enzyme the AHBN consists of the active residues, functional residues, and conservative water molecules, which are connected by hydrogen bonds, forming a three dimensional network. In the catalysis hydrolytic reactions of hydrolytic enzymes AHBN is responsible for the transportation of protons and water molecules, and maintaining the active and dynamic structures of enzymes. The AHBN of pullulanase BNPulA324 from Bacillus naganoensis was constructed based on a homologous model structure using Swiss Model Protein-modeling Server according to the template structure of pullulanase BAPulA (2WAN). The pullulanase BNPulA324 are mutated at the mutation sites selected by means of the AHBN method. Both thermal stability and pH-sensitivity of pullulanase BNPulA324 were successfully improved. The mutations at the residues located at the out edge of AHBN may yield positive effects. On the other hand the mutations at the residues inside the AHBN may deprive the bioactivity of enzymes. The AHBN method, proposed in this study, may provide an assistant and alternate tool for protein rational design and protein engineering.

Structure-activity relationships on the study of β-galactosidase folding/unfolding due to interactions with immobilization additives: Triton X-100 and ethanol

Improving the enzyme stability is a challenge for allowing their practical application. The surfactants are stabilizing agents, however, there are still questions about their influence on enzyme properties. The structure-activity/stability relationship for β-galactosidase from Bacillus circulans is studied here by Circular Dichroism and activity measurements, as a function of temperature and pH. The tendency of preserving the β-sheet and α-helix structures at temperatures below 65°C and different pH is the result of the balance between the large- and short-range effects, respecting to the active site. This information is fundamental for explaining the structural changes of this enzyme in the presence of Triton X-100 surfactant and ethanol. The enzyme thermal stabilization in the presence of this surfactant responds to the rearrangement of the secondary structure for having optimal activity/stability. The effect of ethanol is more related with changes in the dielectric properties of the aqueous solution than with protein structural transformations. These results contribute to understand the effects of surfactant-enzyme interactions on the enzyme behavior, from the structural point of view and to rationalize the surfactant-based stabilizing strategies for β-galactosidades.

De novo design of a four-fold symmetric TIM-barrel protein with atomic-level accuracy

Despite efforts for over 25 years, de novo protein design has not succeeded in achieving the TIM-barrel fold. Here we describe the computational design of four-fold symmetrical (β/α)8 barrels guided by geometrical and chemical principles. Experimental characterization of 33 designs revealed the importance of side chain-backbone hydrogen bonds for defining the strand register between repeat units. The X-ray crystal structure of a designed thermostable 184-residue protein is nearly identical to that of the designed TIM-barrel model. PSI-BLAST searches do not identify sequence similarities to known TIM-barrel proteins, and sensitive profile-profile searches indicate that the design sequence is distant from other naturally occurring TIM-barrel superfamilies, suggesting that Nature has sampled only a subset of the sequence space available to the TIM-barrel fold. The ability to design TIM barrels de novo opens new possibilities for custom-made enzymes.

Growth of Chitinophaga pinensis on Plant Cell Wall Glycans and Characterisation of a Glycoside Hydrolase Family 27 β-l-Arabinopyranosidase Implicated in Arabinogalactan Utilisation

The genome of the soil bacterium Chitinophaga pinensis encodes a diverse array of carbohydrate active enzymes, including nearly 200 representatives from over 50 glycoside hydrolase (GH) families, the enzymology of which is essentially unexplored. In light of this genetic potential, we reveal that C. pinensis has a broader saprophytic capacity to thrive on plant cell wall polysaccharides than previously reported, and specifically that secretion of β-l-arabinopyranosidase activity is induced during growth on arabinogalactan. We subsequently correlated this activity with the product of the Cpin_5740 gene, which encodes the sole member of glycoside hydrolase family 27 (GH27) in C. pinensis, CpArap27. Historically, GH27 is most commonly associated with α-d-galactopyranosidase and α-d-N-acetylgalactosaminidase activity. A new phylogenetic analysis of GH27 highlighted the likely importance of several conserved secondary structural features in determining substrate specificity and provides a predictive framework for identifying enzymes with the less common β-l-arabinopyranosidase activity.

Highly diverse protein library based on the ubiquitous (β/α)₈ enzyme fold yields well-structured proteins through in vitro folding selection

Proper protein folding is a prerequisite for protein stability and enzymatic activity. Although directed evolution can be a powerful tool to investigate enzymatic function and to isolate novel activities, well-designed libraries of folded proteins are essential. In vitro selection methods are particularly capable of searching for enzymatic activities in libraries of trillions of protein variants, yet high-quality libraries of well-folded enzymes with such high diversity are lacking. We describe the construction and detailed characterization of a folding-enriched protein library based on the ubiquitous (β/α)₈ barrel fold, which is found in five of the six enzyme classes. We introduced seven randomized loops on the catalytic face of the monomeric, thermostable (β/α)₈ barrel of glycerophosphodiester phosphodiesterase (GDPD) from Thermotoga maritima. We employed in vitro folding selection based on protease digestion to enrich intermediate libraries containing three to four randomized loops for folded variants, and then combined them to assemble the final library (10¹⁴ DNA sequences). The resulting library was analyzed by using the in vitro protease assay and an in vivo GFP-folding assay; it contains ∼10¹² soluble monomeric protein variants. We isolated six library members and demonstrated that these proteins are soluble, monomeric and show (β/α)₈-barrel fold-like secondary and tertiary structure. The quality of the folding-enriched library improved up to 50-fold compared to a control library that was assembled without the folding selection. To the best of our knowledge, this work is the first example of combining the ultra-high throughput mRNA display method with selection for folding. The resulting (β/α)₈ barrel libraries provide a valuable starting point to study the unique catalytic capabilities of the (β/α)₈ fold, and to isolate novel enzymes.

Alternative splice variants in TIM barrel proteins from human genome correlate with the structural and evolutionary modularity of this versatile protein fold

After the surprisingly low number of genes identified in the human genome, alternative splicing emerged as a major mechanism to generate protein diversity in higher eukaryotes. However, it is still not known if its prevalence along the genome evolution has contributed to the overall functional protein diversity or if it simply reflects splicing noise. The (βα)₈ barrel or TIM barrel is one of the most frequent, versatile, and ancient fold encountered among enzymes. Here, we analyze the structural modifications present in TIM barrel proteins from the human genome product of alternative splicing events. We found that 87% of all splicing events involved deletions; most of these events resulted in protein fragments that corresponded to the (βα)₂, (βα)₄, (βα)₅, (βα)₆, and (βα)₇ subdomains of TIM barrels. Because approximately 7% of all the splicing events involved internal β-strand substitutions, we decided, based on the genomic data, to design β-strand and α-helix substitutions in a well-studied TIM barrel enzyme. The biochemical characterization of one of the chimeric variants suggests that some of the splice variants in the human genome with β-strand substitutions may be evolving novel functions via either the oligomeric state or substrate specificity. We provide results of how the splice variants represent subdomains that correlate with the independently folding and evolving structural units previously reported. This work is the first to observe a link between the structural features of the barrel and a recurrent genetic mechanism. Our results suggest that it is reasonable to expect that a sizeable fraction of splice variants found in the human genome represent structurally viable functional proteins. Our data provide additional support for the hypothesis of the origin of the TIM barrel fold through the assembly of smaller subdomains. We suggest a model of how nature explores new proteins through alternative splicing as a mechanism to diversify the proteins encoded in the human genome.

Evolution of a new function in an esterase: simple amino acid substitutions enable the activity present in the larger paralog, BioH

Gene duplication and divergence are essential processes for the evolution of new activities. Divergence may be gradual, involving simple amino acid residue substitutions, or drastic, such that larger structural elements are inserted, deleted or rearranged. Vast protein sequence comparisons, supported by some experimental evidence, argue that large structural modifications have been necessary for certain catalytic activities to evolve. However, it is not clear whether these activities could not have been attained by gradual changes. Interestingly, catalytic promiscuity could play a fundamental evolutionary role: a preexistent secondary activity could be increased by simple amino acid residue substitutions that do not affect the enzyme's primary activity. The promiscuous profile of the enzyme may be modified gradually by genetic drift, making a pool of potentially useful activities that can be selected before duplication. In this work, we used random mutagenesis and in vivo selection to evolve the Pseudomonas aeruginosa PAO1 carboxylesterase PA3859, a small protein, to attain the function of BioH, a much larger paralog involved in biotin biosynthesis. BioH was chosen as a target activity because it provides a highly sensitive selection for evolved enzymatic activities by auxotrophy complementation. After only two cycles of directed evolution, mutants with the ability to efficiently complement biotin auxotrophy were selected. The in vivo and in vitro characterization showed that the activity of one of our mutant proteins was similar to that of the wild-type BioH enzyme. Our results demonstrate that it is possible to evolve enzymatic activities present in larger proteins by discrete amino acid substitutions.

Exploring the Structure-Function Loop Adaptability of a (β/α)(8)-Barrel Enzyme through Loop Swapping and Hinge Variability

Evolution of proteins involves sequence changes that are frequently localized at loop regions, revealing their important role in natural evolution. However, the development of strategies to understand and imitate such events constitutes a challenge to design novel enzymes in the laboratory. In this study, we show how to adapt loop swapping as semiautonomous units of functional groups in an enzyme with the (β/α)(8)-barrel and how this functional adaptation can be measured in vivo. To mimic the natural mechanism providing loop variability in antibodies, we developed an overlap PCR strategy. This includes introduction of sequence diversity at two hinge residues, which connect the new loops with the rest of the protein scaffold, and we demonstrate that this is necessary for a successful exploration of functional sequence space. This design allowed us to explore the sequence requirements to functional adaptation of each loop replacement that may not be sampled otherwise. Libraries generated following this strategy were evaluated in terms of their folding competence and their functional proficiency, an observation that was formalized as a Structure-Function Loop Adaptability value. Molecular details about the function and structure of some variants were obtained by enzyme kinetics and circular dichroism. This strategy yields functional variants that retain the original activity at higher frequencies, suggesting a new strategy for protein engineering that incorporates a more divergent sequence exploration beyond that limited to point mutations. We discuss how this approach may provide insights into the mechanism of enzyme evolution and function.

Identification and analysis of residues contained on beta --> alpha loops of the dual-substrate (beta alpha)8 phosphoribosyl isomerase A specific for its phosphoribosyl anthranilate isomerase activity

A good model to experimentally explore evolutionary hypothesis related to enzyme function is the ancient-like dual-substrate (beta alpha)(8) phosphoribosyl isomerase A (PriA), which takes part in both histidine and tryptophan biosynthesis in Streptomyces coelicolor and related organisms. In this study, we determined the Michaelis-Menten enzyme kinetics for both isomerase activities in wild-type PriA from S. coelicolor and in selected single-residue monofunctional mutants, identified after Escherichia coli in vivo complementation experiments. Structural and functional analyses of a hitherto unnoticed residue contained on the functionally important beta --> alpha loop 5, namely, Arg(139), which was postulated on structural grounds to be important for the dual-substrate specificity of PriA, is presented for the first time. Indeed, enzyme kinetics analyses done on the mutant variants PriA_Ser(81)Thr and PriA_Arg(139)Asn showed that these residues, which are contained on beta --> alpha loops and in close proximity to the N-terminal phosphate-binding site, are essential solely for the phosphoribosyl anthranilate isomerase activity of PriA. Moreover, analysis of the X-ray crystallographic structure of PriA_Arg(139)Asn elucidated at 1.95 A herein strongly implicates the occurrence of conformational changes in this beta --> alpha loop as a major structural feature related to the evolution of the dual-substrate specificity of PriA. It is suggested that PriA has evolved by tuning a fine energetic balance that allows the sufficient degree of structural flexibility needed for accommodating two topologically dissimilar substrates--within a bifunctional and thus highly constrained active site--without compromising its structural stability.

Directed evolution of a quorum-quenching lactonase from Mycobacterium avium subsp. paratuberculosis K-10 in the amidohydrolase superfamily

The PLL(PTE-like lactonase)-group of enzymes within the amidohydrolase superfamily hydrolyze N-acyl-homoserine lactones (AHLs) that are involved in bacterial quorum-sensing pathways. These enzymes possess the (beta/alpha)(8)-barrel fold and serve as attractive templates for in vitro evolution and engineering of quorum-quenching biological molecules that can serve as antivirulence therapeutic agents. Using a quorum-quenching lactonase from Mycobacterium avium subsp. paratuberculosis K-10 (GI: 41409766) as the initial template for in vitro evolution experiments, we enhanced the catalytic efficiency and increased the substrate range of the wild-type enzyme through a single point mutation on the loop at the C-terminal end of the eighth beta-strand. This N266Y mutant had an increased value of k(cat)/K(M) of 30- and 32-fold toward 3-oxo-N-octanoyl-l-homoserine lactone and N-hexanoyl-l-homoserine lactone, respectively; the evolved mutant also exhibited lactonase activity toward 3-oxo-N-hexanoyl-l-homoserine lactone and N-butyryl-l-homoserine lactone, AHLs that were previously not hydrolyzed by the wild-type enzyme. This article reinforces the evolutionary potential of the (beta/alpha)(8)-barrel fold and highlights the possibility of using quorum-quenching lactonases in the amidohydrolase superfamily as templates for engineering biomolecules of therapeutic use.