Basic amino acids in a distinct subset of signal peptides promote interaction with the signal recognition particle

Previous studies have demonstrated that signal peptides bind to the signal recognition particle (SRP) primarily via hydrophobic interactions with the 54-kDa protein subunit. The crystal structure of the conserved SRP ribonucleoprotein core, however, raised the surprising possibility that electrostatic interactions between basic amino acids in signal peptides and the phosphate backbone of SRP RNA may also play a role in signal sequence recognition. To test this possibility we examined the degree to which basic amino acids in a signal peptide influence the targeting of two Escherichia coli proteins, maltose binding protein and OmpA. Whereas both proteins are normally targeted to the inner membrane by SecB, we found that replacement of their native signal peptides with another moderately hydrophobic but unusually basic signal peptide (DeltaEspP) rerouted them into the SRP pathway. Reduction in either the net positive charge or the hydrophobicity of the DeltaEspP signal peptide decreased the effectiveness of SRP recognition. A high degree of hydrophobicity, however, compensated for the loss of basic residues and restored SRP binding. Taken together, the data suggest that the formation of salt bridges between SRP RNA and basic amino acids facilitates the binding of a distinct subset of signal peptides whose hydrophobicity falls slightly below a threshold level.     

Amino acid conformational preferences and solvation of polar backbone atoms in peptides and proteins

Amino acids in peptides and proteins display distinct preferences for alpha-helical, beta-strand, and other conformational states. Various physicochemical reasons for these preferences have been suggested: conformational entropy, steric factors, hydrophobic effect, and backbone electrostatics; however, the issue remains controversial. It has been proposed recently that the side-chain-dependent solvent screening of the local and non-local backbone electrostatic interactions primarily determines the preferences not only for the alpha-helical but also for all other main-chain conformational states. Side-chains modulate the electrostatic screening of backbone interactions by excluding the solvent from the vicinity of main-chain polar atoms. The deficiency of this electrostatic screening model of amino acid preferences is that the relationships between the main-chain electrostatics and the amino acid preferences have been demonstrated for a limited set of six non-polar amino acid types in proteins only. Here, these relationships are determined for all amino acid types in tripeptides, dekapeptides, and proteins. The solvation free energies of polar backbone atoms are approximated by the electrostatic contributions calculated by the finite difference Poisson-Boltzmann and the Langevin dipoles methods. The results show that the average solvation free energy of main-chain polar atoms depends strongly on backbone conformation, shape of side-chains, and exposure to solvent. The equilibrium between the low-energy beta-strand conformation of an amino acid (anti-parallel alignment of backbone dipole moments) and the high-energy alpha conformation (parallel alignment of backbone dipole moments) is strongly influenced by the solvation of backbone polar atoms. The free energy cost of reaching the alpha conformation is by approximately 1.5 kcal/mol smaller for residues with short side-chains than it is for the large beta-branched amino acid residues. This free energy difference is comparable to those obtained experimentally by mutation studies and is thus large enough to account for the distinct preferences of amino acid residues. The screening coefficients gamma(local)(r) and gamma(non-local)(r) correlate with the solvation effects for 19 amino acid types with the coefficients between 0.698 to 0.851, depending on the type of calculation and on the set of point atomic charges used. The screening coefficients gamma(local)(r) increase with the level of burial of amino acids in proteins, converging to 1.0 for the completely buried amino acid residues. The backbone solvation free energies of amino acid residues involved in strong hydrogen bonding (for example: in the middle of an alpha-helix) are small. The hydrogen bonded backbone is thus more hydrophobic than the peptide groups in random coil. The alpha-helix forming preference of alanine is attributed to the relatively small free energy cost of reaching the high-energy alpha-helix conformation. These results confirm that the side-chain-dependent solvent screening of the backbone electrostatic interactions is the dominant factor in determining amino acid conformational preferences.     

Optimization of the electrostatic interactions in proteins of different functional and folding type

The 3-dimensional optimization of the electrostatic interactions between the charged amino acid residues was studied by Monte Carlo simulations on an extended representative set of 141 protein structures with known atomic coordinates. The proteins were classified by different functional and structural criteria, and the optimization of the electrostatic interactions was analyzed. The optimization parameters were obtained by comparison of the contribution of charge-charge interactions to the free energy of the native protein structures and for a large number of randomly distributed charge constellations obtained by the Monte Carlo technique. On the basis of the results obtained, one can conclude that the charge-charge interactions are better optimized in the enzymes than in the proteins without enzymatic functions. Proteins that belong to the mixed αβ folding type are electrostatically better optimized than pure α-helical or β-strand structures. Proteins that are stabilized by disulfide bonds show a lower degree of electrostatic optimization. The electrostatic interactions in a native protein are effectively optimized by rejection of the conformers that lead to repulsive charge-charge interactions. Particularly, the rejection of the repulsive contacts seems to be a major goal in the protein folding process. The dependence of the optimization parameters on the choice of the potential function was tested. The majority of the potential functions gave practically identical results.     

Dominant role of local dipoles in stabilizing uncompensated charges on a sulfate sequestered in a periplasmic active transport protein.

Electrostatic interactions are among the key factors determining the structure and function of proteins. Here we report experimental results that illuminate the functional importance of local dipoles to these interactions. The refined 1.7-A X-ray structure of the liganded form of the sulfate-binding protein, a primary sulfate active transport receptor of Salmonella typhimurium, shows that the sulfate dianion is completely buried and bound by hydrogen bonds (mostly main-chain peptide NH groups) and van der Waals forces. The sulfate is also closely linked, via one of these peptide units, to a His residue. It is also adjacent to the N-termini of three alpha-helices, of which the two shortest have their C-termini "capped" by Arg residues. Site-directed mutagenesis of the recombinant Escherichia coli sulfate receptor had no effect on sulfate-binding activity when an Asn residue was substituted for the positively charged His and the two Arg (changed singly and together) residues. These results, combined with other observations, further solidify the idea that stabilization of uncompensated charges in a protein is a highly localized process that involves a collection of local dipoles, including those of peptide units confined to the first turns of helices. The contribution of helix macrodipoles appears insignificant.     

Calculation of protein backbone geometry from alpha-carbon coordinates based on peptide-group dipole alignment.

An algorithm is proposed for the conversion of a virtual-bond polypeptide chain (connected C alpha atoms) to an all-atom backbone, based on determining the most extensive hydrogen-bond network between the peptide groups of the backbone, while maintaining all of the backbone atoms in energetically feasible conformations. Hydrogen bonding is represented by aligning the peptide-group dipoles. These peptide groups are not contiguous in the amino acid sequence. The first dipoles to be aligned are those that are both sufficiently close in space to be arranged in approximately linear arrays termed dipole paths. The criteria used in the construction of dipole paths are: to assure good alignment of the greatest possible number of dipoles that are close in space; to optimize the electrostatic interactions between the dipoles that belong to different paths close in space; and to avoid locally unfavorable amino acid residue conformations. The equations for dipole alignment are solved separately for each path, and then the remaining single dipoles are aligned optimally with the electrostatic field from the dipoles that belong to the dipole-path network. A least-squares minimizer is used to keep the geometry of the alpha-carbon trace of the resulting backbone close to that of the input virtual-bond chain. This procedure is sufficient to convert the virtual-bond chain to a real chain; in applications to real systems, however, the final structure is obtained by minimizing the total ECEPP/2 (empirical conformational energy program for peptides) energy of the system, starting from the geometry resulting from the solution of the alignment equations. When applied to model alpha-helical and beta-sheet structures, the algorithm, followed by the ECEPP/2 energy minimization, resulted in an energy and backbone geometry characteristic of these alpha-helical and beta-sheet structures. Application to the alpha-carbon trace of the backbone of the crystallographic 5PTI structure of bovine pancreatic trypsin inhibitor, followed by ECEPP/2 energy minimization with C alpha-distance constraints, led to a structure with almost as low energy and root mean square deviation as the ECEPP/2 geometry analog of 5PTI, the best agreement between the crystal and reconstructed backbone being observed for the residues involved in the dipole-path network.     

Peptide tags for enhanced cellular and protein adhesion to single-crystalline sapphire

In order to facilitate a novel means for coupling proteins to metal oxides, peptides were identified from a dodecamer peptide yeast surface display library that bound a model metal oxide material, the C, A, and R crystalline faces of synthetic sapphire (alpha-Al(2)O(3)). Seven rounds of screening yielded peptides enriched in basic amino acids compared to the naive library. While the C-face had a high background of endogenous yeast cell binding, the A- and R faces displayed clear peptide-mediated cell adhesion. Cell detachment assays showed that cell adhesion strength correlated positively with increasing basicity of expressed peptides. Cell adhesion was also shown to be sensitive to buffer ionic strength as well as incubation with soluble peptide (with half maximal inhibition of cell binding at approximately 5 microM peptide). Next, dodecamer peptides cloned into yeast showed that lysine led to stronger interactions than arginine, and that charge distribution affected adhesion strength. We postulate binding to arise from peptide geometries that permit conformation alignment of the basic amino acids towards the surface so that the charged groups can undergo local electrostatic interactions with the surface oxide. Lastly, peptide K1 (-(GK)(6)) was cloned onto the c-terminus of maltose binding protein (MBP) and the resultant mutant protein showed a half-maximal binding at approximately 10(-7)-10(-6) M, which marked a approximately 500- to 1,000-fold binding improvement to sapphire's A-face as compared with wild-type MBP. Targeting proteins to metal oxide surfaces with peptide tags may provide a facile one-step alternative coupling chemistry for the formation of protein bioassays and biosensors.     

Mapping the electrostatic potential within the ribosomal exit tunnel

Electrostatic potentials influence interactions among proteins and nucleic acids, the orientation of dipoles and quadrupoles, and the distribution of mobile charges. Consequently, electrostatic potentials can modulate macromolecular folding and conformational stability, as well as rates of catalysis and substrate binding. The ribosomal exit tunnel, along with its resident nascent peptide, is no less susceptible to these consequences. Yet, the electrostatics inside the tunnel have never been measured. Here we map both the electrostatic potential and accessibilities along the length of the tunnel and determine the electrostatic consequences of introducing a charged amino acid into the nascent peptide. To do this we developed novel probes and strategies. Our findings provide new insights regarding the dielectric of the tunnel and the dynamics of its local electric fields.     

Designed armadillo repeat proteins as general peptide-binding scaffolds: consensus design and computational optimization of the hydrophobic core

Armadillo repeat proteins are abundant eukaryotic proteins involved in several cellular processes, including signaling, transport, and cytoskeletal regulation. They are characterized by an armadillo domain, composed of tandem armadillo repeats of approximately 42 amino acids, which mediates interactions with peptides or parts of proteins in extended conformation. The conserved binding mode of the peptide in extended form, observed for different targets, makes armadillo repeat proteins attractive candidates for the generation of modular peptide-binding scaffolds. Taking advantage of the large number of repeat sequences available, a consensus-based approach combined with a force field-based optimization of the hydrophobic core was used to derive soluble, highly expressed, stable, monomeric designed proteins with improved characteristics compared to natural armadillo proteins. These sequences constitute the starting point for the generation of designed armadillo repeat protein libraries for the selection of peptide binders, exploiting their modular structure and their conserved binding mode.     

Spin-label electron spin resonance studies on the interactions of lysine peptides with phospholipid membranes.

The interactions of lysine oligopeptides with dimyristoyl phosphatidylglycerol (DMPG) bilayer membranes were studied using spin-labeled lipids and electron spin resonance spectroscopy. Tetralysine and pentalysine were chosen as models for the basic amino acid clusters found in a variety of cytoplasmic membrane-associating proteins, and polylysine was chosen as representative of highly basic peripherally bound proteins. A greater motional restriction of the lipid chains was found with increasing length of the peptide, while the saturation ratio of lipids per peptide was lower for the shorter peptides. In DMPG and dimyristoylphosphatidylserine host membranes, the perturbation of the lipid chain mobility by polylysine was greater for negatively charged spin-labeled lipids than for zwitterionic lipids, but for the shorter lysine peptides these differences were smaller. In mixed bilayers composed of DMPG and dimyristoylphosphatidylcholine, little difference was found in selectivity between spin-labeled phospholipid species on binding pentalysine. Surface binding of the basic lysine peptides strongly reduced the interfacial pK of spin-labeled fatty acid incorporated into the DMPG bilayers, to a greater extent for polylysine than for tetralysine or pentalysine at saturation. The results are consistent with a predominantly electrostatic interaction with the shorter lysine peptides, but with a closer surface association with the longer polylysine peptide.     

Electrostatic interactions in proteins: Calculations of the electrostatic term of free energy and the electrostatic potential field

The pH-dependence of the electrostatic energy of interactions between titratable groups is calculated for some well studied globular proteins: basic pancreatic trypsin inhibitor, sperm whale myoglobin and tuna cytochrome c. The calculations are carried out using a semi-empirical appraach in terms of the macroscopic model based on the Kirkwood-Tanford theory. The results are discussed in the light of their physicochemical and biological properties. It was found that the pH-dependence of the electrostatic energy correlates with the III–IV transition of cytochrome c. The electrostatic field of the cysteine proteinase inhibitor, cystatin, was calculated in two ways. In the first one, the electrostatic field created by the pH dependent charges of the ionizable groups and peptide dipoles was calculated using the approach proposed. In the second one, the finite-difference method was used. The results obtained by the two methods are in overall agreement. The calculated field was discussed in terms of the binding of cystatin to papain.     

The optimization of protein-solvent interactions: thermostability and the role of hydrophobic and electrostatic interactions.

Protein-solvent interactions were analyzed using an optimization parameter based on the ratio of the solvent-accessible area in the native and the unfolded protein structure. The calculations were performed for a set of 183 nonhomologous proteins with known three-dimensional structure available in the Protein Data Bank. The dependence of the total solvent-accessible surface area on the protein molecular mass was analyzed. It was shown that there is no difference between the monomeric and oligomeric proteins with respect to the solvent-accessible area. The results also suggested that for proteins with molecular mass above some critical mass, which is about 28 kDa, a formation of domain structure or subunit aggregation into oligomers is preferred rather than a further enlargement of a single domain structure. An analysis of the optimization of both protein-solvent and charge-charge interactions was performed for 14 proteins from thermophilic organisms. The comparison of the optimization parameters calculated for proteins from thermophiles and mesophiles showed that the former are generally characterized by a high degree of optimization of the hydrophobic interactions or, in cases where the optimization of the hydrophobic interactions is not sufficiently high, by highly optimized charge-charge interactions.     

Peptide inhibitors of fibronectin, laminin, and other adhesion molecules: unique and shared features

Synthetic peptides can specifically inhibit the function of certain adhesive glycoproteins in vitro and in vivo. We have compared the relative activities of a set of six variant synthetic peptides based on the sequence of fibronectin in terms of their ability to inhibit the interactions of fibroblasts with fibronectin, spreading factor/vitronectin, laminin, and native collagen gels. BHK (baby hamster kidney) and chick embryo fibroblasts spreading on these adhesive molecules displayed distinctive patterns of sensitivity to inhibition by this panel of peptides, which depended on the adhesive molecule rather than the cell type. For fibronectin, Gly-Arg-Gly-Asp-Ser was considerably more active than Arg-Gly-Asp-Ser, whereas these two peptides displayed little difference in activity in inhibiting cell adhesion to spreading factor. For both proteins, the inverted peptide sequence Ser-Asp-Gly-Arg was also moderately active, whereas closely related peptides containing a transposition, a deletion, or a single, conserved amino acid substitution were much less active. For inhibiting interactions with laminin or native type I collagen gels, Gly-Arg-Gly-Asp-Ser was only weakly active, but the inverted peptide Ser-Asp-Gly-Arg unexpectedly continued to display inhibitory activity for both attachment proteins in both cell types. Our results indicate that different adhesive processes depend on distinct peptide recognition events by a cell. However, there may be a possible common denominator among attachment proteins in a moderate sensitivity to Ser-Asp-Gly-Arg. Our study also underscores the importance of examining a full set of peptide analogs when these novel inhibitors are used to characterize biological processes.     

The HP1 chromo shadow domain binds a consensus peptide pentamer

Heterochromatin-associated protein 1 (HP1) is thought to affect chromatin structure through interactions with other proteins in heterochromatin. Chromo domains located near the amino (amino chromo) and carboxy (chromo shadow) termini of HP1 may mediate such interactions, as suggested by domain swapping, in vitro binding and 3D structural studies . Several HP1-associated proteins have been reported, providing candidates that might specifically complex with the chromo domains of HP1. However, such association studies provide little mechanistic insight and explore only a limited set of potential interactions in a largely non-competitive setting. To determine how chromo domains can selectively interact with other proteins, we probed random peptide phage display libraries using chromo domains from HP1. Our results demonstrate that a consensus pentapeptide is suffident for specific interaction with the HP1 chromo shadow domain. The pentapeptide is found in the amino acid sequence of reported HP1-associated proteins, including the shadow domain itself. Peptides that bind the shadow domain also disrupt shadow domain dimers. Our results suggest that HP1 dimerization, which is thought to mediate heterochromatin compaction and cohesion, occurs via pentapeptide binding. In general, chromo domains may function by avidly binding short peptides at the surface of chromatin-associated proteins.     

Additive method for the prediction of protein-peptide binding affinity. Application to the MHC class I molecule HLA-A*0201

A method has been developed for prediction of binding affinities between proteins and peptides. We exemplify the method through its application to binding predictions of peptides with affinity to major histocompatibility complex class I molecule HLA-A*0201. The method is named "additive" because it is based on the assumption that the binding affinity of a peptide could be presented as a sum of the contributions of the amino acids at each position and the interactions between them. The amino acid contributions and the contributions of the interactions between adjacent side chains and every second side chain were derived using a partial least squares (PLS) statistical methodology using a training set of 420 experimental IC50 values. The predictive power of the method was assessed using rigorous cross-validation and using an independent test set of 89 peptides. The mean value of the residuals between the experimental and predicted pIC50 values was 0.508 for this test set. The additive method was implemented in a program for rapid T-cell epitope search. It is universal and can be applied to any peptide-protein interaction where binding data is known.     

Beta-sheet capping: signals that initiate and terminate beta-sheet formation

In the present work, we address the question of whether different amino acids have different beta-sheet initiating and terminating characteristics. Using a large scale analysis of parallel and antiparallel beta-sheets in a non-redundant dataset of proteins, we observed that most of the amino acids show significant under- or over-representation in at least one of the positions at the two ends of beta-sheets, which are denoted as N-cap and C-cap. In addition, based on statistical data and structural comparison, we found that certain amino acids, especially Asp, Asn, Gly and Pro have strong tendencies to block beta-sheet continuation. Hence, we can consider these residues as beta-sheet terminators. It was also proposed that the dipole moments in parallel beta-sheets, whose direction is from C-terminal (partially negative) to N-terminal (partially positive), are much stronger than has previously been suggested. In fact, enhancement of dipole moments in parallel beta-sheets is a result of the positioning of positively charged residues at N-cap and negatively charged residues at C-cap. This enhancement in dipole moment magnitude leads to strengthened dipolar interactions between parallel beta-sheets dipoles and other partners especially alpha-helices dipoles. The results provide an explanation for the antiparallel alignment of parallel beta-sheets with alpha-helices.     

15-Deoxyspergualin hinders physical interaction between basic residues of transit peptide in PfENR and Hsp70-1

The apicoplast of Plasmodium harbors several metabolic pathways. The enzymes required to perform these reactions are all nuclearly encoded and apicoplast targeted (NEAT) proteins. Plasmodium falciparum Enoyl-ACP Reductase (PfENR) is one such NEAT protein. The NEAT proteins have a transit peptide which is required for crossing the membranes of apicoplast. We studied the importance of basic residues like Arginine and Lysine within the transit peptide. Previous studies have suggested that all basic residues are essential for apicoplast trafficking. In this study, we demonstrate that only some of these residues are essential (K44, R48, K51, and R52), whereas others are dispensable (R40, K42, and K49). On mutating these specific residues, PfENR is not imported into the apicoplast and is mislocalized to the cytoplasm. We also demonstrate that these residues are also crucial for interaction with Hsp70-1, implying that interactions of Lysine 44, Arginine 48, Lysine 51, and Arginine 52 of the transit peptide with PfHsp70-1 are required for apicoplast trafficking. 15-Deoxyspergualin, which has earlier been proposed to interact with EEVD motif of PfHsp70-1 hinders the physical interaction between these cationic residues of PfENR and Hsp70-1. Hence, we propose that in the transport competent state of NEAT proteins some specific positively charged amino acids in the transit peptide interact with PfHsp70-1, and this interaction is essential for apicoplast targeting.     

Ionic Interactions Promote Transmembrane Helix-Helix Association Depending on Sequence Context

Folding and oligomerization of integral membrane proteins frequently depend on specific interactions of transmembrane helices. Interacting amino acids of helix-helix interfaces may form complex motifs and exert different types of molecular forces. Here, a set of strongly self-interacting transmembrane domains (TMDs), as isolated from a combinatorial library, was found to contain basic and acidic residues, in combination with polar nonionizable amino acids and C-terminal GxxxG motifs. Mutational analyses of selected sequences and reconstruction of high-affinity interfaces confirmed the cooperation of these residues in homotypic interactions. Probing heterotypic interaction indicated the presence of interhelical charge-charge interactions. Furthermore, simple motifs of an ionizable residue and GxxxG are significantly overrepresented in natural TMDs, and a specific combination of these motifs exhibits high-affinity heterotypic interaction. We conclude that intramembrane charge-charge interactions depend on sequence context. Moreover, they appear important for homotypic and heterotypic interactions of numerous natural TMDs.     

Accurate Prediction of Peptide Binding Sites on Protein Surfaces

Many important protein–protein interactions are mediated by the binding of a short peptide stretch in one protein to a large globular segment in another. Recent efforts have provided hundreds of examples of new peptides binding to proteins for which a three-dimensional structure is available (either known experimentally or readily modeled) but where no structure of the protein–peptide complex is known. To address this gap, we present an approach that can accurately predict peptide binding sites on protein surfaces. For peptides known to bind a particular protein, the method predicts binding sites with great accuracy, and the specificity of the approach means that it can also be used to predict whether or not a putative or predicted peptide partner will bind. We used known protein–peptide complexes to derive preferences, in the form of spatial position specific scoring matrices, which describe the binding-site environment in globular proteins for each type of amino acid in bound peptides. We then scan the surface of a putative binding protein for sites for each of the amino acids present in a peptide partner and search for combinations of high-scoring amino acid sites that satisfy constraints deduced from the peptide sequence. The method performed well in a benchmark and largely agreed with experimental data mapping binding sites for several recently discovered interactions mediated by peptides, including RG-rich proteins with SMN domains, Epstein-Barr virus LMP1 with TRADD domains, DBC1 with Sir2, and the Ago hook with Argonaute PIWI domain. The method, and associated statistics, is an excellent tool for predicting and studying binding sites for newly discovered peptides mediating critical events in biology.