Stabilizing Salt Bridge Research Articles

Electrostatic Solvation Energy for Two Oppositely Charged Ions in a Solvated Protein System: Salt Bridges Can Stabilize Proteins

Born-type electrostatic continuum methods have been an indispensable ingredient in a variety of implicit-solvent methods that reduce computational effort by orders of magnitude compared to explicit-solvent MD simulations and thus enable treatment using larger systems and/or longer times. An analysis of the limitations and failures of the Born approaches serves as a guide for fundamental improvements without diminishing the importance of prior works. One of the major limitations of the Born theory is the lack of a liquidlike description of the response of solvent dipoles to the electrostatic field of the solute and the changes therein, a feature contained in the continuum Langevin-Debye (LD) model applied here to investigate how Coulombic interactions depend on the location of charges relative to the protein/water boundary. This physically more realistic LD model is applied to study the stability of salt bridges. When compared head to head using the same (independently measurable) physical parameters (radii, dielectric constants, etc.), the LD model is in good agreement with observations, whereas the Born model is grossly in error. Our calculations also suggest that a salt bridge on the protein's surface can be stabilizing when the charge separation is ≤4 Å.

How Stable is an Enzyme from a Thermophilic Organism? Denaturation Studies with the Esterase from Pyrococcus Furiosus - The Role of Charge-Charge Interactions

We have expressed the gene of the hyperthermophilic esterase (Pf2001Δ60) from Pyrococcus furiosus. This esterase showed to be active after boiling and has a half-life of 120 min at 75 °C. We decided to study the unfolding of this enzyme by fluorescence spectroscopy induced by urea, guanidinium hydrochloride (GndHCl) and high hydrostatic pressure (HHP). Ours results pointed out that urea is about three times more efficient than GndHCl in promoting structural perturbations in PF2001Δ60 ([D1/2]=7.1 and 2.3 M, respectively), suggesting that ion pairs are important stabilizing factors on its structure. There was almost no change in the tryptophan center of mass of PF2001Δ60 (342 nm) under HHP up to 3.1 kbar, even by staying 120 min under 2,500 bar. Interestingly, the combination of HHP with a subdenaturating concentration of urea (1 M) displaced the center of mass by ≈ 7 nm after 40 min at 2,500 bar. Since HHP enhances the electrostriction, this result reinforces the crucial contribution of salt bridges in esterase's stability. The binding of the bis(8-anilinonaphthalene-1-sulfonate) to PF2001Δ60 was increased by 2.5- and 3-fold after treatment with 2 M GndHCl or under HHP combined with urea, suggesting that these treatments convert the enzyme into a partially folded intermediate with exposed hydrophobic regions. Altogether, ours results may be an indication that the optimization of charge-charge interactions on the protein surface is a key factor for its stability. To our knowledge, this is the first time that HHP is used to access the ion pair contributions to the stability of a hyperthermophilic esterase.

Molecular Dynamics Investigation of Anesthetic Halothane Interactions with the Proton-Gated Ion Channel GLIC

Experimental studies found that general anesthetics inhibited the function of GLIC. To investigate the underlying inhibition mechanism, we studied the interaction of halothane with GLIC and determined the potential of mean force (PMF) for transporting a sodium ion across the GLIC channel in the presence and absence of halothane. Multiple halothane binding sites were identified through docking and multi-nanosecond molecular dynamics (MD) simulations; three important ones are: A) inter-subunit site near TM2-3 loop; B) inter-subunit site near D86; and C) intra-subunit site near W160. MD simulations demonstrated that halothane altered GLIC dynamics and electrostatic interactions between critical residues. Halothane near W160 reduced the stability of salt-bridges between D32 and R192, whose homologous electrostatic interaction was suggested to be important for channel gating in the ligand-gated ion channels. To assess potential impact of halothane binding to channel conductance, we calculated PMF using adaptive biased force method. Protonation state of E222 is a determinant for the PMF. PMF profile with five deprotonated E222 residues in the pentamer showed a significantly deeper energy well than that with three deprotonated E222, suggesting that five deprotonated E222 may trap an ion and that partial protonation of E222 is necessary for ion leaving the trap. While halothane near sites A and B introduced comparatively small effects on PMF, a profound PMF change near E222 was observed in the presence of halothane at W160. This change may likely affect the single channel conductance. Taken together, halothane may modulate GLIC function through altering salt bridges crucial for gating, coupled with a breaking of the fivefold symmetry of the pore-lining helices, leading to a change in the energy profile for ion passage through the channel. Supported by NIH (R01GM66358 and R01GM56257) and NCSA through PSC.

Molecular Modeling of Two Distinct Triangular Oligomers in Amyloid β-protein

Amyloid-beta (Abeta) peptides exhibit many distinct structural morphology at the early aggregate stage, some of which are biological relevant to the pathogenesis of Alzheimer's disease (AD). Atomic-resolution structures of the early Abeta aggregates and their conformational changes in amyloid aggregation remain elusive. Here, we perform all-atom molecular modeling and dynamics simulations to obtain two stable triangular-like Abeta structures with the lowest packing energy, one corresponding to the Tycko's model (Paravastu, A.; Leapman, R.; Yau, W.; Tycko, R. Proc. Nat. Acad. Soc. U.S.A. 2008, 105, 18349-18354) (referred to C-WT model) and the other corresponding to computational model (N-WT model). Both models have the same 3-fold symmetry but distinct beta-sheet organizations in which three Abeta hexamers pack together via either C-terminal beta-strand residues or N-terminal beta-strand residues forming distinct hydrophobic cross section. Structural and energetic comparisons of two 3-fold Abeta oligomers, coupled with structural changes upon the mutations occurring at the interacting interfaces, reveal that although hydrophobic interactions are still dominant forces, electrostatic interactions are more favorable in the N-WT model due to the formation of more and stable intersheet salt bridges, while solvation energy is more favorable in the C-WT model due to more exposed hydrophilic residues to solvent. Both models display many common features similar to other amyloid oligomers and therefore are likely to be biologically relevant.

Nonproteolytic Induction of Catalytic Activity into the Single-Chain Form of Urokinase-Type Plasminogen Activator by Dipeptides

Serine proteases are initially synthesized as single-chain proenzymes with activities that are many orders of magnitude lower than those of the mature enzyme. Proteolytic cleavage of an exposed loop liberates a new amino terminus that inserts into a hydrophobic pocket and forms a stabilizing salt bridge with a ubiquitously conserved aspartate residue, resulting in a conformational change organizing the mature oxyanion hole. In a decisive 1976 work, Huber and Bode [Bode, W., and Huber, R. (1976) FEBS Lett. 68, 231-236] demonstrated that peptides sequentially similar to the new amino terminus in combination with a catalytic site inhibitor could specifically induce a trypsin-like conformation in trypsinogen. We now demonstrate that an Ile-Ile or Ile-Val dipeptide can induce limited enzyme activity in the single-chain zymogen form of urokinase-type plasminogen activator (uPA) or its K158A variant, which cannot be activated proteolytically. Furthermore, the slow formation of a covalent serpin-protease complex between single-chain uPA and PAI-1 is significantly accelerated in the presence of specific dipeptide sequences. The technique of using a dipeptide mimic as a surrogate for the liberated amino terminus further provides a novel means by which to covalently label the immature active site of single-chain uPA with a fluorescent probe, permitting fluorescence approaches for direct observations of conformational changes within the protease domain during zymogen activation. These data demonstrate the structural plasticity of the protease domain, reinforce the notion of "molecular sexuality", and provide a novel way of studying conformational changes of zymogens during proteolytic activation.

Role of Conserved Salt Bridges in Homeodomain Stability and DNA Binding

The sequence information available for homeodomains reveals that salt bridges connecting pairs 19/30, 31/42, and 17/52 are frequent, whereas aliphatic residues at these sites are rare and mainly restricted to proteins from homeotherms. We have analyzed the influence of salt and hydrophobic bridges at these sites on the stability and DNA binding properties of human Hesx-1 homeodomain. Regarding the protein stability, our analysis shows that hydrophobic side chains are clearly preferred at positions 19/30 and 31/42. This stabilizing influence results from the more favorable packing of the aliphatic side chains with the protein core, as illustrated by the three-dimensional solution structure of a thermostable variant, herein reported. In contrast only polar side chains seem to be tolerated at positions 17/52. Interestingly, despite the significant influence of pairs 19/30 and 31/42 on the stability of the homeodomain, their effect on DNA binding ranges from modest to negligible. The observed lack of correlation between binding strength and conformational stability in the analyzed variants suggests that salt/hydrophobic bridges at these specific positions might have been employed by evolution to independently modulate both properties.

Structures of Protonated Dipeptides: The Role of Arginine in Stabilizing Salt Bridges

Structures of protonated dipeptides containing N-terminal Gly, Val, Pro, Lys, His, or Arg and C-terminal Arg are investigated with infrared multiple photon dissociation (IRMPD) spectroscopy between 900 and 1850 cm(-1) and theory. The IRMPD spectra clearly indicate that, when Gly, Val, Pro, Lys, or His are N-terminal to Arg, these protonated dipeptides adopt gas-phase structures with a single formal charge site (SCS), whereas ArgArg x H(+) has a salt-bridge (SB) structure in which the C-terminus is deprotonated and two basic sites are protonated. There are only subtle differences in the IRMPD spectra for dipeptides containing Gly, Val, Pro, and Lys. A sharp, intense peak at 1080 cm(-1) is observed for HisArg x H(+) that is attributed to the neutral histidine side chain, an assignment that is confirmed by comparison to the IRMPD spectrum of (HisArg x H(2))(2+). Lowest-energy B3LYP/6-31+G(d,p) structures and energies for the SCS and SB forms of these protonated dipeptides indicate that stability of the SB form relative to the SCS form generally increases with increasing gas-phase basicity of the N-terminal amino acid, but only ArgArg x H(+) is calculated to have a SB ground state at 298 K, in agreement with the results from IRMPD spectroscopy. This is the first direct experimental evidence for a salt-bridge structure in a gaseous protonated peptide, and ArgArg x H(+) is the smallest protonated peptide for which a SB structure has been reported. These results suggest that SB structures should be common for protonated peptides containing at least two arginine residues and may also occur for large protonated peptides or proteins with at least one arginine residue and other basic residues, such as lysine or histidine.

Anionic Lipid and Cholesterol Interactions with α4β2 nAChR: Insights from MD Simulations

Anionic lipids and cholesterols (CHOL) are critical to the function of nicotinic acetylcholine receptors (nAChR). We investigated their interactions with an open- and closed-channel alpha4beta2 nAChR by over 10 ns molecular dynamics simulations in a ternary lipid mixture of 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl phosphatidic acid (POPA), and CHOL with a ratio of 3:1:1 (Haddadian et al., J. Phys. Chem. B 2008, 112, 13981). On average there were 65 and 74 interfacial lipids around the closed- and open-channel alpha4beta2 nAChR, respectively, in the equilibrated simulation systems. In the open-channel system, 42% of the interfacial POPA had acyl chains partially inserted into intra- or intersubunit cavities, as compared to only 7% in the closed-channel alpha4beta2. No CHOL was found in cavities within single subunits, though some CHOL infiltrated into the gaps between subunits. Because of its smaller headgroup, POPA could access some nonannular sites where POPC could not easily reach due to steric exclusion. Furthermore, POPA acted not only as an acceptor for hydrogen bonding (H bonding) as POPC did, but also as a donor through its hydroxyl group for H bonding with the backbone of the protein. The charged headgroup of POPA allowed the lipid to form stable salt bridges with conserved Arg and Lys residues at the interfaces of the transmembrane (TM) and extracellular (EC) or intracellular (IC) domains of the alpha4beta2. A higher number of salt bridges and hydrogen bonds (H bonds) between POPA and the alpha4beta2 nAChR were found in the open system than in the closed system, suggesting a potential role of POPA in the equilibrium between different channel states. Most interfacial POPA molecules showed lower order parameters than the bulk POPA due to the mixed effect of gauche defects, hydrophobic mismatch, and the lipid orientations near the magic angle. These unique properties enable the interfacial POPA to achieve what POPC cannot with regard to specific interactions with the protein, thereby making POPA essential for the function of nAChR.

Ionic Surfactant Aggregates in Saline Solutions: Sodium Dodecyl Sulfate (SDS) in the Presence of Excess Sodium Chloride (NaCl) or Calcium Chloride (CaCl2)

The properties of sodium dodecyl sulfate (SDS) aggregates in saline solutions of excess sodium chloride (NaCl) or calcium chloride (CaCl(2)) ions were studied through extensive molecular dynamics simulations with explicit solvent. We find that the ionic strength of the solution affects not only the aggregate size of the resulting anionic micelles but also their structure. Specifically, the presence of CaCl(2) induces more compact and densely packed micelles with a significant reduction in gauche defects in the SDS hydrocarbon chains in comparison with NaCl. Furthermore, we observe significantly more stable salt bridges between the charged SDS head groups mediated by Ca(2+) than Na(+). The presence of these salt bridges helps stabilize the more densely packed micelles.

Conformational Changes at the Tetramerization Site of Erythroid α-Spectrin upon Binding β-Spectrin: A Spin Label EPR Study

We used cysteine scanning, isothermal titration calorimetry (ITC) and spin label EPR methods to study the two regions that flank the partial domain Helix C' of the N-terminal end of alpha-spectrin (residues 14-20 and residues 44-54) in the absence and presence of a model protein of the beta-spectrin C-terminal end. In the absence of beta-spectrin, residues 14-20 and 46-52 were known to be unstructured. The EPR spectral values of the inverse line width (Delta H (-1)) and of the width between the low field peak and the central peak ( aZ) of residues in part of the first unstructured region (residues 17-20) and of most residues in the second unstructured junction region (residues 46-52) changed dramatically upon association with beta-spectrin, suggesting that the two regions undergo a conformational change, becoming more rigid and likely becoming helical. ITC results showed that three of the seven residues in the junction region (residues 46-52) were very important in its association with beta-spectrin, in the following order: L49 > G46 > K48. In general, our results suggest that any mutations that affect the propensity of helical formation in the region spanning residues 17-52 in alpha-spectrin, or that affect hydrophobic clustering and/or salt-bridge stabilization of the bundled helices, would affect spectrin tetramer formation, and may lead to blood disorders.

Insights into thermal stability of thermophilic nitrile hydratases by molecular dynamics simulation

Thermal stability is of great importance for industrial enzymes. Here we explored the thermal-stable mechanism of thermophilic nitrile hydratases (NHases) utilizing a molecular dynamic simulation. At a nanosecond timescale, profiles of root mean square fluctuation (RMSF) of two thermophilic NHases, 1UGQ and 1V29, under enhancing thermal stress were carried out at 300 K, 320 K, 350 K and 370 K, respectively. Results showed that the region A1 (211–231 aa) and A2 (305–316 aa) in 1UGQ, region B1 (186–192 aa) in 1V29, and most of terminal ends in both enzymes are hyper-sensitive. Salt-bridge analyses revealed that in one hand, salt-bridges contributed to maintaining the rigid structure and stable performance of the thermophilic 1UGQ and 1V29; in the other hand, salt-bridges involved in thermal sensitive regions are relatively weak and prone to be broken at elevated temperature, thereby cannot hold the stable conformation of the spatial neighborhood. In 1V29, region A1 was stabilized by a well-organized hook–hook like cluster with multiple salt-bridge interactions, region A2 was stabilized by two strong salt-bridge interactions of GLU52-ARG332 and GLU334-ARG332. In 1UGQ, the absence of a charged residue decreased its thermal sensitivity of region B1, and the formation of a small β-sheet containing a stable salt-bridge in C-β-terminal significantly enhanced its thermal stability. By radius of gyration calculation containing or eliminating the thermal sensitive regions, we quantified the contribution of thermal sensitive regions for thermal sensitivity of 1UGQ and 1V29. Consequently, we presented strategies to improve thermal stability of the industrialized mesophilic NHase by introducing stable salt-bridge interactions into its thermal sensitive regions.

Molecular dynamics guided study of salt bridge length dependence in both fluorinated and non-fluorinated parallel dimeric coiled-coils.

The alpha-helical coiled-coil is one of the most common oligomerization motifs found in both native and engineered proteins. To better understand the stability and dynamics of the coiled-coil motifs, including those modified by fluorination, several fluorinated and nonfluorinated parallel dimeric coiled-coil protein structures were designed and modeled. We also attempt to investigate how changing the length and geometry of the important stabilizing salt bridges influences the coiled-coil protein structure. Molecular dynamics (MD) and free energy simulations with AMBER used a particle mesh Ewald treatment of the electrostatics in explicit TIP3P solvent with balanced force field treatments. Preliminary studies with legacy force fields (ff94, ff96, and ff99) show a profound instability of the coiled-coil structures in short MD simulation. Significantly, better behavior is evident with the more balanced ff99SB and ff03 protein force fields. Overall, the results suggest that the coiled-coil structures can readily accommodate the larger acidic arginine or S-2,7-diaminoheptanedoic acid mutants in the salt bridge, whereas substitution of the smaller L-ornithine residue leads to rapid disruption of the coiled-coil structure on the MD simulation time scale. This structural distortion of the secondary structure allows both the formation of large hydration pockets proximal to the charged groups and within the hydrophobic core. Moreover, the increased structural fluctuations and movement lead to a decrease in the water occupancy lifetimes in the hydration pockets. In contrast, analysis of the hydration in the stable dimeric coiled-coils shows high occupancy water sites along the backbone residues with no water occupancy in the hydrophobic core, although transitory water interactions with the salt bridge residues are evident. The simulations of the fluorinated coiled-coils suggest that in some cases fluorination electrostatically stabilizes the intermolecular coiled-coil salt bridges. Structural analyses also reveal different side chain rotamer preferences for leucine when compared with 5,5,5,5',5',5'-hexafluoroleucine mutants. These observed differences in the side chain rotamer populations suggest differential changes in the side chain conformational entropy upon coiled-coil formation when the protein is fluorinated. The free energy of hydration of the isolated 5,5,5,5',5',5'-hexafluoroleucine amino acid is calculated to be 1.1 kcal/mol less stable than leucine; this hydrophobic penalty in the monomer may provide a driving force for coiled-coil dimer formation. Estimation of the ellipticity at 222 nm from a series of snapshots from the MD simulations with DicroCalc shows distinct increases in the ellipticity when the coiled-coil is fluorinated, which suggests that the helicity in the folded coiled-coils is greater when fluorinated.

The Atomistic Mechanism of Conformational Transition in Adenylate Kinase: A TEE-REX Molecular Dynamics Study

We report on an atomistic molecular dynamics simulation of the complete conformational transition of Escherichia coli adenylate kinase (ADK) using the recently developed TEE-REX algorithm. Two phases characterize the transition pathway of ADK, which folds into the domains CORE and LID and the AMP binding domain AMPbd. Starting from the closed conformation, half-opening of the AMPbd precedes a partially correlated opening of the LID and AMPbd, defining the second phase. A highly stable salt bridge D118-K136 at the LID-CORE interface, contributing substantially to the total nonbonded LID-CORE interactions, was identified as a major factor that stabilizes the open conformation. Alternative transition pathways, such as AMPbd opening following LID opening, seem unlikely, as full transition events were not observed along this pathway. The simulation data indicate a high enthalpic penalty, possibly obstructing transitions along this route.

Investigation of double‐stranded DNA/drug interaction by ESI/FT ICR: orientation of dissociations relates to stabilizing salt bridges

Noncovalent complexes of DNA and Hoechst 33258 were investigated by ESI-FT/ICR MS in various activation modes (collision-induced dissociation (CID), sustained off-resonance irradiation collision-induced dissociation (SORI-CID), infrared multiphoton dissociation (IRMPD) and electron detachment dissociation (EDD)). The binding selectivity of Hoechst 33258 was confirmed by the comparative study of its noncovalent association with different DNA sequences. The CID spectra of [ds + HO - 5H](5-) obtained with a linear hexapole ion trap resulted in unzipping of the strands. This outcome is a clue to the drug-binding mode, shading light on the localization of the binding sites of Hoechst 33258 to the DNA sequence. The IRMPD and SORI-CID experiments mainly gave DNA backbone cleavages and internal fragment ions. From this result, information on the localization of the binding sites of Hoechst 33258 in the DNA sequence was obtained. No sodium cationization was observed on the DNA sequence ions although they were present on fragmentation of the duplex, indicating that the backbone cleavages were generated from the single strand associated with the Hoechst 33258 where the number of alkali cation is restricted. Under electron detachment (ED) conditions, multiple EDs were achieved for the [ds + HO - 5H](5-) ion without any significant dissociation. The presence of drug appears to enhance the stability of the multiply charged system. It was proposed that the studied noncovalent complex involved the formation of zwitterions and consequently strong salt-bridge interactions between DNA and drug.

Effects of replacing active site residues in a cold‐active alkaline phosphatase with those found in its mesophilic counterpart fromEscherichia coli

Alkaline phosphatase (AP) from a North Atlantic marine Vibrio bacterium was previously characterized as being kinetically cold-adapted. It is still unknown whether its characteristics originate locally in the active site or are linked to more general structural factors. There are three metal-binding sites in the active site of APs, and all three metal ions participate in catalysis. The amino acid residues that bind the two zinc ions most commonly present are conserved in all known APs. In contrast, two of the residues that bind the third metal ion (numbered 153 and 328 in Escherichia coli AP) are different in various APs. This may explain their different catalytic efficiencies, as the Mg2+ most often present there is important for both structural stability and the reaction mechanism. We have mutated these key residues to the corresponding residues in E. coli AP to obtain the double mutant Asp116/Lys274, and both single mutants. All these mutants displayed reduced substrate affinity and lower overall reaction rates. The Lys274 and Asp116/Lys274 mutants also displayed an increase in global heat stability, which may be due to the formation of a stabilizing salt bridge. Overall, the results show that a single amino acid substitution in the active site is sufficient to alter the structural stability of the cold-active Vibrio AP both locally and globally, and this influences kinetic properties.

Reversal of Agonist Selectivity by Mutations of Conserved Amino Acids in the Binding Site of Nicotinic Acetylcholine Receptors

Homomeric alpha7 and heteromeric alpha4beta2 nicotinic acetylcholine receptors (nAChR) can be distinguished by their pharmacological properties, including agonist specificity. We introduced point mutations of conserved amino acids within the C loop, a region of the receptor critical for agonist binding, and we examined the expression of the mutant receptors in Xenopus oocytes. Mutation of either a conserved C loop tyrosine (188) to phenylalanine or a nearby conserved aspartate (197) to alanine resulted in alpha7 receptors for which the alpha7-selective agonist 3-(4-hydroxy, 2-methoxybenzylidene) anabaseine (4OH-GTS-21) had roughly the same potency as for wild-type receptors, whereas the physiologic agonist acetylcholine (ACh) showed drastically reduced potency for these mutant receptors. Corresponding mutations in alpha4 receptors co-expressed with beta2 resulted in alpha4beta2 receptors for which ACh potency was relatively unchanged, although the efficacy of the alpha7-selective agonist 4OH-GTS-21 was increased greatly relative to that of ACh. We also investigated the significance of a conserved lysine (145 in alpha7), proposed to form a stable salt bridge with Asp-197 in the resting state of the receptor. Mutations of this residue in both alpha7 and alpha4 resulted in receptors that were largely unresponsive to both ACh and 4OH-GTS-21. Our results suggest that initiation of gating depends both on specific interactions between residues in the C loop domain and, depending on receptor subtype, the physiochemical properties of the agonist, so that in the altered environment of the alpha4Y190F-binding site, large hydrophobic benzylidene anabaseines may close the C loop and initiate channel gating more effectively than the polar agonist ACh.

Thrombin-activable Fibrinolysis Inhibitor (TAFI) Zymogen Is an Active Carboxypeptidase

Thrombin-activable fibrinolysis inhibitor (TAFI) is a carboxypeptidase found in human plasma, presumably as an inactive zymogen. The current dogma is that proteolytic activation by thrombin/thrombomodulin generates the active enzyme (TAFIa), which down-regulates fibrinolysis by removing C-terminal lysine residues from partially degraded fibrin. In this study, we have shown that the zymogen exhibits continuous and stable carboxypeptidase activity against large peptide substrates, and we suggest that the activity down-regulates fibrinolysis in vivo.

Relationship between Salt-Bridge Identity and 14-Helix Stability of β3-Peptides in Aqueous Buffer

We report a systematic analysis of the relationship between salt bridge composition and 14-helix structure within a family of model beta-peptides in aqueous buffer. We find an inverse relationship between side-chain length and the extent of 14-helix structure as judged by CD. Introduction of a stabilizing salt bridge pair within a previously reported beta-peptide ligand for hDM2 led to changes in structure that were detectable by NMR.

On the Influence of Charged Side Chains on the Folding–Unfolding Equilibrium of β‐Peptides: A Molecular Dynamics Simulation Study

The influence of charged side chains on the folding-unfolding equilibrium of beta-peptides was investigated by means of molecular dynamics simulations. Four different peptides containing only negatively charged side chains, positively charged side chains, both types of charged side chains (with the ability to form stabilizing salt bridges) or no charged side chains were studied under various conditions (different simulation temperatures, starting structures and solvent environment). The NMR solution structure in methanol of one of the peptides (A) has already been published; the synthesis and NMR analysis of another peptide (B) is described here. The other peptides (C and D) studied herein have hitherto not been synthesized. All four peptides A-D are expected to adopt a left-handed 3(14)-helix in solution as well as in the simulations. The resulting ensembles of structures were analyzed in terms of conformational space sampled by the peptides, folding behavior, structural properties such as hydrogen bonding, side chain-side chain and side chain-backbone interactions and in terms of the level of agreement with the NMR data available for two of the peptides. It was found that the presence of charged side chains significantly slows down the folding process in methanol solution due to the stabilization of intermediate conformers with side chain-backbone interactions. In water, where the solvent competes with the solute-solute polar interactions, the folding process to the 3(14)-helix is faster in the simulations.

Investigation of Salt Bridge Stability in a Generalized Born Solvent Model.

Potentials of mean force (PMFs) of salt bridge formation between oppositely charged amino acid side chains were calculated both in explicit solvent and in a Generalized Born (GB) continuum solvent model to quantify the potential overstabilization of side chain ion pairs in GB relative to explicit solvation. These show that salt bridges are too stable by as much as 3-4 kcal/mol in the GB solvent models that we tested, consistent with previously reported observations of significantly different structural ensembles in GB models and explicit solvent for proteins containing ionizable groups. We thus investigated a simple empirical correction, wherein the intrinsic GB radii of hydrogen atoms bound to charged nitrogen atoms are reduced, effectively increasing the desolvation penalty of the positively charged groups. The thermodynamics of salt bridge formation were considerably improved, as exemplified by the close match of the corrected GB PMF to the reference explicit solvent PMF, and more significantly by our ability to closely reproduce the experimental temperature melting profile of the TC5b Trp-cage miniprotein, which is otherwise highly distorted by prevalent non-native salt bridges when using standard GB parameters.