Effect of sequence-specific interactions on the stability of helical conformations in polypeptides

A modification of the Zimm–Bragg two-state model for the helix–coil transition in polypeptides, which considers the effect of charge–dipole, charge–charge, and other specific interactions on helix stability, is presented. The new model introduces a series of adjustable parameters whose values are estimated by fitting to recent spectroscopic results on medium-sized peptides. This formalism, based on traditional two-state helix–coil transition models, provides a framework in which data on the helix contents of peptides of specific sequence can be rationalized by a statistical mechanical theory.     

Helix-coil transition theory including long-range electrostatic interactions: application to globular proteins

An extension of the Zimm–Bragg two-state theory for the helix–coil transition in polypeptides, which takes into account the effect of peptide charge–dipole interactions on helix stability, is presented. This new theory incorporates these interactions in an expression that is parameterized on recently obtained experimental data on polypeptides for which electrostatic effects are known to influence helix content. Unlike previous two-state or multistate models, which are parameterized on protein x-ray data, the present theoretical treatment in independent of such protein data. The theoretical model is applied to a series of peptides derived from the C-peptide of ribonuclease A, which have been the object of recent spectroscopic studies. The new theoretical approach can account for most of the structural information derived from studies of these C-peptides, and for overall average helix probabilities that are close in magnitude to those observed for these polypeptides in solution. An application of this new formulation for the prediction of the locations of α-helices in globular proteins from their amino acid sequence is also presented.     

Examine the characterization of biofilm formation and inhibition by targeting SrtA mechanism in Bacillus subtilis: a combined experimental and theoretical study

A unified coarse-grained model of three major classes of biological molecules—proteins, nucleic acids, and polysaccharides—has been developed. It is based on the observations that the repeated units of biopolymers (peptide groups, nucleic acid bases, sugar rings) are highly polar and their charge distributions can be represented crudely as point multipoles. The model is an extension of the united residue (UNRES) coarse-grained model of proteins developed previously in our laboratory. The respective force fields are defined as the potentials of mean force of biomacromolecules immersed in water, where all degrees of freedom not considered in the model have been averaged out. Reducing the representation to one center per polar interaction site leads to the representation of average site–site interactions as mean-field dipole–dipole interactions. Further expansion of the potentials of mean force of biopolymer chains into Kubo’s cluster-cumulant series leads to the appearance of mean-field dipole–dipole interactions, averaged in the context of local interactions within a biopolymer unit. These mean-field interactions account for the formation of regular structures encountered in biomacromolecules, e.g., α-helices and β-sheets in proteins, double helices in nucleic acids, and helicoidally packed structures in polysaccharides, which enables us to use a greatly reduced number of interacting sites without sacrificing the ability to reproduce the correct architecture. This reduction results in an extension of the simulation timescale by more than four orders of magnitude compared to the all-atom representation. Examples of the performance of the model are presented.

Helical peptides with three pairs of Asp-Arg and Glu-Arg residues in different orientations and spacings.

The helix-stabilizing effects of repeating pairs of Asp-Arg and Glu-Arg residues have been characterized using a peptide system of the same design used earlier to study Glu-Lys (Marqusee, S. & Baldwin, R.L., 1987, Proc. Natl. Acad. Sci. USA 84, 8898-8902) and Asp-Lys ion pairs (Marqusee, S. & Baldwin, R.L., 1990, In Protein Folding [Gierasch, L.M. & King, J., Eds.], pp. 85-94, AAAS, Washington, D.C.). The consequences of breaking ion pair and charge-helix dipole interactions by titration to pH 2 have been compared with the results of screening these interactions with NaCl at pH 7.0 and pH 2.5. The four peptides in each set contain three pairs of acidic (A) and basic (B) residues spaced either i, i + 4 or i, i + 3 apart. In one peptide of each kind the pairwise order of residues is AB, with the charges oriented favorably to the helix macrodipole, and in the other peptide the order is BA. The results are as follows: (1) Remarkably, both Asp-Arg and Glu-Arg peptides show the same pattern of helix stabilization at pH 7.0 found earlier for Glu-Lys and Asp-Lys peptides: i + 4 AB > i + 4 BA approximately i + 3 AB > i + 3 BA. (2) The ion pairs and charge-helix dipole interactions cannot be cleanly separated, but the results suggest that both interactions make important contributions to helix stability.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structure and Function of Membrane-Lytic Peptides

This article reviews the membrane interactions of a variety of peptides including alamethicin, melittin, cecropins, magainins, and defensins. The biological activities of the peptides are discussed and correlated to results from biophysical and structural studies. A picture emerges that allows one to understand the mechanisms of lysis and the regulation of the peptides' activities. Specific peptide–lipid interactions are particularly important in the case of antibiotic peptides, which affect the functionality of bacterial membranes, fungal membranes, or both but leave the bilayers of higher organisms, including those of the host cells, intact. Several models are presented and discussed in view of the ensemble of experimental data. These include the barrel stave, the wormhole, the carpet, and the detergent-like model.     

Positional dependence of the effects of negatively charged Glu side chains on the stability of two-stranded α-helical coiled-coils

The effects on protein stability of negatively charged Glu side chains at different positions along the length of the α-helix were investigated in the two-stranded α-helical coiled-coil. A native coiled-coil has been designed which consists of two identical 35 residue polypeptide chains with a heptad repeat Q_gV_aG_bA_cL_dQ_eK_f and a Cys residue at position 2 to allow the formation of an interchain 2-2′ disulphide bridge. This coiled-coil contains no intra- or interchain electrostatic interactions and served as a control for peptides in which Glu was substituted for Gln in the e or g heptad positions. The effect of the substitutions on stability was determined by urea denaturation at 20°C with the degree of unfolding monitored by circular dichroism spectroscopy. A Glu substituted for Gln near the N-terminus in each chain of the coiled-coil stabilizes the coiled-coil at pH 7, consistent with the charge–helix dipole interaction model. This stability increase is modulated by pH change and the addition of salt (KCl or guanidine hydrochloride), confirming the electrostatic nature of the effect. In contrast, Glu substitution in the middle of the helix destabilizes the coiled-coil because of the lower helical propensity and hydrophobicity of Glu compared with Gln at pH 7. Taking the intrinsic differences into account, the apparent charge–helix dipole interaction at the N-terminus is approximately 0.35 kcal/mol per Glu substitution. A Glu substitution at the C-terminus destabilizes the coiled-coil more than in the middle owing to the combined effects of intrinsic destabilization and unfavourable charge–helix dipole interaction with the negative pole of the helix dipole. The estimated destabilizing charge–helix dipole interaction of 0.08 kcal/mol is smaller than the stabilizing interaction at the N-terminus. The presence of a 2-2′disulphide bridge appears to have little influence on the magnitude of the charge–helix dipole interactions at either end of the coiled-coil. © 1997 European Peptide Society and John Wiley & Sons, Ltd.     

Side chain-backbone hydrogen bonding contributes to helix stability in peptides derived from an alpha-helical region of carboxypeptidase A

Recently, Presta and Rose proposed that a necessary condition for helix formation is the presence of residues at the N- and C-termini (called NTBs and CTBs) whose side chains can form hydrogen bonds with the initial four amides and the last four carbonyls of the helix, which otherwise lack intrahelical hydrogen bonding partners. We have tested this hypothesis by conformational analysis by circular dichroism (CD) of a synthetic peptide corresponding to a region (171-188) of the protein carboxypeptidase A; in the protein, residues 174 to 186 are helical and are flanked by NTBs and CTBs. Since helix formation in this peptide may also be stabilized by electrostatic interactions, we have compared the helical content of the native peptide with that of several modified peptides designed to enable dissection of different contributions to helix stability. As expected, helix dipole interactions appear to contribute substantially, but we conclude that hydrogen bonding interactions as proposed by Presta and Rose also stabilize helix formation. To assist in comparison of different peptides, we have introduced two concentration-independent CD parameters which are sensitive probes of helix formation.     

Folding of a peptide corresponding to the alpha-helix in bovine pancreatic trypsin inhibitor

A short peptide corresponding to the alpha-helical region of BPTI shows partial folding in aqueous solution (pH 7) as judged by circular dichroism (CD). Folding is temperature and denaturant sensitive, and the peptide is monomeric. The difference CD spectrum, obtained from spectra at two temperatures, indicates that the peptide folds as an alpha-helix. Difference CD spectroscopy provides a sensitive assay for helix formation in peptides exhibiting small amounts of structure. Helix stability in this peptide shows a marked pH dependence which is consistent with stabilizing charged side-chain interactions with the helix dipole and/or salt bridge formation.     

Circular dichroism studies of α-aminoisobutyric acid-containing peptides: Chain length and solvent effects in alternating Aib-L-Ala and Aib-L-Val sequences

The CD spectra of the peptides Boc-X-(Aib-X)_n-OMe (n = 1, 2, 3) and Boc-(Aib-X)₅-OMe, where X = L -Ala or L -Val have been examined in several solvents. The X = Ala and Val peptides behave similarly in all solvents, suggesting that the Aib residues dominate the folding preferences of these peptides. The decapeptides adopt helical conformations in methanol and trifluoroethanol, with characteristic negative CD bands at 222 and 205 nm. In the heptapeptides, similar spectra with reduced intensities are observed. Comparison with nmr studies suggest that estimates of helical content in oligopeptides by CD methods may lead to erroneous conclusions. The pentapeptides yield solvent-dependent spectra indicative of conformational perturbations. Peptide association in dioxane results in an unusual spectrum with a single negative band at 210 nm for the decapeptides. Disaggregation is induced by the addition of methanol or water to dioxane solutions. Aggregation of the heptapeptides is less pronounced in dioxane, suggesting that a critical helix length may be necessary to promote association stabilized by helix dipole–dipole interactions.     

Myocardial electrodynamics during electrical systole]

A dipole equivalent generator of heart is considered, which represents a circuit containing an active, capacitive and inductive resistance. The model makes it possible to solve a direct diagnostic task for the electric systole. A differential equation for the integral electrical vector of the heart in the dipole approximation was solved. A model vector electrocardiogram and a linear electrocardiogram for the electrical heart parameters are presented.     

Liposome-mediated labeling of adrenocorticotropin fragments parallels their biological activity

To test our hypothesis that specific interactions of ACTH peptides with model lipid membranes reflect the biological importance of similar interactions on target cells, we investigated the liposome-mediated labeling of ACTH fragments with the extremely hydrophobic photolabel, 3-trifluoromethyl-3-(m-[125I]iodophenyl)diazirine. Correlations were found between the labeling rates and the agonistic and antagonistic potencies of the peptides for in vitro steroidogenesis and inhibition of a synaptosomal protein kinase. A model for the cross-reactivity between ACTH and opioid peptides is discussed.     

The RH 421 styryl dye induced,pore model-dependent modulation of antimicrobial peptides activity in reconstituted planar membranes

Background

Antimicrobial agents, with different pore-formation mechanisms, may be differently influenced by alteration of the dipolar electric field of a lipid membrane.

Methods

By using electrophysiological measurements on reconstituted lipid membranes, we used alamethicin, melittin and magainin to report on how controlled manipulation of the membrane dipole potential by the styrylpyridinium dye RH 421 affects the kinetic and transport features of peptides within membranes.

Results

Our data demonstrate that the increase of the membrane dipole potential caused by RH 421 decreases the activity and single-channel conductance of alamethicin. Surprisingly, we found that RH 421 increases the activity of melittin and magainin, suggesting that RH 421 may contribute via electrostatic repulsions, among others, to an increase in the monolayer spontaneous curvature of the membrane. We propose that RH 421-induced dipole potential and membrane elasticity changes alter the peptide-induced channel dynamics, and the prevalence of one mechanism over the other for particular classes of peptides is dictated by the electrical and mechanical interactions which rule the pore-formation mechanism of such peptides.

General significance

These results point to a novel paradigm in which electrical and mechanical effects promoted by chemicals which preferentially alter the electrostatics of the membrane, may be employed to help distinguish among various pore-formation mechanisms of membrane-permeabilizing peptides.     

Electrostatic interactions in gramicidin channels

A model based on the solution of the electrostatic potential for a geometry of three dielectric regions associated with a gramicidin A channel (GA) is presented. The model includes a cylindrical dielectric layer to represent the peptide backbone and dipole rings to account for dipolar side chains. Image potential and dipolar contributions for different orientations and positions along the channel are analyzed. The conductance of GA and two analogues obtained by substituting the amino acid at position 1 are studied. The numerical simulation reproduces experimental results (Barrett et al. 1986, Biophys J 49, 673–686) and supports the idea that electrostatic dipole-ion interactions are of primary importance in gramicidin channel function. Correspondence to: G. Martinez     

A theoretical model for the association of amphiphilic transmembrane peptides in lipid bilayers

A theoretical model is proposed for the association of trans-bilayer peptides in lipid bilayers. The model is based on a lattice model for the pure lipid bilayer, which accounts accurately for the most important conformational states of the lipids and their mutual interactions and statistics. Within the lattice formulation the bilayer is formed by two independent monolayers, each represented by a triangular lattice, on which sites the lipid chains are arrayed. The peptides are represented by regular objects, with no internal flexibility, and with a projected area on the bilayer plane corresponding to a hexagon with seven lattice sites. In addition, it is assumed that each peptide surface at the interface with the lipid chains is partially hydrophilic, and therefore interacts with the surrounding lipid matrix via selective anisotropic forces. The peptides would therefore assemble in order to shield their hydrophilic residues from the hydrophobic surroundings. The model describes the self-association of peptides in lipid bilayers via lateral and rotational diffusion, anisotropic lipid-peptide interactions, and peptide-peptide interactions involving the peptide hydrophilic regions. The intent of this model study is to analyse the conditions under which the association of trans-bilayer and partially hydrophilic peptides (or their dispersion in the lipid matrix) is lipid-mediated, and to what extent it is induced by direct interactions between the hydrophilic regions of the peptides. The model properties are calculated by a Monte Carlo computer simulation technique within the canonical ensemble. The results from the model study indicate that direct interactions between the hydrophilic regions of the peptides are necessary to induce peptide association in the lipid bilayer in the fluid phase. Furthermore, peptides within each aggregate are oriented in such a way as to shield their hydrophilic regions from the hydrophobic environment. The average number of peptides present in the aggregates formed depends on the degree of mismatch between the peptide hydrophobic length and the lipid bilayer hydrophobic thickness: The lower the degree of mismatch is the higher this number is. Received: 30 December 1996 / Accepted: 9 May 1997     

A model of the electric field of the brain at EEG and microwave frequencies

A simple calculation of the current dipole moment of the extracellular electric field of the cortex is proposed; it is based on the dipole layer model. The model is extended to the range of microwave frequencies. Arguments in favor of emission of microwave radiation by the dendritic membranes of pyramidal neurons are presented and the strength of the radiative electric field at a distance from the head is calculated and discussed.     

The computer modelling of human TRH receptor, TRH and TRH-like peptides

The aim of this work was to verify the possibility of interactions between the human TRH receptor (an integral membrane protein which belongs to family 1 of G-protein coupled receptors) and TRH-like peptides presented in the prostate gland. These peptides are characterized by substitution of basic amino acid histidine (related to authentic TRH) for neutral or acidic amino acid, such as glutamic acid, phenylalanine, glutamine or tyrosine. The physiological function of TRH-like peptides in peripheral tissues is not precisely known. However, according to our recent experiments, we assume the existence of a local hormonal network formed by TRH-like peptides and TSH in the prostate gland. The network can be associated with circulating thyroid and steroid hormones, and may represent a new regulatory mechanism influencing the proliferative ability of prostatic tissue. A similar network of authentic TRH and TSH was already found in the gastrointestinal tract. The experimentally determined 3D-structures of human TRH receptor (hTRHr) and TRH-like peptides are not available. From this point of view we used de novo modeling procedures of G-protein coupled receptors on an automated protein modeling server used at the Glaxo Wellcome Experimental Research (Geneva, Switzerland). 3D-structures of TRH-like peptides were determined with a computer program CORINA (written by the team of J. Gasteiger, Computer-Chemie-Centrum and Institute for Organic Chemistry, University of Erlangen-Nurenberg, Germany). The generated PDB files with 3D-coordinates were visualized with Swiss-Pdb Viewer Release 3.51 (Glaxo Wellcome). From recent results it is evident that polar amino acids belonging to the extracellular terminus of hTRHr transmembrane regions can participate in interactions between TRH and hTRHr. There is no direct evidence that TRH-like peptides interact with the presented hTRHr model. On the contrary, with respect to the similar 3D-shape and the identity of terminal amino acids, it appears that these interactions are highly probable as well as the nearly 100 % cross-reactions between TRH or TRH-like peptides and antibody specific against authentic TRH. Closed terminal amino acids (pyroglutamic acid and proline-amide) of TRH or TRH-like peptides are important for these interactions. Desamido-TRH or glutamyl metabolites will be repelled by the negative potential of ASP195 (E: D93) and GLU298 (G: E137).     

Studies of synthetic helical peptides using circular dichroism and nuclear magnetic resonance

We have designed a set of 17-residue synthetic peptides to be monomeric helices in aqueous solution. Circular dichrosim experiments indicate the presence of helical structure in aqueous solution at low temperature and low pH. The two-dimensional nuclear magnetic resonance results for one of the peptides show a segment of ten residues which clearly meets all of the criteria for the existence of helical structure at both 5 degrees C and 15 degrees C. The first four residues of the peptide are in a largely extended conformation. Calculations suggest that residues 5 through 14 are significantly helical at 5 degrees C. When the temperature is increased, circular dichroism spectra indicate that the helical content decreases. At 15 degrees C, the 3JN alpha coupling constants increase in the helical region, indicating an increase in motion or conformational averaging in the helical segment. None of the peptides has pH titration behavior consistent with salt bridge stabilization of helical conformation. Our data lend themselves to interpretation with the helix dipole model and specific side-chain interactions. When the N and C termini charges are removed the helical content of the peptides increases. The amount of helicity increases as the pH is lowered, due to the ionization of His16. Much of the helical stabilization appears to be due to a specific side-chain interaction between His16 and Tyr12.     

Correlation between computed conformational properties of cytochrome c peptides and their antigenicity in a T-lymphocyte proliferation assay

By means of conformational energy calculations, we previously showed that the antigenic strength of a series of oligopeptides (derived from the carboxyl terminal sequence of cytochrome c) in a T-lymphocyte proliferation assay depends on their ability to adopt the α-helix conformation. Using experimentally determined statistical weights (within the framework of the Zimm–Bragg theory for the helix–coil transition), here we present a simple free energy analysis of the ability of these peptides to adopt the α-helix conformation in water. The experimental statistical weights have been modified to include the effect of long-range charge–dipole interactions on helix stability. We find that there is a close correlation between the tendency of a peptide to adopt the α-helix conformation and its ability to stimulate antigen-primed T cells. The shortest peptide with a tendency to adopt the α-helix conformation is also the shortest one that exhibits antigenic activity. The rapid and simple method presented here can thus be used to predict relative antigenicities for different peptides derived from cytochrome c.