A phage display system for studying the sequence determinants of protein folding

Phage display and human disease detection

Phage display is a significant and active molecular method and has continued crucial for investigative sector meanwhile its unearthing in 1985. This practice has numerous benefits: the association among physiology and genome, the massive variety of variant proteins showed in sole collection and the elasticity of collection that can be achieved. It suggests a diversity of stages for manipulating antigen attachment; yet, variety and steadiness of exhibited library are an alarm. Additional improvements, like accumulation of non-canonical amino acids, resulting in extension of ligands that can be recognized through collection, will support in expansion of the probable uses and possibilities of technology. Epidemic of COVID-19 had taken countless lives, and while indicative prescriptions were provided to diseased individuals, still no prevention was observed for the contamination. Phage demonstration has presented an in-depth understanding into protein connections included in pathogenesis. Phage display knowledge is developing as an influential, inexpensive, quick, and effectual method to grow novel mediators for the molecular imaging and analysis of cancer.

Direct Observation of the Mechanical Role of Bacterial Chaperones in Protein Folding

Protein folding under force is an integral source of generating mechanical energy in various cellular processes, ranging from protein translation to degradation. Although chaperones are well known to interact with proteins under mechanical force, how they respond to force and control cellular energetics remains unknown. To address this question, we introduce a real-time magnetic tweezer technology herein to mimic the physiological force environment on client proteins, keeping the chaperones unperturbed. We studied two structurally distinct client proteins--protein L and talin with seven different chaperones─independently and in combination and proposed a novel mechanical activity of chaperones. We found that chaperones behave differently, while these client proteins are under force, than their previously known functions. For instance, tunnel-associated chaperones (DsbA and trigger factor), otherwise working as holdase without force, assist folding under force. This process generates an additional mechanical energy up to ∼147 zJ to facilitate translation or translocation. However, well-known cytoplasmic foldase chaperones (PDI, thioredoxin, or DnaKJE) do not possess the mechanical folding ability under force. Notably, the transferring chaperones (DnaK, DnaJ, and SecB) act as holdase and slow down the folding process, both in the presence and absence of force, to prevent misfolding of the client proteins. This provides an emerging insight of mechanical roles of chaperones: they can generate or consume energy by shifting the energy landscape of the client proteins toward a folded or an unfolded state, suggesting an evolutionary mechanism to minimize energy consumption in various biological processes.

Trigger factor chaperone acts as a mechanical foldase

Proteins fold under mechanical forces in a number of biological processes, ranging from muscle contraction to co-translational folding. As force hinders the folding transition, chaperones must play a role in this scenario, although their influence on protein folding under force has not been directly monitored yet. Here, we introduce single-molecule magnetic tweezers to study the folding dynamics of protein L in presence of the prototypical molecular chaperone trigger factor over the range of physiological forces (4–10 pN). Our results show that trigger factor increases prominently the probability of folding against force and accelerates the refolding kinetics. Moreover, we find that trigger factor catalyzes the folding reaction in a force-dependent manner; as the force increases, higher concentrations of trigger factor are needed to rescue folding. We propose that chaperones such as trigger factor can work as foldases under force, a mechanism which could be of relevance for several physiological processes.

Proteins fold under mechanical force during co-translational folding at the ribosome. Here, the authors use single molecule magnetic tweezers to study the influence of chaperones on protein folding and show that the ribosomal chaperone trigger factor acts as a mechanical foldase by promoting protein folding under force.

Centenary Award and Sir Frederick Gowland Hopkins Memorial Lecture. Protein folding, structure prediction and design

I describe how experimental studies of protein folding have led to advances in protein structure prediction and protein design. I describe the finding that protein sequences are not optimized for rapid folding, the contact order-protein folding rate correlation, the incorporation of experimental insights into protein folding into the Rosetta protein structure production methodology and the use of this methodology to determine structures from sparse experimental data. I then describe the inverse problem (protein design) and give an overview of recent work on designing proteins with new structures and functions. I also describe the contributions of the general public to these efforts through the Rosetta@home distributed computing project and the FoldIt interactive protein folding and design game.

Evolution of the potential energy landscape with static pulling force for two model proteins

The energy landscape is analyzed for off-lattice bead models of protein L and protein G as a function of a static pulling force. Two different pairs of attachment points (pulling directions) are compared in each case, namely, residues 1/56 and 10/32. For the terminal residue pulling direction 1/56, the distinct global minimum structures are all extended, aside from the compact geometry that correlates with zero force. The helical turns finally disappear at the highest pulling forces considered. For the 10/32 pulling direction, the changes are more complicated, with a variety of competing arrangements for beads outside the region where the force is directly applied. These alternatives produce frustrated energy landscapes, with low-lying minima separated by high barriers. The calculated folding pathways in the absence of force are in good agreement with previous work. The N-terminal hairpin folds first for protein L and the C-terminal hairpin for protein G, which exhibits an intermediate. However, for a relatively low static force, where the global minimum retains its structure, the folding mechanisms change, sometimes dramatically, depending on the protein and the attachment points. The scaling relations predicted by catastrophe theory are found to hold in the limit of short path lengths.

Small-angle X-ray scattering and single-molecule FRET spectroscopy produce highly divergent views of the low-denaturant unfolded state

The results of more than a dozen single-molecule Förster resonance energy transfer (smFRET) experiments suggest that chemically unfolded polypeptides invariably collapse from an expanded random coil to more compact dimensions as the denaturant concentration is reduced. In sharp contrast, small-angle X-ray scattering (SAXS) studies suggest that, at least for single-domain proteins at non-zero denaturant concentrations, such compaction may be rare. Here, we explore this discrepancy by studying protein L, a protein previously studied by SAXS (at 5 °C), which suggested fixed unfolded-state dimensions from 1.4 to 5 M guanidine hydrochloride (GuHCl), and by smFRET (at 25 °C), which suggested that, in contrast, the chain contracts by 15-30% over this same denaturant range. Repeating the earlier SAXS study under the same conditions employed in the smFRET studies, we observe little, if any, evidence that the unfolded state of protein L contracts as the concentration of GuHCl is reduced. For example, scattering profiles (and thus the shape and dimensions) collected within ∼4 ms after dilution to as low as 0.67 M GuHCl are effectively indistinguishable from those observed at equilibrium at higher denaturant. Our results thus argue that the disagreement between SAXS and smFRET is statistically significant and that the experimental evidence in favor of obligate polypeptide collapse at low denaturant cannot be considered conclusive yet.

IgG-binding proteins of bacteria

Proteins capable of non-immune binding of immunoglobulins G (IgG) of various mammalian species, i.e. without the involvement of the antigen-binding sites of the immunoglobulins, are widespread in bacteria. These proteins are located on the surface of bacterial cells and help them to evade the host's immune response due to protection against the action of complement and to decrease in phagocytosis. This review summarizes data on the structure of immunoglobulin-binding proteins (IBP) and their complexes with IgG. Common and distinctive structural features of IBPs of gram-positive bacteria (staphylococci, streptococci, peptostreptococci) are discussed. Conditions for IBP expression by bacteria and their functional heterogeneity are considered. Data on IBPs of gram-negative bacteria are presented.

Identification of a mechanical rheostat in the hydrophobic core of protein L

The ability of proteins and their complexes to withstand or respond to mechanical stimuli is vital for cells to maintain their structural organisation, to relay external signals and to facilitate unfolding and remodelling. Force spectroscopy using the atomic force microscope allows the behaviour of single protein molecules under an applied extension to be investigated and their mechanical strength to be quantified. protein L, a simple model protein, displays moderate mechanical strength and is thought to unfold by the shearing of two mechanical sub-domains. Here, we investigate the importance of side-chain packing for the mechanical strength of protein L by measuring the mechanical strength of a series of protein L variants containing single conservative hydrophobic volume deletion mutants. Of the five thermodynamically destabilised variants characterised, only one residue (I60V) close to the interface between two mechanical sub-domains was found to differ in mechanical properties to wild type (ΔF_I60V–WT = − 36 pN at 447 nm s^− 1, Δx_uI60V–WT = 0.2 nm). Φ-value analysis of the unfolding data revealed a highly native transition state. To test whether the number of hydrophobic contacts across the mechanical interface does affect the mechanical strength of protein L, we measured the mechanical properties of two further variants. protein L L10F, which increases core packing but does not enhance interfacial contacts, increased mechanical strength by 13 ± 11 pN at 447 nm s^− 1. By contrast, protein L I60F, which increases both core and cross-interface contacts, increased mechanical strength by 72 ± 13 pN at 447 nm s^− 1. These data suggest a method by which nature can evolve a varied mechanical response from a limited number of topologies and demonstrate a generic but facile method by which the mechanical strength of proteins can be rationally modified.

Biomimicry as a route to new materials: what kinds of lessons are useful?

We consider the attributes of a successful engineered material, acknowledging the contributions of composition and processing to properties and performance. We recognize the potential for relevant lessons to be learned from nature, at the same time conceding both the limitations of such lessons and our need to be selective. We then give some detailed attention to the molecular biomimicry of filamentous phage, the process biomimicry of silk and the structure biomimicry of hippopotamus 'sweat', in each case noting that the type of lesson now being learned is not the same as the potential lesson that originally motivated the study.

The use of phage display in neurobiology

Phage display is a technique that involves the coupling of phenotype to genotype in a selectable format. It has been extensively used in molecular biology to study protein-protein interactions, receptor and antibody binding sites, and immune responses; to modify protein properties; and to select antibodies against a wide range of different antigens. In the format most often used, a polypeptide is displayed on the surface of a filamentous phage by genetic fusion to one of the coat proteins, creating a chimeric coat protein. As the gene encoding the chimeric coat protein is packaged within the phage, selection of the phage on the basis of the binding properties of the polypeptide displayed on the surface simultaneously results in the isolation of the gene encoding the polypeptide. This unit describes the background of the technique and illustrates how it has been applied to a number of different problems, each of which has its neurobiological counterparts.

Exploiting Elements of Transcriptional Machinery to Enhance Protein Stability

The correlation between protein structure and function is well established, yet the role stability/flexibility plays in protein function is being explored. Here, we describe an in vivo screen in which the thermal stability of a test protein is correlated directly to the transcriptional regulation of a reporter gene. The screen readout is independent of the function of the test protein, proteolytic resistance, solubility or propensity to aggregate indiscriminately, and is thus dependent solely on the overall stability of the test protein. The system entails the use of an engineered chimeric construct that consists of three covalently linked domains; a constant N-terminal DNA-binding domain, a variable central test protein, and a constant C-terminal transcriptional activation domain. The test proteins are mutant variants of the beta1 domain of streptococcal protein G that span fairly evenly a thermal stability range from as low as 38 degrees C to above 100 degrees C. When the chimeric construct contains a test variant of low thermal stability, the reporter gene is up-regulated to a greater extent relative to that of more stable/less flexible variants. A panel of nine Gbeta1 mutant variants was used to benchmark the screen, and spectroscopic methods were employed to characterize the thermal and structural properties of each variant accurately. The screen was combined with in silico methods to interrogate a library of randomized variants for selection of mutants of greater structural integrity.

The Highly Cooperative Folding of Small Naturally Occurring Proteins Is Likely the Result of Natural Selection

To illuminate the evolutionary pressure acting on the folding free energy landscapes of naturally occurring proteins, we have systematically characterized the folding free energy landscape of Top7, a computationally designed protein lacking an evolutionary history. Stopped-flow kinetics, circular dichroism, and NMR experiments reveal that there are at least three distinct phases in the folding of Top7, that a nonnative conformation is stable at equilibrium, and that multiple fragments of Top7 are stable in isolation. These results indicate that the folding of Top7 is significantly less cooperative than the folding of similarly sized naturally occurring proteins, suggesting that the cooperative folding and smooth free energy landscapes observed for small naturally occurring proteins are not general properties of polypeptide chains that fold to unique stable structures but are instead a product of natural selection.

Kinetic studies of protein L aggregation and disaggregation

We have investigated the aggregation of protein L in 25% (vol/vol) TFE and 10 mM HCl. Under both conditions, aggregates adopt a fibrillar structure and bind dyes Congo Red and Thioflavin T consistent with the presence of amyloid fibrils. The kinetics of aggregation in 25% TFE suggest a linear-elongation mechanism with critical nucleus size of either two or three monomers. Aggregation kinetics in 10 mM HCl show a prolonged lag phase prior to a rapid increase in aggregation. The lag phase is time-dependent, but the time dependence can be eliminated by the addition of pre-formed seeds. Disaggregation studies show that for aggregates formed in TFE, aggregate stability is a strong function of aggregate age. For example, after 200 min of aggregation, 40% of the aggregation reaction is irreversible, while after 3 days over 60% is irreversible. When the final concentration of the denaturant, TFE, is reduced from 5% to 0, the amount of reversible aggregation doubles. Disaggregation studies of aggregates formed in TFE and 10 mM HCl reveal a complicated effect of pH on aggregate stability.

Coil-globule transition in the denatured state of a small protein

Upon transfer from strongly denaturing to native conditions, proteins undergo a collapse that either precedes folding or occurs simultaneously with it. This collapse is similar to the well known coil-globule transition of polymers. Here we employ single-molecule fluorescence methods to fully characterize the equilibrium coil-globule transition in the denatured state of the IgG-binding domain of protein L. By using FRET measurements on freely diffusing individual molecules, we determine the radius of gyration of the protein, which shows a gradual expansion as the concentration of the denaturant, guanidinium hydrochloride, is increased all the way up to 7 M. This expansion is observed also in fluorescence correlation spectroscopy measurements of the hydrodynamic radius of the protein. We analyze the radius of gyration measurements using the theory of the coil-globule transition of Sanchez [Sanchez, I. C. (1979) Macromolecules 12, 980-988], which balances the excluded volume entropy of the chain with the average interresidue interaction energy. In particular, we calculate the solvation energy of the denatured protein, a property that is not readily accessible in other experiments. The dependence of this energy on denaturant concentration is nonlinear, contrasting with the common linear extrapolation method used to describe denaturation energy. Interestingly, a fit to the binding model of chemical denaturation suggests a single denaturant binding site per protein residue. The size of the denatured protein under native conditions can be extrapolated from the data as well, showing that the fully collapsed state of protein is only approximately 10% larger than the folded state.

Total chemical synthesis of the B1 domain of protein L from Peptostreptococcus magnus

Reported here is a native chemical ligation strategy for the total chemical synthesis of the B1 domain of protein L. A synthetic construct of this 76 amino acid protein domain was prepared by the chemoselective ligation of two unprotected polypeptide fragments, one containing an N-terminal cysteine residue and one containing a C-terminal thioester moiety. The polypeptide fragments utilized in the ligation reaction were readily prepared by stepwise solid phase peptide synthesis (SPPS) methods for Boc-chemistry. The milligram quantities of protein required for conventional biophysical studies were readily accessible using the synthetic protocol described here. The folding properties of the synthetic protein L construct were also determined and found to be very similar to those of a similar wild-type protein L constructs prepared by recombinant-DNA methods. This work facilitates future unnatural amino acid mutagenesis experiments on this model protein system to further dissect the molecular basis of its folding and stability.

Phage display systems and their applications

Screening phage display libraries of proteins and peptides has, for almost two decades, proven to be a powerful technology for selecting polypeptides with desired biological and physicochemical properties from huge molecular libraries. The scope of phage display applications continues to expand. Recent applications and technical improvements driving further developments in the field of phage display are discussed.

Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.

Intermediates and the folding of proteins L and G

We use a minimalist protein model, in combination with a sequence design strategy, to determine differences in primary structure for proteins L and G, which are responsible for the two proteins folding through distinctly different folding mechanisms. We find that the folding of proteins L and G are consistent with a nucleation-condensation mechanism, each of which is described as helix-assisted beta-1 and beta-2 hairpin formation, respectively. We determine that the model for protein G exhibits an early intermediate that precedes the rate-limiting barrier of folding, and which draws together misaligned secondary structure elements that are stabilized by hydrophobic core contacts involving the third beta-strand, and presages the later transition state in which the correct strand alignment of these same secondary structure elements is restored. Finally, the validity of the targeted intermediate ensemble for protein G was analyzed by fitting the kinetic data to a two-step first-order reversible reaction, proving that protein G folding involves an on-pathway early intermediate, and should be populated and therefore observable by experiment.

Searching for folded proteins in vitro and in silico

Understanding the sequence determinants of protein structure, stability and folding is critical for understanding how natural proteins have evolved and how proteins can be engineered to perform novel functions. The complexity of the protein folding problem requires the ability to search large volumes of sequence space for proteins with specific structural or functional characteristics. Here we describe our efforts to identify novel proteins using a phage-display selection strategy from a 'mini-exon' shuffling library generated from the yeast genome and from completely random sequence libraries, and compare the results to recent successes in generating novel proteins using in silico protein design.

Phage-display as a tool for quantifying protein stability determinants

To address questions of protein stability, researchers have increasingly turned to combinatorial approaches that permit the rapid analysis of libraries of protein variants. Phage-display has proved to be a powerful tool for analyzing protein stability due to the large library size and the robustness of the phage particle to a variety of denaturing conditions. With the B1 domain of protein G (GB1) and a camelid heavy chain antibody as model systems, we are using phage-display libraries to experimentally address questions that have generally been addressed in silico, either through computational studies or statistical analysis of known protein structures. One effort has focused on identifying novel solutions to repacking the hydrophobic core of GB1, while maintaining stability comparable to the wild type protein. In a second study, a small set of substitutions in complimentarity-determining region 3 was found to stabilize the framework of the camelid antibody. Another major focus has been to obtain quantitative data on beta-sheet stability determinants. We have successfully adapted a phage-display method for quantitating affinities of protein variants (shotgun alanine scanning) to analysis of GB1 stability. Using this method, we have analyzed the energetic contributions of cross-strand side chain-side chain interactions. Finally, we discuss parameters to consider in using phage-display to discriminate subtle stability differences among fully folded variants. Overall, this method provides a fast approach for quantitatively addressing biophysical questions.

Combinatorial approaches to protein stability and structure

Why do proteins adopt the conformations that they do, and what determines their stabilities? While we have come to some understanding of the forces that underlie protein architecture, a precise, predictive, physicochemical explanation is still elusive. Two obstacles to addressing these questions are the unfathomable vastness of protein sequence space, and the difficulty in making direct physical measurements on large numbers of protein variants. Here, we review combinatorial methods that have been applied to problems in protein biophysics over the last 15 years. The effects of hydrophobic core composition, the most important determinant of structure and stability, are still poorly understood. Particular attention is given to core composition as addressed by library methods. Increasingly useful screens and selections, in combination with modern high-throughput approaches borrowed from genomics and proteomics efforts, are making the empirical, statistical correlation between sequence and structure a tractable problem for the coming years.

Dual surface selection methodology for the identification of thrombin binding epitopes from hotspot biased phage-display libraries

Protein libraries biased towards amino-acid residues found at so-called 'hotspots' were incorporated into the beta-sheet region of the thermostable variant (HTB1) of the B1 domain of the immunoglobulin (IgG) binding protein G and expressed as gene 3 fusions on M13 bacteriophage. The HTB1 library (2.2 x 10(9)) variants with a minimal 12 amino acid basis set were selected for binding IgG, to ensure structural conservation, and subsequently to thrombin to evolve a thrombin-binding function. We believe that this dual surface selection strategy will have great utility in evolving new bi-functional proteins without compromising structure. Furthermore the discrete beta-sheet epitopes identified by our methodology will lend itself to small-molecule mimicry of beta-sheets.

Functional epitope on human neutrophil flavocytochrome b558

mAb NL7 was raised against purified flavocytochrome b(558), important in host defense and inflammation. NL7 recognized the gp91(phox) flavocytochrome b(558) subunit by immunoblot and bound to permeabilized neutrophils and neutrophil membranes. Epitope mapping by phage display analysis indicated that NL7 binds the (498)EKDVITGLK(506) region of gp91(phox). In a cell-free assay, NL7 inhibited in vitro activation of the NADPH oxidase in a concentration-dependent manner, and had marginal effects on the oxidase substrate Michaelis constant (K(m)). mAb NL7 did not inhibit translocation of p47(phox), p67(phox), or Rac to the plasma membrane, and bound its epitope on gp91(phox) independently of cytosolic factor translocation. However, after assembly of the NADPH oxidase complex, mAb NL7 bound the epitope but did not inhibit the generation of superoxide. Three-dimensional modeling of the C-terminal domain of gp91(phox) on a corn nitrate reductase template suggests close proximity of the NL7 epitope to the proposed NADPH binding site, but significant separation from the proposed p47(phox) binding sites. We conclude that the (498)EKDVITGLK(506) segment resides on the cytosolic surface of gp91(phox) and represents a region important for oxidase function, but not substrate or cytosolic component binding.

Low free energy cost of very long loop insertions in proteins

Long insertions into a loop of a folded host protein are expected to have destabilizing effects because of the entropic cost associated with loop closure unless the inserted sequence adopts a folded structure with amino- and carboxy-termini in close proximity. A loop entropy reduction screen based on this concept was used in an attempt to retrieve folded sequences from random sequence libraries. A library of long random sequences was inserted into a loop of the SH2 domain, displayed on the surface of M13 phage, and the inserted sequences that did not disrupt SH2 function were retrieved by panning using beads coated with a phosphotyrosine containing SH2 peptide ligand. Two sequences of a library of 2 x 10(8) sequences were isolated after multiple rounds of panning, and were found to have recovery levels similar to the wild-type SH2 domain and to be relatively intolerant to further mutation in PCR mutagenesis experiments. Surprisingly, although these inserted sequences exhibited little nonrandom structure, they do not significantly destabilize the host SH2 domain. Additional insertion variants recovered at lower levels in the panning experiments were also found to have a minimal effect on the stability and peptide-binding function of the SH2 domain. The additional level of selection present in the panning experiments is likely to involve in vivo folding and assembly, as there was a rough correlation between recovery levels in the phage-panning experiments and protein solubility. The finding that loop insertions of 60-80 amino acids have minimal effects on SH2 domain stability suggests that the free energy cost of inserting long loops may be considerably less than polymer theory estimates based on the entropic cost of loop closure, and, hence, that loop insertion may have provided an evolutionary route to multidomain protein structures.

Designing gene libraries from protein profiles for combinatorial protein experiments

Protein combinatorial libraries provide new ways to probe the determinants of folding and to discover novel proteins. Such libraries are often constructed by expressing an ensemble of partially random gene sequences. Given the intractably large number of possible sequences, some limitation on diversity must be imposed. A non-uniform distribution of nucleotides can be used to reduce the number of possible sequences and encode peptide sequences having a predetermined set of amino acid probabilities at each residue position, i.e., the amino acid sequence profile. Such profiles can be determined by inspection, multiple sequence alignment or physically-based computational methods. Here we present a computational method that takes as input a desired sequence profile and calculates the individual nucleotide probabilities among partially random genes. The calculated gene library can be readily used in the context of standard DNA synthesis to generate a protein library with essentially the desired profile. The fidelity between the desired profile and the calculated one coded by these partially random genes is quantitatively evaluated using the linear correlation coefficient and a relative entropy, each of which provides a measure of profile agreement at each position of the sequence. On average, this method of identifying such codon frequencies performs as well or better than other methods with regard to fidelity to the original profile. Importantly, the method presented here provides much better yields of complete sequences that do not contain stop codons, a feature that is particularly important when all or large fractions of a gene are subject to combinatorial mutation.

Directed evolution of barnase stability using proteolytic selection

We report the construction of a phage-displayed repertoire of mutants of the ribonuclease barnase from Bacillus amyloliquefaciens. The construction was guided by the natural variability between two closely related ribonucleases, barnase and binase from Bacillus intermedius. This repertoire was selected using a proteolytic selection method, allowing sorting of the library according to the resistance of the mutants toward proteolysis. Susceptibility toward proteolysis has been correlated with flexibility and unfolding, and is thus expected to yield mutants with increased thermal stability. Enrichment of barnase mutants with specific combinations of amino acid residues at four of the randomised positions was observed. Three of these enriched amino acid residues are present in neither barnase nor binase. For some of the mutations, the improvement in proteolytic stability does not lead to a pronounced improvement in thermodynamic stability, indicating that the factors governing the proteolytic stability in some cases may be different from those governing the thermodynamic stability, e.g. propensity to local unfolding.The results obtained add important knowledge to a novel use of phage display technology for selection of thermodynamically stable proteins. Only by carefully establishing the parameters that can be adjusted, and recognising the influence this will have on the outcome of selection, will it be possible to realise the powerful technique of proteolytic selection.

Accurate computer-based design of a new backbone conformation in the second turn of protein L

The rational design of loops and turns is a key step towards creating proteins with new functions. We used a computational design procedure to create new backbone conformations in the second turn of protein L. The Protein Data Bank was searched for alternative turn conformations, and sequences optimal for these turns in the context of protein L were identified using a Monte Carlo search procedure and an energy function that favors close packing. Two variants containing 12 and 14 mutations were found to be as stable as wild-type protein L. The crystal structure of one of the variants has been solved at a resolution of 1.9 A, and the backbone conformation in the second turn is remarkably close to that of the in silico model (1.1 A RMSD) while it differs significantly from that of wild-type protein L (the turn residues are displaced by an average of 7.2 A). The folding rates of the redesigned proteins are greater than that of the wild-type protein and in contrast to wild-type protein L the second beta-turn appears to be formed at the rate limiting step in folding.

Genomic biodiversity, phylogenetics and coevolution in proteins

Comprehensive sampling of genomic biodiversity is fast becoming a reality for some genomic regions and complete organelle genomes. Genomic biodiversity is defined as large genomic sequences from many species, and here some recent work is reviewed that demonstrates the potential benefits of genomic biodiversity for molecular evolutionary analysis and phylogenetic reconstruction. This work shows that using likelihood-based approaches, taxon addition can dramatically improve phylogenetic reconstruction. Features or dynamics of the evolutionary process are much more easily inferred with large numbers of taxa, and large numbers are essential for discriminating differences in evolutionary patterns between sites. Accurate prediction of site-specific patterns can improve phylogenetic reconstruction by an amount equivalent to quadrupling sequence length. Genomic biodiversity is particularly central to research relating patterns of evolution, adaptation and coevolution to structural and functional features of proteins. Research on detecting coevolution between amino acid residues in proteins demonstrates a clear need for much greater numbers of closely related taxa to better discriminate site-specific patterns of interaction, and to allow more detailed analysis of coevolutionary interactions between subunits in protein complexes. It is argued that parsing out coevolutionary and other context-dependent substitution probabilities is essential for discriminating between coevolution and adaptation, and for more realistically modelling the evolution of proteins. Also reviewed is research that argues for increasing the efficiency of acquiring genomic biodiversity, and suggests that this might be done by simultaneously shotgun cloning and sequencing genomic mixtures from many species. Increased efficiency is a prerequisite if genomic biodiversity levels are to rapidly increase by orders of magnitude, and thus lead to dramatically improved understanding of interactions between protein structure, function and sequence evolution.

A substrate phage enzyme-linked immunosorbent assay to profile panels of proteases

It is estimated that proteases comprise nearly 2% of the human genome. Given that the primary structure of all known proteases will soon be available, an important challenge is to define the structure-activity relationships that govern substrate hydrolysis. Ideally this would be accomplished on a genome-wide scale. To this end, we have developed a one-pot phage selection system that yields the substrate recognition profile of multiple proteases from a single round of selection. The system meets five key criteria: (i) multiple proteases can be analyzed simultaneously, (ii) prior knowledge of substrate preference is not required, (iii) information regarding substrate preferences on both side of the scissile bond is obtained, (iv) the system yields selective substrates that distinguish closely related proteases, and (v) semiquantitative information on substrate hydrolysis is obtained, allowing for the assignment of initial rank-order preferences. As an illustration, a phage selection with a mixture of thrombin and factor Xa (serine proteases) along with matrix-metalloproteinase-9 and atrolysin C (metalloproteinases) was performed. Peptide substrates were identified that (i) have high k(cat)/K(m) ratios, (ii) are selective for individual proteases, and (iii) match the sequences of known physiological substrates.

A high-throughput screen for identifying transmembrane pore-forming peptides

We have developed a visual microwell plate assay for rapid, high-throughput screening for membrane-disrupting molecules such as de novo designed pore formers, antibiotic peptides, bacterial toxins, and lipases. The detectability is based on the strong fluorescence emission of the lanthanide metal terbium(III) (Tb(3+)) when it interacts with the aromatic chelator dipicolinic acid (DPA). While Tb(3+) is not strongly fluorescent alone, the binary complex emits bright green fluorescence when irradiated with uv light. For the microwell plate assay, we prepared unilamellar phospholipid vesicles that had either Tb(3+) or DPA entrapped and the opposite molecule in the external solution. Disruption of the membranes allows the Tb(3+)/DPA complex to form, giving rise to a visibly fluorescent solution. In plates with 20-microl wells, the lower limit of visual detectability of the Tb(3+)/DPA complex in solution was about 2.5 microM. The lower limit of detectability using vesicles with entrapped Tb(3+) or DPA was about 50 microM phospholipid. We show that the membrane-disrupting effect of as little as 0.25 microM or 5 pmol of the pore-forming, antibiotic peptide alamethicin can be detected visually with this system. This sensitive, high-throughput assay is readily automatable and makes possible the visual screening of combinatorial peptide libraries for members that permeabilize lipid bilayer membranes.

Mutational analysis of protein substrate presentation in the post-translational attachment of biotin to biotin domains

Biotinylation in vivo is an extremely selective post-translational event where the enzyme biotin protein ligase (BPL) catalyzes the covalent attachment of biotin to one specific and conserved lysine residue of biotin-dependent enzymes. The biotin-accepting lysine, present in a conserved Met-Lys-Met motif, resides in a structured domain that functions as the BPL substrate. We have employed phage display coupled with a genetic selection to identify determinants of the biotin domain (yPC-104) of yeast pyruvate carboxylase 1 (residues 1075-1178) required for interaction with BPL. Mutants isolated using this strategy were analyzed by in vivo biotinylation assays performed at both 30 degrees C and 37 degrees C. The temperature-sensitive substrates were reasoned to have structural mutations, leading to compromised conformations at the higher temperature. This interpretation was supplemented by molecular modeling of yPC-104, since these mutants mapped to residues involved in defining the structure of the biotin domain. In contrast, substitution of the Met residue N-terminal to the target lysine with either Val or Thr produced mutations that were temperature-insensitive in the in vivo assay. Furthermore, these two mutant proteins and wild-type yPC-104 showed identical susceptibility to trypsin, consistent with these substitutions having no structural effect. Kinetic analysis of enzymatic biotinylation using purified Met --> Thr/Val mutant proteins with both yeast and Escherichia coli BPLs revealed that these substitutions had a strong effect upon K(m) values but not k(cat). The Met --> Thr mutant was a poor substrate for both BPLs, whereas the Met --> Val substitution was a poor substrate for bacterial BPL but had only a 2-fold lower affinity for yeast BPL than the wild-type peptide. Our data suggest that substitution of Thr or Val for the Met N-terminal of the biotinyl-Lys results in mutants specifically compromised in their interaction with BPL.

Studies on a single immunoglobulin-binding domain of protein L from Peptostreptococcus magnus: the role of tyrosine-53 in the reaction with human IgG

Chemical modification experiments with tetranitromethane (TNM) have been used to investigate the role of tyrosine residues in the formation of the complex between PpL (the single Ig-binding domain of protein L, isolated from P. magnus strain 3316) and the kappa light chain (kappa-chain). Reaction of PpL with TNM causes the modification of 1.9 equiv. of tyrosine (Tyr(51) and Tyr(53)) and results in an approx. 140-fold decrease in affinity for human IgG. Similar experiments with mutated PpL proteins suggest that nitration predominantly inactivates the protein by modification of Tyr(53). Reduction of the nitrotyrosine groups to aminotyrosine by incubation with sodium hydrosulphite does not restore high affinity for IgG. Modification of kappa-chain by TNM resulted in the nitration of 3.1+/-0.09 tyrosine residues. When the PpL-kappa-chain complex was incubated with TNM, 4.1+/-0.04 tyrosine residues were nitrated, indicating that one tyrosine residue previously modified by the reagent was protected from TNM when the proteins are in complex with each other. The K(d) for the equilibrium between PpL, human IgG and their complex has been shown by ELISA to be 112+/-20 nM. A similar value (153+/-33 nM) was obtained for the complex formed between IgG and the Tyr(64)-->Trp mutant (Y64W). However, the K(d) values for the equilibria involving the PpL mutants Y53F and Y53F,Y64W were found to be 3.2+/-0.2 and 4.6+/-1 microM respectively. These suggest that the phenol group of Tyr(53) in PpL is important to the stability of the PpL-kappa-chain complex.

Hydrophobic core malleability of a de novo designed three-helix bundle protein

De novo protein design provides a tool for testing the principles that stabilize the structures of proteins. Recently, we described the design and structure determination of alpha(3)D, a three-helix bundle protein with a well-packed hydrophobic core. Here, we test the malleability and adaptability of this protein's structure by mutating a small, Ala residue (A60) in its core to larger, hydrophobic side-chains, Leu and Ile. Such changes introduce strain into the structures of natural proteins, and therefore generally destabilize the native state. By contrast, these mutations were slightly stabilizing ( approximately 1.5 kcal mol(-1)) to the tertiary structure of alpha(3)D. The value of DeltaC(p) for unfolding of these mutants was not greatly affected relative to wild-type, indicating that the change in solvent accessibility for unfolding was similar. However, two-dimensional heteronuclear single quantum coherence spectra indicate that the protein adjusts to the introduction of steric bulk in different ways. A60L-alpha(3)D showed serious erosion in the dispersion of both the amide backbone as well as the side-chain methyl chemical shifts. By contrast, A60I-alpha(3)D showed excellent dispersion of the backbone resonances, and selective changes in dispersion of the aliphatic side-chains proximal to the site of mutation. Together, these data suggest that alpha(3)D, although folded into a unique three-dimensional structure, is nevertheless more malleable and flexible than most natural, native proteins.

A "loop entropy reduction" phage-display selection for folded amino acid sequences

As a step toward selecting folded proteins from libraries of randomized sequences, we have designed a 'loop entropy reduction'-based phage-display method. The basic premise is that insertion of a long disordered sequence into a loop of a host protein will substantially destabilize the host because of the entropic cost of closing a loop in a disordered chain. If the inserted sequence spontaneously folds into a stable structure with the N and C termini close in space, however, this entropic cost is diminished. The host protein function can, therefore, be used to select folded inserted sequences without relying on specific properties of the inserted sequence. This principle is tested using the IgG binding domain of protein L and the lck SH2 domain as host proteins. The results indicate that the loop entropy reduction screen is capable of discriminating folded from unfolded sequences when the proper host protein and insertion point are chosen.

NMR characterization of residual structure in the denatured state of protein L

Triple-resonance NMR experiments were used to assign the (13)C(alpha), (13)C(beta), (15)N and NH resonances for all the residues in the denatured state of a destabilized protein L variant in 2 M guanidine. The chemical shifts of most resonances were very close to their random coil values. Significant deviations were observed for G22, L38 and K39; increasing the denaturant concentration shifted the chemical shifts of these residues towards theory random coil values. Medium-range nuclear Overhauser enhancements were detected in segments corresponding to the turn between the first two strands, the end of the second strand through the turn between the second strand and the helix, and the turn between the helix and the third strand in 3D H(1), N(15)-HSQC-NOESY-HSQC experiments on perdeuterated samples. Longer-range interactions were probed by measuring the paramagnetic relaxation enhancement produced by nitroxide spin labels introduced via cysteine residues at five sites around the molecule. Damped oscillations in the magnitude of the paramagnetic relaxation enhancement as a function of distance along the sequence suggested native-like chain reversals in the same three turn regions. The more extensive interactions within the region corresponding to the first beta-turn than in the region corresponding to the second beta-turn suggests that the asymmetry in the folding reaction evident in previous studies of the protein L folding transition state is already established in the denatured state.

A breakdown of symmetry in the folding transition state of protein L

The 62 residue IgG binding domain of protein L consists of a central alpha-helix packed on a four-stranded beta-sheet formed by N and C-terminal beta-hairpins. The overall topology of the protein is quite symmetric: the beta-hairpins have similar lengths and make very similar interactions with the central helix. Characterization of the effects of 70 point mutations distributed throughout the protein on the kinetics of folding and unfolding reveals that this symmetry is completely broken during folding; the first beta-hairpin is largely structured while the second beta-hairpin and helix are largely disrupted in the folding transition state ensemble. The results are not consistent with a "hydrophobic core first" picture of protein folding; the first beta-hairpin appears to be at least as ordered at the rate limiting step in folding as the hydrophobic core.

Robustness of protein folding kinetics to surface hydrophobic substitutions

We use both combinatorial and site-directed mutagenesis to explore the consequences of surface hydrophobic substitutions for the folding of two small single domain proteins, the src SH3 domain, and the IgG binding domain of Peptostreptococcal protein L. We find that in almost every case, destabilizing surface hydrophobic substitutions have much larger effects on the rate of unfolding than on the rate of folding, suggesting that nonnative hydrophobic interactions do not significantly interfere with the rate of core assembly.

Selection for improved protein stability by phage display

A library of mutants of a single-chain Fv fragment (scFv) was generated by a combination of directed and random mutagenesis, using oligonucleotides randomized at defined positions and two rounds of DNA shuffling. The library was based on the already well folding and stable scFv fragment 4D5Flu. In order to further improve this framework and test the efficiency of various selection strategies, phage display selection was carried out under different selective pressures for higher thermodynamic stability. Incubation of the display phages at elevated temperatures was compared to exposure of the phages to high concentrations of guanidinium chloride. Temperature stress-guided selection yielded the most stable scFv mutant after two rounds of mutagenesis and selection, due to the irreversibility of the unfolding process. It possessed only two mutations (His(L27d)Asn and Phe(L55)Val) and showed a thermodynamic stability improved by roughly 4 kcal/mol, threefold better expression yields in Escherichia coli as well as a 20-fold better binding constant than the 4D5Flu wild-type. The selection results obtained in this study delineate the advantages, disadvantages and limitations of different stability stress selection methods in phage display.

The single helix in protein L is largely disrupted at the rate-limiting step in folding

To investigate the role of helix formation in the folding of protein L, a 62 residue alpha/beta protein, we studied the consequences of both single and multiple mutations in the helix on the kinetics of folding. A triple mutant with 11 additional carbon atoms in core residues in the amino-terminal portion of the helix folded substantially faster than wild type, suggesting that hydrophobic association with residues elsewhere in the protein occurs at the rate-limiting step in folding. However, helix-destabilizing mutations had little effect on the rate of folding; in particular, a triple glycine substitution on the solvent-exposed side of the helix increased the unfolding rate 56-fold while reducing the folding rate less than threefold. Thus, in contrast to the predictions of models of folding involving the coalescence of well-formed secondary structure elements, the single helix in protein L appears to be largely disrupted at the rate-limiting step in folding and unfolding.

Proteolytic selection for protein folding using filamentous bacteriophages

BACKGROUND

Filamentous bacteriophages have been used for the selection of folded peptide and protein 'ligands' by binding the phage to 'receptor'-coated solid phase. Here, using proteolysis, we have developed a technique for the selection of folded and stable proteins that is independent of their binding activities.

RESULTS

When a 21-residue peptide comprising a protease cleavage site was introduced into the flexible linker between the second and third domains of the minor coat protein p3 of filamentous bacteriophage, the phages could be cleaved by trypsin and were rendered non-infective. By contrast, phages displaying mutant barnases at this site were resistant to proteolysis, but were cleaved and their infectivity was destroyed as the temperature was raised. By mixing phages bearing two barnase mutants of differing stability, and adding protease at a temperature at which one mutant was resistant and the other was sensitive, we were able to enrich by 1.6 x 10(4)-fold for phages bearing the more stable barnase.

CONCLUSIONS

The approach provides a means for the selection of folded and stable proteins, and may be applicable to the selection of de novo proteins.

Stabilizing the subtilisin BPN' pro-domain by phage display selection: how restrictive is the amino acid code for maximum protein stability?

We have devised a procedure using monovalent phage display to select for stable mutants in the pro-domain of the serine protease, subtilisin BPN'. In complex with subtilisin, the pro-domain assumes a compact structure with a four-stranded antiparallel beta-sheet and two three-turn alpha-helices. When isolated, however, the pro-domain is 97% unfolded. These experiments use combinatorial mutagenesis to select for stabilizing amino acid combinations at a particular structural locus and determine how many combinations are close to the maximum protein stability. The selection for stability is based on the fact that the independent stability of the pro-domain is very low and that binding to subtilisin is thermodynamically linked to folding. Two libraries of mutant pro-domains were constructed and analyzed to determine how many combinations of amino acids at a particular structural locus result in the maximum stability. A library comprises all combinations of four amino acids at a structural locus. Previous studies using combinatorial genetics have shown that many different combinations of amino acids can be accommodated in a selected locus without destroying function. The present results indicate that the number of sequence combinations at a structural locus, which are close to the maximum stability, is small. The most striking example is a selection at an interior locus of the pro-domain. After two rounds of phagemid selection, one amino acid combination is found in 40% of sequenced mutants. The most frequently selected mutant has a deltaG(unfolding) = 4 kcal/mol at 25 degrees C, an increase of 6 kcal/mol relative to the naturally occurring sequence. Some implications of these results on the amount of sequence information needed to specify a unique tertiary fold are discussed. Apart from possible implications on the folding code, the phage display selection described here should be useful in optimizing the stability of other proteins, which can be displayed on the phage surface.

Reproducing the natural evolution of protein structural features with the selectively infective phage (SIP) technology. The kink in the first strand of antibody kappa domains

The beta-sandwich structure of immunoglobulin variable domains is characterized by a typical kink in the first strand, which allows the first part of the strand to hydrogen bond to the outer beta-sheet (away from the VH-VL interface) and the second part to the inner beta-sheet. This kink differs in length and sequence between the Vkappa, Vlambda and VH domains and yet is involved in several almost perfectly conserved interactions with framework residues. We have used the selectively infective phage (SIP) system to select the optimal kink region from several defined libraries, using an anti-hemagglutinin single-chain Fv (scFv) fragment as a model system. Both for the kink with the Vkappa domain length and that with the Vlambda length, a sequence distribution was selected that coincides remarkably well with the sequence distribution of natural antibodies. The selected scFv fragments were purified and characterized, and thermodynamic stability was found to be the prime factor responsible for selection. These data show that the SIP technology can be used for optimizing protein structural features by evolutionary approaches.

The sequences of small proteins are not extensively optimized for rapid folding by natural selection

The thermodynamic stabilities of small protein domains are clearly subject to natural selection, but it is less clear whether the rapid folding rates typically observed for such proteins are consequences of direct evolutionary optimization or reflect intrinsic physical properties of the polypeptide chain. This issue can be investigated by comparing the folding rates of laboratory-generated protein sequences to those of naturally occurring sequences provided that the method by which the sequences are generated has no kinetic bias. Herein we report the folding thermodynamics and kinetics of 12 heavily mutated variants of the small IgG binding domain of protein L retrieved from high-complexity combinatorial libraries by using a phage-display selection for proper folding that does not discriminate between rapidly and slowly folding proteins. Although the stabilities of all variants were decreased, many of the variants fold faster than wild type. Taken together with similar results for the src homology 3 domain, this observation suggests that the sequences of small proteins have not been extensively optimized for rapid folding; instead, rapid folding appears to be a consequence of selection for stability.

Simplified proteins: minimalist solutions to the 'protein folding problem'

Recent research has suggested that stable, native proteins may be encoded by simple sequences of fewer than the full set of 20 proteogenic amino acids. Studies of the ability of simple amino acid sequences to encode stable, topologically complex, native conformations and to fold to these conformations in a biologically relevant time frame have provided insights into the sequence determinants of protein structure and folding kinetics. They may also have important implications for protein design and for theories of the origins of protein synthesis itself.

Contrasting roles for symmetrically disposed beta-turns in the folding of a small protein

To investigate the role of turns in protein folding, we have characterized the effects of combinatorial and site-directed mutations in the two beta-turns of peptostreptococcal protein L on folding thermodynamics and kinetics. Sequences of folded variants recovered from combinatorial libraries using a phase display selection method were considerably more variable in the second turn than in the first turn. These combinatorial mutants as well as strategically placed point mutants in the two turns had a similar range of thermodynamic stabilities, but strikingly different folding kinetics. A glycine to alanine substitution in the second beta-turn increased the rate of unfolding more than tenfold but had little effect on the rate of folding, while mutation of a symmetrically disposed glycine residue in the first turn had little effect on unfolding but slowed the rate of folding nearly tenfold. These results demonstrate that the role of beta-turns in protein folding is strongly context-dependent, and suggests that the first turn is formed and the second turn disrupted in the folding transition state.