Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “electrostatics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12Linked to original sources
PubMed1981

Electrostatic interaction of cytochrome c with cytochrome c1 and cytochrome oxidase.

The reactions of horse heart cytochrome c with succinate-cytochrome c reductase and cytochrome oxidase were studied as a function of ionic strength using both spectrophotometric and oxygen electrode assay techniques. The kinetic parameter Vmax/Km for both reactions decreased very rapidly as the ionic strength was increased, indicating that electrostatic interactions were important to the reactions. A new semiempirical relationship for the electrostatic energy of interaction between cytochrome c and its oxidation-reduction partners was developed, in which specific complementary charge-pair interactions between lysine amino groups on cytochrome c and negatively charged carboxylate groups on the other protein are assumed to dominate the interaction. The contribution of individual cytochrome c lysine amino groups to the electrostatic interaction was estimated from the decrease in reaction rate caused by specific modification of the lysine amino groups by reagents that change the charge to 0 or -1. These estimates range from -0.9 kcal/mol for lysines immediately surrounding the heme crevice of cytochrome c to 0 kcal/mol for lysines well removed from the heme crevice region. The semiempirical relationship for the total electrostatic energy of interaction was in quantitative agreement with the experimental ionic strength dependence of the reaction rates when the parameters were based on the specific lysine modification results. The electrostatic energies of interaction between cytochrome c and its reductase and oxidase were nearly the same, providing additional evidence that the two reactions take place at similar sites on cytochrome c.

Animals↗
PubMed1996

Electrostatic properties deduced from refined structures of NADH-cytochrome b5 reductase and the other flavin-dependent reductases: pyridine nucleotide-binding and interaction with an electron-transfer partner.

Electrostatic properties on the protein surface were examined on the basis of the crystal structure of NADH-cytochrome b5 reductase refined to a crystallographic R factor of 0.223 at 2.1 A resolution and of the other three flavin-dependent reductases. A structural comparison of NADH-cytochrome b5 reductase with the other flavin-dependent reductases, ferredoxin-NADP+ reductase, phthalate dioxygenase reductase, and nitrate reductase, showed that the alpha/beta structure is the common motif for binding pyridine nucleotide. Although the amino acid residues associated with pyridine nucleotide-binding are not conserved, the electrostatic properties and the location of the pyridine nucleotide-binding pockets of NADH-requiring reductases were similar to each other. The electrostatic potential of the surface near the flavin-protruding side (dimethylbenzene end of the flavin ring) of NADH-cytochrome b5 reductase was positive over a wide area while that of the surface near the heme-binding site of cytochrome b5 was negative. This implied that the flavin-protruding side of NADH-cytochrome b5 reductase is suitable for interacting with its electron-transfer partner, cytochrome b5. This positive potential area is conserved among four flavin-dependent reductases. A comparison of the electron-transfer partners of four flavin-dependent reductases showed that there are significant differences in the distribution of electrostatic potential between inter-molecular and inter-domain electron-transfer reactions.

Binding Sites↗
PubMed1997

Simulation of electrostatic and hydrodynamic properties of Serratia endonuclease.

We analyze the electrostatic and hydrodynamic properties of a nuclease from the pathogenic gram-negative bacterium Serratia marcescens using finite-difference Poisson-Boltzmann methods for electrostatic calculations and a bead-model approach for diffusion coefficient calculations. Electrostatic properties are analyzed for the enzyme in monomeric and dimeric forms and also in the context of DNA binding by the nuclease. Our preliminary results show that binding of a double-stranded DNA dodecamer by nuclease causes an overall shift in the charge of the protein by approximately three units of elementary charge per monomer, resulting in a positively charged protein at physiologic pH. In these calculations, the free enzyme was found to have a negative (-1 e) charge per monomer at pH 7. The most dramatic shift in pKa involves His 89 whose pKa increases by three pH units upon DNA binding. This shift leads to a protonated residue at pH 7, in contrast to the unprotonated form in the free enzyme. DNA binding also leads to a decrease in the energetic distances between the most stable protonation states of the enzyme. Dimerization has no significant effect on the electrostatic properties of each of the monomers for both free enzyme and that bound to DNA. Results of hydrodynamic calculations are consistent with the dimeric form of the enzyme in solution. The computed translational diffusion coefficient for the dimer model of the enzyme is in very good agreement with measurements from light scattering experiments. Preliminary electrooptical calculations indicate that the dimer should possess a large dipole moment (approximately 600 Debye units) as well as substantial optical anisotropy (limiting reduced linear electric dichroism of about 0.3). Therefore, this system may serve as a good model for investigation of electric and hydrodynamic properties by relaxation electrooptical experiments.

Endonucleases↗
PubMed2006

Electrostatic properties of cowpea chlorotic mottle virus and cucumber mosaic virus capsids.

Electrostatic properties of cowpea chlorotic mottle virus (CCMV) and cucumber mosaic virus (CMV) were investigated using numerical solutions to the Poisson-Boltzmann equation. Experimentally, it has been shown that CCMV particles swell in the absence of divalent cations when the pH is raised from 5 to 7. CMV, although structurally homologous, does not undergo this transition. An analysis of the calculated electrostatic potential confirms that a strong electrostatic repulsion at the calcium-binding sites in the CCMV capsid is most likely the driving force for the capsid swelling process during the release of calcium. The binding interaction between the encapsulated genome material (RNA) inside of the capsid and the inner capsid shell is weakened during the swelling transition. This probably aids in the RNA release process, but it is unlikely that the RNA is released through capsid openings due to unfavorable electrostatic interaction between the RNA and capsid inner shell residues at these openings. Calculations of the calcium binding energies show that Ca(2+) can bind both to the native and swollen forms of the CCMV virion. Favorable binding to the swollen form suggests that Ca(2+) ions can induce the capsid contraction and stabilize the native form.

Binding Sites↗
PubMed2004

Use of protein charge ladders to study electrostatic interactions during protein ultrafiltration.

Recent studies have demonstrated the importance of electrostatic interactions in membrane systems, but there is still controversy about the underlying phenomena. Protein charge ladders, consisting of a set of chemical derivatives of a given protein that differ by single charge groups, were used to quantify the electrostatic interactions during protein ultrafiltration. Myoglobin charge ladders were generated by acylation, with the different derivatives analyzed simultaneously by capillary electrophoresis. Filtration experiments were performed using polyethersulfone and composite regenerated cellulose membranes, with the membrane charge determined from the streaming potential. As expected, the rejection increased as the protein became more heavily charged due to the increase in electrostatic repulsion. However, the transmission of the weakly charged myoglobin species increased dramatically at very low ionic strength. This increase in transmission was attributed to a shift in pH within the pore caused by hydrogen ion partitioning into the charged membrane. The sieving data were in good agreement with theoretical calculations accounting for the effects of this pH shift on the electrostatic interactions.

Computer Simulation↗
PubMed2004

A mixed mechanistic-electrostatic model to explain pH dependence of glycosyl hydrolase enzyme activity.

Glycosyl hydrolases are a vast group of enzymes that share a common topology at their active site with two acidic residues that are responsible for activity. In spite of their similarity, they exhibit a wide range of pH optima that must depend on other factors. Using structural and mechanistic knowledge about glycosyl hydrolases from families 7, 10, and 16, we have formulated a new mathematical model that can include not only the ionization behavior of the catalytic residues but also as many ionizable residues as desired in the active site. In addition, the model can incorporate electrostatic influences via acid dissociation equilibrium constants and chemical relationships such as hydrogen bonds. The results of the simulations indicate a clear shift in the pH dependence of activity for the enzymes only when a close interrelation (hydrogen bond) between the catalytic and auxiliary residues in the active site is taken into account. This explains the observations from mutagenesis studies that show this type of shift and cannot be explained by a purely electrostatic interaction theory. Moreover, the presence of the kind of chemical interaction proposed could provide stabilization of the activity in the presence of environmental, structural, pH and electrostatic variations. These findings and the implications for the design of new mutagenesis strategies are discussed. The results suggest a way to modify, via site-directed mutagenesis, the acid dissociation of one of the catalytic residues in the active site independently of the other, which could have clear advantages over the purely electrostatic modifications that usually affect both residues simultaneously.

Catalysis↗
PubMed2003

Electrostatic calculation of the substituent effect: an efficient test on isolated molecules.

The energy of a disubstituted molecule has often been approximated by simple electrostatic formulas that represent the substituents as poles or dipoles. Herein, we test this approach on a new model system that is more direct and more efficient than testing on acid-base properties. The energies of 27 1,4-derivatives of bicyclo[2.2.2]octane were calculated within the framework of the density functional theory at the B3LYP/6-311+G(d,p) level; interaction of the two substituents was evaluated in terms of isodesmic homodesmotic reactions. This interaction energy, checked previously on some experimental gas-phase acidities, was considered to be accurate and served as reference to test the electrostatic approximation. This approximation works well in the qualitative sense as far as the sign and the order of magnitude are concerned: beginning with the strongest interaction between two poles, a weaker interaction between pole and dipole, and the weakest between two dipoles. However, all the electrostatic calculations yield energies that are too small, particularly for weak interaction, and this fundamental defect is not remedied by some possible improvements. In particular, variation of the effective permittivity would require a physically impossible value less than unity. The explanation must lie in a more complex distribution of electron density than anticipated in the electrostatic model. It also follows that possible conclusions about the transmission of substituent effects "through space" have little validity.

Journal Article↗
PubMed1993

Analysis of molecular recognition: steric electrostatic and hydrophobic complementarity.

We discuss three important aspects of molecular recognition: steric, electrostatic and hydrophobic. Steric fit means that interacting atoms may not approach each other beyond their van der Waals radii and, simultaneously, crevices should be filled as densely as possible. Electrostatic fit requires the maximum ionic and polar (hydrogen bond or other) interaction between host and guest atoms while the hydrophobic fit corresponds to the association trend between apolar groups in an aqueous medium. Space-filling models, obtained by molecular graphics, illustrate steric complementarity while we use molecular electrostatic potentials (MEPs) and fields (MEFs) to investigate electrostatic and hydrophobic matching. Molecular regions with negative and positive MEPs attract and repel a positive probe charge, respectively, so we consider them as attracting each other. Furthermore we postulate that regions with MEFs of similar magnitude tend to associate more strongly than those with very different fields (similis simili gaudet principle). We apply the above rules to the study of complementarity in the trypsin-BPTI complex and in a crystalline association between styrene epoxide as guest and a camphor-based anthracene derivative as host. We discuss molecular similarity on the same footing as complementarity and give some examples on the application of the concept to the rationalization of relative strengths of trypsin inhibitors.

Allosteric Site↗
PubMed1993

Electrostatic stabilization in four-helix bundle proteins.

Charge substitutions generated by site-directed mutagenesis at the termini of adjacent anti-parallel alpha-helices in a four-helix bundle protein were used to determine a precise value for the contribution of indirect charge-charge interactions to overall protein stability, and to simulate the electrostatic effects of alpha-helix macrodipoles. Thermodynamic double mutant cycles were constructed to measure the interaction energy between such charges on adjacent anti-parallel helices in the four-helix bundle cytochrome b562 from Escherichia coli. Previously, theoretical calculations of helix macrodipole interactions using modeled four-helix bundle proteins have predicted values ranging over an order of magnitude from 0.2 to 2.5 kcal/mol. Our system represents the first experimental evidence for electrostatic interactions such as those between partial charges due to helix macrodipole charges. At the positions mutated, we have measured a favorable interaction energy of 0.6 kcal/mol between opposite charges simulating an anti-parallel helix pair. Pairs of negative or positive charges simulating a parallel orientation of helices produce an unfavorable interaction of similar magnitude. The interaction energies show a strong dependence upon ionic strength, consistent with an electrostatic effect. Indirect electrostatic contacts do appear to confer a limited stabilization upon the association of anti-parallel packing of helices, favoring this orientation by as much as 1 kcal/mol at 20 mM K phosphate.

Apoproteins↗
PubMed1994

Receptor binding properties of four-helix-bundle growth factors deduced from electrostatic analysis.

Hormones of the hematopoietin class mediate signal transduction by binding to specific transmembrane receptors. Structural data show that the human growth hormone (hGH) forms a complex with a homodimeric receptor and that hGH is a member of a class of hematopoietins possessing an antiparallel 4-alpha-helix bundle fold. Mutagenesis experiments suggest that electrostatic interactions may have an important influence on hormone-receptor recognition. In order to examine the specificity of hormone-receptor complexation, an analysis was made of the electrostatic potentials of hGH, interleukin-2 (IL-2), interleukin-4 (IL-4), granulocyte colony-stimulating factor (G-CSF), granulocyte-macrophage colony-stimulating factor (GM-CSF), and the hGH and IL-4 receptors. The binding surfaces of hGH and its receptor, and of IL-4 and its receptor, show complementary electrostatic potentials. The potentials of the hGH and its receptor display approximately 2-fold rotational symmetry because the receptor subunits are identical. In contrast, the potentials of GM-CSF and IL-2 lack such symmetry, consistent with their known high affinity for hetero-oligomeric receptors. Analysis of the electrostatic potentials supports a recently proposed hetero-oligomeric model for a high-affinity IL-4 receptor and suggests a possible new receptor binding mode for G-CSF; it also provides valuable information for guiding structural and mutagenesis studies of signal-transducing proteins and their receptors.

Electrochemistry↗
PubMed1995

Electrostatic potential distribution of the gene V protein from Ff phage facilitates cooperative DNA binding: a model of the GVP-ssDNA complex.

The crystal structure of the gene V protein (GVP) from the Ff filamentous phages (M13, fl, fd) has been solved for the wild-type and two mutant (Y41F and Y41H) proteins at high resolution. The Y41H mutant crystal structure revealed crystal packing interactions, which suggested a plausible scheme for constructing the polymeric protein shell of the GVP-single-stranded DNA (ssDNA) complex (Guan Y, et al., 1994, Biochemistry 33:7768-7778). The electrostatic potentials of the isolated and the cooperatively formed protein shell have been calculated using the program GRASP and they revealed a highly asymmetric pattern of the electrostatic charge distribution. The inner surface of the putative DNA-binding channel is positively charged, whereas the opposite outer surface is nearly neutral. The electrostatic calculation further demonstrated that the formation of the helical protein shell enhanced the asymmetry of the electrostatic distribution. A model of the GVP-ssDNA complex with the n = 4 DNA-binding mode could be built with only minor conformational perturbation to the GVP protein shell. The model is consistent with existing biochemical and biophysical data and provides clues to the properties of GVP, including the high cooperatively of the protein binding to ssDNA. The two antiparallel ssDNA strands form a helical ribbon with the sugar-phosphate backbones at the middle and the bases pointing away from each other. The bases are stacked and the Phe 73 residue is intercalated between two bases. The optimum binding to a tetranucleotide unit requires the participation of four GVP dimers, which may explain the cooperativity of the GVP binding to DNA.

Amino Acid Sequence↗
PubMed1998

Role of the electrostatic loop charged residues in Cu,Zn superoxide dismutase.

We have expressed and characterized a mutant of Xenopus laevis Cu,Zn superoxide dismutase in which four highly conserved charged residues belonging to the electrostatic loop have been replaced by neutral side chains: Lys120 --> Leu, Asp130 --> Gln, Glu131 --> Gln, and Lys134 --> Thr. At low ionic strength, the mutant enzyme is one of the fastest superoxide dismutases ever assayed (k = 6.7 x 10(9) M(-1) s(-1), at pH 7 and mu = 0.02 M). Brownian dynamics simulations give rise to identical enzyme-substrate association rates for both wild-type and mutant enzymes, ruling out the possibility that enhancement of the activity is due to pure electrostatic factors. Comparative analysis of the experimental catalytic rate of the quadruple and single mutants reveals the nonadditivity of the mutation effects, indicating that the hyperefficiency of the mutant is due to a decrease of the energy barrier and/or to an alternative pathway for the diffusion of superoxide within the active site channel. At physiological ionic strength the catalytic rate of the mutant at neutral pH is similar to that of the wild-type enzyme as it is to the catalytic rate pH dependence. Moreover, mutation effects are additive. These results show that, at physiological salt conditions, electrostatic loop charged residues do not influence the diffusion pathway of the substrate and, if concomitantly neutralized, are not essential for high catalytic efficiency of the enzyme, pointing out the role of the metal cluster and of the invariant Arg141 in determining the local electrostatic forces facilitating the diffusion of the substrate towards the active site.

Animals↗
PubMed2003

Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials.

A continuum electrostatics approach for molecular dynamics (MD) simulations of macromolecules is presented and analyzed for its performance on a peptide and a globular protein. The approach incorporates the screened Coulomb potential (SCP) continuum model of electrostatics, which was reported earlier. The model was validated in a broad set of tests some of which were based on Monte Carlo simulations that included single amino acids, peptides, and proteins. The implementation for large-scale MD simulations presented in this article is based on a pairwise potential that makes the electrostatic model suitable for fast analytical calculation of forces. To assess the suitability of the approach, a preliminary validation is conducted, which consists of (i) a 3-ns MD simulation of the immunoglobulin-binding domain of streptococcal protein G, a 56-residue globular protein and (ii) a 3-ns simulation of Dynorphin, a biological peptide of 17 amino acids. In both cases, the results are compared with those obtained from MD simulations using explicit water (EW) molecules in an all-atom representation. The initial structure of Dynorphin was assumed to be an alpha-helix between residues 1 and 9 as suggested from NMR measurements in micelles. The results obtained in the MD simulations show that the helical structure collapses early in the simulation, a behavior observed in the EW simulation and consistent with spectroscopic data that suggest that the peptide may adopt mainly an extended conformation in water. The dynamics of protein G calculated with the SCP implicit solvent model (SCP-ISM) reveals a stable structure that conserves all the elements of secondary structure throughout the entire simulation time. The average structures calculated from the trajectories with the implicit and explicit solvent models had a cRMSD of 1.1 A, whereas each average structure had a cRMSD of about 0.8A with respect to the X-ray structure. The main conformational differences of the average structures with respect to the crystal structure occur in the loop involving residues 8-14. Despite the overall similarity of the simulated dynamics with EW and SCP models, fluctuations of side-chains are larger when the implicit solvent is used, especially in solvent exposed side-chains. The MD simulation of Dynorphin was extended to 40 ns to study its behavior in an aqueous environment. This long simulation showed that the peptide has a tendency to form an alpha-helical structure in water, but the stabilization free energy is too weak, resulting in frequent interconversions between random and helical conformations during the simulation time. The results reported here suggest that the SCP implicit solvent model is adequate to describe electrostatic effects in MD simulation of both peptides and proteins using the same set of parameters. It is suggested that the present approach could form the basis for the development of a reliable and general continuum approach for use in molecular biology, and directions are outlined for attaining this long-term goal.

Amino Acid Sequence↗
PubMed2004

An electrostatic basis for the stability of thermophilic proteins.

Two factors provide key contributions to the stability of thermophilic proteins relative to their mesophilic homologues: electrostatic interactions of charged residues in the folded state and the dielectric response of the folded protein. The dielectric response for proteins in a "thermophilic series" globally modulates the thermal stability of its members, with the calculated dielectric constant for the protein increasing from mesophiles to hyperthermophiles. This variability results from differences in the distribution of charged residues on the surface of the protein, in agreement with structural and genetic observations. Furthermore, the contribution of electrostatic interactions to the stability of the folded state is more favorable for thermophilic proteins than for their mesophilic homologues. This leads to the conclusion that electrostatic interactions play an important role in determining the stability of proteins at high temperatures. The interplay between electrostatic interactions and dielectric response also provides further rationalization for the enhanced stability of thermophilic proteins with respect to cold-denaturation. Taken together, the distribution of charged residues and their fluctuations have been shown to be factors in modulating protein stability over the entire range of biologically relevant temperatures.

Amino Acid Sequence↗
PubMed2005

Reversal of some viral IL-6 electrostatic properties compared to IL-6 contributes to a loss of alpha receptor component recruitment.

Human interleukin-6 (hIL-6) is a pleiotropic mediator of activation and proliferation across a large number of different cell types. Human herpesvirus-8 (HHV-8) has been associated with classical and AIDS-related Kaposi's sarcoma (KS). HHV-8 encodes viral IL-6 (vIL-6), a functional homolog of human interleukin-6, that promotes the growth of KS and of some lymphoma cells. Signaling induced by human IL-6 requires recruitment of the glycoprotein gp130, which acts as the signal transducing chain, and of IL-6Ralpha, which is necessary for cognate recognition and high affinity receptor complex formation. In contrast, the formation of a functional complex between vIL-6 and gp130 does not require the presence of IL-6Ralpha. The physico-chemical properties of vIL-6 have been analyzed and compared to those of hIL-6 and of the receptor chains, gp130 and IL-6Ralpha. Interaction sites on vIL-6 involve more hydrophobic residues than those of hIL-6. The electrostatic fields induced by vIL-6 and IL-6Ralpha are repulsive and prevent interaction between vIL-6 and IL-6Ralpha, whereas the electrostatic field induced by hIL-6 steers the complex formation with IL-6Ralpha. Subsequently, electrostatic binding free energy in the vIL-6/IL-6Ralpha complex is destabilizing, whereas it is stabilizing in the complex comprising hIL-6. These properties result from charge reversals between viral and human IL-6, an unusual phenomenon of amino acid substitutions within a homologous protein family. This suggests a selection pressure for vIL-6 to by-pass the IL-6Ralpha control of host defense against virus infection. This selection pressure has yielded the reversal of electrostatic properties of vIL-6 when compared to hIL-6.

Amino Acid Sequence↗
PubMed1988

Calculation of the total electrostatic energy of a macromolecular system: solvation energies, binding energies, and conformational analysis.

In this report we describe an accurate numerical method for calculating the total electrostatic energy of molecules of arbitrary shape and charge distribution, accounting for both Coulombic and solvent polarization terms. In addition to the solvation energies of individual molecules, the method can be used to calculate the electrostatic energy associated with conformational changes in proteins as well as changes in solvation energy that accompany the binding of charged substrates. The validity of the method is examined by calculating the hydration energies of acetate, methyl ammonium, ammonium, and methanol. The method is then used to study the relationship between the depth of a charge within a protein and its interaction with the solvent. Calculations of the relative electrostatic energies of crystal and misfolded conformations of Themiste dyscritum hemerythrin and the VL domain of an antibody are also presented. The results indicate that electrostatic charge-solvent interactions strongly favor the crystal structures. More generally, it is found that charge-solvent interactions, which are frequently neglected in protein structure analysis, can make large contributions to the total energy of a macromolecular system.

Binding, Competitive↗
PubMed1989

Molecular dynamics effects on protein electrostatics.

Electrostatic calculations have been carried out on a number of structural conformers of tuna cytochrome c. Conformers were generated using molecular dynamics simulations with a range of solvent simulating, macroscopic dielectric formalisms, and one solvent model that explicitly included solvent water molecules. Structures generated using the lowest dielectric models were relatively tight, with side chains collapsed on the surface, while those from the higher dielectric models had more internal and external fluidity, with surface side chains exploring a fuller range of conformational space. The average structure generated with the explicitly solvated model corresponded most closely with the crystal structure. Individual pK values, overall titration curves, and electrostatic potential surfaces were calculated for average structures and structures along each simulation. Differences between structural conformers within each simulation give rise to substantial changes in calculated local electrostatic interactions, resulting in pK value fluctuations for individual sites in the protein that vary by 0.3-2.0 pK units from the calculated time average. These variations are due to the thermal side chain reorientations that produce fluctuations in charge site separations. Properties like overall titration curves and pH dependent stability are not as sensitive to side chain fluctuations within a simulation, but there are substantial effects between simulations due to marked differences in average side chain behavior. These findings underscore the importance of proper dielectric formalism in molecular dynamics simulations when used to generate alternate solution structures from a crystal structure, and suggest that conformers significantly removed from the average structure have altered electrostatic properties that may prove important in episodic protein properties such as catalysis.

Animals↗
PubMed1992

Continuum electrostatics of the C-peptide: anatomy of the problem.

A computational study of the role of all ionizable groups of the C-peptide in its helix-coil transition is performed within the framework of continuum electrostatics. The method employed in our computations involves a numeric solution of the Poisson equation with the Boundary Element Method. Our calculations correctly predict the experimentally observed trends in the helix-coil equilibrium of the C-peptide, and suggest that the mechanisms involved are more complex than usually presumed in the literature. Our results suggest that electrostatic interactions in the unfolded conformation are often more important than in the helix, total electrostatic contribution to the helix-coil transition due to the side chains of the C-peptide destabilizes the helix, changes in the helix stability produced by the changes in the ionization state of the side chains are dominated by side chain effects, the effect of the helix dipole on the energetics of the helix-coil transition of the C-peptide is either minor or similar to other contributions in magnitude; while the formation of a salt bridge is electrostatically favorable, formation of the hydrogen bond between a charged and a polar side chains is not. Factors limiting the accuracy of the computations are discussed.

C-Peptide↗