The structure of interfaces between subunits of dimeric and tetrameric proteins

The structures of the interfaces of nine dimeric and nine tetramericproteins have been analyzed and have been seen to follow generalprinciples. These interfaces are combinations of four structuralmotifs, which resemble features of monomeric proteins. Theseare: (i) extended beta sheet; (ii) helix–helix packing;(iii) sheet–sheet packing; and (iv) loop interactions.Other common structural features in the interfaces studied aretwo-fold symmetry, charged hydrogen bonds and channel formation(found only in tetramers). Monomer–monomer interfacesare intermediate in hydrophobicity and charge between the interfacesbetween secondary structures of monomeric proteins and the exteriorsof monomeric proteins. A typical interface has one of the firstthree of the structural motifs at its centre and loop interactionsaround the outside, where most of the charge resides.     

Polarity as a criterion in protein design

Hypothetical proteins can be tested computationally by determiningwhether or not the designed sequence-structure pair has thecharacteristics of a typical globular protein. We have developedsuch a test by deriving quantities with approximately constantvalue for all globular proteins, based on empirical analysisof the exposed and buried surfaces of 128 structurally knownproteins. The characteristic quantities that best appear tosegregate badly designed or deliberately misfolded proteinsfrom their properly folded natural relatives are the polar fractionof side chains on the protein surface and, independently, inthe protein interior. Three of the seven hypothetical structurestested here can be rejected as having too many polar side-chaingroups in the interior or too few on the protein surface. Inaddition, a recently designed nutritional protein is identifiedas being very much unlike globular proteins. These database-derivedcharacteristic quantities are useful in screening designed proteinsprior to experiment and may be useful in screening experimentallydetermined (X-ray, NMR) protein structures for possible errors.     

A novel method for predicting transmembrane segments in proteins based on a statistical analysis of the SwissProt database: the PRED-TMR algorithm

We present a novel method that predicts transmembrane domainsin proteins using solely information contained in the sequenceitself. The PRED-TMR algorithm described, refines a standardhydrophobicity analysis with a detection of potential termini(`edges', starts and ends) of transmembrane regions. This allowsone both to discard highly hydrophobic regions not delimitedby clear start and end configurations and to confirm putativetransmembrane segments not distinguishable by their hydrophobiccomposition. The accuracy obtained on a test set of 101 non-homologoustransmembrane proteins with reliable topologies compares wellwith that of other popular existing methods. Only a slight decreasein prediction accuracy was observed when the algorithm was appliedto all transmembrane proteins of the SwissProt database (release35). A WWW server running the PRED-TMR algorithm is availableat http://o2.db.uoa.gr/PRED-TMR/     

The role of heat-shock and chaperone proteins in protein folding: possible molecular mechanisms

Recently some heat-shock proteins have been linked to functionsof ‘chaperoning’ protein folding in vivo. Here currentexperimental evidence is reviewed and possible requirementsfor such an activity are discussed. It is proposed that onemode of chaperone action is to actively unfold misfolded orbadly aggregated proteins to a conformation from whkh they couldrefold spontaneously; that improperly folded proteins are recognizedby excessive stretches of solvent-exposed backbone, rather thanby exposed hydrophobic patches; and that the molecular mechanismfor unfolding is either repeated binding and dissociation (‘plucking’)or translocation of the protein backbone through a binding cleft(‘threading’), allowing the threaded chain to refoldspontaneously. The observed hydrolysis of ATP would providethe energy for active unfolding. These hypotheses can be appliedto both monomeric folding and oligomeric assembly and are sufficientlydetailed to be open to directed experimental verification.     

Electrostatic effects in proteins: comparison of dielectric and charge models

Two approaches for calculating electrostatic effects in proteinsare compared and an analysis is presented of the dependenceof calculated properties on the model used to define the chargedistribution. Changes in electrostatic free energy have beencalculated using a screened Coulomb potential (SCP) with a distance-dependenteffective dielectric permittivity to model bulk solvent effectsand a finite difference approach to solve the Poisson-Boltzmann(FDPB) equation. The properties calculated include shifts indissociation constants of ionizable groups, the effect of annihilatingsurface charges on the binding of metals, and shifts in redoxpotentials due to changes in the charge of ionizable groups.In the proteins considered the charged sites are separated by3.5–12 Å. It is shown that for the systems studiedin this distance range the SCP yields calculated values whichare at least as accurate as those obtained from solution ofthe FDPB equation. In addition, in the distance range 3–5Å the SCP gives substantially better results than theFDPB equation. Possible sources of this difference between thetwo methods are discussed. Shifts in binding constants and redoxpotentials were calculated with several standard charge sets,and the resulting values show a variation of 20–40% betweenthe best and worst cases. From this study it is concluded thatin most applications, changes in electrostatic free energiescan be calculated economically and reliably using an SCP approachwith a single functional form of the screening function.     

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.     

A new approach to artificial and modified proteins: theory-based design, synthesis in a cell-free system and fast testing of structural properties by radiolabels

A novel approach to the creation of artificial and modifiedproteins has been elaborated. The approach includes a sequencedesign based on the molecular theory of protein secondary structureand folding patterns, gene expression in a cell-free systemand testing of structural properties of the synthesized polypeptidesat a nanogram level using radiolabelled chains. The approachhas been applied to a new synthetic protein albebetin whichhas been designed to form a 3-D fold which does not contradictany structural rule but has been never observed up to now innatural proteins. Using size-exclusion chromatography, urea-gradientelectrophoresis and limited proteolysis of a radiolabelled chain,it has been shown that the artificial protein is nearly as compactas natural proteins, cooperatively unfolds at high urea concentrationsand has some structural features of a definite structure consistentwith the designed one. As albebetin has been designed as consistingof two structural repeats, a ‘halfalbebetin’ (oneof these repeats) has also been synthesized and studied. Itwas shown that ‘half-albebetin’ is also compact     

Can a simple function for the dielectric response model electrostatic effects in globular proteins?

The relationship between the effective dielectric constant thatmodels the electrostatic effect from a charged side chain ina protein was evaluated both experimentally and theoretically.Experimental values were obtained from the shifts in pK_a thatresulted from point mutations of side chains in subtilisin.Theoretical values were obtained from an iterative solutionto Poisson's equation that considers the dielectric responseof the protein and the solvent together with charge positions.There is no simple relationship between the effective dielectricconstant and the distance from the charge responsible for theinteractions. For some charge positions a linear but not a directproportional relationship of the effective dielectric with distanceof separation was observed. Thus, simple models such as a lineardistance-dependence for the dielectric response are not suitableto evaluate electrostatic effects in proteins.     

Structure activity relationships of monocyte chemoattractant proteins in complex with a blocking antibody

Monocyte chemoattractant proteins (MCPs) are cytokines that direct immune cells bearing appropriate receptors to sites of inflammation or injury and are therefore attractive therapeutic targets for inhibitory molecules. 11K2 is a blocking mouse monoclonal antibody active against several human and murine MCPs. A 2.5 A structure of the Fab fragment of this antibody in complex with human MCP-1 has been solved. The Fab blocks CCR2 receptor binding to MCP-1 through an adjacent but distinct binding site. The orientation of the Fab indicates that a single MCP-1 dimer will bind two 11K2 antibodies. Several key residues on the antibody and on human MCPs were predicted to be involved in antibody selectivity. Mutational analysis of these residues confirms their involvement in the antibody-chemokine interaction. In addition to mutations that decreased or disrupted binding, one antibody mutation resulted in a 70-fold increase in affinity for human MCP-2. A key residue missing in human MCP-3, a chemokine not recognized by the antibody, was identified and engineering the preferred residue into the chemokine conferred binding to the antibody.     

GRAFTER: a computational aid for the design of novel proteins

The GRAFTER suite of programs provides geometric search andevaluation functions that simplify and automate the processof identifying the best scaffolds for a particular structuralmotif. Three applications of the GRAFTER suite are presented.Potential grafts between repressor and 434 repressor were identifiedthat should change the DNA binding specificity of these repressors.These results are compared with site-directed mutagenesis experimentsthat have been shown to alter repressor-DNA binding specificity.Next, 26 loops from antibody structures were grouped into familiesof similar structure. Grafts of antibody loops onto a pre-existingscaffold are an essential component of antibody humanization.Finally, interleukin (lL)-4 was searched as a scaffold thatmight accept the graft of a five residue epitope from humangrowth hormone (hGH). The existence of a crystal structure ofthe hGH-hGH receptor complex, extensive mutagenesis studiesof the hGH residues that contribute to the energetics of ligand-receptorinteractions and the gross structural homology between hGH andIL-4 make this an appealing computational target. The approachpresented here could aid the development of novel enzymes andbinding proteins     

Engineering a de novo-designed coiled-coil heterodimerization domain for the rapid detection, purification and characterization of recombinantly expressed peptides and proteins

Using the techniques of genetic engineering and the principlesof protein de novo design, we have developed a unique affinitymatrix protein tag system as a rapid, convenient and sensitivemethod to detect, purify and characterize newly expressed recombinantpeptides or proteins from cell extracts. The method utilizestwo de novo-designed linear peptide sequences that can selectivelydimerize to form the stable protein motif, the two-stranded-helical coiled-coil. In this method, a recombinant bacterialexpression vector pRLDE has been engineered so that one of thedimerization strands (E-coil) is expressed as a C-terminal fusiontag on newly expressed peptides or proteins, while the other(K-coil) is either biotin-labeled for detection in a Westernblot-type format or immobilized on an insoluble silica supportfor selective dimerization affinity chromatography. Recombinantlyexpressed peptides from Escherichia coli containing the dimerizationtag have been produced, detected and purified using this method.The recombinant peptides were easily and clearly identifiedusing the biotin-labeled coil, while the single-step affinitypurification results indicated the purity of the affinity purifiedexpressed peptides to be >95%, as assessed by reversed-phasechromatography. The stability of the dimerization domain alsoallows for the purified peptide to be left attached to the matrix,thus creating a new peptide-bound column that can be used tostudy peptide–protein or peptide–ligand interactions.Therefore this system offers a new alternative to existing peptideor protein fusion tags and demonstrates the utility of a denovo-designed system.     

Convex constraint analysis: a natural deconvolution of circular dichroism curves of proteins

A new algorithm, called convex constraint analysis, has beendeveloped to deduce the chiral contribution of the common secondarystructures directly from experimental CD curves of a large numberof proteins. The analysis is based on CD data reported by YangJ.T., Wu,C.-S.C. and Martinez,H.M. [Methods Enzymol., 130, 208–269(1986)]. Application of the decomposition algorithm for simulatedprotein data sets resulted in component spectra [B(, i)] identicalto the originals and weights [C(i, k)] with excellent Pearsoncorrelation coefficients (R) [Chang,C.T., Wu,C.-S.C. and Yang,J.T.(1978) Anal. Biochem., 91,12–31]. Test runs were performedon sets of simulated protein spectra created by the Monte Carlotechnique using poly-L-lysine-based pure component spectra.The significant correlational coefficients (R >0.9) demonstratedthe high power of the algorithm. The algorithm, applied to globularprotein data, independent of X-ray data, revealed that the CDspectrum of a given protein is composed of at least four independentsources of chirality. Three of the computed component curvesshow remarkable resemblance to the CD spectra of known proteinsecondary structures. This approach yields a significant improvementin secondary structural evaluations when compared with previousmethods, as compared with X-ray data, and yields a realisticset of pure component spectra. The new method is a useful toolnot only in analyzing CD spectra of globular proteins but alsohas the potential for the analysis of integral membrane proteins.     

Use of a quantitative structure-property relationship to design larger model proteins that fold rapidly

A quantitative structure–property relationship (QSPR)was used to design model protein sequences that fold repeatedlyand relatively rapidly to stable target structures. The specificmodel was a 125-residue heteropolymer chain subject to MonteCarlo dynamics on a simple cubic lattice. The QSPR was derivedfrom an analysis of a database of 200 sequences by a statisticalmethod that uses a genetic algorithm to select the sequenceattributes that are most important for folding and a neuralnetwork to determine the corresponding functional dependenceof folding ability on the chosen attributes. The QSPR dependson the number of anti-parallel sheet contacts, the energy gapbetween the native state and quasi-continuous part of the spectrumand the total energy of the contacts between surface residues.Two Monte Carlo procedures were used in series to optimize boththe target structures and the sequences. We generated 20 fullyoptimized sequences and 60 partially optimized control sequencesand tested each for its ability to fold in dynamic MC simulations.Although sequences in which either the number of anti-parallelsheet contacts or the energy of the surface residues is non-optimalare capable of folding almost as well as fully optimized ones,sequences in which only the energy gap is optimized fold markedlymore slowly. Implications of the results for the design of proteinsare discussed.     

Mutational and structural analysis of the lectin activity in binding domain 2 of ricin B chain

The study of the lectin binding sites of ricin B chain and ofother homologous members of the small gene family that makeup ricin-like molecules has revealed a number of key contactresidues involved in sugar binding. In particular, on the basisof data generated by the X-ray crystallographic structure ofricin, comparisons of sequence homologies to other ricin-likemolecules and substrate binding studies with these molecules,it has been proposed that His248 of Ricinus communis agglutinin(RCA) B chain may interfere with galactose binding in the secondbinding domain of that lectin. To test that hypothesis, singlebinding domain 2 (SBD2) of ricin B chain was expressed as agene 3 fusion protein on the surface of fd phage to measuredirectly the effect of mutational changes on this binding site.Replacement of tyrosine with histidine at amino acid position248 of SBD2 of ricin B chain was shown to reduce lectin activity.The sequences of RCA and ricin B chains were aligned and comparedwith the tertiary structure of ricin B chain to select variousmutations that were introduced as controls in the study. Oneof these controls, Leu247 to Val247, displayed increased affinityfor galactosides. The role of sequence changes is discussedin relation to the structural and functional divergence in thesemolecules.     

Construction of a dictionary of sequence motifs that characterize groups of related proteins

An automatic procedure is proposed to identify, from the proteinsequence database, conserved amino acid patterns (or sequencemotifs) that are exclusive to a group of functionally relatedproteins. This procedure is applied to the PIR database anda dictionary of sequence motifs that relate to specific superfamiliesconstructed. The motifs have a practical relevance in identifyingthe membership of specific superfamilies without the need toperform sequence database searches in 20% of newly determinedsequences. The sequence motifs identified represent functionallyimportant sites on protein molecules. When multiple blocks existin a single motif they are often close together in the 3-D structure.Furthermore, occasionally these motif blocks were found to besplit by introns when the correlation with exon structures wasexamined.     

Distance distributions in proteins: a six-parameter representation

We present a statistical analysis of protein structures basedon interatomic C_a distances. The overall distance distributionsreflect in detail the contents of sequence-specific substructuresmaintained by local interactions (such as -helixes) and longerrange interactions (such as disulfide bridges and ß-sheets).We also show that a volume scaling of the distances makes distancedistributions for protein chains of different length superimposable.Distance distributions were also calculated specifically foramino acids separated by a given number of residues. Specificfeatures in these distributions are visible for sequence separationsof up to 20 amino acid residues. A simple representation, whichpreserves most of the information in the distance distributions,was obtained using six parameters only. The parameters giverise to canonical distance intervals and when predicting coarse-graineddistance constraints by methods such as data-driven artificialneural networks, these should preferably be selected from theseintervals. We discuss the use of the six parameters for determiningor reconstructing 3-D protein structures.     

Identification of important functional environs in protein tertiary structures from the analysis of residue variation in 3-D: application to cytochromes c and carboxypeptidases A and B

A simple methodology is described to apply to aligned proteinsequence sets for which at least one representative 3-D C structureis known. The evolutionary variation observed at each residueposition in the sequence alignment is qualified by taking intoaccount the residue variation that has occurred at other positionslocated within 7 A (according to the probable chain fold). Thisexpresses the evolutionary behaviour of any residue positionin the more appropriate context of its immediate surroundingsand distinguishes between invariant residues on the basis ofthe variation of their environment. The highest mechanisticsignificance is attached to conserved residues in conservedsurroundings, but the quantitative nature of the analysis meansthat all residue vicinities can be ranked and merged accordingto the degree of conservation that they exhibit and the residuepositions that comprise them. Therefore, with the aid of thechain fold, contour maps can be constructed that show gradedfoci of evolutionary conservation in the underlying superstructureof the protein type, and the irregular shapes and extents oflarge conserved areas. To test the methodology, it was appliedto cytochromes c and the carboxypeptidases A and B.     

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.     

Physicochemical factors for discriminating between soluble and membrane proteins: hydrophobicity of helical segments and protein length

The average hydrophobicity of a polypeptide segment is consideredto be the most important factor in the formation of transmembranehelices, and the partitioning of the most hydrophobic (MH) segmentinto the alternative nonpolar environment, a membrane or hydrophobiccore of a globular protein may determine the type of proteinproduced. In order to elucidate the importance of the MH segmentin determining which of the two types of protein results froma given amino acid sequence, we statistically studied the characteristicsof MH helices, longer than 19 residues in length, in 97 membraneproteins whose three-dimensional structure or topology is known,as well as 397 soluble proteins selected from the Protein DataBank. The average hydrophobicity of MH helices in membrane proteinshad a characteristic relationship with the length of the protein.All MH helices in membrane proteins that were longer than 500residues had a hydrophobicity greater than 1.75 (Kyte and Doolittlescale), while the MH helices in membrane proteins smaller than100 residues could be as hydrophilic as 0.1. The possibilityof developing a method to discriminate membrane proteins fromsoluble ones, based on the effect of size on the type of proteinproduced, is discussed.