A new substitution matrix for protein sequence searches based on contact frequencies in protein structures

The instabilities of the native structures of mutant proteinswith an amino acid exchange are estimated by using the contactenergy and the number of contacts for each type of amino acidpair, which were estimated from 18 192 residue–residuecontacts observed in 42 crystals of globular proteins. Theywere then used to evaluate a transition probability matrix ofcodon substitutions and a log relatedness odds matrix, whichis used as a scoring matrix to measure the similarity betweenprotein sequences. To consider amino acid substitutions in homologousproteins, base mutation rates and the effects of the geneticcode are also taken into account. The average fitness of anamino acid exchange is approximated to be proportional to thestructural stability of the mutant protein, which is then approximatedby the average energy change of the protein native structureexpected for the ammo acid exchange with neglect of the energychange of the denatured state. In global and local homologysearches, this scoring matrix tends to yield significantly higheralignment scores than either the unitary matrix or the geneticcode matrix, and also may yield higher alignment scores fordistantly related protein pairs than MDM78. One of advantagesof this scoring matrix is that the equilibrium frequencies ofcodons and also base mutation rates can be adjusted.     

Modulating D-amino acid oxidase substrate specificity: production of an enzyme for analytical determination of all D-amino acids by directed evolution

Recent research on the flavoenzyme D-amino acid oxidase from Rhodotorula gracilis (RgDAAO) has revealed new, intriguing properties of this catalyst and offers novel biotechnological applications. Among them, the reaction of RgDAAO has been exploited in the analytical determination of the D-amino acid content in biological samples. However, because the enzyme does not oxidize acidic D-amino acids, it cannot be used to detect the total amount of D-amino acids. We now present the results obtained using a random mutagenesis approach to produce RgDAAO mutants with a broader substrate specificity. The libraries of RgDAAO mutants were generated by error-prone PCR, expressed in BL21(DE3)pLysS Escherichia coli cells and screened for their ability to oxidize different substrates by means of an activity assay. Five random mutants that have a 'modified' substrate specificity, more useful for the analytical determination of the entire content of D-amino acids than wild-type RgDAAO, have been isolated. With the only exception of Y223 and G199, none of the effective amino acid substitutions lie in segments predicted to interact directly with the bound substrate. The substitutions appear to cluster on the protein surface: it would not have been possible to predict that these substitutions would enhance DAAO activity. We can only conclude that these substitutions synergistically generate small structural changes that affect the dynamics and/or stability of the protein in a way that enhances substrate binding or subsequently catalytic turnover.     

Directed evolution of a type I antifreeze protein expressed in Escherichia coli with sodium chloride as selective pressure and its effect on antifreeze tolerance

Both freezing tolerance and NaCI tolerance are improved whenantifreeze proteins are expressed as fusion proteins with twodomains of staphylococcal protein A (SPA) in Escherichia coli.To characterize these properties further we created a randomlymutated expression library in E.coli, based on the winter flounderantifreeze protein HPLC-8 component gene. Low-fidelity PCR productsof this gene were fused to the spa gene encoding two domainsof the SPA. The library was screened for enhanced NaCl toleranceand four clones were selected. The freezing tolerance of eachof the selected clones was enhanced to varying extents. DNAsequencing of the isolated mutants revealed that the amphiphilicproperties of the native antifreeze protein were essentiallyconserved. Furthermore, by studying the primary sequence ofthe randomly mutated clones, in comparison with the degree offreezing tolerance, we have identified clues which help in understandingthe relationship between salt and freezing tolerance.     

Structural identity between the iron- and manganese-containing superoxide dismutases

We have recently reported the first complete amino acid sequenceof an iron-containing superoxide dismutase. The iron enzymeis thought to be closely homologous to the manganese-containingsuperoxide dismutases. The availability of complete amino acidsequence information for four manganese superoxide dismutasesand the crystal structures for two iron and two manganese superoxidedismutases prompted us to investigate the degree of homologybetween the two proteins at various levels. We report that itis not possible to clearly distinguish the two proteins on thebasis of their secondary or tertiary structures. It would appearthat a small number of single site substitutions are responsiblefor conferring distinguishing properties between the two proteins.Substitution of glyclne 77 and glutamine 154 by a glutamineand an alanine respectively in Photobacterium leiognathi ironsuperoxide dismutase may distinguish the kinetic and other particularproperties of this protein from the manganese protein (and otheriron superoxide dismutases). Furthermore the primary structureof both the iron and manganese proteins does not appear to haveany homology with any other known amino acid sequence.     

Position-dependent protein mutant profile based on mean force field calculation

The application of the mean force field in protein mutant stabilityprediction is explored. Based on protein main chain characteristics,including polar fraction, accessibility and dihedral angles,the mean force field was constructed to evaluate the compatibilitybetween an amino acid residue and its environment, from whicha position-dependent protein mutant profile was constructed.At each position along a protein sequence, the native residuewas replaced by the other 19 types of amino acid residues. Thematches were evaluated by energies from mean force field calculation,from which a mutant profile along the protein sequence was derived.General characteristics of such a profile were analyzed. Mutantstabilities for two sets of mutants in two proteins were foundto be reasonable compared with experimental data, which indicatesthat the present method can act as a guide in protein engineeringand as an effective scoring matrix in protein sequence–structurealignment studies.     

Analysis of protein conformational characteristics related to thermostability

The thermal stability of proteins was studied, 195 single aminoacid residue replacements reported elsewhere being analysedfor several protein conformational characteristics: type ofresidue replacement; conservative versus nonconservative substitution;replacement being in a homologous stretch of amino acid residues;change in hydrogen bond, van der Waals and secondary structurepropensities; solvent-accessible versus inaccessible replacement;type of secondary structure involved in the substitution; thephysico-chemical characteristics to which the thermostabilityenhancement can be attributed; and the relationship of the replacementsite to the folding intermediates of the protein, when known.From the above analyses, some general rules arise which suggestwhere amino acid substitutions can be made to enhance proteinthermostability: substitutions are conservative according tothe Dayhoff matrix; mainly occur on conserved stretches of residues;preferentially occur on solvent-accessible residues; maintainor enhance the secondary structure propensity upon substitution;contribute to neutralize the dipole moment of the caps of helicesand strands; and tend to increase the number of potential hydrogenbonding or van der Waals contacts or improve hydrophobic packing.     

A simple approach for protein structure discrimination based on the network pattern of conserved hydrophobic residues

Evolutionarily conserved hydrophobic residues at the core of protein structures are generally assumed to play a structural role in protein folding and stability. Recent studies have implicated that their importance to protein structures is uneven, with a few of them being crucial and the rest of them being secondary. In this work, we explored the possibility of employing this feature of native structures for discriminating non-native structures from native ones. First, we developed a network tool to quantitatively measure the structural contributions of individual amino acid residues. We systematically applied this method to diverse fold-type sets of native proteins. It was confirmed that this method could grasp the essential structural features of native proteins. Next, we applied it to a number of decoy sets of proteins. The results indicate that such an approach indeed identified non-native structures in most test cases. This finding should be of help for the investigation of the fundamental problem of protein structure prediction.     

Synthesis and sequence optimization of GFP mutants containing aromatic non-natural amino acids at the Tyr66 position

In order to alter the fluorescence properties of green fluorescent protein (GFP), aromatic non-natural amino acids were introduced into the Tyr66 position of GFP in a cell-free translation system using a four-base codon method. Two non-natural mutants (O-methyltyrosine and p-aminophenylalanine mutants) out of 18 mutants showed blue-shifted but weak fluorescence compared with wild-type GFP. Then the aminophenylalanine mutant was sequence optimized by introducing random mutations around the Tyr66 site. For this purpose, a method for random mutation of non-natural proteins in a cell-free system was developed. Three aminophenylalanine mutants with Y145F, Y145L and Y145 M mutations were obtained, which exhibited increased fluorescence by 1.5-, 3- and 4-fold, respectively. These results indicate that random mutation around non-natural amino acids is useful strategy in order to improve protein functions that are reduced by non-natural amino acid incorporation. The method described here will be applicable to other non-natural mutant proteins in a high-throughput manner.     

Coordinated amino acid changes in homologous protein families

In the tobamovirus coat protein family, amino acid residuesat some spatially close positions are found to be substitutedin a coordinated manner [Altschuh et al. (1987) J. Mol. Biol.,193,693]. Therefore, these positions show an identical patternof amino acid substitutions when amino acid sequences of thesehomologous proteins are aligned. Based on this principle, coordinatedsubstitutions have been searched for in three additional proteinfamilies: serine proteases, cysteine proteases and the haemoglobins.Coordinated changes have been found in all three protein familiesmostly within structurally constrained regions. This methodworks with a varying degree of success depending on the functionof the proteins, the range of sequence similarities and thenumber of sequences considered. By relaxing the criteria forresidue selection, the method was adapted to cover a broaderrange of protein families and to study regions of the proteinshaving weaker structural constraints. The information derivedby these methods provides a general guide for engineering ofa large variety of proteins to analyse structure–functionrelationships.     

Negatively charged purification tags for selective anion-exchange recovery

A novel strategy for the highly selective purification of recombinant fusion proteins using negatively charged protein domains, which were constructed by protein design, is described. A triple alpha-helical domain of 58 amino acids was used as scaffold. Far-ultraviolet circular dichroism measurements showed that the designed domains had very low alpha-helicity in a low-conductivity environment in contrast to the scaffold. The secondary structure could be induced by adding salt, giving a structure comparable to the parental molecule. Further studies showed that the new domains were able to bind to an anion exchanger even at pH values down to 5 and 6. Gene fusions between one of the designed domains and different target proteins, such as green fluorescent protein (GFP), maltose binding protein (MBP) and firefly luciferase, were also constructed. These gene products could be efficiently purified from whole cell lysates at pH 6 using anion-exchange chromatography.     

Contribution of amino acid substitutions at two different interior positions to the conformational stability of human lysozyme

To elucidate correlative relationships between structural changeand thermodynamic stability in proteins, a series of mutanthuman lysozymes modified at two buried positions (Ile56 andIle59) were examined. Their thermodynamic parameters of denaturationand crystal structures were studied by calorimetry and X-raycrystallography. The mutants at positions 56 and 59 exhibiteddifferent responses to a series of amino acid substitutions.The changes in stability due to substitutions showed a linearcorrelation with changes in hydrophobicity of substituted residues,having different slopes at each mutation site. However, thestability of each mutant was found to be represented by a uniqueequation involving physical properties calculated from mutantstructures. By fitting present and previous stability data formutant human lysozymes substituted at various positions to theequation, the magnitudes of the hydrophobicity of a carbon atomand the hydrophobicity of nitrogen and neutral oxygen atomswere found to be 0.178 and –0.013 kJ/mol.Ų, respectively.It was also found that the contribution of a hydrogen bond witha length of 3.0 Å to protein stability was 5.1 kJ/moland the entropy loss of newly introduction of a water moleculeswas 7.8 kJ/mol.     

The prediction and orientation of {alpha}-helices from sequence alignments: the combined use of environment-dependent substitution tables, Fourier transform methods and helix capping rules

Amino acid substitution tables are used to estimate the extentto which amino acids in families of homologous proteins areexposed to the solvent. The approach depends on the comparisonof difference environment-dependent tables for solvent accessible/inaccessibleresidues with amino acid substitutions at each position in analigned set of sequences. The periodicity in the predicted accessible/inaccessibleresidues is calculated using a Fourier transform procedure modifiedfrom that used to calculate hydrophobic moments. a-Helices areidentified from the characteristic periodicities and the solventaccessible face of the helix is defined. The initial helix predictionsare refined using rules for identifying the N- and C-terminiof helices from sequence alignments. These rules have been definedfrom a study of protein structures. The combined method correctlypredicts 79% of the residues in helices and incorrectly predictsonly 12% of the nonhelical residues as helical. In addition,since the method is reliable at predicting the correct numberof helices in the correct position in the sequence and sinceit also predicts the internal face of each helix, the resultscan be used to postulate 3-D arrangements of the secondary structureelements.     

A combinatorial library of an {alpha}-helical bacterial receptor domain

The construction and characterization of a combinatorial libraryof a solvent-exposed surface of an -helical domain derived froma bacterial receptor is described. Using a novel solid-phaseapproach, the library was assembled in a directed and successivemanner utilizing single-stranded oligonucleotides containingmultiple random substitutions for the variegated segments ofthe gene fragment The simultaneous substitution of 13 residuesto all 20 possible amino acids was carried out in a region spanning81 nucleotides. The randomization was made in codons for aminoacids that were modelled to be solvent accessible at a surfacemade up from two of the three a-helices of a monovalent Fc-bindingdomain of staphylococcal protein A. After cloning of the PCR-amplifiedlibrary into a phagemid vector adapted for phage display ofthe mutants, DNA sequencing analysis suggested a random distributionof codons in the mutagenized positions. Four members of thelibrary with multiple substitutions were produced in Escherichiacoli as fusions to an albumin-binding affinity tag derived fromstreptococcal protein G. The fusion proteins were purified byhuman serum albumin affinity chromatography and subsequentlycharacterized by SDSelectrophoresis, CD spectroscopy and biosensoranalysis. The analyses showed that the mutant protein A derivativescould all be secreted as soluble full-length proteins. Furthermore,the CD analysis showed that all mutants, except one with a prolineintroduced into helix 2, have secondary structures in closeagreement with the wild-type domain. These results proved thatmembers of this -helical receptor library with multiple substitutionsin the solvent-exposed surface remain stable and soluble inE.coli. The possibility of using this library for a phenotypicselection strategy to obtain artificial antibodies with novelfunctions is discussed.     

Engineering surface charges in a subtilisin: the effects on electrophoretic and ion-exchange behaviour

The introduction or removal of multiple charged amino acid residuesin the subtilisin Savinase^TM by genetic engineering allowedus to modify the electrostatic properties of this enzyme ina systematic way. The effects of these charge changes were investigatedexperimentally using ion-exchange chromatography and electrophoreticmobility in native gels all under identical conditions. Theexperiments clearly demonstrated that the overall charge ofa given protein is not the only factor determining electrophoreticmobility at low or moderate ionic strengths. For a series ofvariants having identical overall positive charge a linear relationwas observed between mobility towards the cathode and the totalnumber of charged residues present. This effect was found todepend on the type of (chloride) salt used: calcium ions giverise to complete screening of all negative charges, whereasonly partial screening is found for magnesium and sodium ions.In contrast, in the presence of sodium phosphate the overallcharge of the enzyme becomes slightly negative. These data indicatethat cations as well as anions may strongly perturb the overallcharge of proteins depending on the type of salt and on thenumber of charged amino acid residues present. The ion-exchangebehaviour demonstrated similar results, i.e. showing strongerenzyme adsorption with increasing numbers of surface chargeson a cation-exchange column run below the isoelectric pointof the proteins. However, the apparent sign reversal noted abovefor electrophoresis with sodium phosphate did not appear inthe ion-exchange experiments. This work provides further insightinto the adsorption of proteins to surfaces and the role playedby small ions, particularly when electrostatic forces dominatethe adsorption process.     

Does a backwardly read protein sequence have a unique native state?

Amino acid sequences of native proteins are generally not palindromic.Nevertheless, the protein molecule obtained as a result of readingthe sequence backwards, i.e. a retro-protein, obviously hasthe same amino acid composition and the same hydrophobicityprofile as the native sequence. The important questions whicharise in the context of retro-proteins are: does a retro-proteinfold to a well defined native-like structure as natural proteinsdo and, if the answer is positive, does a retro-protein foldto a structure similar to the native conformation of the originalprotein? In this work, the fold of retro-protein A, originatedfrom the retro-sequence of the B domain of Staphylococcal proteinA, was studied. As a result of lattice model simulations, itis conjectured that the retro-protein A also forms a three-helixbundle structure in solution. It is also predicted that thetopology of the retro-protein A three-helix bundle is that ofthe native protein A, rather than that corresponding to themirror image of native protein A. Secondary structure elementsin the retro-protein do not exactly match their counterpartsin the original protein structure; however, the amino acid sidechain contact pattern of the hydrophobic core is partly conserved.     

Mapping of residues involved in the interaction between the Bacillus subtilis xylanase A and proteinaceous wheat xylanase inhibitors

The Bacillus subtilis xylanase A was subjected to site-directed mutagenesis, aimed at changing the interaction with Triticum aestivum xylanase inhibitor, the only wheat endogenous proteinaceous xylanase inhibitor interacting with this xylanase. The published structure of Bacillus circulans XynA was used to target amino acids surrounding the active site cleft of B.subtilis XynA for mutation. Twenty-two residues were mutated, resulting in 62 different variants. The catalytic activity of active mutants ranged from 563 to 5635 XU/mg and the interaction with T.aestivum xylanase inhibitor showed a similar variation. The results indicate that T.aestivum xylanase inhibitor interacts with several amino acid residues surrounding the active site of the enzyme. Three different amino acid substitutions in one particular residue (D11) completely abolished the interaction between T.aestivum xylanase inhibitor and B.subtilis xylanase A.     

Improvement in the alkaline stability of subtilisin using an efficient random mutagenesis and screening procedure

An efficient random mutagenesis procedure coupled to a replicaplate screen facilitated the isolation of mutant subtilisinsfrom Bacillus amyloliquefaciens that had altered autolytic stabilityunder alkaline conditions. Out of about 4000 clones screened,approximately 70 produced subtilisins with reduced stability(negatives). Two dones produced a more stable subtilisin (positives)and were identified as having a single mutation, either IIe107Valor Lys2l3Arg (the wild-type amino acid is followed by the codonposition and the mutant amino acid). One of the negative mutants,Met50Val, was at a site where other homologous subtilisins containeda Phe. When the Met50Phe mutation was introduced into the B.amyloliquefaciens gene, the mutant subtilisin was more alkalinestable. The double mutant IIe107Val/Lys2l3Arg) was more stablethan the isolated single mutant parents. The triple mutant (Met50Phe/IIel07Val/Lys2l3Arg)was even more stable than IIe107Val/Lys2l3Arg (up to two timesthe autolytic half-time of wild-type at pH 12). These studiesdemonstrate the feasibility for improving the alkaline stabilityof proteins by random mutagenesis and identifying potentialsites where substitutions from homologous proteins can improvealkaline stability.     

Red fluorescent protein variants with incorporated non-natural amino acid analogues

Fluorescent proteins are important tools in biotechnology applications and biosensing. DsRed, a red fluorescent protein, has expanded the colors of fluorescent proteins beyond the more commonly used green fluorescent protein. Many genetic modifications have been performed on DsRed to overcome some of its drawbacks. These primarily focused on overcoming the oligomerization detrimental to DsRed activity, and the parasitic green fluorescence caused by the immature chromophore. One such variant, DsRed-monomer, has minimal green fluorescence and no oligomerization. A few traditional mutagenesis studies have been done with DsRed and its mutants to shift the fluorescence wavelengths creating additions to the pallet of fluorescent protein colors. We have explored incorporation of non-natural amino acid analogues into DsRed-Monomer, obtaining variants with differing emission properties. In this work, two such analogues of tyrosine have been incorporated into DsRed-Monomer: 3-amino-l-tyrosine and 3-fluoro-l-tyrosine. Tyrosine analogues were chosen due to the role of tyrosine in the formation and structure of the protein's chromophore. The variants obtained in our study showed altered emission wavelengths and spectral characteristics. Our study demonstrates that incorporation of non-natural analogues into DsRed-Monomer is a viable approach to alter the spectral characteristics of the protein. We envision that this study will open up the door to non-natural mutagenesis studies with red fluorescent proteins and its mutants.     

Downregulation of eRF1 by RNA interference increases mis-acylated tRNA suppression efficiency in human cells

The site-specific incorporation of non-natural amino acids into proteins by nonsense suppression has been widely used to investigate protein structure and function. Usually this technique exhibits low incorporation efficiencies of non-natural amino acids into proteins. We describe for the first time an approach for achieving an increased level of nonsense codon suppression with synthetic suppressor tRNAs in cultured human cells. We find that the intracellular concentration of the eukaryotic release factor 1 (eRF1) is a critical parameter influencing the efficiency of amino acid incorporation by nonsense suppression. Using RNA interference we were able to lower eRF1 gene expression significantly. We achieved a five times higher level of amino acid incorporation as compared with non-treated control cells, as demonstrated by enhanced green fluorescent protein (EGFP) fluorescence recovery after importing a mutated reporter mRNA together with an artificial amber suppressor tRNA.     

Local polarity analysis: a sensitive method that discriminates between native proteins and incorrectly folded models

The evaluation of calculated protein structures is an importantstep in the protein design cycle. Known criteria for this assessmentof proteins are the polar and apolar, accessible and buriedsurface area, electrostatic interactions and other interactionsbetween the protein atoms (e.g. HO, S-S),atomic packing, analysisof amino acid environment and surface charge distribution. Weshow that a powerful test of accuracy of protein structure canbe derived by analysing the water contact of atoms and additionallytaking into account their polarity. On the basis of estimatedreference values of the polar fraction of typical globular proteinswith known structure (mean, SD and distribution), the evaluationof misfolded structures can be improved significantly. The referencevalues are derived by moving windows of different length (3–99amino acid residues) over the amino acid sequence. Model proteins,which are included in the Brookhaven protein structure databank,deliberately misfolded proteins, hypothetical proteins and predictedprotein structures are diagnosed as at least partially incorrectlyfolded. The local fault, mostly observed, is that polar groupsare buried too frequently in the interior of the protein. Thedatabase-derived quantities are useful in screening the designedproteins prior to experimentation and may also be useful inthe assessment of errors in the experimentally determined proteinstructures.