A Pareto-Optimal Refinement Method for Protein Design Scaffolds

Computational design of protein function involves a search for amino acids with the lowest energy subject to a set of constraints specifying function. In many cases a set of natural protein backbone structures, or “scaffolds”, are searched to find regions where functional sites (an enzyme active site, ligand binding pocket, protein – protein interaction region, etc.) can be placed, and the identities of the surrounding amino acids are optimized to satisfy functional constraints. Input native protein structures almost invariably have regions that score very poorly with the design force field, and any design based on these unmodified structures may result in mutations away from the native sequence solely as a result of the energetic strain. Because the input structure is already a stable protein, it is desirable to keep the total number of mutations to a minimum and to avoid mutations resulting from poorly-scoring input structures. Here we describe a protocol using cycles of minimization with combined backbone/sidechain restraints that is Pareto-optimal with respect to RMSD to the native structure and energetic strain reduction. The protocol should be broadly useful in the preparation of scaffold libraries for functional site design.     

An active enzyme constructed from a 9-amino acid alphabet

Nature employs a set of 20 amino acids to produce a repertoire of protein structures endowed with sophisticated functions. Here, we combined design and selection to create an enzyme composed entirely from a set of only 9 amino acids that can rescue auxotrophic cells lacking chorismate mutase. The simplified protein captures key structural features of its natural counterpart but appears to be somewhat less stable and more flexible. The potential of a dramatically reduced amino acid alphabet to produce an active catalyst supports the notion that primordial enzymes may have possessed low amino acid diversity and suggests that combinatorial engineering strategies, such as the one used here, may be generally applied to create enzymes with novel structures and functions.     

Genomic analysis of carbon source metabolism of Shewanella oneidensis MR-1: Predictions versus experiments

Genomic sequences have been used to find the genetic foundation for carbon source metabolism in Shewanella oneidensis MR-1. Annotated S. oneidensis MR-1 gene products were examined for their sequence similarity to enzymes participating in pathways for utilization of carbon and energy as described in the BioCyc database (http://www.biocyc.org/) or in the primary literature. A picture emerges that relegates five- and six-carbon sugars to minor roles as carbon sources, whereas multiple pathways for utilization of up to three-carbon carbohydrates seem to be present. Capacity to utilize amino acids for carbon and energy is also present. A few contradictions emerged in which enzymes appear to be present by annotations but are not active in the cell according to physiological experiments. Annotations are based on close sequence similarity and will not reveal inactivity due to deleterious mutations or due to lack of coordination of regulation and transport. Genes for a few enzymes known by experiment to be active are not found in the genome. This may be due to extensive divergence after duplication or convergence of function in separate lines in evolution rendering activities undetectable by sequence similarity. To minimize false predictions from protein sequences, we have been conservative in predicting pathways. We did not predict any pathway when, although a partial pathway was seen it was composed largely of enzymes already accounted for in any other complete pathway. This is an example of how a biochemically oriented sequence analysis can generate questions and direct further experimental investigation.     

Relationships between the osmoadaptation strategy, amino acid composition of total cellular protein, and properties of certain enzymes of haloalkaliphilic bacteria

Haloalkaliphilic microorganisms isolated from soda lakes were compared in terms of the amino acid composition of total cellular protein and the reaction of a number of key enzymes to salts and pH of the medium. In the extremely halophilic bacterium Natroniella acetigena (salt-inside osmoadaptation strategy), acidic amino acids (glutamic and aspartic) made up 30.91 mol % of the total of cellular protein amino acids. In the moderate haloalkaliphiles Tindallia magadiensis, Halomonas campisalis, and Halomonas sp. AIR-1 (compatible-solutes osmoadaptation strategy), the proportion of acidic amino acids (24.36, 23.15, and 23.58 mol %, respectively) was lower than in N. acetigena but higher than in the freshwater Acetobacterium paludosum (20.77 mol %). The excess of acidic amino acids over basic amino acids (lysine and arginine) increased with the degree of halophily. The enzymes of haloalkaliphiles proved to be tolerant to salts and high pH values, although the degree of tolerance varied. The activity of N. acetigena CO dehydrogenase was maximum in the presence of 0.7 M NaCl, but it was virtually independent of the NaHCO3 concentration. The hydrogenase and CO dehydrogenase of T. magadiensis exhibited maximum activity in the absence of NaCl; the CO dehydrogenase was most active at 0.25 M NaHCO3, and hydrogenase activity was only weakly dependent on NaHCO3 in the concentration range of 0-1.2 M. The nitrate reductases of H. campisalis and Halomonas sp. AIR-2 were active in broad ranges of NaCl and KCl concentrations; the activity maxima were recorded at moderate concentrations of these salts. The pH optima of most of the studied enzymes of haloalkaliphiles were in the alkaline zone. Thus, it was shown that the amino acid composition of total cellular protein is determined by the osmoadaptation strategy employed by the bacterium. A correlation was found between the salt tolerance of enzymes and the proportion of acidic amino acids in the total cellular protein. The ability of enzymes to function at high pH values is one of the mechanisms of adaptation of microorganisms to high pH values.     

Distance‐dependent classification of amino acids by information theory

Reduced amino acid alphabets are useful to understand molecular evolution as they reveal basal, shared properties of amino acids, which the structures and functions of proteins rely on. Several previous studies derived such reduced alphabets and linked them to the origin of life and biotechnological applications. However, all this previous work presupposes that only direct contacts of amino acids in native protein structures are relevant. We show in this work, using information–theoretical measures, that an appropriate alphabet reduction scheme is in fact a function of the maximum distance amino acids interact at. Although for small distances our results agree with previous ones, we show how long‐range interactions change the overall picture and prompt for a revised understanding of the protein design process. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Triosephosphate isomerase by consensus design: dramatic differences in physical properties and activity of related variants

Consensus design, the selection of mutations based on the most common amino acid in each position of a multiple sequence alignment, has proven to be an efficient way to engineer stabilized mutants and even to design entire proteins. However, its application has been limited to small motifs or small families of highly related proteins. Also, we have little idea of how information that specifies a protein's properties is distributed between positional effects (consensus) and interactions between positions (correlated occurrences of amino acids). Here, we designed several consensus variants of triosephosphate isomerase (TIM), a large, diverse family of complex enzymes. The first variant was only weakly active, had molten globular characteristics, and was monomeric at 25 °C despite being based on nearly all dimeric enzymes. A closely related variant from curation of the sequence database resulted in a native-like dimeric TIM with near-diffusion-controlled kinetics. Both enzymes vary substantially (30–40%) from any natural TIM, but they differ from each other in only a relatively small number of unconserved positions. We demonstrate that consensus design is sufficient to engineer a sophisticated protein that requires precise substrate positioning and coordinated loop motion. The difference in oligomeric states and native-like properties for the two consensus variants is not a result of defects in the dimerization interface but rather disparate global properties of the proteins. These results have important implications for the role of correlated amino acids, the ability of TIM to function as a monomer, and the ability of molten globular proteins to carry out complex reactions.     

A functional core of the mitochondrial genome maintenance protein Mgm101p in Saccharomyces cerevisiae determined with a temperature-conditional allele

Analysis of Mgm101p isolated from mitochondria shows that the mature protein of 27.6 kDa lacks 22 amino acids from the N-terminus. This mitochondrial targeting sequence has been incorporated in the design of oligonucleotides used to determine a functional core of Mgm101p. Progressive deletions, although retaining the targeting sequence, reveal that 76 N-terminal and six C-terminal amino acids of Mgm101p can be removed without altering the ability to complement an mgm101-1(ts) temperature-sensitive mutant. However, this active core is unable to complement mgm101 null mutants, suggesting that the Mgm101p might need to form a dimer or multimer to be functional in vivo. The active core, enriched in basic residues, contains 165 amino acids with a pI of 9.2. Alignment with 22 Mgm101p sequences from other lower eukaryotes shows that a number of amino acids are highly conserved in this region. Random mutagenesis confirms that certain critical amino acids required for function are invariant across the 23 proteins. Searches in the PFAM database revealed a low level of structural similarity between the active core and the Rad52 protein family.     

A viable amino acid editing activity in the leucyl-tRNA synthetase CP1-splicing domain is not required in the yeast mitochondria

Aminoacyl-tRNA synthetases are a family of enzymes that are responsible for translating the genetic code in the first step of protein synthesis. Some aminoacyl-tRNA synthetases have editing activities to clear their mistakes and enhance fidelity. Leucyl-tRNA synthetases have a hydrolytic active site that resides in a discrete amino acid editing domain called CP1. Mutational analysis within yeast mitochondrial leucyl-tRNA synthetase showed that the enzyme has maintained an editing active site that is competent for post-transfer editing of mischarged tRNA similar to other leucyl-tRNA synthetases. These mutations that altered or abolished leucyl-tRNA synthetase editing were introduced into complementation assays. Cell viability and mitochondrial function were largely unaffected in the presence of high levels of non-leucine amino acids. In contrast, these editing-defective mutations limited cell viability in Escherichia coli. It is possible that the yeast mitochondria have evolved to tolerate lower levels of fidelity in protein synthesis or have developed alternate mechanisms to enhance discrimination of leucine from non-cognate amino acids that can be misactivated by leucyl-tRNA synthetase.     

Structural and functional comparisons between vanadium haloperoxidase and acid phosphatase enzymes

The crystallographic structures of both the vanadium chloroperoxidase and bromoperoxidase enzymes have been determined with either vanadium or phosphate bound at their active site. The amino acids that are involved in phosphate binding in the acid phosphatase enzymes and those that are coordinated to vanadium in the haloperoxidases appear to be conserved between the two classes of enzyme. The detailed active site architecture for enzymes that recognize and use either vanadium or phosphate will be discussed in relation to their proposed enzymatic mechanism.     

Mutagenesis study on amino acids around the molybdenum centre of the periplasmic nitrate reductase from Ralstonia eutropha

Molybdenum enzymes containing the pterin cofactor are a diverse group of enzymes that catalyse in general oxygen atom transfer reactions. Aiming at studying the amino acid residues, which are important for the enzymatic specificity, we used nitrate reductase from Ralstonia eutropha (R.e.NAP) as a model system for mutational studies at the active site. We mutated amino acids at the Mo active site (Cys181 and Arg421) as well as amino acids in the funnel leading to it (Met182, Asp196, Glu197, and the double mutant Glu197-Asp196). The mutations were made on the basis of the structural comparison of nitrate reductases with formate dehydrogenases (FDH), which show very similar three-dimensional structures, but clear differences in amino acids surrounding the active site. For mutations Arg421Lys and Glu197Ala we found a reduced nitrate activity while the other mutations resulted in complete loss of activity. In spite of the partial of total loss of nitrate reductase activity, these mutants do not, however, display FDH activity.     

Broad‐substrate screen as a tool to identify substrates for bacterial Gcn5‐related N‐acetyltransferases with unknown substrate specificity

Due to a combination of efforts from individual laboratories and structural genomics centers, there has been a surge in the number of members of the Gcn5‐related acetyltransferasesuperfamily that have been structurally determined within the past decade. Although the number of three‐dimensional structures is increasing steadily, we know little about the individual functions of these enzymes. Part of the difficulty in assigning functions for members of this superfamily is the lack of information regarding how substrates bind to the active site of the protein. The majority of the structures do not show ligand bound in the active site, and since the substrate‐binding domain is not strictly conserved, it is difficult to predict the function based on structure alone. Additionally, the enzymes are capable of acetylating a wide variety of metabolites and many may exhibit promiscuity regarding their ability to acetylate multiple classes of substrates, possibly having multiple functions for the same enzyme. Herein, we present an approach to identify potential substrates for previously uncharacterized members of the Gcn5‐related acetyltransferase superfamily using a variety of metabolites including polyamines, amino acids, antibiotics, peptides, vitamins, catecholamines, and other metabolites. We have identified potential substrates for eight bacterial enzymes of this superfamily. This information will be used to further structurally and functionally characterize them.     

Relationships between the osmoadaptation strategy, amino acid composition of bulk protein, and properties of certain enzymes of haloalkaliphilic bacteria

Haloalkaliphilic microorganisms isolated from soda lakes were compared in terms of the amino acid composition of bulk protein and the reaction of a number of key enzymes to salts and pH of the medium. In the extremely haloalkaliphilic bacterium Natroniella acetigena (selt-in osmoadaptation strategy), acidic amino acids (glutamic and aspartic) made up 30.91 mol % of the total of bulk protein amino acids. In the moderate haloalkaliphiles Tindallia magadiensis, Halomonas campisalis, and Halomonas sp. AIR-2 (compatible-solutes osmoadaptation strategy), the proportion of acidic amino acids (24.36, 23.15, and 23.58 mol %, respectively) was lower than in N. acetigena but higher than in the freshwater Acetobacterium paludosum (20.77 mol %). The excess of acidic amino acids over basic amino acids (lysine and arginine) increased with the degree of halophily. The enzymes of haloalkaliphiles proved to be tolerant to salts and high pH values, although the degree of tolerance varied. The activity of N. acetigena CO dehydrogenase was maximum in the presence of 0.7 M NaCl, but it was virtually independent of the NaHCO₃ concentration. The hydrogenase and CO dehydrogenase of T. magadiensis exhibited maximum activity in the absence of NaCl; the Co dehydrogenase was most active at 0.25 M NaHCO₃, and hydrogenase activity was only weakly dependent on NaHCO₃ in the concentration range of 0–1.2 M. The nitrate reductases of H. campisalis and Halomonas sp. AIR-2 were active in broad ranges of NaCl and KCl concentrations; the activity maxima were recorded at moderate concentrations of these salts. The pH optima of most of the studied enzymes of haloalkaliphiles were in the alkaline zone. Thus, it was shown that the amino acid composition of bulk protein is determined by the osmoadaptation strategy employed by the bacterium. A correlation was found between the salt tolerance of enzymes and the proportion of acidic amino acids in the bulk protein. The ability of enzymes to function at high pH values is one of the mechanisms of adaptation of microorganisms to high pH values.     

A comprehensive analysis of the geranylgeranylglyceryl phosphate synthase enzyme family identifies novel members and reveals mechanisms of substrate specificity and quaternary structure organization

Geranylgeranylglyceryl phosphate synthase (GGGPS) family enzymes catalyse the formation of an ether bond between glycerol‐1‐phosphate and polyprenyl diphosphates. They are essential for the biosynthesis of archaeal membrane lipids, but also occur in bacterial species, albeit with unknown physiological function. It has been known that there exist two phylogenetic groups (I and II) of GGGPS family enzymes, but a comprehensive study has been missing. We therefore visualized the variability within the family by applying a sequence similarity network, and biochemically characterized 17 representative GGGPS family enzymes regarding their catalytic activities and substrate specificities. Moreover, we present the first crystal structures of group II archaeal and bacterial enzymes. Our analysis revealed that the previously uncharacterized bacterial enzymes from group II have GGGPS activity like the archaeal enzymes and differ from the bacterial group I enzymes that are heptaprenylglyceryl phosphate synthases. The length of the isoprenoid substrate is determined in group II GGGPS enzymes by ‘limiter residues’ that are different from those in group I enzymes, as shown by site‐directed mutagenesis. Most of the group II enzymes form hexamers. We could disrupt these hexamers to stable and catalytically active dimers by mutating a single amino acid that acts as an ‘aromatic anchor’.     

Analysis of the reaction mechanism and substrate specificity of haloalkane dehalogenases by sequential and structural comparisons

Haloalkane dehalogenases catalyse environmentally important dehalogenation reactions. These microbial enzymes represent objects of interest for protein engineering studies, attempting to improve their catalytic efficiency or broaden their substrate specificity towards environmental pollutants. This paper presents the results of a comparative study of haloalkane dehalogenases originating from different organisms. Protein sequences and the models of tertiary structures of haloalkane dehalogenases were compared to investigate the protein fold, reaction mechanism and substrate specificity of these enzymes. Haloalkane dehalogenases contain the structural motifs of alpha/beta-hydrolases and epoxidases within their sequences. They contain a catalytic triad with two different topological arrangements. The presence of a structurally conserved oxyanion hole suggests the two-step reaction mechanism previously described for haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The differences in substrate specificity of haloalkane dehalogenases originating from different species might be related to the size and geometry of an active site and its entrance and the efficiency of the transition state and halide ion stabilization by active site residues. Structurally conserved motifs identified within the sequences can be used for the design of specific primers for the experimental screening of haloalkane dehalogenases. Those amino acids which were predicted to be functionally important represent possible targets for future site-directed mutagenesis experiments.     

Computational exploration of the network of sequence flow between protein structures

We investigate small sequence adjustments (of one or a few amino acids) that induce large conformational transitions between distinct and stable folds of proteins. Such transitions are intriguing from evolutionary and protein‐design perspectives. They make it possible to search for ancient protein structures or to design protein switches that flip between folds and functions. A network of sequence flow between protein folds is computed for representative structures of the Protein Data Bank. The computed network is dense, on an average each structure is connected to tens of other folds. Proteins that attract sequences from a higher than expected number of neighboring folds are more likely to be enzymes and alpha/beta fold. The large number of connections between folds may reflect the need of enzymes to adjust their structures for alternative substrates. The network of the Cro family is discussed, and we speculate that capacity is an important factor (but not the only one) that determines protein evolution. The experimentally observed flip from all alpha to alpha + beta fold is examined by the network tools. A kinetic model for the transition of sequences between the folds (with only protein stability in mind) is proposed. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Suppression of isotope scrambling in cell-free protein synthesis by broadband inhibition of PLP enymes for selective 15N-labelling and production of perdeuterated proteins in H2O

Selectively isotope labelled protein samples can be prepared in vivo or in vitro from selectively labelled amino acids but, in many cases, metabolic conversions between different amino acids result in isotope scrambling. The best results are obtained by cell-free protein synthesis, where metabolic enzymes are generally less active, but isotope scrambling can never be suppressed completely. We show that reduction of E. coli S30 extracts with NaBH₄ presents a simple and inexpensive way to achieve cleaner selective isotope labelling in cell-free protein synthesis reactions. The purpose of the NaBH₄ is to inactivate all pyridoxal-phosphate (PLP) dependent enzymes by irreversible reduction of the Schiff bases formed between PLP and lysine side chains of the enzymes or amino groups of free amino acids. The reduced S30 extracts retain their activity of protein synthesis, can be stored as well as conventional S30 extracts and effectively suppress conversions between different amino acids. In addition, inactivation of PLP-dependent enzymes greatly stabilizes hydrogens bound to α-carbons against exchange with water, minimizing the loss of α-deuterons during cell-free production of proteins from perdeuterated amino acids in H₂O solution. This allows the production of highly perdeuterated proteins that contain protons at all exchangeable positions, without having to back-exchange labile deuterons for protons as required for proteins that have been synthesized in D₂O.     

Opioid peptide biosynthesis: enzymatic selectivity and regulatory mechanisms

Certain general principles determine the biosynthesis of most biologically active peptides, including the opioid peptides, from large protein precursors. In almost all instances, the active peptide is embedded in the precursor flanked on both sides by pairs of basic amino acids. The first step in processing involves a trypsinlike enzyme, cleaving to the carboxyl terminus of basic amino acids, and leaving the active peptide with a basic amino acid on the carboxyl terminus. A carboxy-peptidase peptidase B-like enzyme then removes the remaining basic amino acid. It has been unclear whether any endopeptidases with trypsinlike activity are selective for one or another basic amino acid. Recently a soluble endopeptidase has been identified that can cleave to both the carboxyl and amino termini of basic amino acids. Enkephalin convertase (carboxypeptidase E, H) (EC 3.4.17.10) has considerable selectivity, and appears to be physiologically associated with the biosynthesis of enkephalin as well as a limited number of other neuropeptides. The turnover of opioid peptides and other neuropeptides is most effectively ascertained by measuring levels of mRNA either biochemically or by in situ hybridization. Striking dynamic alterations include a pronounced increase in levels of proenkephalin mRNA in the corpus striatum after blockade of dopamine receptors, but changes in opioid peptide mRNA after opiate addiction are less clear.     

Active site specificity of plasmepsin II.

Members of the aspartic proteinase family of enzymes have very similar three-dimensional structures and catalytic mechanisms. Each, however, has unique substrate specificity. These distinctions arise from variations in amino acid residues that line the active site subsites and interact with the side chains of the amino acids of the peptides that bind to the active site. To understand the unique binding preferences of plasmepsin II, an enzyme of the aspartic proteinase class from the malaria parasite, Plasmodium falciparum, chromogenic octapeptides having systematic substitutions at various positions in the sequence were analyzed. This enabled the design of new, improved substrates for this enzyme (Lys-Pro-Ile-Leu-Phe*Nph-Ala/Glu-Leu-Lys, where * indicates the cleavage point). Additionally, the crystal structure of plasmepsin II was analyzed to explain the binding characteristics. Specific amino acids (Met13, Ser77, and Ile287) that were suspected of contributing to active site binding and specificity were chosen for site-directed mutagenesis experiments. The Met13Glu and Ile287Glu single mutants and the Met13Glu/Ile287Glu double mutant gain the ability to cleave substrates containing Lys residues.     

The Universal Epitope of Influenza A Viral Neuraminidase Fundamentally Contributes to Enzyme Activity and Viral Replication

The only universally conserved sequence among all influenza A viral neuraminidases is located between amino acids 222 and 230. However, the potential roles of these amino acids remain largely unknown. Through an array of experimental approaches including mutagenesis, reverse genetics, and growth kinetics, we found that this sequence could markedly affect viral replication. Additional experiments revealed that enzymes with mutations in this region demonstrated substantially decreased catalytic activity, substrate binding, and thermostability. Consistent with viral replication analyses and enzymatic studies, protein modeling suggests that these amino acids could either directly bind to the substrate or contribute to the formation of the active site in the enzyme. Collectively, these findings reveal the essential role of this unique region in enzyme function and viral growth, which provides the basis for evaluating the validity of this sequence as a potential target for antiviral intervention and vaccine development.