Disruption and formation of surface salt bridges are coupled to DNA binding by the integration host factor: a computational analysis

Revealing the thermodynamic driving force of protein-DNA interactions is crucial to the understanding of factors that dictate the properties and function of protein-DNA complexes. For the binding of DNA to DNA-wrapping proteins, such as the integration host factor (IHF), Record and co-workers proposed that the disruption of a large number of preexisting salt bridges is coupled with the binding process [Holbrook, J. A., et al. (2001) J. Mol. Biol. 310, 379]. To test this proposal, we have conducted explicit solvent MD simulations (multiple ～25-50 ns trajectories for each salt concentration) to examine the behavior of charged residues in IHF, especially concerning their ability to form salt bridges at different salt concentrations. Of the 17 cationic residues noted by Record and co-workers, most are engaged in salt bridge interactions for a significant portion of the trajectories, especially in the absence of salt. This observation suggests that, from a structural point of view, their proposal is plausible. However, the complex behaviors of charged residues observed in the MD simulations also suggest that the unusual thermodynamic characteristics of IHF-DNA binding likely arise from the interplay between complex dynamics of charged residues both in and beyond the DNA binding site. Moreover, a comparison of MD simulations at different salt concentrations suggests that the strong dependence of the IHF-DNA binding enthalpy on salt concentration may not be due to a significant decrease in the number of stable salt bridges in apo IHF at high salt concentrations. In addition to the Hofmeister effects quantified in more recent studies of IHF-DNA binding, we recommend consideration of the variation of the enthalpy change of salt bridge disruption at different salt concentrations. Finally, the simulation study presented here explicitly highlights the fact that the electrostatic properties of DNA-binding proteins can be rather different in the apo and DNA-bound states, which has important implications for the design of robust methods for predicting DNA binding sites in proteins.     

Role of salt bridges in homeodomains investigated by structural analyses and molecular dynamics simulations

Homeodomains are a class of helix-turn-helix DNA-binding protein motifs that play an important role in the control of cellular development in eukaryotes. They fold in a three alpha-helix structural module, where the third helix is the recognition helix that fits into the major groove of DNA. Structural analysis of the members of the homeodomain family led to the identification of interactions likely to stabilize the protein domains. Linking the helices pairwise, three salt bridges were found to be well preserved within the family. Also well conserved were two cation-pi interactions between aromatic and positively charged side chains. To analyze the structural role of the salt bridges, molecular dynamics simulations (MD) were carried out on the wild-type homeodomain from the Drosophila paired protein (1fjl) and on three mutants, which lack one or two salt bridges and mimic natural mutations in other homeodomains. Analysis of the trajectories revealed only small structural rearrangements of the three helices in all MD simulations, thereby suggesting that the salt bridges have no essential stabilizing role at room temperature, but rather might be important for improving thermostability. The latter hypothesis is supported by a good correlation between the melting midpoint temperatures of several homeodomains and the number of salt bridges and cation-pi interactions that connect secondary structures.     

Differences in electrostatic properties at antibody-antigen binding sites: implications for specificity and cross-reactivity

Antibodies HyHEL8, HyHEL10, and HyHEL26 (HH8, HH10, and HH26, respectively) recognize highly overlapping epitopes on hen egg-white lysozyme (HEL) with similar affinities, but with different specificities. HH8 binding to HEL is least sensitive toward mutations in the epitope and thus is most cross-reactive, HH26 is most sensitive, whereas the sensitivity of HH10 lies in between HH8 and HH26. Here we have investigated intra- and intermolecular interactions in three antibody-protein complexes: theoretical models of HH8-HEL and HH26-HEL complexes, and the x-ray crystal structure of HH10-HEL complex. Our results show that HH8-HEL has the lowest number and HH26-HEL has the highest number of intra- and intermolecular hydrogen bonds. The number of salt bridges is lowest in HH8-HEL and highest in HH26-HEL. The binding site salt bridges in HH8-HEL are not networked, and are weak, whereas, in HH26-HEL, an intramolecular salt-bridge triad at the binding site is networked to an intermolecular triad to form a pentad. The pentad and each salt bridge of this pentad are exceptionally stabilizing. The number of binding-site salt bridges and their strengths are intermediate in HH10-HEL, with an intramolecular triad. Our further calculations show that the electrostatic component contributes the most to binding energy of HH26-HEL, whereas the hydrophobic component contributes the most in the case of HH8-HEL. A "hot-spot" epitope residue Lys-97 forms an intermolecular salt bridge in HH8-HEL, and participates in the intermolecular pentad in the HH26-HEL complex. Mutant modeling and surface plasmon resonance (SPR) studies show that this hot-spot epitope residue contributes significantly more to the binding than an adjacent epitope residue, Lys-96, which does not form a salt bridge in any of the three HH-HEL complexes. Furthermore, the effect of mutating Lys-97 is most severe in HH26-HEL. Lys-96, being a charged residue, also contributes the most in HH26-HEL among the three complexes. The SPR results on these mutants also highlight that the apparent "electrostatic steering" on net on rates actually act at post-collision level stabilization of the complex. The significance of this work is the observed variations in electrostatic interactions among the three complexes. Our work demonstrates that higher electrostatics, both as a number of short-range electrostatic interactions and their contributions, leads to higher binding specificity. Strong salt bridges, their networking, and electrostatically driven binding, limit flexibilities through geometric constrains. In contrast, hydrophobic driven binding and low levels of electrostatic interactions are associated with conformational flexibility and cross-reactivity.     

Salt enhances calmodulin-target interaction

Calmodulin (CaM) operates as a Ca(2+) sensor and is known to interact with and regulate hundreds of proteins involved in a great many aspects of cellular function. It is of considerable interest to understand the balance of forces in complex formation of CaM with its target proteins. Here we have studied the importance of electrostatic interactions in the complex between CaM and a peptide derived from smooth-muscle myosin light-chain kinase by experimental methods and Monte Carlo simulations of electrostatic interactions. We show by Monte Carlo simulations that, in agreement with experimental data, the binding affinity between CaM and highly charged peptides is surprisingly insensitive to changes in the net charge of both the protein and peptide. We observe an increase in the binding affinity between oppositely charged partners with increasing salt concentration from zero to 100 mM, showing that formation of globular CaM-kinase type complexes is facilitated at physiological ionic strength. We conclude that ionic interactions in complex formation are optimized at pH and saline similar to the cell environment, which probably overrules the electrostatic repulsion between the negatively charged Ca(2+)-binding domains of CaM. We propose a conceivable rationalization of CaM electrostatics associated with interdomain repulsion.     

Contribution of arginine-glutamate salt bridges to helix stability

Peptide side chain interactions were studied by molecular dynamics simulation using explicit solvent on a peptide with the sequence AAARAAAAEAAEAAAARA. Three different protonation states of the glutamic acid side chains were simulated for four 20 ns runs each, a total simulation time of 240 ns. Two different salt bridge geometries were observed and the preferred geometry was found to depend on Glu — Arg residue spacing. Stable charge clusters were also observed, particularly in the fully charged peptide. Salt bridges were selectively interrupted upon protonation, with concomitant changes in secondary structure. The fully charged peptide was highly helical between residues 9 and 13, although protonation increased helicity near the N-terminus. The contribution of salt bridges to helix stability therefore depends on both position and relative position of charged residues within a sequence.     

Do salt bridges stabilize proteins? A continuum electrostatic analysis

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).     

Kinetic analysis of sequence-specific recognition of ssDNA by an autoantibody

11F8 is a pathogenic monoclonal anti-ssDNA autoantibody isolated from a lupus prone mouse. Previous studies established that 11F8 is sequence-specific and identified the thermodynamic and kinetic basis for the specific recognition of ssDNA, and binding site mutations of a single-chain construct reveal that (Y32)LCDR1, (R31)HCDR1, (W33)HCDR1, (R98)HCDR3, (L97)HCDR3, and (Y100)HCDR3 are responsible for approximately 80% of the binding free energy. Here we evaluate the role of these residues along with a group of basic residues (K62, K64, R24, K52) within the context of the binding mechanism. Binding of 11F8 takes place in two steps. In the first step, the overall positive charge of the antigen binding site attracts the negatively charged DNA to form an encounter complex that is stabilized by two salt bridges and a hydrogen bond. The second step is a slow process in which minor conformational changes occur. During this step, aromatic side chains become desolvated, presumably through stacking interactions involving two thymine bases within the DNA recognition epitope. Although the stability of the complex arises primarily from interactions formed in the second step, sequence specificity results from interactions with residues involved in both steps. These studies also show that the way in which 11F8 achieves high affinity sequence-specific binding is more closely related to RNA binding proteins than those that bind DNA and point to strategies for disrupting DNA binding that could prove to be therapeutically useful.     

Importance of context in protein folding: secondary structural propensities versus tertiary contact-assisted secondary structure formation

Molecular dynamics simulations can be used to reveal the detailed conformational behaviors of peptides and proteins. By comparing fragment and full-length protein simulations, we can investigate the role of each peptide segment in the folding process. Here, we take advantage of information regarding the helix formation process from our previous simulations of barnase and protein A as well as new simulations of four helical fragments from these proteins at three different temperatures, starting with both helical and extended structures. Segments with high helical propensity began the folding process by tethering the chain through side chain interactions involving either polar interactions, such as salt bridges, or hydrophobic staples. These tethers were frequently nonnative (i.e., not i --> i + 4 spacing) and provided a scaffold for other residues, thereby limiting the conformational search. The helical structure then propagated on both sides of the tether. Segments with low stability and propensity formed later in the folding process and utilized contacts with other portions of the protein when folding. These helices formed via a tertiary contact-assisted mechanism, primarily via hydrophobic contacts between residues distant in sequence. Thus, segments with different helical propensities appear to play different roles during protein folding. Furthermore, the active role of nonlocal side chains in helix formation highlights why we must move beyond simple hierarchical models of protein folding.     

Nonnative interactions in coupled folding and binding processes of intrinsically disordered proteins

Proteins function by interacting with other molecules, where both native and nonnative interactions play important roles. Native interactions contribute to the stability and specificity of a complex, whereas nonnative interactions mainly perturb the binding kinetics. For intrinsically disordered proteins (IDPs), which do not adopt rigid structures when being free in solution, the role of nonnative interactions may be more prominent in binding processes due to their high flexibilities. In this work, we investigated the effect of nonnative hydrophobic interactions on the coupled folding and binding processes of IDPs and its interplay with chain flexibility by conducting molecular dynamics simulations. Our results showed that the free-energy profiles became rugged, and intermediate states occurred when nonnative hydrophobic interactions were introduced. The binding rate was initially accelerated and subsequently dramatically decreased as the strength of the nonnative hydrophobic interactions increased. Both thermodynamic and kinetic analysis showed that disordered systems were more readily affected by nonnative interactions than ordered systems. Furthermore, it was demonstrated that the kinetic advantage of IDPs ("fly-casting" mechanism) was enhanced by nonnative hydrophobic interactions. The relationship between chain flexibility and protein aggregation is also discussed.     

Segregation of negatively charged phospholipids by the polycationic and farnesylated membrane anchor of Kras

The Kras protein, a member of the Ras family of bio-switches that are frequently mutated in cancer and developmental disorders, becomes functional when anchored to the inner surface of the plasma membrane. It is well known that membrane attachment involves the farnesylated and poylcationic C-terminus of the protein. However, little is known about the structure of the complex and the specific protein-lipid interactions that are responsible for the binding. On the basis of data from extensive (>0.55 μs) molecular dynamics simulations of multiple Kras anchors in bilayers of POPC/POPG lipids (4:1 ratio), we show that, as expected, Kras is tethered to the bilayer surface by specific lysine-POPG salt bridges and by nonspecific farnesyl-phospholipid van der Waals interactions. Unexpectedly, however, only the C-terminal five of the eight Kras Lys side chains were found to directly interact with the bilayer, with the N-terminal ones staying in water. Furthermore, the positively charged Kras anchors pull the negatively charged POPG lipids together, leading to the clustering of the POPG lipids around the proteins. This selective Kras-POPG interaction is directly related to the specific geometry of the backbone, which exists in two major conformational states: 1), a stable native-like ensemble of structures characterized by an extended geometry with a pseudohelical turn; and 2), less stable nonnative ensembles of conformers characterized by severely bent geometries. Finally, although the interface-bound anchor has little effect on the overall structure of the bilayer, it induces local thinning within a persistence length of ∼12 Å. Our results thus go beyond documenting how Kras attaches to a mixed bilayer of charged and neutral lipids; they highlight a fascinating process of protein-induced lipid sorting coupled with the (re)shaping of a surface-bound protein by the host lipids.     

Composition-dependent energetic contribution of complex salt bridges to collagen stability

Complex salt bridges, on which three or more charged residues interplay simultaneously, cannot be considered as addition of individual salt bridges. This is still an intriguing problem in protein folding and stability. Here, we used an obligated ABC-type collagen heterotrimer as a platform to study the relationship between energetic contributions and conformational details of three-body complex salt bridges anchored by positively charged residues, K and R. Eight complex salt bridges were constructed by engineering point mutations in the heterotrimer. The circular dichroism measurements showed that the K-anchored complex salt bridges were stronger than the R-anchored ones. The molecular dynamics simulation revealed that both types of salt bridges had distinct dynamic features. The energetic contribution of K-anchored salt bridges was mainly determined by strong single bridges. In the R-anchored complex salt bridges, both side-chain electrostatic interactions and side-chain-backbone hydrogen bonding were involved. An empirical equation was proposed to predict the energetic contributions with high accuracy (R² = 0.93). This work could help us take insights into the mechanisms of composition-dependent behaviors of the complex salt bridges on protein surface.     

Protein-protein interactions in actin-myosin binding and structural effects of R405Q mutation: a molecular dynamics study

Detailed residue-wise interactions involved in the binding of myosin to actin in the rigor conformation without nucleotides have been examined using molecular dynamics simulations of the chicken skeletal myosin head complexed with two actin monomers, based on the cryo-microscopic model of Holmes et al. (Nature 2003;425:423-427). The overall interaction is largely electrostatic in nature, because of the charged residues in the four loops surrounding the central primary binding site. The 50k/20k loop, disordered in crystal structures and in simulations of free myosin in solution, was found to be in a conformation stabilized with 1 - 2 internal salt bridges. The cardiomyopathy loop forms 2 - 3 interprotein salt bridges with actin monomers upon binding, whereas its Arg405 residue, the mutation site associated with the hypertrophic cardiomyopathy, forms a strong salt bridge with Glu605 in the neighboring helix away from actin in the actin-bound myosin. The myopathy loop of the R405Q mutant maintains a high degree of two-strand beta-sheet character when bound to actin with the corresponding salt bridges broken.     

Retention of Conformational Entropy upon Calmodulin Binding to Target Peptides Is Driven by Transient Salt Bridges

Calmodulin (CaM) is a highly flexible calcium-binding protein that mediates signal transduction through an ability to differentially bind to highly variable binding sequences in target proteins. To identify how binding affects CaM motions, and its relationship to conformational entropy and target peptide sequence, we have employed fully atomistic, explicit solvent molecular dynamics simulations of unbound CaM and CaM bound to five different target peptides. The calculated CaM conformational binding entropies correlate with experimentally derived conformational entropies with a correlation coefficient R² of 0.95. Selected side-chain interactions with target peptides restrain interhelical loop motions, acting to tune the conformational entropy of the bound complex via widely distributed CaM motions. In the complex with the most conformational entropy retention (CaM in complex with the neuronal nitric oxide synthase binding sequence), Lys-148 at the C-terminus of CaM forms transient salt bridges alternating between Glu side chains in the N-domain, the central linker, and the binding target. Additional analyses of CaM structures, fluctuations, and CaM-target interactions illuminate the interplay between electrostatic, side chain, and backbone properties in the ability of CaM to recognize and discriminate against targets by tuning its conformational entropy, and suggest a need to consider conformational dynamics in optimizing binding affinities.     

Cooperative helix stabilization by complex Arg-Glu salt bridges

Among the interactions that stabilize the native state of proteins, the role of electrostatic interactions has been difficult to quantify precisely. Surface salt bridges or ion pairs between acidic and basic side chains have only a modest stabilizing effect on the stability of helical peptides or proteins: estimates are roughly 0.5 kcal/mol or less. On the other hand, theoretical arguments and the occurrence of salt bridge networks in thermophilic proteins suggest that multiple salt bridges may exert a stronger stabilizing effect. We show here that triads of charged side chains, Arg(+)-Glu(-)-Arg(+) spaced at i,i+4 or i,i+3 intervals in a helical peptide stabilize alpha helix by more than the additive contribution of two single salt bridges. The free energy of the triad is more than 1 kcal/mol in excess of the sum of the individual pairs, measured in low salt concentration (10 mM). The effect of spacing the three groups is severe; placing the charges at i,i+4 or i,i+3 sites has a strong effect on stability relative to single bridges; other combinations are weaker. A conservative calculation suggests that interactions of this kind between salt bridges can account for much of the stabilization of certain thermophilic proteins.     

Insights into the Structure,Correlated Motions,and Electrostatic Properties of Two HIV-1 gp120 V3 Loops

The V3 loop of the glycoprotein 120 (gp120) is a contact point for cell entry of HIV-1 leading to infection. Despite sequence variability and lack of specific structure, the highly flexible V3 loop possesses a well-defined role in recognizing and selecting cell-bound coreceptors CCR5 and CXCR4 through a mechanism of charge complementarity. We have performed two independent molecular dynamics (MD) simulations to gain insights into the dynamic character of two V3 loops with slightly different sequences, but significantly different starting crystallographic structures. We have identified highly populated trajectory-specific salt bridges between oppositely charged stem residues Arg9 and Glu25 or Asp29. The two trajectories share nearly identical correlated motions within the simulations, despite their different overall structures. High occupancy salt bridges play a key role in the major cross-correlated motions in both trajectories, and may be responsible for transient structural stability in preparation for coreceptor binding. In addition, the two V3 loops visit conformations with similarities in spatial distributions of electrostatic potentials, despite their inherent flexibility, which may play a role in coreceptor recognition. It is plausible that cooperativity between overall electrostatic potential, charged residue interactions, and correlated motions could be associated with a coreceptor selection and binding.     

Experimental assignment of the structure of the transition state for the association of barnase and barstar

Association of a protein complex follows a two step reaction mechanism, with the first step being the formation of an encounter complex which evolves into the final complex. Here we present new experimental data for the association of the bacterial ribonuclease barnase and its polypeptide inhibitor barstar which shed light on the thermodynamics and structure of the transition state and preceding encounter complex of association at diminishing electrostatic attraction. We show that the activation entropy at the transition state is close to zero, with the activation enthalpy being equal to the free energy of binding. This observation was independent of the magnitude of the mutual electrostatic attraction, which were altered by mutagenesis or by addition of salt. The low activation entropy implies that the transition state is mostly solvated at all ionic strengths. The structure of the transition state was probed by measuring pairwise interaction energies using double-mutant-cycles. While at low ionic strength all proximal charge-pairs form contacts, at high salt only a subset of these interactions are maintained. More specifically, charge-charge interactions between partially buried residues are lost, while exposed charged residues maintain their ability to form specific interactions even at the highest salt concentration. Uncharged residues do not interact at any ionic strength. The results presented here suggest that the barnase-barstar binding sites are correctly aligned during the transition state even at diminishing electrostatic attraction, although specific short range interactions of uncharged residues are not yet formed. Furthermore, most of the interface desolvation (which contributes to the entropy of the system) has not yet occurred. This picture seems to be valid at low and high salt. However, at high salt, interactions of the activated complex are limited to a more restricted set of residues which are easier approached during diffusion, prior to final docking. This suggest that the steering region at high salt is more limited, albeit maintaining its specificity.     

Salt effects on hydrophobic‐core formation in folding of a helical miniprotein studied by molecular dynamics simulations

We have investigated effects of salt ions on folding events of a helical miniprotein chicken villin headpiece subdomain HP36. Low concentrations of ions alter electrostatic interactions between charged groups of a protein and can change the populations of conformers. Here, we compare two data sets of folding simulations of HP36 in explicit water solvent with or without ions. For efficient sampling of the conformational space of HP36, the multicanonical replica‐exchange molecular dynamics method was employed. Our analyses suggest that salt alters salt‐bridging nature of the protein at later stages of folding at room temperature. Especially, more nonnative, nonlocal salt bridges are formed at near‐native conformations in pure water. Our analyses also show that such salt‐bridge formation hinders the fully native hydrophobic‐core packing at the final stages of folding. Proteins 2014; 82:933–943. © 2013 Wiley Periodicals, Inc.     

Penetratin-membrane association: W48/R52/W56 shield the peptide from the aqueous phase

Using molecular dynamics simulations, we studied the mode of association of the cell-penetrating peptide penetratin with both a neutral and a charged bilayer. The results show that the initial peptide-lipid association is a fast process driven by electrostatic interactions. The homogeneous distribution of positively charged residues along the axis of the helical peptide, and especially residues K46, R53, and K57, contribute to the association of the peptide with lipids. The bilayer enhances the stability of the penetratin helix. Oriented parallel to the lipid-water interface, the subsequent insertion of the peptide through the bilayer headgroups is significantly slower. The presence of negatively charged lipids considerably enhances peptide binding. Lateral side-chain motion creates an opening for the helix into the hydrophobic core of the membrane. The peptide aromatic residues form a pi-stacking cluster through W48/R52/W56 and F49/R53, protecting the peptide from the water phase. Interaction with the penetratin peptide has only limited effect on the overall membrane structure, as it affects mainly the conformation of the lipids which interact directly with the peptide. Charge matching locally increases the concentration of negatively charged lipids, lateral lipid diffusion locally decreases. Lipid disorder increases, through decreased order parameters of the lipids interacting with the penetratin side chains. Penetratin molecules at the membrane surface do not seem to aggregate.     

Molecular dynamics simulations of human prion protein: importance of correct treatment of electrostatic interactions.

Molecular dynamics simulations have been used to investigate the dynamical and structural behavior of a homology model of human prion protein HuPrP(90-230) generated from the NMR structure of the Syrian hamster prion protein ShPrP(90-231) and of ShPrP(<90-231) itself. These PrPs have a large number of charged residues on the protein surface. At the simulation pH 7, HuPrP(90-230) has a net charge of -1 eu from 15 positively and 14 negatively charged residues. Simulations for both PrPs, using the AMBER94 force field in a periodic box model with explicit water molecules, showed high sensitivity to the correct treatment of the electrostatic interactions. Highly unstable behavior of the structured region of the PrPs (127-230) was found using the truncation method, and stable trajectories could be achieved only by including all the long-range electrostatic interactions using the particle mesh Ewald (PME) method. The instability using the truncation method could not be reduced by adding sodium and chloride ions nor by replacing some of the sodium ions with calcium ions. The PME simulations showed, in accordance with NMR experiments with ShPrP and mouse PrP, a flexibly disordered N-terminal part, PrP(90-126), and a structured C-terminal part, PrP(127-230), which includes three alpha-helices and a short antiparallel beta-strand. The simulations showed some tendency for the highly conserved hydrophobic segment PrP(112-131) to adopt an alpha-helical conformation and for helix C to split at residues 212-213, a known disease-associated mutation site (Q212P). Three highly occupied salt bridges could be identified (E146/D144<-->R208, R164<-->D178, and R156<-->E196) which appear to be important for the stability of PrP by linking the stable main structured core (helices B and C) with the more flexible structured part (helix A and strands A and B). Two of these salt bridges involve disease-associated mutations (R208H and D178N). Decreased PrP stability shown by protein unfolding experiments on mutants of these residues and guanidinium chloride or temperature-induced unfolding studies indicating reduced stability at low pH are consistent with stabilization by salt bridges. The fact that electrostatic interactions, in general, and salt bridges, in particular, appear to play an important role in PrP stability has implications for PrP structure and stability at different pHs it may encounter physiologically during normal or abnormal recycling from the pH neutral membrane surface into endosomes or lysomes (acidic pHs) or in NMR experiments (5.2 for ShPrP and 4.5 for mouse PrP).     

Molecular dynamics simulations of peptides from the central domain of smooth muscle caldesmon

The central domain of smooth muscle caldesmon contains a highly charged region consisting of ten 13-residue repeats. Experimental evidence obtained from the intact protein and fragments thereof suggests that this entire region forms a single stretch of stable alpha-helix. We have carried out molecular dynamics simulations on peptides consisting of one, two and three repeats to examine the mechanism of alpha-helical stability of the central domain at the atomic level. All three peptides show high helical stability on the timescale of the MD simulations. Deviations from alpha-helical structure in all the simulations arise mainly from the formation of long stretches of pi-helix. Interconversion between alpha-helical and pi-helical conformations occurs through insertion of water molecules into alpha-helical hydrogen bonds and subsequent formation of reverse turns. The alpha-helical structure is stabilized by electrostatic interactions (salt bridges) between oppositely charged sidechains with i,i+4 spacings, while the pi-helix is stabilized by i,i+5 salt bridge interactions. Possible i,i+3 salt bridges are of minor importance. There is a strong preference for salt bridges with a Glu residue N-terminal to a basic sidechain as compared to the opposite orientation. In the double and triple repeat peptides, strong i,i+4 salt bridges exist between the last Glu residue of one repeat and the first Lys residue of the next. This demonstrates a relationship between the repetitive nature of the central domain sequence and its ability to form very long stretches of alpha-helical structure.