Position-dependent interactions between cysteine residues and the helix dipole

A protein model was developed for studying the interaction between cysteine residues and the helix dipole. Site-directed mutagenesis was used to introduce cysteine residues at the N-terminus of helix H in recombinant sperm whale myoglobin. Based on the difference in thiol pK(a) between folded proteins and an unfolded peptide, the energy of interaction between the thiolate and the helix dipole was determined. Thiolates at the N1 and N2 positions of the helix were stabilized by 0.3 kcal/mole and 0.7 kcal/mole, respectively. A thiolate at the Ncap position was stabilized by 2.8 kcal/mole, and may involve a hydrogen bond. In context with other studies, an experimentally observed helix dipole effect may be defined in terms of two distinct components. A charge-dipole component involves electrostatic interactions with peptide bond dipoles in the first two turns of the helix and affects residues at all positions of the terminus; a hydrogen bond component involves one or more backbone amide groups and is only possible at the capping position due to conformational restraints elsewhere. The nature and magnitude of the helix dipole effect is, therefore, position-dependent. Results from this model system were used to interpret cysteine reactivity in rodent hemoglobins and the thioredoxin family.     

Electrostatic field effects on membrane domain segregation and on lateral diffusion

Natural membranes are organized structures of neutral and charged molecules bearing dipole moments which generate local non-homogeneous electric fields. When subjected to such fields, the molecules experience net forces that can modify the lipid and protein organization, thus modulating cell activities and influencing (or even dominating) the biological functions. The energetics of electrostatic interactions in membranes is a long-range effect which can vary over distance within r⁻¹ to r⁻³. In the case of a dipole interacting with a plane of dipoles, e.g. a protein interacting with a lipid domain, the interaction is stronger than two punctual dipoles and depends on the size of the domain. In this article, we review several contributions on how electrostatic interactions in the membrane plane can modulate the phase behavior, surface topography and mechanical properties in monolayers and bilayers.     

Helix-capping interaction in λ cro protein: A free energy simulation analysis

The stability mutant Tyr-26 → Asp was studied in the Cro protein from bacteriophage λ using free energy molecular dynamics simulations. The mutant was calculated to be more stable than the wild type by 3.0 ± 1.7 kcal/mol/monomer, in reasonable agreement with experiment (1.4 kcal/mol/monomer). Moreover, the aspartic acid in the mutant was found to form a capping interation with the amino terminus of the third α-helix of Cro. The simulations were analyzed to understand better the source of the stability of this helix-capping interaction and to examine the results in light of previous explanations of stabilizing helix caps-namely, a model of local unsatisfied hydrogen bonds at the helix termini and the helix macro dipole model. Analysis of the simulations shows that the stabilizing effect of this charged helical cap is due both to favorable hydrogen bonds with backbone NH groups at the helix terminus and to favorable electrostatic interactions (but not hydrogen bonds) with their carbonyls (effectively the next row of local dipoles in the helix). However, electrostatic interactions are weak or negligible with backbone dipolar groups in the helix further away from the terminus. Moreover, the importance of other local electrostatic interactions with polar side chains near the helix terminus, which are neglected in most treatments of this effect, are shown to be important. Thus, the results support a model that is intermediate between the two previous explanations: both unsatisfied hydrogen bonds at the helix terminus and other, local preoriented dipolar groups stabilize the helix cap. These findings suggest that similar interactions with preoriented dipolar groups may be important for cooperativity in other charge–dipole interactions and may be employed to advantage for molecular design. © 1994 Wiley-Liss, Inc.     

Solution NMR structures of the antimicrobial peptides phylloseptin-1, -2, and -3 and biological activity: the role of charges and hydrogen bonding interactions in stabilizing helix conformations

Phylloseptins are antimicrobial peptides of 19-20 residues which are found in the skin secretions of the Phyllomedusa frogs that inhabit the tropical forests of South and Central Americas. The peptide sequences of PS-1, -2, and -3 carry an amidated C-terminus and they exhibit 74% sequence homology with major variations of only four residues close to the C-terminus. Here we investigated and compared the structures of the three phylloseptins in detail by CD- and two-dimensional NMR spectroscopies in the presence of phospholipid vesicles or in membrane-mimetic environments. Both CD and NMR spectroscopies reveal a high degree of helicity in the order PS-2> or =PS-1>PS-3, where the differences accumulate at the C-terminus. The conformational variations can be explained by taking into consideration electrostatic interactions of the negative ends of the helix dipoles with potentially cationic residues at positions 17 and 18. Whereas two are present in the sequence of PS-1 and -2 only one is present in PS-3. In conclusion, the antimicrobial phylloseptin peptides adopt alpha-helical conformations in membrane environments which are stabilized by electrostatic interactions of the helix dipole as well as other contributions such hydrophobic and capping interactions.     

The pore helix dipole has a minor role in inward rectifier channel function

Ion channels lower the energetic barrier for ion passage across cell membranes and enable the generation of bioelectricity. Electrostatic interactions between permeant ions and channel pore helix dipoles have been proposed as a general mechanism for facilitating ion passage. Here, using genetic selections to probe interactions of an exemplar potassium channel blocker, barium, with the inward rectifier Kir2.1, we identify mutants bearing positively charged residues in the potassium channel signature sequence at the pore helix C terminus. We show that these channels are functional, selective, resistant to barium block, and have minimally altered conductance properties. Both the experimental data and model calculations indicate that barium resistance originates from electrostatics. We demonstrate that potassium channel function is remarkably unperturbed when positive charges occur near the permeant ions at a location that should counteract pore helix electrostatic effects. Thus, contrary to accepted models, the pore helix dipole seems to be a minor factor in potassium channel permeation.     

The alpha helix dipole: screened out?

Aligned alpha helix peptide dipoles sum to a "macroscopic" dipole parallel to the helix axis that has been implicated in protein folding and function. However, in aqueous solution the dipole is counteracted by an electrostatic reaction field generated by the solvent, and the strength of the helix dipole may reduce drastically from its value in vacuum. Here, using atomic-detail helix models and Poisson-Boltzmann continuum electrostatics calculations, the net effective dipole moment, mu(eff), is calculated. Some initially surprising results are found. Whereas in vacuum mu(eff) increases with helix length, the opposite is found to be the case for transmembrane helices. In soluble proteins, mu(eff) is found to vary strongly with the orientation and position of the helix relative to the aqueous medium. A set of rules is established to estimate of the strength of mu(eff) from graphical inspection of protein structures.     

Analysis of the interaction between charged side chains and the alpha-helix dipole using designed thermostable mutants of phage T4 lysozyme

It was shown previously that the introduction of a negatively charged amino acid at the N-terminus of an alpha-helix could increase the thermostability of phage T4 lysozyme via an electrostatic interaction with the "helix dipole" [Nicholson, H., Becktel, W. J., & Matthews, B. W. (1988) Nature 336, 651-656]. The prior report focused on the two stabilizing substitutions Ser 38----Asp (S38D) and Asn 144----Asp (N144D). Two additional examples of stabilizing mutants, T109D and N116D, are presented here. Both show the pH-dependent increase in thermal stability expected for the interaction of an aspartic acid with an alpha-helix dipole. Control mutants were also constructed to further characterize the nature of the interaction with the alpha-helix dipole. High-resolution crystal structure analysis was used to determine the nature of the interaction of the substituted amino acids with the end of the alpha-helix in both the primary and the control mutants. Control mutant S38N has stability essentially the same as that of wild-type lysozyme but hydrogen bonding similar to that of the stabilizing mutant S38D. This confirms that it is the electrostatic interaction between Asp 38 and the helix dipole, rather than a change in hydrogen-bonding geometry, that gives enhanced stability. Structural and thermodynamic analysis of mutant T109N provide a similar control for the stabilizing replacement T109D. In the case of mutant N116D, there was concern that the enhanced stability might be due to a favorable salt-bridge interaction between the introduced aspartate and Arg 119, rather than an interaction with the alpha-helix dipole. The additivity of the stabilities of N116D and R119M seen in the double mutant N116D/R119M indicates that favorable interactions are largely independent of residue 119. As a further control, Asp 92, a presumed helix-stabilizing residue in wild-type lysozyme, was replaced with Asn. This decreased the stability of the protein in the manner expected for the loss of a favorable helix dipole interaction. In total, five mutations have been identified that increase the thermostability of T4 lysozyme and appear to do so by favorable interactions with alpha-helix dipoles. As measured by the pH dependence of stability, the strength of the electrostatic interaction between the charged groups studied here and the helix dipole ranges from 0.6 to 1.3 kcal/mol in 150 mM KCl. In the case of mutants S38D and N144H, NMR titration was used to measure the pKa's of Asp 38 and His 144 in the folded structures.(ABSTRACT TRUNCATED AT 400 WORDS)     

Dipole lattice membrane model for protein calculations

A dipole lattice model for lipid membranes and their interactions with peptides is presented. It uses the Langevin dipole method to calculate electrostatic interactions in the heterogeneous membrane environment. A series of test cases are presented, including spherical charges, dipoles, side chain analogs, and helical peptides. The model consistently produces qualitatively correct results.     

Electrostatics significantly affect the stability of designed homeodomain variants

The role of electrostatic interactions in determining the stability of designed proteins was studied by constructing and analyzing a set of designed variants of the Drosophila engrailed homeodomain. Computational redesign of 29 surface positions results in a 25-fold mutant with moderate stability, similar to the wild-type protein. Incorporating helix dipole and N-capping considerations into the design algorithm by restricting amino acid composition at the helix termini and N-capping positions yields a ninefold mutant of the initial design (a 23-fold mutant of wild-type) that is over 3 kcal mol(-1) more stable than the protein resulting from the unbiased design. Four additional proteins were constructed and analyzed to isolate the effects of helix dipole and N-capping interactions in each helix. Based on the results of urea-denaturation experiments and calculations using the finite difference Poisson-Boltzmann method, both classes of interaction are found to increase the stability of the designed proteins significantly. The simple electrostatic model used in the optimization of rotamers by iterative techniques (ORBIT) force-field, which is similar to the electrostatic models used in other protein design force-fields, is unable to predict the experimentally determined stabilities of the designed variants. The helix dipole and N-capping restrictions provide a simple but effective method to incorporate two types of electrostatic interactions that impact protein stability significantly.     

Effects of aligned α‐helix peptide dipoles on experimental electrostatic potentials

Aligned protein α‐helix dipoles have been implicated in protein function and structure. The recent breakthroughs in high‐resolution electron microscopy (EM) of macromolecules makes it possible to explore fundamental aspects of structural biology at the detailed molecular level. The electrostatic potential (ESP) generated by aligned protein α‐helix dipole should be observable in high‐resolution EM maps despite the fact that the effect may be partially screened by induced electric fields. Here, we show that aligned backbone dipoles in protein α‐helices account for long‐range features in the protein ESP functions. Our results are consistent with experimental EM maps and density functional theory calculations, including direct Fourier summation for proper calculation of the ESP due to the nonlocal nature of the ESP function from aligned dipoles and other partial atomic charges.     

Electrostatic potentials of the alpha helix dipole and of elastase

Molecular graphics has been used to display the electrostatic potentials of the α-helix dipole and that of elastase calculated using atomic charges obtained by a new, simple method^1–3. Calculations on the α-helix dipole support the simple dipole model in which the helix is represented by single, half integral charges at the helix termini. The potentials of elastase show some interesting features which may be related to the binding processes.     

Orientational oscillations of the peptide groups of an α-helix

The low-energy orientational oscillations of the peptide groups of an -helix are considered and the value of the frequency is estimated to be in agreement with experiments. Approximate formulae are derived for the projection of a dipole moment on the helix axis and for the helix parameters. Within the framework of a three-chain model, the asymptotics of the soliton solution is obtained using a discrete approach.The analysis of -helix geometry exhibits two types of low-frequency oscillations of the -helix. The first one is connected with atom movements along the helix axis with the peptide groups twisting around the helix axis. Accordingly, it changes the hydrogen bond lengths between neighbouring peptide groups. In the second case, the slopes of the peptide groups to the helix axis oscillate without the helix parameters changing. Here, the energy of interactions between peptide-group dipoles is changed and, as a result, the oscillations have an optical nature. The frequency of the optical orientational oscillations is approximately 100 cm^-1.     

Further studies of the helix dipole model: effects of a free alpha-NH3+ or alpha-COO- group on helix stability

Interactions between the alpha-helix peptide dipoles and charged groups close to the ends of the helix were found to be an important determinant of alpha-helix stability in a previous study. The charge on the N-terminal residue of the C-peptide from ribonuclease A was varied chiefly by changing the alpha-NH2 blocking group, and the correlation of helix stability with N-terminal charge was demonstrated. An alternative explanation for some of those results is that the succinyl and acetyl blocking groups stabilize the helix by hydrogen bonding to an unsatisfied main-chain NH group. The helix dipole model is tested here with peptides that contain either a free alpha-NH3+ or alpha-COO- group, and no other charged groups that would titrate with similar pKa's. This model predicts that alpha-NH3+ and alpha-COO- groups are helix-destabilizing and that the destabilizing interactions are electrostatic in origin. The hydrogen bonding model predicts that alpha-NH3+ and alpha-COO- groups are not themselves helix-destabilizing, but that an acetyl or amide blocking group at the N- or C-terminus, respectively, stabilizes the helix by hydrogen bonding to an unsatisfied main-chain NH or CO group. The results are as follows: (1) Removal of the charge from alpha-NH3+ and alpha-COO- groups by pH titration stabilizes an alpha-helix. (2) The increase in helix stability on pH titration of these groups is close to the increase produced by adding an acetyl or amide blocking group.(ABSTRACT TRUNCATED AT 250 WORDS)     

Helix-helix interactions in lipid bilayers.

Using a continuum model, we calculated the electrostatic interaction free energy between two alpha-helices in three environments: the aqueous phase, a low dielectric alkane phase, and a simple representation of a lipid bilayer. As was found in previous work, helix-helix interactions in the aqueous phase are quite weak, because of solvent screening, and slightly repulsive, because of desolvation effects that accompany helix assembly. In contrast, the interactions can be quite strong in a hypothetical alkane phase because desolvation effects are essentially nonexistent and because helix-helix interactions are not well screened. In this type of environment, the antiparallel helix orientation is strongly favored over the parallel orientation. In previous work we found that the free energy penalty associated with burying helix termini in a bilayer is quite high, which is why the termini tend to protrude into the solvent. Under these conditions the electrostatic interaction is strongly screened by solvent; indeed, it is sufficient for the termini to protrude a few angstroms from the two surfaces of the bilayer for their interaction to diminish almost completely. The effect is consistent with the classical model of the helix dipole in which the dipole moment is represented by point charges located at either terminus. Our results suggest, in agreement with previous models, that there is no significant nonspecific driving force for helix aggregation and, hence, that membrane protein folding must be driven by specific interactions such as close packing and salt-bridge and hydrogen bond formation.     

Molecular dipole effects on tuning electron transfer in a porphine–quinone complex: a DFT and TDDFT study

The effect of a strong electric field generated by molecular dipoles on the ground state electronic structure and the Q and B states as well as the lowest charge transfer (CT) excited state of porphine–2,5-dimethyl-1,4-benzoquinone (PQ) complex has been investigated theoretically. Density functional theory DFT and time-dependent DFT (TDDFT) with the BH&HLYP hybrid functional have been applied in these calculations. The molecular dipole effect was generated by imposing one or two helical homopeptides consisting of eight α-aminoisobutyric acid residues (Aib₈) close to the PQ complex. The molecular dipoles in a close proximity to the PQ complex expose it to an electric field of the order of magnitude of 10⁹ V/m. The presence of the ambient molecular dipoles affects mainly the energy of the lowest CT state and barely the energies of the Q and B states. The molecular dipoles affect the energies of the excited states in a similar way as an external electrostatic field. Hence, the electric field induced by the molecular dipoles of the helical peptides could be used analogously to the external electrostatic field to control electron transfer (ET) in the PQ complex.     

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.     

Interactions between tricyclic antidepressants and phospholipid bilayer membranes

Participation of electrostatic and other noncovalent interactions in the binding of tricyclic antidepressants (TCAs) to the lipid bilayers was estimated from pH-dependencies of imipramine, desipramine, amitriptyline and nortriptyline binding to the lipid bilayers prepared from different phospholipids, both electroneutral and acidic. The binding was studied using a radioligand binding assay. It was found that the membrane phospholipid composition and methylation of the acyl side chain of TCA has a decisive effect on participation of particular noncovalent interactions in the binding. Apparent high-affinity binding of TCAs to the phosphatidylcholine or phosphatidylethanolamine membranes are achieved mainly by incorporation of uncharged drug molecules into the hydrophobic core of the bilayers. Van der Waals forces and hydrophobic effect are responsible for this binding. Both charged and uncharged drug molecules bind to phosphatidylserine membranes, therefore coulomb- or ion-induced dipole interactions play a role in these binding. Different spatial distribution of charged residues within the interface causes different electrostatic interactions between charged TCAs and vesicles formed from phosphatidylserine and phosphatidylinositol. The data supports the hypothesis under which TCAs could have effect on affective disorders partially via binding to the lipid part of the membrane and following changes of lipid-protein interactions.     

Molecular dynamics simulations of water within models of ion channels.

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.     

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.