Comprehensive mutational analysis of the M13 major coat protein: improved scaffolds for C-terminal phage display

A peptide was fused to the C terminus of the M13 bacteriophage major coat protein (P8), and libraries of P8 mutants were screened to select for variants that displayed the peptide with high efficiency. Over 600 variants were sequenced to compile a comprehensive database of P8 sequence diversity compatible with assembly into the wild-type phage coat. The database reveals that, while the alpha-helical P8 molecule was highly tolerant to mutations, certain functional epitopes were required for efficient incorporation. Three hydrophobic epitopes were located approximately equidistantly along the length of the alpha-helix. In addition, a positively charged epitope was required directly opposite the most C-terminal hydrophobic epitope and on the same side as the other two epitopes. Both ends of the protein were highly tolerant to mutations, consistent with the use of P8 as a scaffold for both N and C-terminal phage display. Further rounds of selection were used to enrich for P8 variants that supported higher levels of C-terminal peptide display. The largest improvements in display resulted from mutations around the junction between P8 and the C-terminal linker, and additional mutations in the N-terminal region were selected for further improvements in display. The best P8 variants improved C-terminal display more than 100-fold relative to the wild-type, and these variants could support the simultaneous display of N and C-terminal fusions. These finding provide information on the requirements for filamentous phage coat assembly, and provide improved scaffolds for phage display technology.     

A minimized M13 coat protein defines the requirements for assembly into the bacteriophage particle

The M13 filamentous bacteriophage coat is a symmetric array of several thousand alpha-helical major coat proteins (P8) that surround the DNA core. P8 molecules initially reside in the host membrane and subsequently transition into their role as coat proteins during the phage assembly process. A comprehensive mutational analysis of the 50-residue P8 sequence revealed that only a small subset of the side-chains were necessary for efficient incorporation into a wild-type (wt) coat. In the three-dimensional structure of P8, these side-chains cluster into three functional epitopes: a hydrophobic epitope located near the N terminus and two epitopes (one hydrophobic and the other basic) located near the C terminus on opposite faces of the helix. The results support a model for assembly in which the incorporation of P8 is mediated by intermolecular interactions involving these functional epitopes. In this model, the N-terminal hydrophobic epitope docks with P8 molecules already assembled into the phage particle in the periplasm, and the basic epitope interacts with the acidic DNA backbone in the cytoplasm. These interactions could facilitate the transition of P8 from the membrane into the assembling phage, and the incorporation of a single P8 would be completed by the docking of additional P8 molecules with the second hydrophobic epitope at the C terminus. We constructed a minimized P8 that contained only nine non-Ala side-chains yet retained all three functional epitopes. The minimized P8 assembled into the wt coat almost as efficiently as wt P8, thus defining the minimum requirements for protein incorporation into the filamentous phage coat. The results suggest possible mechanisms of natural viral evolution and establish guidelines for the artificial evolution of improved coat proteins for phage display technology.     

Comprehensive mutagenesis of the C-terminal domain of the M13 gene-3 minor coat protein: the requirements for assembly into the bacteriophage particle

Filamentous bacteriophage assemble at the host membrane in a non-lytic process; the gene-3 minor coat protein (P3) is required for release from the membrane and subsequently, for recognition and infection of a new host. P3 contains at least three distinct domains: two N-terminal domains that mediate host recognition and infection, and a C-terminal domain (P3-C) that is required for release from the host cell following phage assembly and contributes to the structural stability of the phage particle. A comprehensive mutational analysis of the 150 residue P3-C revealed that only 24 side-chains, located within the last 70 residues of sequence, were necessary for efficient incorporation into a wild-type coat. The results reveal that the requirements for the assembly of P3 into the phage particle are quite lax and involve only a few key side-chains. These findings shed light on the functional and structural requirements for filamentous phage assembly, and they may provide guidelines for the engineering of improved coat proteins as scaffolds for phage display technology.     

Mutational analysis of the major coat protein of M13 identifies residues that control protein display

We have reported variants of the M13 bacteriophage major coat protein (P8) that enable high copy display of monomeric and oligomeric proteins, such as human growth hormone and steptavidin, on the surface of phage particles (Sidhu SS, Weiss GA, Wells JA. 2000. High copy display of large proteins on phage for functional selections. J Mol Biol 296:487-495). Here, we explore how an optimized P8 variant (opti-P8) could evolve the ability to efficiently display a protein fused to its N-terminus. Reversion of individual opti-P8 residues back to the wild-type P8 residue identifies a limited set of hydrophobic residues responsible for the high copy protein display. These hydrophobic amino acids bracket a conserved hydrophobic face on the P8 alpha helix thought to be in contact with the phage coat. Mutations additively combine to promote high copy protein display, which was further enhanced by optimization of the linker between the phage coat and the fusion protein. These data are consistent with a model in which protein display-enhancing mutations allow for better packing of the fusion protein into the phage coat. The high tolerance for phage coat protein mutations observed here suggests that filamentous phage coat proteins could readily evolve new capabilities.     

A small protein in model membranes: a time-resolved fluorescence and ESR study on the interaction of M13 coat protein with lipid bilayers

Model membranes with unsaturated lipid chains containing various amounts of M13 coat protein in the -helical form were studied using time-resolved fluorescence and ESR spectroscopy. The lipid-to-protein (L/P) ratios used were > 12 to avoid protein-protein contacts and irreversible aggregation leading to -polymeric coat protein. In the ESR spectra of the 12-SASL probe in dioleoyl phosphatidylcholine (DOPC) bilayers no second protein induced component is observed upon incorporation of M13 coat protein. However, strong effects are detected on the ESR lineshapes upon changing the protein concentration. The ESR lineshapes are simulated by assuming a fixed ratio between the parallel (D) and perpendicular (D) diffusion coefficients of 4, and an order parameter equal to zero. It is found that increasing the protein concentration from L/P to L/P 15 results in a decrease of the rotational diffusion coefficient D from 3.4 × 10⁷ to 1.9 × 10⁷ s^–1. In the time-resolved fluorescence experiments with DPH-propionic acid as a probe, it is observed that increasing the M13 coat protein concentration causes an increase of the two fluorescent lifetimes, indicating an increase in bilayer order. Analysis of the time-resolved fluorescence anisotropy decay allows one to quantitatively determine the order parameters P₂ and P₄, and the rotational diffusion coefficient D of the fluorescent probe. The order parameters P₂ and P₄ increase from 0.34 to 0.55 and from 0.59 to 0.77, respectively, upon adding M13 coat protein to DOPC bilayers with an L/P ratio of 35. The rotational diffusion coefficient D of the DPH-propionic acid probe decreases on incorporating M13 coat protein, in accordance with the ESR results. It is concluded that M13 coat protein in the -monomeric state is not able to produce a long living lipid boundary shell and consequently an immobilization of the lipids. An overall effect on the lipids is induced, resulting in a reduction in the dynamics and an increase in average lipid order. The hydrophobic region of M13 coat protein is proposed to perfectly match the lipid bilayer, resulting in a relatively small distortion of the bilayer structure of the lipid system.     

Major coat proteins of bacteriophage Pf3 and M13 as model systems for Sec-independent protein transport

Abstract: The membrane insertion of bacteriophage coat proteins occurs independent of the Sec-translocase of Escherichia coli . Detailed study of the Pf3 and M13 coat proteins has elucidated two fundamental mechanisms of how proteins invade the membrane, most likely by direct interaction with the lipid bilayer. The Sec-independent translocation of amino-terminal regions across the inner membrane is limited to a short length and a small number of charged residues. Protein regions that contain several charged residues are efficiently translocated across the membrane when these regions are flanked by two adjacent hydrophobic segments interacting synergistically. The relevance of these findings for the membrane insertion mechanism of multispanning membrane proteins is discussed.     

Crystal structure of the coat protein from the GA bacteriophage: model of the unassembled dimer.

There are four groups of RNA bacteriophages with distinct antigenic and physicochemical properties due to differences in surface residues of the viral coat proteins. Coat proteins also play a role as translational repressor during the viral life cycle, binding an RNA hairpin within the genome. In this study, the first crystal structure of the coat protein from a Group II phage GA is reported and compared to the Group I MS2 coat protein. The structure of the GA dimer was determined at 2.8 A resolution (R-factor = 0.20). The overall folding pattern of the coat protein is similar to the Group I MS2 coat protein in the intact virus (Golmohammadi R, Valegård K, Fridborg K, Liljas L. 1993, J Mol Biol 234:620-639) or as an unassembled dimer (Ni Cz, Syed R, Kodandapani R. Wickersham J, Peabody DS, Ely KR, 1995, Structure 3:255-263). The structures differ in the FG loops and in the first turn of the alpha A helix. GA and MS2 coat proteins differ in sequence at 49 of 129 amino acid residues. Sequence differences that contribute to distinct immunological and physical properties of the proteins are found at the surface of the intact virus in the AB and FG loops. There are six differences in potential RNA contact residues within the RNA-binding site located in an antiparallel beta-sheet across the dimer interface. Three differences involve residues in the center of this concave site: Lys/Arg 83, Ser/Asn 87, and Asp/Glu 89. Residue 87 was shown by molecular genetics to define RNA-binding specificity by GA or MS2 coat protein (Lim F. Spingola M, Peabody DS, 1994, J Biol Chem 269:9006-9010). This sequence difference reflects recognition of the nucleotide at position -5 in the unpaired loop of the translational operators bound by these coat proteins. In GA, the nucleotide at this position is a purine whereas in MS2, it is a pyrimidine.     

The molecular basis for genetic polymorphism of human deoxyribonuclease II (DNase II): a single nucleotide substitution in the promoter region of human DNase II changes the promoter activity

We report that, contrary to common belief, polypeptides fused to the carboxy-terminus of the M13 gene-3 minor coat protein are functionally displayed on the phage surface. In a phagemid display system, carboxy-terminal fusion through optimized linker sequences resulted in display levels comparable to those achieved with conventional amino-terminal fusions. These findings are of considerable importance to phage display technology because they enable investigations not suited to amino-terminal display, including the study of protein–protein interactions requiring free carboxy-termini, functional cDNA cloning efforts, and the display of intracellular proteins.     

Assessment of chemical-crosslink-assisted protein structure modeling in CASP13

With the advance of experimental procedures obtaining chemical crosslinking information is becoming a fast and routine practice. Information on crosslinks can greatly enhance the accuracy of protein structure modeling. Here, we review the current state of the art in modeling protein structures with the assistance of experimentally determined chemical crosslinks within the framework of the 13th meeting of Critical Assessment of Structure Prediction approaches. This largest-to-date blind assessment reveals benefits of using data assistance in difficult to model protein structure prediction cases. However, in a broader context, it also suggests that with the unprecedented advance in accuracy to predict contacts in recent years, experimental crosslinks will be useful only if their specificity and accuracy further improved and they are better integrated into computational workflows.     

C‐terminal juxtamembrane region of full‐length M2 protein forms a membrane surface associated amphipathic helix

The influenza A M2 protein is a 97‐residue integral membrane protein involved in viral budding and proton conductance. Although crystal and NMR structures exist of truncated constructs of the protein, there is disagreement between models and only limited structural data are available for the full‐length protein. Here, the structure of the C‐terminal juxtamembrane region (sites 50–60) is investigated in the full‐length M2 protein using site‐directed spin‐labeling electron paramagnetic resonance (EPR) spectroscopy in lipid bilayers. Sites 50–60 were chosen for study because this region has been shown to be critical to the role the M2 protein plays in viral budding. Continuous wave EPR spectra and power saturation data in the presence of paramagnetic membrane soluble oxygen are consistent with a membrane surface associated amphipathic helix. Comparison between data from the C‐terminal juxtamembrane region in full‐length M2 protein with data from a truncated M2 construct demonstrates that the line shapes and oxygen accessibilities are remarkably similar between the full‐length and truncated form of the protein.     

The in situ aggregational and conformational state of the major coat protein of bacteriophage M13 in phospholipid bilayers mimicking the inner membrane of host Escherichia coli

The major coat protein of bacteriophage M13 has been reconstituted into phospholipids with a composition comparable to that found in the host (Escherichia coli) inner membrane. Reconstitution experiments have revealed conditions in which the alpha-oligomeric state is favored over the beta-polymeric state. Discrimination between the two states of the membrane-bound coat protein (alpha-oligomeric and beta-polymeric states) has been achieved using high-performance size-exclusion chromatography and circular dichroism. Interprotein electrostatic interactions, probably induced by head-tail binding, are initiated and facilitating the aggregation-related conformational change process, in which alpha-oligomeric coat protein is converted into beta-polymeric coat protein. A model for this beta-polymerization process of the coat protein is presented. The alpha-helical protein has been studied by the in situ Trp fluorescence quantum yield. This shows that the average distances between coat proteins decrease upon lowering the L/P ratio. In situ cross-linking reactions of the coat protein at high L/P ratios reveal a monomeric state, thus excluding specific aggregation of the coat protein. A monomeric state of detergent-solubilized coat protein is also observed using SDS-PAGE and SDS-HPSEC. On the basis of these results, the smallest in situ aggregational entity of the coat protein is proposed to be a monomer. This finding is discussed in relation to the functional state of the M13 coat protein in the membrane-bound assembly and disassembly processes during infection.     

Membrane protein dynamics versus environment: simulations of OmpA in a micelle and in a bilayer

The bacterial outer membrane protein OmpA is one of the few membrane proteins whose structure has been solved both by X-ray crystallography and by NMR. Crystals were obtained in the presence of detergent, and the NMR structure is of the protein in a detergent micelle. We have used 10 ns duration molecular dynamics simulations to compare the behaviour of OmpA in a detergent micelle and in a phospholipid bilayer. The dynamic fluctuations of the protein structure seem to be ca 1.5 times greater in the micelle environment than in the lipid bilayer. There are subtle differences between the nature of OmpA-detergent and OmpA-lipid interactions. As a consequence of the enhanced flexibility of the OmpA protein in the micellar environment, side-chain torsion angle changes are such as to lead to formation of a continuous pore through the centre of the OmpA molecule. This may explain the experimentally observed channel formation by OmpA.     

Investigating energy‐based pool structure selection in the structure ensemble modeling with experimental distance constraints: The example from a multidomain protein Pub1

The structural variations of multidomain proteins with flexible parts mediate many biological processes, and a structure ensemble can be determined by selecting a weighted combination of representative structures from a simulated structure pool, producing the best fit to experimental constraints such as interatomic distance. In this study, a hybrid structure‐based and physics‐based atomistic force field with an efficient sampling strategy is adopted to simulate a model di‐domain protein against experimental paramagnetic relaxation enhancement (PRE) data that correspond to distance constraints. The molecular dynamics simulations produce a wide range of conformations depicted on a protein energy landscape. Subsequently, a conformational ensemble recovered with low‐energy structures and the minimum‐size restraint is identified in good agreement with experimental PRE rates, and the result is also supported by chemical shift perturbations and small‐angle X‐ray scattering data. It is illustrated that the regularizations of energy and ensemble‐size prevent an arbitrary interpretation of protein conformations. Moreover, energy is found to serve as a critical control to refine the structure pool and prevent data overfitting, because the absence of energy regularization exposes ensemble construction to the noise from high‐energy structures and causes a more ambiguous representation of protein conformations. Finally, we perform structure‐ensemble optimizations with a topology‐based structure pool, to enhance the understanding on the ensemble results from different sources of pool candidates.     

Probing the structure of the Ff bacteriophage major coat protein transmembrane helix dimer by solution NMR

The transmembrane (TM) segment of the major coat protein from Ff bacteriophage has been extensively studied as an example of dimerization in detergent and lipid bilayer systems. However, almost all the information regarding this interaction has been gained through mutagenesis studies, with little direct structural information being available. To this end solution NMR has the potential to provide new insights into structure of the dimer. In order to evaluate the utility of this approach we have studied a selectively ¹⁵N-labeled peptide containing the TM segment of MCP (MCP_TM) by solution NMR. This peptide was found to give rise to detergent concentration-dependent spectra that were assigned to monomeric and dimeric forms. The standard free energy of this interaction in SDS was estimated from these spectra and found to be consistent with weak but specific dimerization. In addition, similar spectra could be obtained in β-octyl glucoside with intermolecular paramagnetic relaxation experiments demonstrating a parallel arrangement of TM helices in the dimer. In both detergents backbone chemical shift differences between monomeric and dimeric forms of MCP_TM showed that the largest changes occur around its GXXXG motif. The resulting structural model is consistent with observations made for MCP mutants previously characterized in biological membranes, opening the door to detailed structural characterization of this form of MCP. These results also have general implications for the study of weakly interacting TM segments by solution NMR since the use of similar sample conditions should allow structural data to be accessed for oligomeric states from a wide range systems that undergo biologically relevant but weak associations in the membrane.     

Competing Lipid-Protein and Protein-Protein Interactions Determine Clustering and Gating Patterns in the Potassium Channel from Streptomyces lividans (KcsA)

There is increasing evidence to support the notion that membrane proteins, instead of being isolated components floating in a fluid lipid environment, can be assembled into supramolecular complexes that take part in a variety of cooperative cellular functions. The interplay between lipid-protein and protein-protein interactions is expected to be a determinant factor in the assembly and dynamics of such membrane complexes. Here we report on a role of anionic phospholipids in determining the extent of clustering of KcsA, a model potassium channel. Assembly/disassembly of channel clusters occurs, at least partly, as a consequence of competing lipid-protein and protein-protein interactions at nonannular lipid binding sites on the channel surface and brings about profound changes in the gating properties of the channel. Our results suggest that these latter effects of anionic lipids are mediated via the Trp⁶⁷–Glu⁷¹–Asp⁸⁰ inactivation triad within the channel structure and its bearing on the selectivity filter.     

The use of soluble protein structures in modeling helical proteins in a layered membrane

Major advances have been made in the prediction of soluble protein structures, led by the knowledge-based modeling methods that extract useful structural trends from known protein structures and incorporate them into scoring functions. The same cannot be reported for the class of transmembrane proteins, primarily due to the lack of high-resolution structural data for transmembrane proteins, which render many of the knowledge-based method unreliable or invalid. We have developed a method that harnesses the vast structural knowledge available in soluble protein data for use in the modeling of transmembrane proteins. At the core of the method, a set of transmembrane protein decoy sets that allow us to filter and train features recognized from soluble proteins for transmembrane protein modeling into a set of scoring functions. We have demonstrated that structures of soluble proteins can provide significant insight into transmembrane protein structures. A complementary novel two-stage modeling/selection process that mimics the two-stage helical membrane protein folding was developed. Combined with the scoring function, the method was successfully applied to model 5 transmembrane proteins. The root mean square deviations of the predicted models ranged from 5.0 to 8.8?Å to the native structures.     

Modeling of the N-terminal Section and the Lumenal Loop of Trimeric Light Harvesting Complex II (LHCII) by Using EPR

The major light harvesting complex II (LHCII) of green plants plays a key role in the absorption of sunlight, the regulation of photosynthesis, and in preventing photodamage by excess light. The latter two functions are thought to involve the lumenal loop and the N-terminal domain. Their structure and mobility in an aqueous environment are only partially known. Electron paramagnetic resonance (EPR) has been used to measure the structure of these hydrophilic protein domains in detergent-solubilized LHCII. A new technique is introduced to prepare LHCII trimers in which only one monomer is spin-labeled. These heterogeneous trimers allow to measure intra-molecular distances within one LHCII monomer in the context of a trimer by using double electron-electron resonance (DEER). These data together with data from electron spin echo envelope modulation (ESEEM) allowed to model the N-terminal protein section, which has not been resolved in current crystal structures, and the lumenal loop domain. The N-terminal domain covers only a restricted area above the superhelix in LHCII, which is consistent with the “Velcro” hypothesis to explain thylakoid grana stacking (Standfuss, J., van Terwisscha Scheltinga, A. C., Lamborghini, M., and Kühlbrandt, W. (2005) EMBO J. 24, 919–928). The conformation of the lumenal loop domain is surprisingly different between LHCII monomers and trimers but not between complexes with and without neoxanthin bound.