Salt-Specific Stability and Denaturation of a Short Salt-Bridge-Forming α-Helix

Understanding of formation, gastrointestinal breakdown, and application of whey protein emulsion gels: Insights from intermolecular interactions

Whey protein emulsion gel is an ideal model food for revealing how the multilength scale food structures affect food digestion, as their structure and mechanical properties can be precisely manipulated by controlling the type and intensity of intermolecular interactions between protein molecules. However, there are still significant understanding gaps among intermolecular interactions, protein aggregation and gelation, emulsion gel formation, gel breakdown in the gastrointestinal tract (GIT), and the practical use of whey protein emulsion gels, which limits their GIT-targeted applications. In this regard, the relationship between the structure and digestion behavior of heat-set whey protein emulsion gels is reviewed and discussed mainly from the following aspects: (1) structural characteristics of whey protein molecules; (2) how different types of intermolecular interactions influence heat-induced aggregation and gelation of whey protein in the aqueous solutions and the oil-in-water emulsions, and the mechanical properties of the final gels; (3) functions of the mouth, the stomach, and the small intestine in processing of solid foods, and how different types of intermolecular interactions influence the breakdown properties of heat-set whey protein emulsion gels in GIT (i.e., their respective role in controlling gel digestion). Finally, the implications of knowledge derived from the formation and gastrointestinal breakdown of heat-set whey protein emulsion gels for developing controlled delivery vehicles, human satiety enhancers, and sensory modifiers are highlighted.

The critical role of aromatic residues in the binding of the SARS-CoV-2 fusion peptide to phospholipid bilayer membranes

Based on the SARS-CoV-2 fusion peptide (FP) structure determined from the NMR experiment, we created six FP models under different environmental conditions to explore the effects of salt and cholesterol on FP-membrane binding. The all-atom molecular dynamics (MD) simulation results indicated that ionic environments notably impact the FP structure as well as the stability of the helical elements within the peptide. Our findings highlighted the unpredictable influence of ions on the secondary structures and dynamics of the FP, emphasizing the complexity and sensitivity of the peptide's conformations to ionic conditions. When exploring the peptide's interaction with a cholesterol-free phospholipid bilayer membrane, we found that the helical elements of the FP remain stable irrespective of the salt type (Na+ or Ca2+). This result emphasizes the crucial role of phospholipid bilayer membranes in supporting the secondary structures of the FP. The MD simulation results showed that Ca2+ ions facilitated deeper membrane penetration than Na+ ions, highlighting the critical role of calcium ions in the FP-membrane binding. Our study indicates the essential role of the aromatic residues (such as Phe833 and Tyr837) in the FP-membrane binding process. Finally, we investigated the FP-membrane binding patterns in the presence of cholesterol. The MD simulation results demonstrated that the coupling of Ca2+ ions and cholesterol would also benefit the FP-membrane binding. Furthermore, our findings reveal that while the type of ion and cholesterol content exert varied and unpredictable influences on FP-membrane binding patterns, aromatic residues like tyrosine (Tyr) and phenylalanine (Phe) play an essential role in FP-membrane binding. In particular, deep mutational scanning (DMS) experiments have confirmed that mutating phenylalanine in the FP significantly decreases viral mutational fitness, emphasizing the pivotal role of phenylalanine residues in membrane fusion. This knowledge can aid in developing more effective therapeutic strategies targeting the viral fusion peptide and its key amino acids, ultimately contributing to developing treatments and vaccines against the virus.

Structural and dynamic properties of urinary human serum albumin fragments: a molecular dynamics study

A microalbuminuria level acts as a good index to screen and monitor diabetes and renal failure. However, the urinary albumin loss after sample preservation and storage is the major bottleneck to obtain the accurate microalbuminuria test. Such loss is due to the rapid albumin fragmentation by urinary proteases. Some fragments were suggested to be bioactive biomarkers of diabetes and renal disease, but no structural and dynamical properties of albumin fragments are available. Thus, in this work, the structural and dynamical properties of reported albumin fragments are revealed using molecular dynamics simulations. The properties of nine fragments (F1-F9) discovered recently were studied at the real pH conditions of urine samples (pH 4.5, 7 and 8). The complete loss of secondary structure is found in short fragments (F1-F6), while large-sized polypeptides (F7-F9) can somehow maintain their folds. Especially, F8 (subdomain IIIB) is the most stable fragment. The difference in histidine protonation states has no impact on the structural stability of albumin fragments. The ability of F8 (subdomain IIIB) to maintain its stability and folds suggests it as an alternative albumin biomarker in urine. An insight obtained here will become the fundamental importance for understanding clinical assays for albumin detection, sample stability and peptidomics analysis of urine.Communicated by Ramaswamy H. Sarma.

A computational study of the R120G mutation in human αB-crystallin: implications for structural stability and functionality

The eye is a vital organ in the visual system, which is composed of transparent vascular tissue. αB-crystallin, a significant protein found in the lens, plays a crucial role in our understanding of lens diseases. Mutations in the αB-crystallin protein can cause lens diseases, such as cataracts and myopathy. However, the molecular mechanism underlying the R120G mutation is not fully understood. In this study, we utilized molecular dynamics simulations to illustrate, in atomic detail, how the R120G mutation leads to the aggregation of αB-crystallin and scattering of light in the lens. Our findings show that the R120G mutation alters the dynamic and structural properties of the αB-crystallin protein. Specifically, this mutation causes the angle of the hairpin at the C-terminal to increase from 80° to 150°, while reducing the distance between the hydrophobic patches around residues 10 and 44-55 from 1.5 nm to 1 nm. In addition, our results showed that the mutation could disrupt the IPI motif - β4/β8 interaction. The disruption of this interaction could affect the αB-crystallin oligomerization and the chaperone activity of αB-crystallin protein. The exposed hydrophobic area at the IPI motif - β4/β8 could become the primary site for interprotein interactions, which are responsible for large-scale aggregation. We have demonstrated that, in wild-type αB-crystallin protein, salt bridges R120 and D109, R107 and D80 are formed. However, in the case of the R120G mutation, the salt bridges R120 and R109 are disrupted, and a new salt bridge with a different pattern is formed. In our study, it has been found that all of the changes associated with the R120G mutation are located at the interface of chains A and B, which could impact the multimerization of the αB-crystallin. Previous research on the K92-E99 residue has shown that a salt bridge in the dimer I can reduce the chaperone activity of the protein. Furthermore, the salt bridges R120 and D109, as well as R107 and D80 in dimer II, induce changes in the hydrophobic envelope of β-sheets in the α-crystallin domain (ACD). These changes could have an impact on the multimerization of the αB-crystallin, leading to disruption of the oligomer structure and aggregation. Moreover, the changes in the αB-crystallin resulting from the R120G mutation can lead to faulty interactions with other proteins, which can cause the aggregation of αB-crystallin with other proteins, such as desmin. These findings may provide new insights into the development of treatments for lens diseases.Communicated by Ramaswamy H. Sarma.

Effects of Charge Sequence Pattern and Lysine-to-Arginine Substitution on the Structural Stability of Bioinspired Polyampholytes

A comprehensive study focusing on the combined influence of the charge sequence pattern and the type of positively charged amino acids on the formation of secondary structures in sequence-specific polyampholytes is presented. The sequences of interest consisting exclusively of ionizable amino acids (lysine, K; arginine, R; and glutamic acid, E) are (EKEK)5, (EKKE)5, (ERER)5, (ERRE)5, and (EKER)5. The stability of the secondary structure was examined at three pH values in the presence of urea and NaCl. The results presented here underscore the combined prominent effects of the charge sequence pattern and the type of positively charged monomers on secondary structure formation. Additionally, (ERRE)5 readily aggregated across a wide range of pH. In contrast, sequences with the same charge pattern, (EKKE)5, as well as the sequences with the equivalent amino acid content, (ERER)5, exhibited no aggregate formation under equivalent pH and concentration conditions.

Disassembly and reassembly of the non-conventional thermophilic C-phycocyanin

C-phycocyanin (CPC), which contains open-chain tetrapyrroles, is a major light-harvesting red-fluorescent protein with an important role in aquatic photosynthesis. Recently, we reported a non-conventional CPC from Thermoleptolyngbya sp. O-77 (CPCO77) that contains two different structures, i.e., a hexameric structure and a non-conventional octameric structure. However, the assembly and disassembly mechanisms of the non-conventional octameric form of CPC remain unclear. To understand this assembly mechanism, we performed an in vitro experiment to study the disassembly and reassembly behaviors of CPC using isolated CPC subunits. The dissociation of the CPCO77 subunit was performed using a Phenyl-Sepharose column in 20 mM potassium phosphate buffer (pH 6.0) containing 7.0 M urea. For the first time, crystals of isolated CPC subunits were obtained and analyzed after separation. After the removal of urea from the purified α and β subunits, we performed an in vitro reassembly experiment for CPC and analyzed the reconstructed CPC using spectrophotometric and X-ray crystal structure analyses. The crystal structure of the reassembled CPC was nearly identical to that of the original CPCO77. The findings of this study indicate that the octameric CPCO77 is a naturally occurring form in the thermophilic cyanobacterium Thermoleptolyngbya sp. O-77.

Effect of Bacterial Amyloid Protein Phenol-Soluble Modulin Alpha 3 on the Aggregation of Amyloid Beta Protein Associated with Alzheimer's Disease

Since the proposal of the brainstem axis theory, increasing research attention has been paid to the interactions between bacterial amyloids produced by intestinal flora and the amyloid β-protein (Aβ) related to Alzheimer's disease (AD), and it has been considered as the possible cause of AD. Therefore, phenol-soluble modulin (PSM) α3, the most virulent protein secreted by Staphylococcus aureus, has attracted much attention. In this work, the effect of PSMα3 with a unique cross-α fibril architecture on the aggregation of pathogenic Aβ40 of AD was studied by extensive biophysical characterizations. The results proposed that the PSMα3 monomer inhibited the aggregation of Aβ40 in a concentration-dependent manner and changed the aggregation pathway to form granular aggregates. However, PSMα3 oligomers promoted the generation of the β-sheet structure, thus shortening the lag phase of Aβ40 aggregation. Moreover, the higher the cross-α content of PSMα3, the stronger the effect of the promotion, indicating that the cross-α structure of PSMα3 plays a crucial role in the aggregation of Aβ40. Further molecular dynamics (MD) simulations have shown that the Met1-Gly20 region in the PSMα3 monomer can be combined with the Asp1-Ala2 and His13-Val36 regions in the Aβ40 monomer by hydrophobic and electrostatic interactions, which prevents the conformational conversion of Aβ40 from the α-helix to β-sheet structure. By contrast, PSMα3 oligomers mainly combined with the central hydrophobic core (CHC) and the C-terminal region of the Aβ40 monomer by weak H-bonding and hydrophobic interactions, which could not inhibit the transition to the β-sheet structure in the aggregation pathway. Thus, the research has unraveled molecular interactions between Aβ40 and PSMα3 of different structures and provided a deeper understanding of the complex interactions between bacterial amyloids and AD-related pathogenic Aβ.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

Solubility and Thermal Stability of Thermotoga maritima MreB

The basis of MreB research is the study of the MreB protein from the Thermotoga maritima species, since it was the first one whose crystal structure was described. Since MreB proteins from different bacterial species show different polymerisation properties in terms of nucleotide and salt dependence, we conducted our research in this direction. For this, we performed measurements based on tryptophan emission, which were supplemented with temperature-dependent and chemical denaturation experiments. The role of nucleotide binding was studied through the fluorescent analogue TNP-ATP. These experiments show that Thermotoga maritima MreB is stabilised in the presence of low salt buffer and ATP. In the course of our work, we developed a new expression and purification procedure that allows us to obtain a large amount of pure, functional protein.

Sequence-Controlled Secondary Structures and Stimuli Responsiveness of Bioinspired Polyampholytes

A comprehensive study focusing on the influence of the sequence charge pattern on the secondary structure preferences of annealed polyampholytes and their responsiveness to external stimuli is presented. Two sequences are designed composed entirely of ionizable amino acids (charge fraction, f = 1) and an equal number of positive and negative charges (f+ = f- = 0.5) with distinct charge patterns consisting of lysine and glutamic acid monomers. The study reveals that the sequence charge pattern has a significant influence on the secondary structure preferences of polyampholytes at physiological pH. Furthermore, it shows that external stimuli such as pH, ionic strength, and solvent dielectric constant can be used to modulate the secondary structure of the two studied sequences. The observed secondary structure transformations for the two sequences are also substantially different from those determined for uniformly charged homo-polypeptides under matching conditions.

Ionic strength-sensitive and pH-insensitive interactions between C-reactive protein (CRP) and an anti-CRP antibody

C-reactive protein (CRP) is an important biomarker of infection and inflammation, as CRP is one of the most prominent acute-phase proteins. CRP is usually detected using anti-CRP antibodies (Abs), where the intermolecular interactions between CRP and the anti-CRP Ab are largely affected by the pH and ionic strength of environmental solutions. Therefore, it is important to understand the environmental effects of CRP–anti-CRP Ab interactions when designing highly sensitive biosensors. Here, we investigated the efficiency of fluorescently labeled CRP–anti-CRP monoclonal antibody (mAb) interactions at different pHs and ionic strengths. Our results indicate that the affinity was insensitive to pH changes in the range of 5.9 to 8.1, while it was significantly sensitive to ionic strength changes. The binding affinity decreased by 55% at an ionic strength of 1.6 mM, when compared to that under a physiological condition (~150 mM). Based on the isoelectric focusing results, both the labeled CRP and anti-CRP mAb were negatively charged in the studied pH range, which rendered the system insensitive to pH changes, but sensitive to ionic strength changes. The decreased ionic strength led to a significant enhancement of the repulsive force between CRP and the anti-CRP mAb. Although the versality of the findings is not fully studied yet, the results provide insights into designing highly sensitive CRP sensors, especially field-effect transistor-based sensors.

α-Helices propagating from stable nucleators exhibit unconventional thermal folding

Although the effect of thermal perturbations on protein structure has long been modeled in helical peptides, several details, such as the relation between the thermal stabilities of the propagating and nucleating segments of helices, remain elusive. We had earlier reported on the helix-nucleating propensities of covalent H-bond surrogate-constrained α-turns. Here, we analyze the thermal stabilities of helices that propagate along peptides appended to these α-helix nucleators using their NMR and far-UV CD spectra. Unconventional thermal folding of these helix models reveals that the helical fold in propagating backbones resists thermal perturbations as long as their nucleating template is intact. The threshold temperature of such resistance is also influenced by the extent of similarity between the natures of helical folds in the nucleating and propagating segments. Correlations between helicities and rigidities of helix-nucleating and helix-propagating segments reveal subtle interdependence, which explains cooperativity and residual helix formation during protein folding.

Interaction strength of osmolytes with the anion of a salt-bridge determines its stability

In order to understand the role of osmolytes in regulating physicochemical behavior of proteins, we investigated the influence of protein destabilizing (urea and guanidinium chloride) and stabilizing osmolytes (TMAO, glycerol, and betaine) on a model salt-bridge (SB) formed between structural analogues of arginine and glutamate/aspartate sidechains in a solvent continuum using first-principles quantum chemical calculations based on DFT and MP2 methods. The binding strength of the osmolyte with the SB is found to be in the order of betaine > TMAO > Gdm+ > glycerol > urea. The osmolytes (TMAO and betaine) that preferentially bind to the SB cation have a marginal influence on SB stability. Also, pure π-π stacking interaction between Gdm+ and the SB cation plays an insignificant role in destabilizing the SB. In fact, the interaction strength of osmolytes with the SB anion mainly determines the stability of SB. For instance, a competition between Gdm+ and the SB cation to bind with the SB anion is responsible for instability and subsequent dissociation of the SB. The competition provided by other osmolytes is too weak to break the SB. Exploiting this information, we designed three structural derivatives of Gdm+, all having a stronger interaction with SB anion, and thereby show a stronger SB dissociation potential. Furthermore, we find an excellent linear anti-correlation between SB interaction energy and the energy of interaction between osmolyte and the SB anion, which suggests that by knowing only the strength of osmolyteacetate interaction, one can predict the influence of osmolytes on the salt-bridge instability. This information is useful in fine-tuning the SB dissociation power of Gdm+, which has a practical significance in obtaining the mechanistic insight into the influence of GdmCl on protein stability. Our results also provide a basis for understanding the chemistry of other ion-pairs formed between a cationic hydrogen donor and an anionic acceptor.

Molecular Insight into Cu2+-Induced Conformational Transitions of Amyloid β-Protein from Fast Kinetic Analysis and Molecular Dynamics Simulations

Cu²⁺-mediated amyloid β-protein (Aβ) aggregation is implicated in the pathogenesis of Alzheimer's disease, so it is of significance to understand Cu²⁺-mediated conformational transitions of Aβ. Herein, four Aβ mutants were created by using the environment-sensitive cyanophenylalanine to respectively substitute F4, Y10, F19, and F20 residues of Aβ₄₀. By using stopped-flow fluorescence spectroscopy and molecular dynamics (MD) simulations, the early stage conformational transitions of the mutants mediated by Cu²⁺ binding were investigated. The fast kinetics unveils that Cu²⁺ has more significant influence on the conformational changes of N-terminal (F4 and Y10) than on the central hydrophobic core (CHC, F19, and F20) under different pH conditions (pH 6.6-8.0), especially Y10. Interestingly, lag periods of the conformational transitions are observed for the F19 and F20 mutants at pH 8.0, indicating the slow response of the two mutation sites on the conformational transitions. More importantly, significantly longer lag periods for F20 than for F19 indicate the conduction of the transition from F19 to F20. The conduction time (difference in lag period) decreases from 4.5 s at Cu²⁺ = 0 to undetectable (<1 ms) at Cu²⁺ = 10 μM. The significant difference in the response time of F19 and F20 and the fast local conformational changes of Y10 imply that the conformational transitions of Aβ start around Y10. MD simulations support the observation of hydrophobicity increase at N-terminal during the conformational transitions of Aβ-Cu²⁺. It also reveals that Y10 is immediately approached by Cu²⁺, supporting the speculation that the starting point of conformational transitions of Aβ is near Y10. The work has provided molecular insight into the early stage conformational transitions of Aβ₄₀ mediated by Cu²⁺.

Bulk Phase Behavior vs Interface Adsorption: Specific Multivalent Cation and Anion Effects on BSA Interactions

Proteins are ubiquitous and play a critical role in many areas from living organisms to protein microchips. In humans, serum albumin has a prominent role in the foreign body response since it is the first protein which will interact with, e.g., an implant or stent. In this study, we focused on the influence of salts (i.e., different cations (Y³⁺, La³⁺) and anions (Cl^-, I^-) on bovine serum albumin (BSA) in terms of its bulk behavior as well as the role of charges for protein adsorption at the solid-liquid interface in order to understand and control the underlying molecular mechanisms and interactions. This is part of our group's effort to gain a deeper understanding of protein-protein and protein-surface interactions in the presence of multivalent ions. In the bulk, we established two new phase diagrams and found not only multivalent cation-triggered phase transitions, but also a dependence of the protein behavior on the type of anion. The attractive interactions between proteins were observed to increase from Cl^- < NO₃^- < I^-, resulting in iodide preventing re-entrant condensation and promoting liquid-liquid phase separation in bulk. Using ellipsometry and a quartz-crystal microbalance with dissipation (QCM-D), we obtained insight into the growth of the protein adsorption layer. Importantly, we found that phase transitions at the substrate can be triggered by certain interface properties, whether they exist in the bulk solution or not. Through the use of a hydrophilic, negatively charged surface (native silica), the direct binding of anions to the interface was prevented. Interestingly, this led to re-entrant adsorption even in the absence of re-entrant condensation in bulk. However, the overall amount of adsorbed protein was enhanced through stronger attractive protein-protein interactions in the presence of iodide salts. These findings illustrate how carefully chosen surface properties and salts can directly steer the binding of anions and cations, which guide protein behavior, thus paving the way for specific/triggered protein-protein, protein-salt, and protein-surface interactions.

Physicochemical Salt Solution Parameters Limit the Survival of Planococcus halocryophilus in Martian Cryobrines

Microorganisms living in sub-zero environments can benefit from the presence of dissolved salts, as they significantly increase the temperature range of liquid water by lowering the freezing point. However, high concentrations of salts can reduce microbial growth and survival, and can evoke a physiological stress response. It remains poorly understood how the physicochemical parameters of brines (e.g. water activity, ionic strength, solubility and hydration shell strength between the ions and the surrounding water molecules) influence the survival of microorganisms. We used the cryo- and halotolerant bacterial strain Planococcus halocryophilus as a model organism to evaluate the degree of stress different salts assert. Cells were incubated in liquid media at -15°C containing single salts at eutectic concentrations (CaCl₂, LiCl, LiI, MgBr₂, MgCl₂, NaBr, NaCl, NaClO₄ and NaI). Four of these salts (LiCl, LiI, MgBr₂ and NaClO₄) were also investigated at concentrations with a low water activity (0.635) and, separately, with a high ionic strength (8 mol/L). Water activity of all solutions was measured at -15°C. This is the first time that water activity has been measured for such a large number of liquid salt solutions at constant sub-zero temperatures (-15°C). Colony-Forming Unit (CFU) counts show that the survival of P. halocryophilus has a negative correlation with the salt concentration, molecular weight of the anion and anion radius; and a positive correlation with the water activity and anions' hydration shell strength. The survival of P. halocryophilus did not show a significant correlation with the ionic strength, the molecular weight of the cation, the hydrated and unhydrated cation and hydrated anion radius, and the cations' hydration bond length. Thus, the water activity, salt concentration and anion parameters play the largest role in the survival of P. halocryophilus in concentrated brines. These findings improve our understanding of the limitations of microbial life in saline environments, which provides a basis for better evaluation of the habitability of extraterrestrial environments such as Martian cryobrines.

Replacing H+ by Na+ or K+ in phosphopeptide anions and cations prevents electron capture dissociation

By successively replacing H+ by Na+ or K+ in phosphopeptide anions and cations, we show that the efficiency of fragmentation into c and z˙ or c˙ and z fragments from N-Cα backbone bond cleavage by negative ion electron capture dissociation (niECD) and electron capture dissociation (ECD) substantially decreases with increasing number of alkali ions attached. In proton-deficient phosphopeptide ions with a net charge of 2-, we observed an exponential decrease in electron capture efficiency with increasing number of Na+ or K+ ions attached, suggesting that electrons are preferentially captured at protonated sites. In proton-abundant phosphopeptide ions with a net charge of 3+, the electron capture efficiency was not affected by replacing up to four H+ ions with Na+ or K+ ions, but the yield of c, z˙ and c˙, z fragments from N-Cα backbone bond cleavage generally decreased next to Na+ or K+ binding sites. We interpret the site-specific decrease in fragmentation efficiency as Na+ or K+ binding to backbone amide oxygen in competition with interactions of protonated sites that would otherwise lead to backbone cleavage into c, z˙ or c˙, z fragments. Our findings seriously challenge the hypothesis that the positive charge responsible for ECD into c, z˙ or c˙, z fragments can generally be a sodium or other metal ion instead of a proton.

Salt Effects on Protein Folding Thermodynamics

Salts differ in their ability to stabilize protein conformations, thereby affecting the thermodynamics and kinetics of protein folding. We developed a coarse-grained protein model that can predict salt-induced changes in protein properties by using the transfer free-energy data of various chemical groups from water to salt solutions. Using this model and molecular dynamics simulations, we probed the effect of seven different salts on the folding thermodynamics of the DNA binding domain of lac repressor protein ( lac-DBD) and N-terminal domain of ribosomal protein (NTL9). We show that a salt can act as a protein stabilizing or destabilizing agent depending on the protein sequence and folded state topology. The computed thermodynamic properties, especially the m values for various salts, which reveal the relative ability of a salt to stabilize the protein folded state, are in quantitative agreement with the experimentally measured values. The computations show that the degree of protein compaction in the denatured ensemble strongly depends on the salt identity, and for the same variation in salt concentration, the compaction in the protein dimensions varies from ∼4% to ∼30% depending on the salt. The transition-state ensemble (TSE) of lac-DBD is homogeneous and polarized, while the TSE of NTL9 is heterogeneous and diffusive. Salts induce subtle structural changes in the TSE that are in agreement with Hammond's postulate. The barrier to protein folding tends to disappear in the presence of moderate concentrations (∼3-4 m) of strongly stabilizing salts.

The Thermodynamics of Anion Complexation to Nonpolar Pockets

The interactions between nonpolar surfaces and polarizable anions lie in a gray area between the hydrophobic and Hofmeister effects. To assess the affinity of these interactions, NMR and ITC were used to probe the thermodynamics of eight anions binding to four different hosts whose pockets each consist primarily of hydrocarbon. Two classes of host were examined: cavitands and cyclodextrins. For all hosts, anion affinity was found to follow the Hofmeister series, with associations ranging from 1.6-5.7 kcal mol-1. Despite the fact that cavitand hosts 1 and 2 possess intrinsic negative electrostatic fields, it was determined that these more enveloping hosts generally bound anions more strongly. The observation that the four hosts each possess specific anion affinities that cannot be readily explained by their structures, points to the importance of counter cations and the solvation of the "empty" hosts, free guests, and host-guest complexes, in defining the affinity.

Ability of the Poisson-Boltzmann equation to capture molecular dynamics predicted ion distribution around polyelectrolytes

Here, we examine polyelectrolyte (PE) and ion chemistry specificity in ion condensation via all-atom molecular dynamics (MD) simulations and assess the ability of the Poisson-Boltzmann (PB) equation to describe the ion distribution predicted by the MD simulations. The PB model enables the extraction of parameters characterizing ion condensation. We find that the modified PB equation which contains the effective PE radius and the energy of the ion-specific interaction as empirical fitting parameters describes ion distribution accurately at large distances but close to the PE, especially when strongly localized charge or specific ion binding sites are present, the mean field description of PB fails. However, the PB model captures the MD predicted ion condensation in terms of the Manning radius and fraction of condensed counterions for all the examined PEs and ion species. We show that the condensed ion layer thickness in our MD simulations collapses on a single master curve for all the examined simple, monovalent ions (Na⁺, Br⁺, Cs⁺, Cl^-, and Br^-) and PEs when plotted against the Manning parameter (and consequently the PE line charge density). The significance of this finding is that, contrary to the Manning radius extracted from the mean field PB model, the condensed layer thickness in the all atom detail MD modelling does not depend on the PE chemistry or counterion type. Furthermore, the fraction of condensed counterions in the MD simulations exceeds the PB theory prediction. The findings contribute toward understanding and modelling ion distribution around PEs and other charged macromolecules in aqueous solutions, such as DNA, functionalized nanotubes, and viruses.

Photoelectron spectroscopy of solvated dicarboxylate and alkali metal ion clusters, M+[O2C(CH2)2CO2]2− [H2O]n (M = Na, K; n = 1–6)

We present results of combined experimental photoelectron spectroscopy and theoretical modeling studies of solvated dicarboxylate species (⁻O₂C(CH₂)₂CO₂⁻) in complex with Na⁺ and K⁺ metal cations.

Ion-induced alterations of the local hydration environment elucidate Hofmeister effect in a simple classical model of Trp-cage miniprotein

Protein stability is known to be influenced by the presence of Hofmeister active ions in the solution. In addition to direct ion-protein interactions, this influence manifests through the local alterations of the interfacial water structure induced by the anions and cations present in this region. In our earlier works it was pointed out that the effects of Hofmeister active salts on the stability of Trp-cage miniprotein can be modeled qualitatively using non-polarizable force fields. These simulations reproduced the structure-stabilization and structure-destabilization effects of selected kosmotropic and chaotropic salts, respectively. In the present study we use the same model system to elucidate atomic processes behind the chaotropic destabilization and kosmotropic stabilization of the miniprotein. We focus on changes of the local hydration environment of the miniprotein upon addition of NaClO₄ and NaF salts to the solution. The process is separated into two parts. In the first, 'promotion' phase, the protein structure is fixed, and the local hydration properties induced by the simultaneous presence of protein and ions are investigated, with a special focus on the interaction of Hofmeister active anions with the charged and polar sites. In the second, 'rearrangement' phase we follow changes of the hydration of ions and the protein, accompanying the conformational relaxation of the protein. We identify significant factors of an enthalpic and entropic nature behind the ion-induced free energy changes of the protein-water system, and also propose a possible atomic mechanism consistent with the Collins's rule, for the chaotropic destabilization and kosmotropic stabilization of protein conformation.

Determinants of Helix Formation for a Kv1.3 Transmembrane Segment inside the Ribosome Exit Tunnel

Proteins begin to fold in the ribosome, and misfolding has pathological consequences. Among the earliest folding events in biogenesis is the formation of a helix, an elementary structure that is ubiquitously present and required for correct protein folding in all proteomes. The determinants underlying helix formation in the confined space of the ribosome exit tunnel are relatively unknown. We chose the second transmembrane segment, S2, of a voltage-gated potassium channel, Kv1.3, as a model to probe this issue. Since the N terminus of S2 is initially in an extended conformation in the folding vestibule of the ribosome yet ultimately emerges at the exit port as a helix, S2 is ideally suited for delineating sequential events and folding determinants of helix formation inside the ribosome. We show that S2's extended N terminus inside the tunnel is converted into a helix by a single, distant mutation in the nascent peptide. This transition depends on nascent peptide sequence at specific tunnel locations. Co-translational secondary folding of nascent chains inside the ribosome has profound physiological consequences that bear on correct membrane insertion, tertiary folding, oligomerization, and biochemical modification of the newborn protein during biogenesis.

Characterization of long and stable de novo single alpha-helix domains provides novel insight into their stability

Naturally-occurring single α-helices (SAHs), are rich in Arg (R), Glu (E) and Lys (K) residues, and stabilized by multiple salt bridges. Understanding how salt bridges promote their stability is challenging as SAHs are long and their sequences highly variable. Thus, we designed and tested simple de novo 98-residue polypeptides containing 7-residue repeats (AEEEXXX, where X is K or R) expected to promote salt-bridge formation between Glu and Lys/Arg. Lys-rich sequences (EK3 (AEEEKKK) and EK2R1 (AEEEKRK)) both form SAHs, of which EK2R1 is more helical and thermo-stable suggesting Arg increases stability. Substituting Lys with Arg (or vice versa) in the naturally-occurring myosin-6 SAH similarly increased (or decreased) its stability. However, Arg-rich de novo sequences (ER3 (AEEERRR) and EK1R2 (AEEEKRR)) aggregated. Combining a PDB analysis with molecular modelling provides a rational explanation, demonstrating that Glu and Arg form salt bridges more commonly, utilize a wider range of rotamer conformations, and are more dynamic than Glu-Lys. This promiscuous nature of Arg helps explain the increased propensity of de novo Arg-rich SAHs to aggregate. Importantly, the specific K:R ratio is likely to be important in determining helical stability in de novo and naturally-occurring polypeptides, giving new insight into how single α-helices are stabilized.

Beyond the Hofmeister Series: Ion-Specific Effects on Proteins and Their Biological Functions

Ions differ in their ability to salt out proteins from solution as expressed in the lyotropic or Hofmeister series of cations and anions. Since its first formulation in 1888, this series has been invoked in a plethora of effects, going beyond the original salting out/salting in idea to include enzyme activities and the crystallization of proteins, as well as to processes not involving proteins like ion exchange, the surface tension of electrolytes, or bubble coalescence. Although it has been clear that the Hofmeister series is intimately connected to ion hydration in homogeneous and heterogeneous environments and to ion pairing, its molecular origin has not been fully understood. This situation could have been summarized as follows: Many chemists used the Hofmeister series as a mantra to put a label on ion-specific behavior in various environments, rather than to reach a molecular level understanding and, consequently, an ability to predict a particular effect of a given salt ion on proteins in solutions. In this Feature Article we show that the cationic and anionic Hofmeister series can now be rationalized primarily in terms of specific interactions of salt ions with the backbone and charged side chain groups at the protein surface in solution. At the same time, we demonstrate the limitations of separating Hofmeister effects into independent cationic and anionic contributions due to the electroneutrality condition, as well as specific ion pairing, leading to interactions of ions of opposite polarity. Finally, we outline the route beyond Hofmeister chemistry in the direction of understanding specific roles of ions in various biological functionalities, where generic Hofmeister-type interactions can be complemented or even overruled by particular steric arrangements in various ion binding sites.

Interaction of human serum albumin with uremic toxins: a thermodynamic study

Interaction of uremic toxins with HSA is studied by ITC and understood in terms of thermodynamic driving forces.

Differential Effects of Hydrophobic Core Packing Residues for Thermodynamic and Mechanical Stability of a Hyperthermophilic Protein

Proteins from organisms that have adapted to environmental extremes provide attractive systems to explore and determine the origins of protein stability. Improved hydrophobic core packing and decreased loop-length flexibility can increase the thermodynamic stability of proteins from hyperthermophilic organisms. However, their impact on protein mechanical stability is not known. Here, we use protein engineering, biophysical characterization, single-molecule force spectroscopy (SMFS), and molecular dynamics (MD) simulations to measure the effect of altering hydrophobic core packing on the stability of the cold shock protein TmCSP from the hyperthermophilic bacterium Thermotoga maritima. We make two variants of TmCSP in which a mutation is made to reduce the size of aliphatic groups from buried hydrophobic side chains. In the first, a mutation is introduced in a long loop (TmCSP L40A); in the other, the mutation is introduced on the C-terminal β-strand (TmCSP V62A). We use MD simulations to confirm that the mutant TmCSP L40A shows the most significant increase in loop flexibility, and mutant TmCSP V62A shows greater disruption to the core packing. We measure the thermodynamic stability (ΔGD-N) of the mutated proteins and show that there is a more significant reduction for TmCSP L40A (ΔΔG = 63%) than TmCSP V62A (ΔΔG = 47%), as might be expected on the basis of the relative reduction in the size of the side chain. By contrast, SMFS measures the mechanical stability (ΔG*) and shows a greater reduction for TmCSP V62A (ΔΔG* = 8.4%) than TmCSP L40A (ΔΔG* = 2.5%). While the impact on the mechanical stability is subtle, the results demonstrate the power of tuning noncovalent interactions to modulate both the thermodynamic and mechanical stability of a protein. Such understanding and control provide the opportunity to design proteins with optimized thermodynamic and mechanical properties.

Tuning protein mechanics through an ionic cluster graft from an extremophilic protein

Proteins from extremophilic organisms provide excellent model systems to determine the role of non-covalent interactions in defining protein stability and dynamics as well as being attractive targets for the development of robust biomaterials. Hyperthermophilic proteins have a prevalence of salt bridges, relative to their mesophilic homologues, which are thought to be important for enhanced thermal stability. However, the impact of salt bridges on the mechanical properties of proteins is far from understood. Here, a combination of protein engineering, biophysical characterisation, single molecule force spectroscopy (SMFS) and molecular dynamics (MD) simulations directly investigates the role of salt bridges in the mechanical stability of two cold shock proteins; BsCSP from the mesophilic organism Bacillus subtilis and TmCSP from the hyperthermophilic organism Thermotoga maritima. Single molecule force spectroscopy shows that at ambient temperatures TmCSP is mechanically stronger yet, counter-intuitively, its native state can withstand greater deformation before unfolding (i.e. it is mechanically soft) compared with BsCSP. MD simulations were used to identify the location and quantify the population of salt bridges, and reveal that TmCSP contains a larger number of highly occupied salt bridges than BsCSP. To test the hypothesis that salt-bridges endow these mechanical properties on the hyperthermophilic CSP, a charged triple mutant (CTM) variant of BsCSP was generated by grafting an ionic cluster from TmCSP into the BsCSP scaffold. As expected CTM is thermodynamically more stable and mechanically softer than BsCSP. We show that a grafted ionic cluster can increase the mechanical softness of a protein and speculate that it could provide a mechanical recovery mechanism and that it may be a design feature applicable to other proteins.

Effect of single point mutations in a form of systemic amyloidosis

Amyloid deposits of light-chain proteins are associated with the most common form of systemic amyloidosis. We have studied the effects of single point mutations on amyloid formation of these proteins using explicit solvent model molecular dynamics simulations. For this purpose, we compare the stability of the wild-type immunoglobulin light-chain protein REI in its native and amyloid forms with that of four mutants: R61N, G68D, D82I, and A84T. We argue that the experimentally observed differences in the propensity for amyloid formation result from two effects. First, the mutant dimers have a lower stability than the wild-type dimer due to increase exposure of certain hydrophobic residues. The second effect is a shift in equilibrium between monomers with amyloid-like structure and such with native structures. Hence, when developing drugs against light-chain associated systemic amyloidosis, one should look for components that either stabilize the dimer by binding to the dimer interface or reduce for the monomers the probability of the amyloid form.

Specific ion binding to carboxylic surface groups and the pH dependence of the Hofmeister series

Ion binding to acidic groups is a central mechanism for ion-specificity of macromolecules and surfaces. Depending on pH, acidic groups are either protonated or deprotonated and thus change not only charge but also chemical structure with crucial implications for their interaction with ions. In a two-step modeling approach, we first determine single-ion surface interaction potentials for a few selected halide and alkali ions at uncharged carboxyl (COOH) and charged carboxylate (COO(-)) surface groups from atomistic MD simulations with explicit water. Care is taken to subtract the bare Coulomb contribution due to the net charge of the carboxylate group and thereby to extract the nonelectrostatic ion-surface potential. Even at this stage, pronounced ion-specific effects are observed and the ion surface affinity strongly depends on whether the carboxyl group is protonated or not. In the second step, the ion surface interaction potentials are used in a Poisson-Boltzmann model to calculate the surface charge and the potential distribution in the solution depending on salt type, salt concentration, and solution pH in a self-consistent manner. Hofmeister phase diagrams are derived on the basis of the long-ranged forces between two carboxyl-functionalized surfaces. For cations we predict direct, reversed, and altered Hofmeister series as a function of the pH, qualitatively similar to recent experimental results for silica surfaces. The Hofmeister series reversal for cations is rationalized by a reversal of the single-cation affinity to the carboxyl group depending on its protonation state: the deprotonated carboxylate (COO(-)) surface group interacts most favorably with small cations such as Li(+) and Na(+), whereas the protonated carboxyl (COOH) surface group interacts most favorably with large cations such as Cs(+) and thus acts similarly to a hydrophobic surface group. Our results provide a general mechanism for the pH-dependent reversal of the Hofmeister series due to the different specific ion binding to protonated and deprotonated surface groups.

Improving the description of interactions between Ca2+ and protein carboxylate groups, including γ-carboxyglutamic acid: revised CHARMM22* parameters

We show that the CHARMM22* force-field over-binds the interaction between aqueous carboxylates and Ca²⁺, and introduce a modification that can recover experimentally-determined binding free energies for these systems.

Salt-specific effects observed in calorimetric studies of alkali and tetraalkylammonium salt solutions of poly(thiophen-3-ylacetic acid)

Differences in hydration of counterions are blamed for strong salt-specific effects produced upon dilution and mixing of poly(thiophene-3-ylacetic acid) salts with simple salts.

Structural conservation of the short α-helix in modified higher and lower polarity water solutions

Which conformation is preferred when the polarity of water molecules is scaled byE_P=E_LJ+S²E_C?

Electric field effects on alanine tripeptide in sodium halide solutions

The electric field effects on conformational properties of trialanine in different halide solutions were explored with long-scale molecular dynamics simulations. NaF, NaCl, NaBr and NaI solutions of low (0.2 M) and high (2 M) concentrations were exposed to a constant electric field of 1000 V/m. Generally, the electric field does not disturb trialanine's structure. Large structural changes appear only in the case of the supersaturated 2.0 M NaF solution containing NaF crystals. Although the electric field affects in a complex way, all the ions-water-peptide interactions, it predominantly affects the electroselectivity effect, which describes specific interactions such as the ion-pair formation.

Mixed-salt effects on the conformation of a short salt-bridge-forming α helix: a simulation study

The structure of a single alanine-based ACE-AEAAAKEAAAKA-NH2 peptide in explicit aqueous solutions with mixed inorganic salts (NaCl and KCl) is investigated by using molecular simulations. The concentration of Na(+), c(Na(+)), varies from 0.0M to 1.0M, whereas the concentration of K(+) is 1-c(Na(+)). The simulated peptide is very sensitive to the change of concentration ratio between Na(+) and K(+). When the concentration ratio between Na^{+} and K^{+} is changed from 0.5/0.5, the structure of the peptide becomes loose or disordered. This specific phenomenon is confirmed via checking the changes of helix parameters and mapping the free energy along different coordinates. The higher normalized probability of forming direct and indirect salt bridges between residues Glu7(+) and Lys11(+) and the smallest probability of forming ringlike structures should be responsible for the stabilized helix structure in the 0.5 Na(+)/0.5 K(+) solution. Furthermore, a noticeable conformational transition from an extended helix to an α helix is found in the 0.5 Na(+)/0.5 K(+) solution, where a local ion cloud shows that some Na(+) ions in the inner shells are still directly binding with the peptide, while K(+) in the outer shells are moving into the inner shells, keeping the peptide in the collapsed state.

Side-chain-side-chain interactions and stability of the helical state

Understanding the driving forces that lead to the stability of the secondary motifs found in proteins, namely α-helix and β-sheet, is a major goal in structural biology. The thermodynamic stability of these repetitive units is a result of a delicate balance between many factors, which in addition to the peptide chain involves also the solvent. Despite the fact that the backbones of all amino acids are the same (except of that of proline), there are large differences in the propensity of the different amino acids to promote the helical structure. In this paper, we investigate by explicit-solvent molecular dynamics simulations the role of the side chains (modeled as coarse-grained single sites) in stabilizing α helices in an aqueous solution. Our model systems include four (six-mer-nine-mer) peptide lengths in which the magnitude of the effective attraction between the side chains is systematically increased. We find that these interactions between the side chains can induce (for the nine-mer almost completely) a transition from a coil to a helical state. This transition is found to be characterized by three states in which the intermediate state is a partially folded α-helical conformation. In the absence of any interactions between the side chains the free energy change for helix formation has a small positive value indicating that favorable contributions from the side chains are necessary to stabilize the helical conformation. Thus, the helix-coil transition is controlled by the effective potentials between the side-chain residues and the magnitude of the required attraction per residue, which is on the order of the thermal energy, reduces with the length of the peptide. Surprisingly, the plots of the population of the helical state (or the change in the free energy for helix formation) as a function of the total effective interactions between the side chains in the helical state for all peptide lengths fall on the same curve.

Role of zinc and copper metal ions in amyloid β-peptides Aβ1–40 and Aβ1–42 aggregation

Conformational structural changes of Aβ_1–40 and Aβ_1–42 monomers during the interaction of Cu²⁺ and Zn²⁺ metal ions.

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.

The stability of cylindrin β-barrel amyloid oligomer models-a molecular dynamics study

Small-soluble amyloid oligomers are believed to play a significant role in the pathology of amyloid diseases. Recently, the atomic structure of a toxic oligomer formed by an 11 residue and its tandem repeat was found to have an out-off register antiparallel β-strands in the shape of a β-barrel. In the present article we investigate the effect of mutations in the hydrophobic cores on the structure and dynamic of the β-barrels using all atom multiple molecular dynamics simulations with an explicit solvent. Extending previous experiments with molecular dynamics simulations we systematically test how stability and formation of cylindrin depends on the interplay between hydrophobicity and steric effects of the core residues. We find that strong hydrophobic interactions between geometrically fitting residues keep the strands (both in register and out-off-register interface) in close proximity, which in turn stabilizes the side-chain and main-chain hydrogen bonds, and the salt bridges on the outer surface along the weak out-of-register interface. Our simulations also indicate presence of water molecules in the hydrophobic interior of the cylindrin β-barrel.Proteins 2013.