Towards novel biocatalysts via protein design: the case of lipases

Abstract: During the past 3 years, the tertiary structures of several lipases have been solved by X-ray analysis. The structures revealed unique features such as hydrophobic 'patches' on the surface, presumably involved in lipid supersubstrate binding, and a lid structure which covers the active site in the absence of substrate. Only very recently the first X-ray structure of a bacterial lipase has been solved, and further structural features different from lipases of eukaryotic origin became apparent. Many lipase genes have been cloned and sequenced recently, and expression systems for the preparation of recombinant enzymes in good yields are available. As an example, the lipase from Rhizopus oryzae has been successfully expressed by us in Escherichia coli , and the resulting inclusion bodies were renatured in high yields. Consequently, the mechanism of action of lipases is now being studied via site-directed mutagenesis, and the rational design of lipases for the selective transformation of substrates is presently addressed in several laboratories.     

Conformational changes and protein stability of the pro-apoptotic protein Bax

Pro-apoptotic Bax is a soluble and monomeric protein under normal physiological conditions. Upon its activation substantial structural rearrangements occur: The protein inserts into the mitochondrial outer membrane and forms higher molecular weight oligomers. Subsequently, the cells can undergo apoptosis. In our studies, we focused on the structural rearrangements of Bax during oligomerization and on the protein stability. Both protein conformations exhibit high stability against thermal denaturation, chemically induced unfolding and proteolytic processing. The oligomeric protein is stable up to 90 °C as well as in solutions of 8 M urea or 6 M guanidinium hydrochloride. Helix 9 appears accessible in the monomer but hidden in the oligomer assessed by proteolysis. Tryptophan fluorescence indicates that the environment of the C-terminal protein half becomes more apolar upon oligomerization, whereas the loop region between helices 1 and 2 gets solvent exposed.     

bPE toolkit: toolkit for computational protein engineering

We present a computational toolkit consisting of five utility tools, for performing basic operations on a protein structure file in PDB format. The toolkit consists of five different programs which can be integrated as part of a pipeline for computational protein structure characterization or as a standalone analysis package. The programs include tools for chirality check for amino acids (ProChiral), contact map generation (CoMa), data redundancy (DaRe), hydrogen bond potential energy (HyPE) and electrostatic interaction energy (EsInE). All programs in the toolkit can be accessed and downloaded through the following link: http://www.iitg.ac.in/bpetoolkit/.     

Structure determination of archaea-specific ribosomal protein L46a reveals a novel protein fold

Three archaea-specific ribosomal proteins recently identified show no sequence homology with other known proteins. Here we determined the structure of L46a, the most conserved one among the three proteins, from Sulfolobus solfataricus P2 using NMR spectroscopy. The structure presents a twisted β-sheet formed by the N-terminal part and two helices at the C-terminus. The L46a structure has a positively charged surface which is conserved in the L46a protein family and is the potential rRNA-binding site. Searching homologous structures in Protein Data Bank revealed that the structure of L46a represents a novel protein fold. The backbone dynamics identified by NMR relaxation experiments reveal significant flexibility at the rRNA binding surface. The potential position of L46a on the ribosome was proposed by fitting the structure into a previous electron microscopy map of the ribosomal 50S subunit, which indicated that L46a contacts to domain I of 23S rRNA near a multifunctional ribosomal protein L7ae.     

Development of a novel Gateway-based vector system for efficient, multiparallel protein expression in Escherichia coli

We describe a cloning and expression system which is based on the Escherichia coli T7 expression system and Gateway recombination technology. We have produced numerous destination vectors with selected fusion tags and an additional set of entry vectors containing the gene of interest and optional labeling tags. This powerful system enables us to transfer a cDNA to several expression vectors in parallel and combine them with various labeling tags. To remove the attached amino terminal tags along with the unwanted attB1 site, we inserted PreScission protease cleavage sites. In contrast to the commercially available destination vectors, our plasmids provide kanamycin resistance, which can be an advantage when expressing toxic proteins in E. coli. Some small-scale protein expression experiments are shown to demonstrate the usefulness of these novel Gateway vectors. In summary, this system has some benefits over the widely used and commercially available Gateway standard system, and it enables many different combinations for expression constructs from a single gene of interest.     

Modular peptide binding: From a comparison of natural binders to designed armadillo repeat proteins

Several binding scaffolds that are not based on immunoglobulins have been designed as alternatives to traditional monoclonal antibodies. Many of them have been developed to bind to folded proteins, yet cellular networks for signaling and protein trafficking often depend on binding to unfolded regions of proteins. This type of binding can thus be well described as a peptide–protein interaction. In this review, we compare different peptide-binding scaffolds, highlighting that armadillo repeat proteins (ArmRP) offer an attractive modular system, as they bind a stretch of extended peptide in a repeat-wise manner. Instead of generating each new binding molecule by an independent selection, preselected repeats – each complementary to a piece of the target peptide – could be designed and assembled on demand into a new protein, which then binds the prescribed complete peptide. Stacked armadillo repeats (ArmR), each typically consisting of 42 amino acids arranged in three α-helices, build an elongated superhelical structure which enables binding of peptides in extended conformation. A consensus-based design approach, complemented with molecular dynamics simulations and rational engineering, resulted in well-expressed monomeric proteins with high stability. Peptide binders were selected and several structures were determined, forming the basis for the future development of modular peptide-binding scaffolds.     

Artificial leucine rich repeats as new scaffolds for protein design

The leucine rich repeat (LRR) motif that participates in many biomolecular recognition events in cells was suggested as a general scaffold for producing artificial receptors. We describe here the design and first total chemical synthesis of small LRR proteins, and their structural analysis. When evaluating the tertiary structure as a function of different number of repeating units (1-3), we were able to find that the 3-repeats sequence, containing 90 amino acids, folds into the expected structure.     

The role of protein conformational fluctuations in allostery, function, and evolution

It is now well-known that proteins exist at equilibrium as ensembles of conformational states rather than as unique static structures. Here we review from an ensemble perspective important biological effects of such spontaneous fluctuations on protein allostery, function, and evolution. However, rather than present a thorough literature review on each subject, we focus instead on connecting these phenomena through the ensemble-based experimental, theoretical, and computational investigations from our laboratory over the past decade. Special emphasis is given to insights that run counter to some of the prevailing ideas that have emerged over the past 40 years of structural biology research. For instance, when proteins are viewed as conformational ensembles rather than as single structures, the commonly held notion of an allosteric pathway as an obligate series of individual structural distortions loses its meaning. Instead, allostery can result from energetic linkage between distal sites as one Boltzmann distribution of states transitions to another. Additionally, the emerging principles from this ensemble view of proteins have proven surprisingly useful in describing the role of intrinsic disorder in inter-domain communication, functional adaptation mediated by mutational control of fluctuations, and evolutionary conservation of the energetics of protein stability.     

C-terminal modification of 6-phosphofructo-1-kinase from Saccharomyces cerevisiae and its influence on enzyme structure and activity

Studies on limited proteolysis of 6-phosphofructo-1-kinase (Pfk-1) from Saccharomyces cerevisiae led to the suggestion that the C-terminal part of the alpha-subunit must contribute to the stabilisation of the octameric enzyme structure. To analyse the role of the C-terminus in vivo, the respective terminus of one of both types of subunits of Pfk-1 was sequentially truncated or extended. These modifications resulted in a decrease of the protein level of the mutated subunit and of the specific enzyme activity in the cell-free extract as well as in changes of the kinetic properties. Size exclusion HPLC demonstrated that the modified subunit is still able to assemble with the native counterpart generating an enzymatically active hetero-octamer. On the basis of our results we assume that the C-termini are important for the three-dimensional structure of the subunits determining their susceptibility to proteolysis and the ability to assembly to an active, oligomeric Pfk-1.     

Isolation of nuclear encoded plastid ribosomal protein cDNAs

Summary A pea leaf cDNA library was constructed in the expression vector gt11 and screened with antisera raised against proteins extracted from 30S and 50S ribosomal subunits and 70S ribosomes prepared from isolated pea chloroplasts. Six recombinant phage were identified that encoded fusion proteins containing plastid ribosomal protein antigenic determinants. Phage-induced cell lysate proteins, containing the fusion proteins, were bound to nitrocellulose membranes and used as affinity matrices to prepare monospecific antibodies. These antibodies were then used to identify by Western blotting which plastid ribosomal protein shared antigenic determinants with the fusion proteins. cDNA inserts from the antigen-producing phage were used to hybrid-select complementary mRNAs. The cell-free translation products of these mRNAs were added to a pea chloroplast in vitro transport system and imported proteins analyzed by two-dimensional gel electrophoresis. The imported proteins comigrated with the plastid ribosomal proteins that were identified as being antigenically related to the fusion proteins produced by the corresponding recombinant phage. The imported proteins were 3,500–5,500 daltons smaller than their precursors.     

Unchanged antigenic properties of a nuclear non-histone protein during phylogenesis

Summary The phyletic distribution of Sm antigen in eukaryotes has been investigated by means of immunochemical techniques. This nuclear antigen, which is an acidic glycoprotein of about 7S with DNA binding ability, is present in all organisms examined. Such evidence suggests a very slow evolutionary rate for at least the antigenic determinant(s) of the Sm antigen.     

Amino acid propensities for secondary structures are influenced by the protein structural class

Amino acid propensities for secondary structures were used since the 1970s, when Chou and Fasman evaluated them within datasets of few tens of proteins and developed a method to predict secondary structure of proteins, still in use despite prediction methods having evolved to very different approaches and higher reliability. Propensity for secondary structures represents an intrinsic property of amino acid, and it is used for generating new algorithms and prediction methods, therefore our work has been aimed to investigate what is the best protein dataset to evaluate the amino acid propensities, either larger but not homogeneous or smaller but homogeneous sets, i.e., all-alpha, all-beta, alpha-beta proteins. As a first analysis, we evaluated amino acid propensities for helix, beta-strand, and coil in more than 2000 proteins from the PDBselect dataset. With these propensities, secondary structure predictions performed with a method very similar to that of Chou and Fasman gave us results better than the original one, based on propensities derived from the few tens of X-ray protein structures available in the 1970s. In a refined analysis, we subdivided the PDBselect dataset of proteins in three secondary structural classes, i.e., all-alpha, all-beta, and alpha-beta proteins. For each class, the amino acid propensities for helix, beta-strand, and coil have been calculated and used to predict secondary structure elements for proteins belonging to the same class by using resubstitution and jackknife tests. This second round of predictions further improved the results of the first round. Therefore, amino acid propensities for secondary structures became more reliable depending on the degree of homogeneity of the protein dataset used to evaluate them. Indeed, our results indicate also that all algorithms using propensities for secondary structure can be still improved to obtain better predictive results.     

Artificial di-iron proteins: solution characterization of four helix bundles containing two distinct types of inter-helical loops

Peptide-based models have an enormous impact for the development of metalloprotein models, as they seem appropriate candidates to mimic both the structural characteristics and reactivity of the natural systems. Through the de novo design of four-helix bundles, we developed the DF (Due Ferri) family of artificial proteins, as models of di-iron and di-manganese metalloproteins. The goal of our research is to elucidate how the electrostatic environment, polarity and solvent accessibility of the metal-binding site, influence the functional properties of di-iron proteins. The first two subsets of the DF protein family, DF1 and DF2, consist of two non-covalently associated helix-loop-helix motifs, which bind the di-metal cofactor near the center of the structure. The DF2 subset was designed to improve the properties of DF1: DF2 and DF2t have several changes in their sequences to improve solubility and metal ion access, as well as a change in the loop connecting the two helices. In order to evaluate how these changes affect the overall structure of the model proteins, we solved the NMR structures of the di-Zn(II) complexes of DF2 and DF2t, and compared these structures with those recently obtained from X-ray crystallography. Further, we examined the thermodynamic consequences associated with the mutations, by measuring the stability of DF2t in the presence of different metal ions, and comparing the results with the data already obtained for DF2. Taken together, analysis of all the data showed the importance of the turn conformation in the design and stability of four-helix bundle.Electronic Supplementary Material Supplementary material is available for this article at     

Organization of functional domains in the docking protein p130Cas

The docking protein p130Cas becomes phosphorylated upon cell adhesion to extracellular matrix proteins, and is thought to play an essential role in cell transformation. Cas transmits signals through interactions with the Src-homology 3 (SH3) and Src-homology 2 domains of FAK or v-Crk signaling molecules, or with 14-3-3 protein, as well as phosphatases PTP1B and PTP-PEST. The large (130kDa), multi-domain Cas molecule contains an SH3 domain, a Src-binding domain, a serine-rich protein interaction region, and a C-terminal region that participates in protein interactions implicated in antiestrogen resistance in breast cancer. In this study, as part of a long-term goal to examine the protein interactions of Cas by X-ray crystallography and nuclear magnetic resonance spectroscopy, molecular constructs were designed to express two adjacent domains, the serine-rich domain and the Src-binding domain, that each participate in intermolecular contacts dependent on protein phosphorylation. The protein products are soluble, homogeneous, monodisperse, and highly suitable for structural studies to define the role of Cas in integrin-mediated cell signaling.     

The pH-dependent unfolding mechanism of P2 myelin protein: an experimental and computational study

The P2 protein is a small, extrinsic protein of the myelin membrane in the peripheral nervous system that structurally belongs to the fatty acid binding proteins (FABPs) family, sharing with them a 10 strands beta-barrel structure. FABPs appear to be involved in cellular fatty acid transport, but very little is known about the role of P2 in the metabolism of peripheral myelin lipids. Study of protein conformation at different pHs is a useful tool for the characterization of the unfolding mechanisms and the intrinsic conformational properties of the protein, and may give insight into factors that guide protein folding pathways. In particular, low pH conditions have been shown to induce partially folded states in several proteins. In this paper, the acidic unfolding of purified P2 protein was studied with both spectroscopic techniques and molecular dynamics simulation. Both experimental and computational results indicate the presence of a partly folded state at low pH, which shows structural changes mainly involving the lid that is formed by the helix-turn-helix domain. The opening of the lid, together with a barrel relaxation, could regulate the ligand exchanges near the cell membrane, supporting the hypothesis that the P2 protein may transport fatty acids between Schwann cells and peripheral myelin.     

Hookworm SCP/TAPS protein structure--A key to understanding host-parasite interactions and developing new interventions

SCP/TAPS proteins are a diverse family of molecules in eukaryotes, including parasites. Despite their abundant occurrence in parasite secretomes, very little is known about their functions in parasitic nematodes, including blood-feeding hookworms. Current information indicates that SCP/TAPS proteins (called Ancylostoma-secreted proteins, ASPs) of the canine hookworm, Ancylostoma caninum, represent at least three distinct groups of proteins. This information, combined with comparative modelling, indicates that all known ASPs have an equatorial groove that binds extended structures, such as peptides or glycans. To elucidate structure-function relationships, we explored the three-dimensional crystal structure of an ASP (called Ac-ASP-7), which is highly up-regulated in expression in the transition of A. caninum larvae from a free-living to a parasitic stage. The topology of the N-terminal domain is consistent with pathogenesis-related proteins, and the C-terminal extension that resembles the fold of the Hinge domain. By anomalous diffraction, we identified a new metal binding site in the C-terminal extension of the protein. Ac-ASP-7 is in a monomer-dimer equilibrium, and crystal-packing analysis identified a dimeric structure which might resemble the homo-dimer in solution. The dimer interaction interface includes a novel binding site for divalent metal ions, and is proposed to serve as a binding site for proteins involved in the parasite-host interplay at the molecular level. Understanding this interplay and the integration of structural and functional data could lead to the design of new approaches for the control of parasitic diseases, with biotechnological outcomes.     

S100 and S100 fused-type protein families in epidermal maturation with special focus on S100A3 in mammalian hair cuticles

Epithelial Ca²⁺-regulation, which governs cornified envelope formation in the skin epidermis and hair follicles, closely coincides with the expression of S100A3, filaggrin and trichohyalin, and the post-translational modification of these proteins by Ca²⁺-dependent peptidylarginine deiminases. This review summarizes the current nomenclature and evolutional aspects of S100 Ca²⁺-binding proteins and S100 fused-type proteins (SFTPs) classified as a separate protein family with special reference to the molecular structure and function of S100A3 dominantly expressed in hair cuticular cells. Both S100 and SFTP family members are identified by two distinct types of Ca²⁺-binding loops in an N-terminal pseudo EF-hand motif followed by a canonical EF-hand motif. Seventeen members of the S100 protein family including S100A3 are clustered with seven related genes encoding SFTPs on human chromosome 1q21, implicating their association with epidermal maturation and diseases. Human S100A3 is characterized by two disulphide bridges and a preformed Zn²⁺-pocket, and may transfer Ca²⁺ ions to peptidylarginine deiminases after its citrullination-mediated tetramerization. Phylogenetic analysis utilizing current genome databases suggests that divergence of the S100A3 gene coincided with the emergence of hair, a defining feature of mammals, and that the involvement of S100A3 in epithelial Ca²⁺-cycling occurred as a result of a skin adaptation in terrestrial mammals.     

Membrane protein prediction methods

We survey computational approaches that tackle membrane protein structure and function prediction. While describing the main ideas that have led to the development of the most relevant and novel methods, we also discuss pitfalls, provide practical hints and highlight the challenges that remain. The methods covered include: sequence alignment, motif search, functional residue identification, transmembrane segment and protein topology predictions, homology and ab initio modeling. In general, predictions of functional and structural features of membrane proteins are improving, although progress is hampered by the limited amount of high-resolution experimental information available. While predictions of transmembrane segments and protein topology rank among the most accurate methods in computational biology, more attention and effort will be required in the future to ameliorate database search, homology and ab initio modeling.