Prediction of signal peptides in archaea

Computational prediction of signal peptides (SPs) and theircleavage sites is of great importance in computational biology;however, currently there is no available method capable of predictingreliably the SPs of archaea, due to the limited amount of experimentallyverified proteins with SPs. We performed an extensive literaturesearch in order to identify archaeal proteins having experimentallyverified SP and managed to find 69 such proteins, the largestnumber ever reported. A detailed analysis of these sequencesrevealed some unique features of the SPs of archaea, such asthe unique amino acid composition of the hydrophobic regionwith a higher than expected occurrence of isoleucine, and acleavage site resembling more the sequences of gram-positiveswith almost equal amounts of alanine and valine at the position-3before the cleavage site and a dominant alanine at position-1,followed in abundance by serine and glycine. Using these proteinsas a training set, we trained a hidden Markov model method thatpredicts the presence of the SPs and their cleavage sites andalso discriminates such proteins from cytoplasmic and transmembraneones. The method performs satisfactorily, yielding a 35-foldcross-validation procedure, a sensitivity of 100% and specificity98.41% with the Matthews’ correlation coefficient beingequal to 0.964. This particular method is currently the onlyavailable method for the prediction of secretory SPs in archaea,and performs consistently and significantly better comparedwith other available predictors that were trained on sequencesof eukaryotic or bacterial origin. Searching 48 completely sequencedarchaeal genomes we identified 9437 putative SPs. The method,PRED-SIGNAL, and the results are freely available for academicusers at http://bioinformatics.biol.uoa.gr/PRED-SIGNAL/ andwe anticipate that it will be a valuable tool for the computationalanalysis of archaeal genomes.     

Temperature-sensitive production of rabbit muscle glycogen phosphorylase in Escherichia coli

In order to understand how allosteric switches regulate boththe catalytic activity and molecular interactions of glycogenphosphorylase, it is necessary to design and analyze variantproteins that test hypotheses about the structural details ofthe allosteric mechanism. Essential to such an investigationis the ability to obtain large amounts of variant proteins.We developed a system for obtaining milligram amounts (>20 mg/1) of rabbit muscle phosphorylase from bacteria. Phosphorylaseaggregates as inactive protein when a strong bacterial promoteris used under full inducing conditions and normal growth conditions.However, when the growth temperature of bacteria expressingphosphorylase is reduced to 22°C we obtain active musclephosphorylase. The degree to which the induced expression ofphosphorylase protein is temperature sensitive depends on thestrain of bacteria used. New assay and purification methodswere developed to allow rapid purification of engineered phosphorylaseproteins from bacterial cultures. The rabbit muscle phosphorylaseobtained from the bacterial expression system is enzymaticallyidentical to the enzyme purified from rabbit muscle. The expressedprotein crystallizes in the same conditions used for growingcrystals of protein from rabbit muscle and the crystal formis isomorphous. Rabbit muscle phosphorylase is one of the largestoligomeric mammalian enzymes successfully expressed in Escherichiacoli. Our results indicate that optimization of a combinationof growth and induction conditions will be important in theexpression of other heterologous proteins in bacteria.     

An evolution-based analysis scheme to identify CO2/O2 specificity-determining factors for ribulose 1,5-bisphosphate carboxylase/oxygenase

Ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCo) catalyzes a rate-limiting step in photosynthetic carbon assimilation (reacting with CO2) and its competitive photo-respiratory carbon oxidation (reacting with O2). RuBisCo enzyme with an enhanced CO2/O2 specificity would boost the ability to make great progress in agricultural production and environmental management. RuBisCos in marine non-green algae, resulting from an earlier endo-symbiotic event, diverge greatly from those in green plants and cyanobacteria and, further, have the highest CO2/O2 specificity whereas RuBisCos in cyanobacteria have the lowest. We assumed that there exist different levels of CO2/O2 specificity-determining factors, corresponding to different evolutionary events and specificity levels. Based on this assumption, we devised a scheme to identify these substrate-determining factors. From this analysis, we are able to discover different categories of the CO2/O2 specificity-determining factors that show which residue substitutions account for (relatively) small specificity changes, as happened in green plants, or a tremendous enhancement, as observed in marine non-green algae. Therefore, the analysis can improve our understanding of molecular mechanisms in the substrate specificity development and prioritize candidate specificity-determining surface residues for site-directed mutagenesis.     

Comparison of conservation within and between the Ser/Thr and Tyr protein kinase family: proposed model for the catalytic domain of the epidermal growth factor receptor

The protein kinase family can be subdivided into two main groupsbased on their ability to phosphorylate Ser/Thr or Tyr substrates.In order to understand the basis of this functional difference,we have carried out a comparative analysis of sequence conservationwithin and between the Ser/Thr and Tyr protein kinases. A multiplesequence alignment of 86 protein kinase sequences was generated.For each position in the alignment we have computed the conservationof residue type in the Ser/Thr, in the Tyr and in both of thekinase subfamilies. To understand the structural and/or functionalbasis for the conservation, we have mapped these conservationproperties onto the backbone of the recently determined structureof the cAMP–dependent Ser/Thr kinase. The results showthat the kinase structure can be roughly segregated, based uponconservation, into three zones. The inner zone contains residueshighly conserved in all the kinase family and describes thehydrophobic core of the enzyme together with residues essentialfor substrate and ATP binding and catalysis. The outer zonecontains residues highly variable in all kinases and representsthe solvent–exposed surface of the protein. The thirdzone is comprised of residues conserved in either the Ser/Thror Tyr kinases or in both, but which are not conserved betweenthem. These are sandwiched between the hydrophobic core andthe solvent-exposed surface. In addition to analyzing overallconservation hi the kinase family, we have also looked at conservationof its substrate and ATP binding sites. The ATP site is highlyconserved throughout the kinases, whereas the substrate bindingsite is more variable. The active site contains several positionswhich differ between the Ser/Thr and Tyr kinases and may beresponsible for discriminating between hydroxyl bearing sidechains. Using this information we propose a model for Tyr substratebinding to the catalytic domain of the epidermal growth factorreceptor (EGFR).     

Distinctive properties of signal sequences from bacterial lipoproteins

We have compared a number of attributes (hydrophobicity, aminoacid size, charge and secondary structure propensities) of signalsequences from bacterial lipoproteins with the same attributesof signal peptides from other prokaryotic proteins (non-lipoproteins).Lipoprotein leader sequences tend to be shorter, more hydrophobicand bulky, and they have stronger conformational preferences,the most conspicuous being a predicted ß-turn comprisingpositions 2 or 3 of the mature protein. Another distinctivefeature is a maximum in the local energy profile between positions–1 and +2. With one exception (ß-lactamase III),the lipoproteins do not have Pro in their signal peptides, andthey tend to have fewer Ser and Thr but more Gly than non-lipoproteins.Lipoproteins also lack a net negative charge in the N-terminalregions of the mature proteins. The signal peptides of the bacteriocinplasmid-coded lysis proteins appear to be unique in that theyhave all the ascribed features of lipoprotein signals; thesecharacteristics can be used to guide signal peptide mutagenesisexperiments and to construct new secretion vehicles.     

Citrate synthase from the hyperthermophilic Archaeon, Pyrococcus furiosus

The gene encoding the enzyme citrate synthase has been clonedand sequenced from the hyperthermophilic Archaeon Pyrococcusfuriosus, and the derived amino acid sequence has been phylogeneticallycompared with citrate synthases from archaeal, bacterial andeukaryal organisms.The gene has been over-expressed in Escherichiacoli to produce an active enzyme that has then been characterizedwith respect to its kinetic, oligomeric and hyperthermost-ableproperties. A structurally-based sequence alignment was madeto the citrate synthase from the thermophilic Archaeon Thermoplasmaacidophilum, the crystal structure of which we have determinedrecently. From this alignment,a homology-modelled structurefor the P.furiosus citratesynthase was generated and analysed.     

Enhancing the stability of microsomal cytochrome b5: a rational approach informed by comparative studies with the outer mitochondrial membrane isoform

The outer mitochondrial membrane isoform of mammalian cytochrome b5 (OM b5) is much less prone to lose heme than the microsomal isoform (Mc b5), with a conserved difference at position 71 (leucine versus serine) playing a major role. We replaced Ser71 in Mc b5 with Leu, with the prediction that it would retard heme loss by diminishing polypeptide expansion accompanying rupture of the histidine to iron bonds. The strategy was partially successful in that it slowed dissociation of heme from its less stable orientation in bMc b5 (B). Heme dissociation from orientation A was accelerated to a similar extent, however, apparently owing to increased binding pocket dynamic mobility related to steric strain. A second mutation (L32I) guided by results of previous comparative studies of Mc and OM b5s diminished the steric strain, but much greater relief was achieved by replacing heme with iron deuteroporphyrin IX (FeDPIX). Indeed, the stability of the Mc(S71L) b5 FeDPIX complex is similar to that of the FeDPIX complex of OM b5. The results suggest that maximizing heme binding pocket compactness in the apo state is a useful general strategy for increasing the stability of engineered or designed proteins.     

Comparative modelling and analysis of amino acid substitutions suggests that the family of pregnancy-associated glycoproteins includes both active and inactive aspartic proteinases

The pregnancy-associated glycoproteins (PAGs) are secretoryproducts synthesized by the outer epithelial cell layer (chorion)of the placentas of various ungulate species. The amino acidsequences of eight PAGs have been inferred from cloned cDNAof cattle and sheep, as well as of the non-ruminant pig andhorse. We compare the PAG sequences and present results of thethree-dimensional models of boPAG-1 and ovPAG-1 that were constructedon the basis of the crystal structures of homologous porcinepepsin and bovine chymosin using a rule-based comparative modellingapproach. Further, we compare peptide binding subsites definedby interactions with pepstatin and a decapeptide inhibitor (CH-66)modelled on the basis of crystal structures of other asparticproteinases. We have extended our analysis of the peptide bindingsubsites to the other PAG molecules of known sequence by aligningthe PAG sequences to the structural template derived from thepepsin family and by making use of the three-dimensional modelsof the boPAG-1 and ovPAG-1. The residues that are likely toaffect peptide binding in the boPAG-1, ovPAG-1 and other PAGmolecules have been identified. Sequence comparisons revealthat all PAG molecules may have evolved from a pepsin-like progenitormolecule with the equine PAG most closely related to the pepsins.The presence of substitutions at the S₁ and other subsites relativeto pepsin make it unlikely that either bovine, ovine or theporcine PAG-1 have catalytic activity. Only two of the eightPAGs examined (porcine PAG-2 and equine PAG-1) retain featuresof active aspartic proteinases with pepsin-like activity. Ourresults indicate that in the PAGs so far characterized the peptidebinding specificities differ significantly from each other andfrom pepsin, despite their high sequence identities. Analysisof the various peptide binding subsites demonstrates why bothbovine and ovine PAG-1 are capable of binding pepstatin. Thestrong negative charge in the binding cleft of boPAG-1 and ovPAG-1indicates a preference for lysine- or arginine-rich peptides.PAGs represent a family where the possible peptide binding functionmay be retained through their binding specificities, but wherethe catalytic activity may be lost in some cases, such as theboPAG-1, ovPAG-1 and the poPAG-1.     

Prediction and analysis of SH2 domain-phosphopeptide interactions

Src homology 2 (SH2) domains are small protein modules of -100amino acids that are found in many proteins involved in intracellularsignal transduction. They mediate protein-protein interactionsand modulate enzyme activity by their ability to bind to specificsequence patterns that contain a phosphorylated tyrosine. Asthe three-dimensional structures of the phosphatidylinositol(PI) 3-kinase, Lck, Src and Abl SH2 domains have been shownto be similar, we have modelled other SH2 domains that showdistinct sequence specificity to allow comparative analysisof SH2-phosphopeptide interactions. The SH2 domains of PLC-Nterm.,Nck, Grb2, GAP and Abl have been model-built with high-affinityphosphopeptides fitted into the putative binding sites. Foreach SH2 domain a detailed analysis of the peptide-protein interactionwas performed. It is apparent that specificity is mainly conferredby three to five residues downstream from the phosphotyrosineresidue (Y*), especially, although not exclusively, peptideposition Y* + 3. The SH2 pocket that binds the Y* + 3 residueis mainly composed of three sections: part of strand (ßEgoing into loop EF, part of B and loop BG. The residues thatconstitute the Y* +3 binding pocket show variability that seemsto determine which amino acid binds preferentially. Residueposition ßE4 seems to play a vital role in the SH2specificity. This study shows that the development of modellingprotocols for SH2 domains whose structure has not been determinedcan prove very useful in predicting which residues are involvedin conferring the affinity and binding specificity of thesedomains towards distinct phosphotyrosine-containing sequences.     

Bacterial lipid modification of proteins for novel protein engineering applications

Functioning of proteins efficiently at the solid-liquid interface is critical to not only biological but also modern man-made systems such as ELISA, liposomes and biosensors. Anchoring hydrophilic proteins poses a major challenge in this regard. Lipid modification, N-acyl-S-diacylglyceryl-Cys, providing an N-terminal hydrophobic membrane anchor is a viable solution that bacteria have successfully evolved but remains unexploited. Based on the current understanding of this ubiquitous and unique bacterial lipid modification it is possible to use Escherichia coli, the popular recombinant protein expression host, for converting a non-lipoprotein to a lipoprotein with a hydrophobic anchor at the N-terminal end. We report two strategies applicable to non-lipoproteins (with or without signal sequences) employing minimal sequence change. Taking periplasmic Shigella apyrase as an example, its signal sequence was engineered to include a lipobox, an essential determinant for lipid modification, or its mature sequence was fused to the signal sequence of abundant outer membrane lipoprotein, Lpp. Lipid modification was proved by membrane localization, electrophoretic mobility shift and mass spectrometric analysis. Substrate specificity and specific activity measurements indicated functional integrity after modification. In conclusion, a convenient protein engineering strategy for converting non-lipoprotein to lipoprotein for commercial application has been devised and tested successfully.     

Improving solubility and refolding efficiency of human V(H)s by a novel mutational approach

The antibody V(H) domains of camelids tend to be soluble and to resist aggregation, in contrast to human V(H) domains. For immunotherapy, attempts have therefore been made to improve the properties of human V(H)s by camelization of a small set of framework residues. Here, we have identified through sequence comparison of well-folded llama V(H) domains an alternative set of residues (not typically camelid) for mutation. Thus, the solubility and thermal refolding efficiency of a typical human V(H), derived from the human antibody BT32/A6, were improved by introduction of two mutations in framework region (FR) 1 and 4 to generate BT32/A6.L1. Three more mutations in FR3 of BT32/A6.L1 further improved the thermal refolding efficiency while retaining solubility and cooperative melting profiles. To demonstrate practical utility, BT32/A6.L1 was used to construct a phage display library from which were isolated human V(H)s with good antigen binding activity and solubility. The engineered human V(H) domains described here may be useful for immunotherapy, due to their expected low immunogenicity, and in applications involving transient high temperatures, due to their efficient refolding after thermal denaturation.     

Structure-guided SCHEMA recombination of distantly related beta-lactamases

We constructed a library of beta-lactamases by recombining three naturally occurring homologs (TEM-1, PSE-4, SED-1) that share 34-42% sequence identity. Most chimeras created by recombining such distantly related proteins are unfolded due to unfavorable side-chain interactions that destabilize the folded structure. To enhance the fraction of properly folded chimeras, we designed the library using SCHEMA, a structure-guided approach to choosing the least disruptive crossover locations. Recombination at seven selected crossover positions generated 6561 chimeric sequences that differ from their closest parent at an average of 66 positions. Of 553 unique characterized chimeras, 111 (20%) retained beta-lactamase activity; the library contains hundreds more novel beta-lactamases. The functional chimeras share as little as 70% sequence identity with any known sequence and are characterized by low SCHEMA disruption (E) compared to the average nonfunctional chimera. Furthermore, many nonfunctional chimeras with low E are readily rescued by low error-rate random mutagenesis or by the introduction of a known stabilizing mutation (TEM-1 M182T). These results show that structure-guided recombination effectively generates a family of diverse, folded proteins even when the parents exhibit only 34% sequence identity. Furthermore, the fraction of sequences that encode folded and functional proteins can be enhanced by utilizing previously stabilized parental sequences.     

Evidence for the involvement of tyrosine-69 in the control of stereospecificity of porcine pancreatic phospholipase A2

We have studied the role of Tyr-69 of porcine pancreatic phospholipaseA₂ in catalysis and substrate binding, using site-directed mutagenesis.A mutant was constructed containing Phe at position 69. Kineticcharacterization revealed that the Phe-69 mutant has retainedenzymatic activity on monomeric and micellar substrates, andthat the mutation has only minor effects on k_cat and K_m. Thisshows that Tyr-69 plays no role in the true catalytic eventsduring substrate hydrolysis. In contrast, the mutation has aprofound influence on the stereospecificity of the enzyme. Whereasthe wild-type phospholipase A₂ is only able to catalyse thedegradation of sn-3 phospholipids, the Phe-69 mutant hydrolysesboth the sn-3 isomers and, at a low (1–2%) rate, the sn-1isomers. Despite the fact that the stereospecificity of themutant phospholipase has been altered, Phe-69 phospholipasestill requires Ca²⁺ ions as a cofactor and also retains itsspecificity for the sn-2 ester bond. Our data suggest that inporcine pancreatic phospholipase A₂ the hydroxyl group of Tyr-69serves to fix and orient the phosphate group of phospholipidmonomers by hydrogen bonding. Because no such interaction canoccur between the Phe-69 side-chain and the phosphate moietyof the substrate monomer, the mutant enzyme loses part of itsstereospecificity but not its positional specificity.     

Pattern-induced multi-sequence alignment (PUMA) algorithm employing secondary structure-dependent gap penalties for use in comparative protein modelling

A multiple sequence alignment algorithm is described that usesa dynamic programming-based pattern construction method to aligna set of homologous sequences based on their common patternof conserved sequence elements. This pattern-induced multi-sequencealignment (PUMA) algorithm can employ secondary-structure dependentgap penalties for use in comparative modelling of new sequenceswhen the three-dimensional structure of one or more membersof the same family is known. We show that the use of secondarystructure information can significantly improve the accuracyof aligning structure boundaries in a set of homologous sequenceseven when the structure of only one member of the family isknown     

A theoretical model of restriction endonuclease NlaIV in complex with DNA, predicted by fold recognition and validated by site-directed mutagenesis and circular dichroism spectroscopy

Restriction enzymes (REases) are commercial reagents commonly used in DNA manipulations and mapping. They are regarded as very attractive models for studying protein-DNA interactions and valuable targets for protein engineering. Their amino acid sequences usually show no similarities to other proteins, with rare exceptions of other REases that recognize identical or very similar sequences. Hence, they are extremely hard targets for structure prediction and modeling. NlaIV is a Type II REase, which recognizes the interrupted palindromic sequence GGNNCC (where N indicates any base) and cleaves it in the middle, leaving blunt ends. NlaIV shows no sequence similarity to other proteins and virtually nothing is known about its sequence-structure-function relationships. Using protein fold recognition, we identified a remote relationship between NlaIV and EcoRV, an extensively studied REase, which recognizes the GATATC sequence and whose crystal structure has been determined. Using the 'FRankenstein's monster' approach we constructed a comparative model of NlaIV based on the EcoRV template and used it to predict the catalytic and DNA-binding residues. The model was validated by site-directed mutagenesis and analysis of the activity of the mutants in vivo and in vitro as well as structural characterization of the wild-type enzyme and two mutants by circular dichroism spectroscopy. The structural model of the NlaIV-DNA complex suggests regions of the protein sequence that may interact with the 'non-specific' bases of the target and thus it provides insight into the evolution of sequence specificity in restriction enzymes and may help engineer REases with novel specificities. Before this analysis was carried out, neither the three-dimensional fold of NlaIV, its evolutionary relationships or its catalytic or DNA-binding residues were known. Hence our analysis may be regarded as a paradigm for studies aiming at reducing 'white spaces' on the evolutionary landscape of sequence-function relationships by combining bioinformatics with simple experimental assays.     

Detection of subunit interfacial modifications by tracing the evolution of clamp-loader complex

The archaeal and eukaryal clamp-loader and clamp proteins were investigated with the evolutionary trace method. The molecular phylogeny of the proteins suggested that the hetero-pentameric complex of the archaeal clamp-loader with two subunits (RFCL and RFCS) was not a preserved ancestral type, but a degenerated version of the eukaryal complex of five subunits (RFC1-5). The evolutionary trace of amino acid replacements during the course of subunit differentiation revealed that the replacements had accumulated preferentially at the subunit interface regions. Some of the interfacial modifications that might be responsible for the specific interaction between the subunits were conserved in the current complex.     

Comparison of the construction of a 3-D model for human thromboxane synthase using P450cam and BM-3 as templates: implications for the substrate binding pocket

A 3-D model of human thromboxane A₂ synthase (TXAS) was constructedusing a homology modeling approach based on information fromthe 2.0 crystal structure of the hemoprotein domains of cytochromeP450BM-3 and P450cam. P450BM-3 is a bacterial fatty acid monooxygenaseresembling eukaryotic microsomal cytochrome P450s in primarystructure and function. TXAS shares 26.4% residue identity and48.4% residue similarity with the P450BM-3 hemoprotein domain.The homology score between TXAS and P450BM-3 is much higherthan that between TXAS and P450cam. Alignment between TXAS andthe P450BM-3 hemoprotein domain or P450cam was determined throughsequence searches. The P450BM-3 or P450cam main-chain coordinateswere spplied to the TXAS main chain in those sements where thetwo sequences were well aligned. These segments were linkedto one another using a fragment search method, and the sidechains were added to produce a 3-D model for TXAS. A TXAS substrate,prostaglandin H₂ (PGH₂) was docked into the TXAS cavity correspondingto the arachidonic acid binding pocket in P450BM-3 or camphorbinding site in P450cam. Regions of the heme and putative PGH2binding cavities in the TXAS model were identified and analyzed.The segments and residues involved in the active-site pocketof the TXAS model provide reasonable candidates for TXAS proteinengineering and inhibitor design. Comparison of the TXAS modelbased on P450BM-3 with another TXAS model based on the P450BM-3with another TXAS model based on the P450cam structure indicatedthat P450BM-3 is a more suitable template for homology modelingof TXAS.     

Enhancing the thermal stability of mitochondrial cytochrome b5 by introducing a structural motif characteristic of the less stable microsomal isoform

Outer mitochondrial membrane cytochrome b5 (OM b5) is the most thermostable cytochrome b5 isoform presently known. Herein, we show that OM b5 thermal stability is substantially enhanced by swapping an apparently invariant motif in its heme-independent folding core with the corresponding motif characteristic of its less stable evolutionary relative, microsomal cytochrome b5 (Mc b5). The motif swap involved replacing two residues, Arg15 with His and Glu20 with Ser, thereby introducing a Glu11-His15-Ser20 H-bonding triad on the protein surface along with a His15/Trp22 pi-stacking interaction. The ferric and ferrous forms of the OM b5 R15H/E20S double mutant have thermal denaturation midpoints (Tm values) of approximately 93 degrees C and approximately 104 degrees C, respectively. A 15 degrees C increase in apoprotein Tm plays a key role in the holoprotein thermal stability enhancement, and is achieved by one of the most common natural mechanisms for stabilization of thermophilic versus mesophilic proteins: raising the unfolding free energy along the entire stability curve.     

Modulation of enzyme activity by antibody binding to an alkaline phosphatase-epitope hybrid protein

An epitope from the HTV-1 gpl20 protein V3 loop has been insertedonto the surface of bacterial alkaline phosphatase at differentpositions in the vicinity of the enzyme active site, creatinghybrid proteins that can bind to an anti-gpl20 monoclonal antibody.One of the hybrid proteins, API1, has a 13 amino acid V3 loopsequence inserted between residues 407 and 408 of alkaline phosphatase.The enzymatic activity of this protein is modulated upon antibodybinding. API1 maintains the full activity of the wild type alkalinephosphatase but in the presence of the anti-gpl20 antibody,the enzyme activity is inhibited by 40–50%. Thus, thehybrid enzyme can be used to detect the presence of antibodyin solution. The concept of signalling proteins may have a wideapplication. Two models for the mechanism of modulation, sterichindrance and allosteric regulation, are discussed.     

Engineering stability into Escherichia coli secreted Fabs leads to increased functional expression

The recombinant expression of immunoglobulin domains, Fabs and scFvs in particular, in Escherichia coli can vary significantly from antibody to antibody. We hypothesized that poor Fab expression is often linked to poor intrinsic stability. To investigate this further, we applied a novel approach for stabilizing a poorly expressing anti-tetanus toxoid human Fab with a predisposition for being misfolded and non-functional. Forty-five residues within the Fab were chosen for saturation mutagenesis based on residue frequency analysis and positional entropy calculations. Using automated screening, we determined the approximate midpoint temperature of thermal denaturation (TM) for over 4000 library members with a maximum theoretical diversity of 855 unique mutations. This dataset led to the identification of 11 residue positions, primarily in the Fv region, which when mutated enhanced Fab stability. By combining these mutations, the TM of the Fab was increased to 92 degrees C. Increases in Fab stability correlated with higher expressed Fab yields and higher levels of properly folded and functional protein. The mutations were selected based on their ability to increase the apparent stability of the Fab and therefore the exact mechanism behind the enhanced expression in E.coli remains undefined. The wild-type and two optimized Fabs were converted to an IgG1 format and expressed in mammalian cells. The optimized IgG1 molecules demonstrated identical gains in thermostability compared to the Fabs; however, the expression levels were unaffected suggesting that the eukaryotic secretion system is capable of correcting potential folding issues prevalent in E.coli. Overall, the results have significant implications for the bacterial expression of functional antibody domains as well as for the production of stable, high affinity therapeutic antibodies in mammalian cells.