Polar or Apolar—The Role of Polarity for Urea-Induced Protein Denaturation

Urea-induced protein denaturation is widely used to study protein folding and stability; however, the molecular mechanism and driving forces of this process are not yet fully understood. In particular, it is unclear whether either hydrophobic or polar interactions between urea molecules and residues at the protein surface drive denaturation. To address this question, here, many molecular dynamics simulations totalling ca. 7 µs of the CI2 protein in aqueous solution served to perform a computational thought experiment, in which we varied the polarity of urea. For apolar driving forces, hypopolar urea should show increased denaturation power; for polar driving forces, hyperpolar urea should be the stronger denaturant. Indeed, protein unfolding was observed in all simulations with decreased urea polarity. Hyperpolar urea, in contrast, turned out to stabilize the native state. Moreover, the differential interaction preferences between urea and the 20 amino acids turned out to be enhanced for hypopolar urea and suppressed (or even inverted) for hyperpolar urea. These results strongly suggest that apolar urea–protein interactions, and not polar interactions, are the dominant driving force for denaturation. Further, the observed interactions provide a detailed picture of the underlying molecular driving forces. Our simulations finally allowed characterization of CI2 unfolding pathways. Unfolding proceeds sequentially with alternating loss of secondary or tertiary structure. After the transition state, unfolding pathways show large structural heterogeneity.     

Swelling of a Lecithin Lamellar Phase Induced by Small Carbohydrate Solutes

In this paper, we consider the effect of adding small carbohydrate solutes (small sugars) to DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) L_α dispersions and the consequences on the force balance at zero osmotic pressure (maximal swelling). We show the importance of long incubations required to obtain samples at thermodynamic equilibrium where molecular diffusion has been completed. The monotonic increase of maximal swelling versus sugar content occurs as a combined effect of the screening of the van der Waals contribution and fluctuations in the lamellar stacks. According to this new approach, it is shown that changes in dielectric properties result in a much less pronounced effect than entropic forces (undulations) generated by the softening of the membranes at high sugar content. However, this sugar-induced swelling cannot be explained quantitatively by adding an entropic contribution to molecular interactions. Quantitative disagreement between the proposed mechanism and our observations is due either to nonadditivity of molecular interactions with entropic forces or to the relation used to account for the entropic contribution.     

Assessing the Intricate Balance of Intermolecular Interactions upon Self-Association of Intrinsically Disordered Proteins

Attractive interactions between intrinsically disordered proteins can be crucial for the functionality or, on the contrary, lead to the formation of harmful aggregates. For obtaining a molecular understanding of intrinsically disordered proteins and their interactions, computer simulations have proven to be a valuable complement to experiments. In this study, we present a coarse-grained model and its applications to a system dominated by attractive interactions, namely, the self-association of the saliva protein Statherin. SAXS experiments show that Statherin self-associates with increased protein concentration, and that both an increased temperature and a lower ionic strength decrease the size of the formed complexes. The model captures the observed trends and provides insight into the size distribution. Hydrophobic interaction is considered to be the major driving force of the self-association, while electrostatic repulsion represses the growth. In addition, the model suggests that the decrease of association number with increased temperature is of entropic origin.     

Dynamics of chromatin decondensation reveals the structural integrity of a mechanically prestressed nucleus

Genome organization within the cell nucleus is a result of chromatin condensation achieved by histone tail-tail interactions and other nuclear proteins that counter the outward entropic pressure of the polymeric DNA. We probed the entropic swelling of chromatin driven by enzymatic disruption of these interactions in isolated mammalian cell nuclei. The large-scale decondensation of chromatin and the eventual rupture of the nuclear membrane and lamin network due to this entropic pressure were observed by fluorescence imaging. This swelling was accompanied by nuclear softening, an effect that we quantified by measuring the fluctuations of an optically trapped bead adhered onto the nucleus. We also measured the pressure at which the nuclear scaffold ruptured using an atomic force microscope cantilever. A simple theory based on a balance of forces in a swelling porous gel quantitatively explains the diffusive dynamics of swelling. Our experiments on decondensation of chromatin in nuclei suggest that its compaction is a critical parameter in controlling nuclear stability.     

Molecular origin of biphasic response of main phase-transition temperature of phospholipid membranes to long-chain alcohols

A statistical mechanical theory is proposed which explains the molecular mechanism of the nonlinear response of the phase-transition temperature of phospholipid vesicle membranes to added 1-alkanols. By assuming that the free energy of transfer of 1-alkanols from the aqueous phase to the membrane and the interaction energy between 1-alkanol molecules are linear functions of alkanol alkyl chain-length, the nonlinear behavior is explained in the Bragg-Williams approximation. For dipalmitoylphosphatidylcholine vesicle membranes, the theory reveals a larger free energy of transfer of 1-alkanols from the aqueous phase to the solid-gel membrane than to the liquid-crystalline membrane when the number of carbon atoms of 1-alkanol exceeds 12. When the intermolecular interaction force between 1-alkanol molecules residing in the gel phase is stronger than the interaction force between those residing in the liquid-crystalline phase, the ligand effect is to tighten the lipid matrix structure, causing the transition temperature to rise. The interaction force is a quadratic function of 1-alkanol concentration; hence, the response of the transition temperature to the 1-alkanol concentration is nonlinear. At low concentrations of the long-chain 1-alkanols that predominantly elevate the transition temperature, this intermolecular interaction force is negligible. In this case, the entropic effect of the incorporated ligand molecules, which loosens the lipid matrix, predominates, and the transition temperature decreases. The biphasic action of long-chain 1-alkanols originates from the balance of these two opposing effects: entropy and intermolecular interaction.     

Non-polar solutes in water and in aqueous solutions of protein denaturants. Modeling of solution and transfer processes

A simple molecular model for the thermodynamic behavior of non-polar solutes in water and in aqueous solutions of protein denaturants is presented. Three contributions are considered: (i) combinatorial arising from the mixing process, (ii) interactional characterizing the molecular interactions occurring in the mixture and (iii) a contribution originating from the structural changes occurring in the first shell of water molecules around the solute. The latter is modeled assuming that water molecules in contact with the solute are involved in a chemical equilibrium between two states. The model describes well the temperature and denaturant concentration dependences of the Gibbs energies of solution and transfer for benzene, toluene and alkanes in water and aqueous solutions of urea and guanidine hydrochloride. Model parameters are physically meaningful, allowing a discussion of the molecular interactions involved. A preferential solvation of the solute by the denaturant is found. However, the non-polar solute-denaturant interaction is not specific, i.e. leading to a distinct chemical entity. Urea and guanidine hydrochloride are non-polar solubilizing agents because their interactions with the solute are less unfavorable than those between water and the solute.     

Effects of organic solvents and substrate binding on trypsin in acetonitrile and hexane media

In this work, we used molecular dynamic (MD) simulation to study trypsin with and without a six-amino-acid peptide bound in three different solvents (water, acetonitrile and hexane) in order to provide molecular information for well understanding the structure and function of enzymes in non-aqueous media. The results show that the enzyme is more compact and less native-like in hexane than in the other two polar solvents. The substrate could stabilize the native protein structure in the two polar media, but not in the non-polar hexane. There are no significant differences in the conformation of the S1 pocket upon the substrate binding in water and acetonitrile media while a reverse behavior is observed in hexane media, implying a possible induced fit binding mechanism in the non-polar media. The substrate binding enhances the stability of catalytic H-bond network since it could expel the solvent molecules from the active site. The enzyme and the substrate appear to be more appropriate to the reactive conformation in the organic solvents compared with aqueous solution. There is much greater substrate binding strength in hexane media than the water and acetonitrile ones since the polar solvent significantly weakens electrostatic interactions, which are observed to be the main driving force to the binding. In addition, some residues of the S1 pocket could remain favorable contribution to the binding despite the solvent change, but with differences in the contribution extent, the number and the type of residues between the three media.

Charge effects on the fibril-forming peptide KTVIIE: a two-dimensional replica exchange simulation study

The assembly of peptides into ordered nanostructures is increasingly recognized as both a bioengineering tool for generating new materials and a critical aspect of aggregation processes that underlie neurological diseases such as Alzheimer's disease, Parkinson's disease, and Huntington's disease. There is a major problem in understanding how extremely subtle sequence changes can lead to profound and often unexpected differences in self-assembly behavior. To better delineate the complex interplay of different microscopic driving forces in such cases, we develop a methodology to quantify and compare the propensity of different peptide sequences to form small oligomers during early self-assembly stages. This umbrella-sampling replica exchange molecular dynamics method performs a replica exchange molecular dynamics simulation along peptide association reaction coordinates using umbrella restraints. With this method, we study a set of sequence-similar peptides that differ in net charge: K(+)TVIIE(-), K(+)TVIIE, and (+)K(+)TVIIE. Interestingly, experiments show that only the monovalent peptide, K(+)TVIIE, forms fibrils, whereas the others do not. We examine dimer, trimer, and tetramer formation processes of these peptides, and compute high-accuracy potential of mean force association curves. The potential of mean forces recapitulate a higher stability and equilibrium constant of the fibril-forming peptide, similar to experiment, but reveal that entropic contributions to association free energies can play a surprisingly significant role. The simulations also show behavior reminiscent of experimental aggregate polymorphism, revealed in multiple stable conformational states and association pathways. Our results suggest that sequence changes can have significant effects on self-assembly through not only direct peptide-peptide interactions but conformational entropies and degeneracies as well.     

Adsorption characteristics of lignosulfonates in salt-free and salt-added aqueous solutions

Five sodium lignosulfonate (SL) fractions with narrow molecular weight distribution and known salt content were used as the polyanion to build up layer-by-layer self-assembly multilayers with poly(diallyldimethylammonium chloride) (PDAC) as polycation. It is interesting to find that the salt-free SL is hardly adsorbed on the PDAC surface, but the SL in salt-added solutions can be self-assembled well with PDAC to form SL/PDAC multilayers. When the five SL fractions dissolved in saline solutions are adsorbed on the PDAC surface by a self-assembly technique, SL with higher M(w) shows a higher adsorption amount than does SL with lower M(w). The driving forces of self-assembly of SL and PDAC are discussed based on the solution behaviors and adsorption characteristics of SL in salt-free and salt-added aqueous solutions. A possible self-assembled mechanism of SL and PDAC is electrostatic or cation-π interactions, but the SL cannot be adsorbed onto the PDAC surface without a hydrophobic interaction. With the addition of enough salt, the Coulomb interaction of SL becomes negligible, but the adsorption amount increases, indicating that the electrostatic interaction is not the main driving force of SL/PDAC self-assembly. For adsorption of SL in saline solution onto the PDAC surface, the cation-π interaction is the main driving force, and the hydrophobic interaction plays an important role in the adsorbed amount.     

Investigation into the mechanism of (-)-epigallocatechin-3-gallate-induced precipitation of insulin

The molecular interactions between EGCG and insulin were investigated to probe the mechanism of EGCG-induced insulin precipitation. The results indicated that 1-5mM EGCG induced insulin into reversible globular precipitates of 185-365 nm. The formation of precipitates was facilitated at high salt concentration and pH values close to insulin's isoelectric point, indicating that hydrophobic interaction was the main driving force. The precipitation was positively related to insulin concentration, but for EGCG, there was a suitable concentration (2 mM at 2 mg/mL of insulin) at which the precipitate content reached maximum. Mass spectroscopy analysis indicated that EGCG formed clusters in the aqueous solution and the clusters correlate with the insulin precipitation. Based on extensive investigation, a physical model was proposed to explain the molecular interactions between EGCG and insulin. Namely, EGCG monomers and clusters first bound to insulin dimers via hydrophobic interaction, leading to the reduction of the thickness of the hydration layer and the partial denaturation of insulin. Then, EGCG clusters acted as bridges to induce the aggregation and precipitation of insulin.     

Insight into a molecular interaction force supporting peptide backbones and its implication to protein loops and folding

Although not being classified as the most fundamental protein structural elements like α-helices and β-strands, the loop segment may play considerable roles for protein stability, flexibility, and dynamic activity. Meanwhile, the protein loop is also quite elusive; i.e. its interactions with the other parts of protein as well as its own shape-maintaining forces have still remained as a puzzle or at least not quite clear yet. Here, we report a molecular force, the so-called polar hydrogen–π interaction (Hp–π), which may play an important role in supporting the backbones of protein loops. By conducting the potential energy surface scanning calculations on the quasi π-plane of peptide bond unit, we have observed the following intriguing phenomena: (1) when the polar hydrogen atom of a peptide unit is perpendicularly pointing to the π-plane of other peptide bond units, a remarkable Hp–π interaction occurs; (2) the interaction is distance and orientation dependent, acting in a broad space, and belonging to the ‘point-to-plane’ one. The molecular force reported here may provide useful interaction concepts and insights into better understanding the loop’s unique stability and flexibility feature, as well as the driving force of the protein global folding.     

Backbone and side-chain contributions in protein denaturation by urea

Urea is a commonly used protein denaturant, and it is of great interest to determine its interaction with various protein groups to elucidate the molecular basis of its effect on protein stability. Using the Trp-cage miniprotein as a model system, we report what we believe to be the first computation of changes in the preferential interaction coefficient of the protein upon urea denaturation from molecular-dynamics simulations and examine the contributions from the backbone and the side-chain groups. The preferential interaction is obtained from reversible folding/unfolding replica exchange molecular-dynamics simulations of Trp-cage in presence of urea, over a wide range of urea concentration. The increase in preferential interaction upon unfolding is dominated by the side-chain contribution, rather than the backbone. Similar trends are observed in simulations using two different force fields, Amber94 and Amber99sb, for the protein. The magnitudes of the side-chain and backbone contributions differ in the two force fields, despite containing identical protein-solvent interaction terms. The differences arise from the unfolded ensembles sampled, with Amber99sb favoring conformations with larger surface area and lower helical content. These results emphasize the importance of the side-chain interactions with urea in protein denaturation, and highlight the dependence of the computed driving forces on the unfolded ensemble sampled.     

Direct force measurements between siRNA and chitosan molecules using force spectroscopy

Information on the interaction strength between small interfering RNA (siRNA) and chitosan can contribute to the understanding of the formation and stability of chitosan/siRNA nanoparticles used as siRNA delivery systems for gene silencing. In this study, we utilize atomic force microscopy to obtain force spectroscopy results of the interaction strengths between siRNA and chitosan measured in physiological phosphate buffered saline buffer at different pH. The force measurements revealed that the adhesive interactions decreased in force strength and force frequency as the pH was increased from 4.1 to 6.1, 7.4, and 9.5, exhibiting distinct multimodal distributions of the interaction forces between siRNA and chitosan molecules at acidic pH and only negligible adhesive forces were observed at neutral or high pH. The strong pH dependence of siRNA-chitosan interactions can provide a convincing rationale for siRNA/chitosan complex formation and nanoparticle stability under low acidic conditions. These findings demonstrate that the use of force spectroscopy for the adhesive force measurements allows an evaluation of the complexing ability between siRNA and chitosan that can be utilized to predict nanoparticle stability.     

Modeling the role of disulfide bonds in protein folding: Entropic barriers and pathways

The role of disulfide bonds in directing protein folding is studied using lattice models. We find that the stability and the specificity of the disulfide bond interactions play quite different roles in the folding process: Under some conditions, the stability decreases the overall rate of folding; the specificity, however, by yielding a simpler connectivity of intermediates, always increases the rate of folding. This conclusion is intimately related to the selection mechanism entailed by entropic driving forces, such as the loop formation probability, and entropic barriers separating the native and the many native-like metastable states. The folding time is found to be a minimum for a certain range of the effective disulfide bond interaction. Examination of a model, which allows for the formation of disulfide bonded intermediates, suggests that folding proceeds via a threestage multiple pathways kinetics. We show that there are pathways to the native state involving only native-like intermediates, as well as those that are mediated by nonnative intermediates. These findings are interpreted in terms of the appropriate energy landscape describing the barriers connecting low energy conformations. The consistency of our conclusions with several experimental studies is also discussed. © 1995 Wiley-Liss, Inc.     

Interpreting protein/DNA interactions: distinguishing specific from non-specific and electrostatic from non-electrostatic components

We discuss the effectiveness of existing methods for understanding the forces driving the formation of specific protein-DNA complexes. Theoretical approaches using the Poisson-Boltzmann (PB) equation to analyse interactions between these highly charged macromolecules to form known structures are contrasted with an empirical approach that analyses the effects of salt on the stability of these complexes and assumes that release of counter-ions associated with the free DNA plays the dominant role in their formation. According to this counter-ion condensation (CC) concept, the salt-dependent part of the Gibbs energy of binding, which is defined as the electrostatic component, is fully entropic and its dependence on the salt concentration represents the number of ionic contacts present in the complex. It is shown that although this electrostatic component provides the majority of the Gibbs energy of complex formation and does not depend on the DNA sequence, the salt-independent part of the Gibbs energy--usually regarded as non-electrostatic--is sequence specific. The CC approach thus has considerable practical value for studying protein/DNA complexes, while practical applications of PB analysis have yet to demonstrate their merit.     

Evaluation of intermolecular forces in a circulating system

Intercellular interactions, which are mediated by a variety of complex intercellular molecules through the processes of formation and dissociation of molecular bonds, play a critical role in regulating cellular functions in biological systems. Various approaches are applied to evaluate intercellular or molecular bonding forces. To quantify the intermolecular interaction forces, flow chamber has become a meaningful technique as it can ultimately mimic the cellular microenvironment in vivo under physiological flow conditions. Hydrodynamic forces are usually used to predict the intercellular forces down to the single molecular level. However, results show that only using hydrodynamic force will overestimate up to 30% of the receptor-ligand strength when the non-specific forces such as Derjaguin-Landau-Verway-Overbeek (DLVO) forces become un-neglected. Due to the nature of high ion concentration in the physiological condition, electrostatic force is largely screened which will cause DLVO force unbalanced. In this study, we propose to take account of the DLVO force, including van der Waals (VDW) force and electrostatic force, to predict the intermolecular forces of a cell doublet and cell-substrate model in a circulating system. Results also show that the DLVO force has a nonlinear effect as the cell-cell or cell-substrate distance changes. In addition, we used the framework of high accuracy hydrodynamic theories proved in colloidal systems. It is concluded that DLVO force could not be ignored in quantitative studies of molecular interaction forces in circulating system. More accurate prediction of intercellular forces needs to take account of both hydrodynamic force and DLVO force.     

The key role of solvent in condensation: Mapping water in liquid-liquid phase-separated FUS

Formation of biomolecular condensates through liquid-liquid phase separation (LLPS) has emerged as a pervasive principle in cell biology, allowing compartmentalization and spatiotemporal regulation of dynamic cellular processes. Proteins that form condensates under physiological conditions often contain intrinsically disordered regions with low-complexity domains. Among them, the RNA-binding proteins FUS and TDP-43 have been a focus of intense investigation because aberrant condensation and aggregation of these proteins is linked to neurodegenerative diseases such as amyotrophic lateral sclerosis and frontotemporal dementia. LLPS occurs when protein-rich condensates form surrounded by a dilute aqueous solution. LLPS is per se entropically unfavorable. Energetically favorable multivalent protein-protein interactions are one important aspect to offset entropic costs. Another proposed aspect is the release of entropically unfavorable preordered hydration water into the bulk. We used attenuated total reflection spectroscopy in the terahertz frequency range to characterize the changes in the hydrogen bonding network accompanying the FUS enrichment in liquid-liquid phase-separated droplets to provide experimental evidence for the key role of the solvent as a thermodynamic driving force. The FUS concentration inside LLPS droplets was determined to be increased to 2.0 mM independent of the initial protein concentration (5 or 10 μM solutions) by fluorescence measurements. With terahertz spectroscopy, we revealed a dewetting of hydrophobic side chains in phase-separated FUS. Thus, the release of entropically unfavorable water populations into the bulk goes hand in hand with enthalpically favorable protein-protein interaction. Both changes are energetically favorable, and our study shows that both contribute to the thermodynamic driving force in phase separation.     

Energetics of membrane protein folding and stability

The critical role of membrane proteins in a myriad of biological and physiological functions has spawned numerous investigations over the past several decades with the long-term goal of identifying the molecular origins and energetic forces that stabilize these proteins within the membrane. Parallel structural and thermodynamics studies on several systems have provided significant insight regarding the driving forces governing folding, assembly, insertion, and translocation of membrane proteins. The present review surveys families of membrane-associated proteins including α-helical and β-barrel structures, viral surface receptors, and pore-forming toxins, citing representative proteins within each of these classes for further scrutiny in terms of structure-function relationships and global conformational stability. This overview presents seminal findings from pioneering studies on the energetics of membrane protein folding and stability to modern techniques that are exploiting the use of molecular genetics and single molecule studies. An overall consensus regarding the molecular origins of membrane protein stability is that a number of intrinsic properties resemble features of soluble proteins, yet there are distinct energetic differences arising from specific intra- and intermolecular interactions within the membrane. The combined efforts from structural, energetics, and dynamics approaches offer unique insights and improve our fundamental understanding of the driving forces dictating membrane protein folding and stability.     

Thermodynamics of site-specific small molecular ion interactions with DNA duplex: a molecular dynamics study

The stability and dynamics of a double-stranded DNA (dsDNA) is affected by the preferential occupancy of small monovalent molecular ions. Small metal and molecular ions such as sodium and alkyl ammonium have crucial biological functions in human body, affect the thermodynamic stability of the duplex DNA and exhibit preferential binding. Here, using atomistic molecular dynamics simulations, we investigate the preferential binding of metal ion such as Na⁺ and molecular ions such as tetramethyl ammonium (TMA⁺) and 2-hydroxy-N,N,N-trimethylethanaminium (CHO⁺) to double-stranded DNA. The thermodynamic driving force for a particular molecular ion-DNA interaction is determined by decomposing the free energy of binding into its entropic and enthalpic contributions. Our simulations show that each of these molecular ions preferentially binds to the minor groove of the DNA and the extent of binding is highest for CHO⁺. The ion binding processes are found to be entropically favourable. In addition, the contribution of hydrophobic effects towards the entropic stabilisation (in case of TMA⁺) and the effect of hydrogen bonding contributing to enthalpic stabilisation (in case of CHO⁺) have also been investigated.     

Aromatic interactions in model systems

A thorough knowledge of noncovalent interactions is crucial to the understanding of biological complexity. One of the less well understood but significant weak interactions in nature is the aromatic interaction. Recent studies have provided new insight into the driving force, stability and selectivity of these interactions. The contribution of solvophobic and electrostatic interactions have been shown to be inextricably linked. Moreover, the influence of electrostatic and solvophobic components on the selectivity of aromatic interactions has been demonstrated.