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.     

The role of tyrosine residues in the DNA-binding site of the Pf1 gene 5 protein

The 144 amino acid gene 5 protein of bacteriophage Pf1 bindstightly and cooperatively to single-stranded DNA during replicationof the phage genome. It has been suggested that aromatic aminoacid side chains are important for this interaction, probablythrough base stacking with the DNA. We have analysed the accessibilityof tyrosine residues in the DNA—protein complex, and theirimportance to the DNA-binding activity of the protein, by chemicalmodification and protection experiments using tetranitromethane.Tyrosines 21, 30 and 55 are surface accessible in the free proteinbut are protected from modification in the complex with phageDNA. Moreover, modification of these residues in the free proteinabolishes the ability to bind to DNA or oligonucleotides, asjudged by fluorescence spectroscopy and gel retardation analysis.Modification of the protein also results in the formation ofan intersubunit covalent cross-link between Tyr55 and Phe76,suggesting that Phe76 is located within the DNA-binding cleftof the protein. It is proposed that residues 17–34 ofthe Pf1 gene 5 protein form a beta-hairpin analogous to the‘DNA-binding wing’ of the fd and Ike gene 5 proteins.We suggest the existence of a single-stranded DNA binding motif,in which Tyr30 of the Pf1 protein is equivalent to the functionallyimportant Tyr26 of the fd gene 5 protein.     

Prediction of the antigenic sites of the cystic fibrosis transmembrane conductance regulator protein by molecular modelling

Antibodies are powerful tools for studying the in situ localizationand physiology of proteins. The prediction of epitopes by molecularmodelling has been used successfully for the papilloma virus,and valuable antibodies have been raised [Muller et aL (1990)J. Gen. Virol, 71, 2709–2717]. We have improved the modellingapproach to allow us to predict epitopes from the primary sequencesof the cystic fibrosis transmembrane conductance regulator.The procedure involves searching for fragments of primary sequenceslikely to make amphipathic secondary structures, which are hydrophilicenough to be at the surface of the folded protein and thus accessibleto antibodies. Amphipathic helices were predicted using themethods of Berzofsky, Eisenberg and Jahnig. Their hydrophobichydrophilicinterface was calculated and drawn, and used to predict theorientation of the helices at the surface of the native protein.Amino acids involved in turns were selected using the algorithmof Eisenberg. Tertiary structures were calculated using ‘FOLDING’,a software developed by R.Brasseur for the prediction of smallprotein structures [Brasseur (1995) J. MoL Graphics, in press].We selected sequences that folded as turns with at least fiveprotruding polar residues. One important property of antibodiesis selectivity. To optimize the selectivity of the raised antibodies,each sequence was screened for similarity (FASTA) to the proteinsequences from several databanks. Ubiquitous sequences werediscarded. This approach led to the identification of 13 potentialepitopes in the cystic fibrosis transmembrane conductance regulator:seven helices and six loops.     

Toxicity-based selection of Escherichia coli mutants for functional recombinant protein production: application to an antibody fragment

We propose a novel approach to the selection of Escherichia coli bacterial strains improved for the production of recombinant functional proteins. This approach is based on aggregation-induced toxicity of recombinant proteins. We show that selection of clones displaying a reduced toxicity is an efficient means of isolating bacteria producing recombinant protein with reduced aggregation in favour of correct folding. For an efficient selection, we found that time of toxicity induction must be precisely determined and recombinant protein must be expressed as a fusion with a protein whose activity is easily detectable on plates, thus allowing elimination of non-productive mutants. Choosing the expression to the periplasmic space of an scFv fragment fused to the N-terminus of alkaline phosphatase as a model, we selected chromosomal mutations that reduce aggregation-induced toxicity and showed that they concomitantly improve production of a functional recombinant hybrid. The effects of the mutations isolated could then be cumulated with those of other strategies used for recombinant scFv production. Thus, we could ensure a 6- to 16-fold increase in production of a functional scFv-PhoA hybrid. This is the first report demonstrating the possibility of directly selecting on agar plates E.coli strains improved for functional recombinant protein production from a large bacterial mutant library.     

Application of scaled particle theory to model the hydrophobic effect: implications for molecular association and protein stability

The energetics of alkane dissolution and partition between waterand organic solvent are described in terms of the energy ofcavity formation and solute-solvent interaction using scaledparticle theory. Thermodynamic arguments are proposed that allowcomparison of experimental measurements of the surface areawith values calculated from an all-atom representation of thesolute. While the surface tension relating to the accessiblesurface is shape dependent, it is found that for the molecularsurface it is not. This model rationalizes the change in surfacetension between the microscopic (20–30 cal/mol/A2) andmacroscopic (70–75 cal/mol/A2) regimes without the needto invoke Flory-Huggins theory or to apply other corrections.The difference in the values arises (i) to a small extent asa result of the curvature dependence of surface tension and(ii) to a large extent due to the difference in the molecularsurface derived from the experiment and that calculated froman extended all-atom model. The model suggests that the primarydriving force for alkane association in water is due to thetendency of water to reduce the solute cavity surface. It isargued that to model the energetics of alkane association, thesurface tension should be related to the molecular surface (ratherthan the accessible surface) with a surface tension near themacroscopic limit for water. This model is compared with resultsfrom theoretical simulations of the hydrophobic effect for twowell-studied systems. The implications for antibody– antigeninteractions and the effect of hydrophobic amino acid deletionon protein stability are discussed. The approach can be usedto model the solute cavity formation energy in solution as afirst step in the continuum modelling of biomolecular interactions     

Spectroscopic investigation of structure in octarellin (a de novo protein designed to adopt the {alpha}/{beta}-barred packing)

We present here a spectroscopic structural characterizationof octarellin, a recently reported de novo protein modelledon /ß-barrel proteins [K. Go raj, A.Renard and J.A.Martial(1990) Protein Engng, 3, 259–266]. Infrared and Ramanspectra analyses of octarellin‘s secondary structure revealthe expected percentage of -helices (30%) and a higher ß-sheetcontent (40%) than predicted from the design. When the Ramanspectra obtained with octarellin and native triosephosphateisomerase (a natural /ß-barrel) are compared, similarpercentages of secondary structures are found. Thermal denaturationof octarellin monitored by CD confirms that its secondary structuresare quite stable, whereas its native-like tertiary fold is not.Tyrosine residues, predicted to be partially hidden from solvent,are actually exposed as revealed by Raman and UV absorptionspectra. We conclude that the attempted /ß-barrelconformation in octarellin may be loosely packed. The criteriaused to design octarellin are discussed and improvements suggested.     

Search for the most stable folds of protein chains: I. Application of a self-consistent molecular field theory to a problem of protein three-dimensional structure prediction

We present a general approach to the prediction of 3-D foldsof protein chains from their amino acid sequences. The approachis based on the use of the self-consistent molecular field theoryfor long-range interactions, the use of 1-D statistical mechanicsfor short-range interactions and on the discovery that thereis and should only be a relatively small discrete set of foldingpatterns. This makes it possible to examine the full varietyof ‘potentially stable’ folds and to determine thethermodynamically stable structure. In this paper, we give thegeneral theoretical background of the approach. The encouragingresults of the application of this approach to ß-domainsare described in another paper.     

Identifying the mechanism of protein loop closure: a molecular dynamics simulation of the Bacillus stearothermophilus LDH loopin solution

The ‘loop’ involving residues 98–110 in Bacillusstearothermophilus lactate dehydrogenase (BSLDH) is of greatinterest as substrate-induced ‘loop’ closure isthought to berate-limiting and essential in catalyzing the reactionand in determining substrate specificity. Consequently, we haveexplored the mechanism underlying ‘loop’ openingin BSLDH through a molecular dynamics simulation at high temperature(1000 K) in the presence of explicit solvent, starting fromthe X-ray structure of BSLDH complexed with the co-enzyme NAD+and oxamate at 2.5 Å. During the simulation, a significantconformational change occurred, as evidenced by sharp dihedralangle transitions, hydrogen bond breaking and formation andlarge root mean square deviations from the starting structure;these changes define the criteria for ‘loop’ opening.The mechanical elements responsible for ‘loop’ opening,i.e. ‘loop’ hinges and flap, are defined througha combination of order parameters, dihedral angle changes andtheir correlations and the dynamical cross-correlation map ofatomic displacements for the ‘loop’ residues. Theresults indicate that the ‘loop’ consists of twoflexible hinge regions on either side of a relatively rigidthree-residue segment that undergoes a significant spatial displacementduring ‘loop’opening. ‘Loop’ openingis made possible through an array of correlated dihedral anglechanges and intra-& ‘loop’ rearrangements ofhydrogen-bond interactions. The presentfindings are comparedto previous work related to ‘loop’ opening and site-directedmutagenesis experiments.     

Modelling the structure of latexin-carboxypeptidase A complex based on chemical cross-linking and molecular docking

We have determined the three-dimensional structure of the protein complex between latexin and carboxypeptidase A using a combination of chemical cross-linking, mass spectrometry and molecular docking. The locations of three intermolecular cross-links were identified using mass spectrometry and these constraints were used in combination with a speed-optimised docking algorithm allowing us to evaluate more than 3 x 10(11) possible conformations. While cross-links represent only limited structural constraints, the combination of only three experimental cross-links with very basic molecular docking was sufficient to determine the complex structure. The crystal structure of the complex between latexin and carboxypeptidase A4 determined recently allowed us to assess the success of this structure determination approach. Our structure was shown to be within 4 A r.m.s. deviation of Calpha atoms of the crystal structure. The study demonstrates that cross-linking in combination with mass spectrometry can lead to efficient and accurate structural modelling of protein complexes.     

A dual-affinity gene fusion system to express small recombinant proteins in a soluble form: expression and characterization of protein A deletion mutants

A novel gene fusion system to express and purify small recombinantproteins in Escherichia coli has been constructed. The conceptallows for affinity purification of soluble gene products bysequential albumin- and Zn²+-affinity chromatography. The dual-affinitysystem is well suited for expression of unstable proteins asonly full-length protein is obtained after purification andproteins gain proteolytic stability in the fusion protein. Herewe show that the dual-affinity approach can be used for theexpression of various unstable derivatives of a single IgG-bindingdomain based on staphylococcal protein A. Analysis of the proteolyticstabilities and the IgG-binding properties of the differentmutant proteins suggest that the model for the structure ofan IgG-binding domain must be re-evaluated.     

DOMPLOT: a program to generate schematic diagrams of the structural domain organization within proteins, annotated by ligand contacts

A program is described for automatically generating schematiclinear representations of protein chains in terms of their structuraldomains. The program requires the co-ordinates of the chain,the domain assignment, PROSITE information and a file listingall intermolecular interactions in the protein structure. Theoutput is a PostScript file in which each protein is representedby a set of linked boxes, each box corresponding to all or partof a structural domain. PROSITE motifs and residues involvedin ligand interactions are highlighted. The diagrams allow immediatevisualization of the domain arrangement within a protein chain,and by providing information on sequence motifs, and metal ion,ligand and DNA binding at the domain level, the program facilitatesdetection of remote evolutionary relationships between proteins.     

Grow to Fit Molecular Dynamics (G2FMD): an ab initio method for protein side-chain assignment and refinement

The rough energy landscapes and tight packing of protein interiors are two of the critical factors that have prevented the wide application of physics-based models in protein side-chain assignment and protein structure prediction in general. Complementing the rotamer-based methods, we propose an ab initio method that utilizes molecular mechanics simulations for protein side-chain assignment and refinement. By reducing the side-chain size, a smooth energy landscape was obtained owing to the increased distances between the side chains. The side chains then gradually grow back during molecular dynamics simulations while adjusting to their surrounding driven by the interaction energies. The method overcomes the barriers due to tight packing that limit conformational sampling of physics-based models. A key feature of this approach is that the resulting structures are free from steric collisions and allow the application of all-atom models in the subsequent refinement. Tests on a small set of proteins showed nearly 100% accuracy on both chi1 and chi2 of buried residues and 94% of them were within 20 degrees from the native conformation, 79% were within 10 degrees and 42% were within 5 degrees . However, the accuracy decreased when exposed side chains were involved. Further improvement and application of the method and the possible reasons that affect the accuracy on the exposed side chains are discussed.     

The stability and unfolding of an IgG binding protein based upon the B domain of protein A from Staphylococcus aureus probed by tryptophan substitution and fluorescence spectroscopy

The stability and unfolding of an immunoglobulin (Ig) G bindingprotein based upon the B domain of protein A (SpAB) from Staphylococcusaureus were studied by substituting tryptophan residues at strategiclocations within each of the three a-helical regions (al-a3)of the domain. The role of the C-terminal helix, a3, was investigatedby generating two protein constructs, one corresponding to thecomplete SpAB, the other lacking a part of ct3; the Trp substitutionswere made in both one-and two-domain versions of each of theseconstructs. The fluorescence properties of each of the single-tryptophanmutants were studied in the native state and as a function ofguanidine-HCl-mediated unfolding, and their IgG binding activitieswere determined by a competitive enzyme-linked immunosorbentassay. The free energies of folding and of binding to IgG foreach mutant were compared with those for the native domains.The effect of each substitution upon the overall structure andupon the IgG binding interface was modelled by molecular graphicsand energy minimization. These studies indicate that (i) 3 contributesto the overall stability of the domain and to the formationof the IgG binding site in l and 2, and (ii) al unfolds first,followed by 2 and 3 together.     

Efficient site specific removal of a C-terminal FLAG fusion from a Fab' using copper(II) ion catalysed protein cleavage

The peptide sequence ^NDKTH^C was investigated as a site for efficient,specific cleavage of a fusion protein by cupric ions using ahumanised 1 Fab' as a model protein. The native upper hinge^NDKTH^C sequence was mutated to create a site resistant to cleavageby cupric ions and a ^NDKTH^C sequence introduced between thehinge and a C-terminal FLAG peptide. Incubation of Fab' withCu²⁺ at 62°C at alkaline pHs resulted in removal of theFLAG peptide with efficiencies of up to 86%. Cleavage conditionsdid not detrimentally affect the Fab' protein. Use of the ^NDKTH^Csequence along with cupric ions may provide a cost-effectivemethod for large scale proteolytic cleavage of fusion proteins.     

Highly effective protease inhibitors from variants of human pancreatic secretory trypsin inhibitor (hPSTI): an assessment of 3-D structure-based protein design

The results of a protein design project are used to comparedifferent predictive strategies with respect to proteinproteininteractions. We have been able to generate variants of humanpancreatic secretory trypsin inhibitor (hPSTI) optimized withrespect to the affinity and specificity for human leukocyteelastase relative to trypsin and chymotrypsin, and in particularchymotrypsin. The extremely strong and specific human leukocyteelastase inhibitors were thus developed in three rounds of mutagenesisand two rounds of 3-D modelling; only 24 variants in total weresynthesized, although variations at seven different amino acidpositions were involved (i.e. from 20⁷ possible variants). Anexcellent elastase inhibitor could be designed with the minimumof two amino acid exchanges. The value of structural modellingand actual structure determination is discussed in the lightof the experimental results of the designed protein variantsand the results of tertiary structure determinations of thefree variant and the inhibitorprotease complex. Particular referenceis given to the strategy to be followed in protein design projectsin general and to the development of protease inhibitors inparticular.     

Molecular cartography of proteins: surface relief analysis of the calf eye lens protein gamma-crystallin

Methods of calculating the protein molecular surface and differentmap representations are described. The maps are obtained byprojection of the space-filling molecular model on the surfaceof the ellipsoid of inertia. A new approach to surface analysisis proposed which is based on the use of three general maps:an identification map with all residues outlined, a surfacerelief map and a coloured map with a specific colour for eachof the surface atoms. Superposition of these maps greatly simplifiesmolecular surface analysis. The usefulness of such an approachhas been demonstrated by the study of the relief of the calfeye lens protein -crystallin II. Protrusions of the relief havebeen shown to be occupied generally by charged residues, butin some cases by the hydrophobic ones. It is interesting tonote that in crystal medium the protruding residues are involved,in the majority of cases, in intermolecular contacts. The protrudingregions have been found to be pseudosymmetrical to each otherin accordance with the two-fold rotation axis of the molecule.However, the colours of these regions, i.e. the atoms of thecorresponding side chains, differ greatly.     

Search for the most stable folds of protein chains: II. Computation of stable architectures of {beta}-proteins using a self-consistent molecular field theory

In a preceding paper we presented a novel approach to computationof 3-D folds of protein chains from their amino acid sequences.This approach is a physically correct generalization of the‘threading’ methods. It is based on a self-consistentmolecular field theory and on a physical theory of protein foldingpatterns, which make it possible to examine all the varietyof ‘potentially stable’ folding patterns and allthe variety of the chain conformations within each of them andto determine the thermodynamically stable structure. In thispaper, we apply this approach to single out stable folding patternsand conformations for the chains of ß-sandwich proteinsand show that the similarity of the calculated and observedstructures is usually rather close.     

Free energy perturbation study on a Trp-binding mutant (Ser88 ->Cys) of the trp-repressor

The Ser⁸⁸Cys mutant of the trp-repressor showed a lower affinityfor the corepressor than the wild-type repressor [G = 1.7 ±0.3 kcal/mol, Chou and Matthews (1989) J. Biol. Chem., 264,18314–18319].A molecular dynamics/free energy cycle perturbation study wasperformed to understand the origin of the decreased affinity.A value (G = 1.58 ± 0.28 kcal/mol) comparable with theexperimental value was obtained by the simulation. Free energycomponent analysis revealed that destabilization of the vander Waals interaction between Ser⁸⁸ and Trp¹⁰⁹ (corepressor)mainly contributed to the decreased affinity of the mutant.The rotational transition of the hydroxyl (sulfhydryl) groupof Ser⁸⁸ (Cys⁸⁸) during the simulations affected the contributionsof Arg⁸⁴ and water to the free energy change in the aporepressorand those of Arg⁸⁴ and Trp ¹⁰⁹ to that in the holorepressor.However, the contributions from different residues compensatedeach other, and the total free energy changes were almost invariablein the various simulations.