Unfolding of hydrophobic polymers in guanidinium chloride solutions

Guanidine Hydrochloride-Induced Hepatitis B Virus Capsid Disassembly Hysteresis

Hepatitis B virus (HBV) displays remarkable self-assembly capabilities that interest the scientific community and biotechnological industries as HBV is leading to an annual mortality of up to 1 million people worldwide (especially in Africa and Southeast Asia). When the ionic strength is increased, hepatitis B virus-like particles (VLPs) can assemble from dimers of the first 149 residues of the HBV capsid protein core assembly domain (Cp149). Using solution small-angle X-ray scattering, we investigated the disassembly of the VLPs by titrating guanidine hydrochloride (GuHCl). Measurements were performed with and without 1 M NaCl, added either before or after titrating GuHCl. Fitting the scattering curves to a linear combination of atomic models of Cp149 dimer (the subunit) and T = 3 and T = 4 icosahedral capsids revealed the mass fraction of the dimer in each structure in all the titration points. Based on the mass fractions, the variation in the dimer-dimer association standard free energy was calculated as a function of added GuHCl, showing a linear relation between the interaction strength and GuHCl concentration. Using the data, we estimated the energy barriers for assembly and disassembly and the critical nucleus size for all of the assembly reactions. Extrapolating the standard free energy to [GuHCl] = 0 showed an evident hysteresis in the assembly process, manifested by differences in the dimer-dimer association standard free energy obtained for the disassembly reactions compared with the equivalent assembly reactions. Similar hysteresis was observed in the energy barriers for assembly and disassembly and the critical nucleus size. The results suggest that above 1.5 M, GuHCl disassembled the capsids by attaching to the protein and adding steric repulsion, thereby weakening the hydrophobic attraction.

Nanopore single-molecule biosensor in protein denaturation analysis

Unclear issues in protein studies include but not limited to the stability and denaturation mechanism in the presence of denaturants. Herein, we report a dynamic monitoring approach based on nanopore single-molecule biosensor, which can detect the protein's folding and unfolding transitions by recording a nanopore ionic current. When gradually increasing the concentration of denaturant guanidine hydrochloride (GdmCl), sensitive responses were observed with lysozyme unfolding. The emergence of the featured biphasic-pulse demonstrated the existence of a stable intermediate. It was the first time to experimentally confirm the dynamic equilibrium between the intermediate and the native states at single molecule level, therefore consolidating the standpoint of lysozyme denaturation process following the three-state model. Additionally, we got more insights into the conformation about the intermediate as globular-like structure, larger gyration radius, and enhanced positive charge density. We considered that the manner of denaturant toward lysozyme adopts the "direct" model based on stronger electrostatic and van der Waals forces. Nanopore biosensor exhibited excellent sensitivity with a low detection concentration of 280 pM and reproducibility in analysing the folding intermediate of lysozyme.

Towards a predictive model for polymer solubility using the noncovalent interaction index: polyethylene as a case study

In this work we present the development of a novel, quantitative solubility descriptor based on the non-covalent interaction index. It is presented as a more insightful alternative to Hansen's solubility parameters and the COSMO model to assess and predict polymer solubility in different solvents. To this end, we studied the solvation behaviour as a function of the chain length of a single chain of arguably the most simple polymer, polyethylene, in anisole (solvent) and methanol (poor solvent) via molecular dynamics simulations. It was found that in anisole the solute maximized its interface with the solvent, whereas in methanol the macromolecule formed rod-like structures by folding on itself once the chain length surpassed a certain barrier. We assessed this behaviour - which can be related to solubility - quantitatively and qualitatively via well-known descriptors, namely the solvation free energy, and the solvent accessible surface area. In addition, we propose the non-covalent interaction (NCI) index as a versatile descriptor, providing information on the strength, as well as the nature, of the solute-solvent interactions, the solute's intramolecular interactions and on the solute's conformation, both qualitatively and quantitatively. Finally, as a quantitative measure for solubility, defined in this context as the solute's tendency to maximize its interactions with the solvent, we propose two new NCI-based descriptors: the relative integrated NCI density and the integrated NCI difference. The former represents the quantitative difference in solute-solvent interactions between a fully extended coil and the actual conformation during simulation and the latter the quantitative difference between the intermolecular (solute-solvent) and the intramolecular (in the solute) non-covalent interactions. The easy interpretation and calculation of these novel quantities open up the possibility of fast, reliable and insightful high-throughput screening of different (anti)solvent and solute combinations.

General Counteraction Exerted by Sugars against Denaturants

The conformational stability of globular proteins is strongly influenced by the addition to water of different co-solutes. Some of the latter destabilize the native state, while others stabilize it. It is emerging that stabilizing agents are able to counteract the action of destabilizing agents. We have already provided experimental evidence that this counteraction is a general phenomenon and offered a rationalization. In the present work, we show that four different sugars, namely fructose, glucose, sucrose, and trehalose, counteract the effect of urea, tetramethylurea, sodium perchlorate, guanidinium chloride, and guanidinium thiocyanate despite the chemical and structural differences of those destabilizing agents. The rationalization we provide is as follows: (a) the solvent-excluded volume effect, a purely entropic effect, stabilizes the native state, whose solvent-accessible surface area is smaller than the one of denatured conformations; (b) the magnitude of the solvent-excluded volume effect increases markedly in ternary solutions because the experimental density of such solutions is larger than that of pure water.

Intercalation-deintercalation of water-in-salt electrolytes in nanoscale hydrophobic confinement

Intercalation-deintercalation of water-in-salt (WIS) electrolytes in nanoscale confinement is an important phenomenon relevant to energy storage and self-assembly applications. In this article, we use molecular simulations to investigate the effects of intersurface separation on the structure and free energy underlying the intercalation-deintercalation of the Li bis(trifluoromethane)sulfonimide ([Li][TFSI]) water-in-salt (WIS) electrolyte confined between nanoscale hydrophobic surfaces. We employ enhanced sampling to estimate the free energy profiles for the intercalation behaviour of WIS in confining sheets at several intersurface separations. We observe that the relative stability of the condensed and vapour phases of WIS in the confinement depends on the separation between the confining surfaces and the WIS concentration. We find that the critical separation at which the condensed and vapour phases are equally stable in confinement depends on the concentration of WIS. The relative height of the free energy barrier also strongly depends on the concentration of [Li][TFSI] inside the confined space, and we find that this concentration dependence can be attributed to changes in line tension. The process of deintercalation passes through vapour tube formation inside the confined space, and this process is initiated by vapour bubble formation. The size of the critical vapour tube required for spontaneous evaporation of WIS from the confinement is also found to depend on the intersurface separation and WIS concentration.

Exploring the role of hydrophilic amino acids in unfolding of protein in aqueous ethanol solution

Hydrophobic association is the key contributor behind the formation of well packed core of a protein which is often believed to be an important step for folding from an unfolded chain to its compact functional form. While most of the protein folding/unfolding studies have evaluated the changes in the hydrophobic interactions during chemical denaturation, the role of hydrophilic amino acids in such processes are not discussed in detail. Here we report the role of the hydrophilic amino acids behind ethanol induced unfolding of protein. Using free energy simulations, we show that chicken villin head piece (HP-36) protein unfolds gradually in presence of water-ethanol binary mixture with increasing composition of ethanol. However, upon mutation of hydrophilic amino acids by glycine while keeping the hydrophobic amino acids intact, the compact state of the protein is found to be stable at all compositions with gradual flattening of the free energy landscape upon increasing compositions. The local environment around the protein in terms of ethanol/water number significantly differs in wild type protein compared to the mutated protein. The calculated Wyman-Tanford preferential binding coefficient of ethanol for wild type protein reveals that a greater number of cosolutes (here ethanol) bind to the unfolded state compared to its folded state. However, no significant increase in binding coefficient of ethanol at the unfolded state is found for mutated protein. Local-bulk partition coefficient calculation also suggests similar scenarios. Our results reveal that the weakening of hydrophobic interactions in aqueous ethanol solution along with larger preferential binding of ethanol to the unfolded state mediated by hydrophilic amino acids combinedly helps unfolding of protein in aqueous ethanol solution.

Distinct and Nonadditive Effects of Urea and Guanidinium Chloride on Peptide Solvation

Using enhanced-sampling replica exchange fully atomistic molecular dynamics simulations, we show that, individually, urea and guanidinium chloride (GdmCl) denature the Trpcage protein, but remarkably, the helical segment ¹NLYIQWL⁷ of the protein is stabilized in mixed denaturant solutions. GdmCl induces protein denaturation via a combination of direct and indirect effects involving dehydration of the protein and destabilization of stabilizing salt bridges. In contrast, urea denatures the protein through favorable protein-urea preferential interactions, with peptide-specific indirect effects of urea on the water structure around the protein. In the case of the helical segment of Trpcage, urea "oversolvates" the peptide backbone by reorganizing water molecules from the peptide side chains to the peptide backbone. An intricate nonadditive thermodynamic balance between GdmCl-induced dehydration of the peptide and the urea-induced changes in solvation structure triggers partial counteraction to urea denaturation and stabilization of the helix.

Counteraction of denaturant-induced protein unfolding is a general property of stabilizing agents

DSC measurements on RNase A at neutral pH show that five stabilizing agents, namely trimethylamine N-oxide, glucose, sucrose, betaine and sodium sulfate, can counteract the destabilizing action of urea, sodium perchlorate, guanidinium chloride and guanidinium thiocyanate. This is an important finding inferring that counteraction has a common physical origin, regardless of the chemical differences among the stabilizing agents and among the destabilizing ones. A rationalization is provided grounded on the following line of reasoning: (a) the decrease in solvent-excluded volume effect is the main stabilizing contribution of the native state; (b) its magnitude increases on increasing the density of the aqueous solution; (c) the density increases significantly in the ternary solutions containing water, a stabilizing agent and a destabilizing one, as indicated by the present experimental data.

Anionic and cationic Hofmeister effects are non-additive for guanidinium salts

To understand specific ion effects on a molecular level we explore the effect of salts on the rotational mobility of a model amide using dielectric spectroscopy. Based on our previous studies on the effect of strong denaturing anions or cations, here we study the additivity of the anionic and cationic effect. Using salts consisting of denaturing spherical anions and spherical cations we find such salts to affect the amide according to what one expects based on the additive activity of the individual ions. The guanidinium (Gdm⁺) cation appears to be a notable exception, as our results suggest that GdmI (and accordingly GdmSCN) is less efficient in hindering the rotation of the amide than KI or GdmCl.

Conformations of a Metastable SH3 Domain Characterized by smFRET and an Excluded-Volume Polymer Model

Conformational states of the metastable drkN SH3 domain were characterized using single-molecule fluorescence techniques. Under nondenaturing conditions, two Förster resonance energy transfer (FRET) populations were observed that corresponded to a folded and an unfolded state. FRET-estimated radii of gyration and hydrodynamic radii estimated by fluorescence correlation spectroscopy of the two coexisting conformations are in agreement with previous ensemble x-ray scattering and NMR measurements. Surprisingly, when exposed to high concentrations of urea and GdmCl denaturants, the protein still exhibits two distinct FRET populations. The dominant conformation is expanded, showing a low FRET efficiency, consistent with the expected behavior of a random chain with excluded volume. However, approximately one-third of the drkN SH3 conformations showed high, nearly 100%, FRET efficiency, which is shown to correspond to denaturation-induced looped conformations that remain stable on a timescale of at least 100 μs. These loops may contain interconverting conformations that are more globally collapsed, hairpin-like, or circular, giving rise to the observed heterogeneous broadening of this population. Although the underlying mechanism of chain looping remains elusive, FRET experiments in formamide and dimethyl sulfoxide suggest that interactions between hydrophobic groups in the distal regions may play a significant role in the formation of the looped state.

Quantifying transient interactions between amide groups and the guanidinium cation

We study the interaction of the guanidinium cation, a widely used protein denaturant, with amide groups, the common structural motif of proteins. Our results provide evidence for direct contact between guanidinium and ∼2 amide groups, but the interaction is transient and weaker than for other cations with high charge-density.

How osmolytes influence hydrophobic polymer conformations: A unified view from experiment and theory

It is currently the consensus belief that protective osmolytes such as trimethylamine N-oxide (TMAO) favor protein folding by being excluded from the vicinity of a protein, whereas denaturing osmolytes such as urea lead to protein unfolding by strongly binding to the surface. Despite there being consensus on how TMAO and urea affect proteins as a whole, very little is known as to their effects on the individual mechanisms responsible for protein structure formation, especially hydrophobic association. In the present study, we use single-molecule atomic force microscopy and molecular dynamics simulations to investigate the effects of TMAO and urea on the unfolding of the hydrophobic homopolymer polystyrene. Incorporated with interfacial energy measurements, our results show that TMAO and urea act on polystyrene as a protectant and a denaturant, respectively, while complying with Tanford-Wyman preferential binding theory. We provide a molecular explanation suggesting that TMAO molecules have a greater thermodynamic binding affinity with the collapsed conformation of polystyrene than with the extended conformation, while the reverse is true for urea molecules. Results presented here from both experiment and simulation are in line with earlier predictions on a model Lennard-Jones polymer while also demonstrating the distinction in the mechanism of osmolyte action between protein and hydrophobic polymer. This marks, to our knowledge, the first experimental observation of TMAO-induced hydrophobic collapse in a ternary aqueous system.

Free energetics of carbon nanotube association in aqueous inorganic NaI salt solutions: Temperature effects using all-atom molecular dynamics simulations

In this study, we examine the temperature dependence of free energetics of nanotube association using graphical processing unit-enabled all-atom molecular dynamics simulations (FEN ZI) with two (10,10) single-walled carbon nanotubes in 3 m NaI aqueous salt solution. Results suggest that the free energy, enthalpy and entropy changes for the association process are all reduced at the high temperature, in agreement with previous investigations using other hydrophobes. Via the decomposition of free energy into individual components, we found that solvent contribution (including water, anion, and cation contributions) is correlated with the spatial distribution of the corresponding species and is influenced distinctly by the temperature. We studied the spatial distribution and the structure of the solvent in different regions: intertube, intratube and the bulk solvent. By calculating the fluctuation of coarse-grained tube-solvent surfaces, we found that tube-water interfacial fluctuation exhibits the strongest temperature dependence. By taking ions to be a solvent-like medium in the absence of water, tube-anion interfacial fluctuation shows similar but weaker dependence on temperature, while tube-cation interfacial fluctuation shows no dependence in general. These characteristics are discussed via the malleability of their corresponding solvation shells relative to the nanotube surface. Hydrogen bonding profiles and tetrahedrality of water arrangement are also computed to compare the structure of solvent in the solvent bulk and intertube region. The hydrophobic confinement induces a relatively lower concentration environment in the intertube region, therefore causing different intertube solvent structures which depend on the tube separation. This study is relevant in the continuing discourse on hydrophobic interactions (as they impact generally a broad class of phenomena in biology, biochemistry, and materials science and soft condensed matter research), and interpretations of hydrophobicity in terms of alternative but parallel signatures such as interfacial fluctuations, dewetting transitions, and enhanced fluctuation probabilities at interfaces.

Hydration of guanidinium depends on its local environment

Hydration of gaseous guanidinium (Gdm+) with up to 100 water molecules attached was investigated using infrared photodissociation spectroscopy in the hydrogen stretch region between 2900 and 3800 cm-1. Comparisons to IR spectra of low-energy computed structures indicate that at small cluster size, water interacts strongly with Gdm+ with three inner shell water molecules each accepting two hydrogen bonds from adjacent NH2 groups in Gdm+. Comparisons to results for tetramethylammonium (TMA+) and Na+ enable structural information for larger clusters to be obtained. The similarity in the bonded OH region for Gdm(H2O)20+vs. Gdm(H2O)100+ and the similarity in the bonded OH regions between Gdm+ and TMA+ but not Na+ for clusters with <50 water molecules indicate that Gdm+ does not significantly affect the hydrogen-bonding network of water molecules at large size. These results indicate that the hydration around Gdm+ changes for clusters with more than about eight water molecules to one in which inner shell water molecules only accept a single H-bond from Gdm+. More effective H-bonding drives this change in inner-shell water molecule binding to other water molecules. These results show that hydration of Gdm+ depends on its local environment, and that Gdm+ will interact with water even more strongly in an environment where water is partially excluded, such as the surface of a protein. This enhanced hydration in a limited solvation environment may provide new insights into the effectiveness of Gdm+ as a protein denaturant.

Role of arginine in mediating protein-carbon nanotube interactions

Arginine-rich proteins (e.g., lysozyme) or poly-L-arginine peptides have been suggested as solvating and dispersing agents for single-wall carbon nanotubes (CNTs) in water. In addition, protein structure-function in porous and hydrophobic materials is of broad interest. The amino acid residue, arginine (Arg(+)), has been implicated as an important mediator of protein/peptide-CNT interactions. To understand the structural and thermodynamic aspects of this interaction at the molecular level, we employ molecular dynamics (MD) simulations of the protein lysozyme in the interior of a CNT, as well as of free solutions of Arg(+) in the presence of a CNT. To dissect the Arg(+)-CNT interaction further, we also perform simulations of aqueous solutions of the guanidinium ion (Gdm(+)) and the norvaline (Nva) residue in the presence of a CNT. We show that the interactions of lysozyme with the CNT are mediated by the surface Arg(+) residues. The strong interaction of Arg(+) residue with the CNT is primarily driven by the favorable interactions of the Gdm(+) group with the CNT wall. The Gdm(+) group is not as well-hydrated on its flat sides, which binds to the CNT wall. This is consistent with a similar binding of Gdm(+) ions to a hydrophobic polymer. In contrast, the Nva residue, which lacks the Gdm(+) group, binds to the CNT weakly. We present details of the free energy of binding, molecular structure, and dynamics of these solutes on the CNT surface. Our results highlight the important role of Arg(+) residues in protein-CNT or protein-carbon-based material interactions. Such interactions could be manipulated precisely through protein engineering, thereby offering control over protein orientation and structure on CNTs, graphene, or other hydrophobic interfaces.

Protein denaturants at aqueous-hydrophobic interfaces: self-consistent correlation between induced interfacial fluctuations and denaturant stability at the interface

The notion of direct interaction between denaturing cosolvent and protein residues has been proposed in dialogue relevant to molecular mechanisms of protein denaturation. Here we consider the correlation between free energetic stability and induced fluctuations of an aqueous-hydrophobic interface between a model hydrophobically associating protein, HFBII, and two common protein denaturants, guanidinium cation (Gdm(+)) and urea. We compute potentials of mean force along an order parameter that brings the solute molecule close to the known hydrophobic region of the protein. We assess potentials of mean force for different relative orientations between the protein and denaturant molecule. We find that in both cases of guanidinium cation and urea relative orientations of the denaturant molecule that are parallel to the local protein-water interface exhibit greater stability compared to edge-on or perpendicular orientations. This behavior has been observed for guanidinium/methylguanidinium cations at the liquid-vapor interface of water, and thus the present results further corroborate earlier findings. Further analysis of the induced fluctuations of the aqueous-hydrophobic interface upon approach of the denaturant molecule indicates that the parallel orientation, displaying a greater stability at the interface, also induces larger fluctuations of the interface compared to the perpendicular orientations. The correlation of interfacial stability and induced interface fluctuation is a recurring theme for interface-stable solutes at hydrophobic interfaces. Moreover, observed correlations between interface stability and induced fluctuations recapitulate connections to local hydration structure and patterns around solutes as evidenced by experiment (Cooper et al., J. Phys. Chem. A 2014, 118, 5657.) and high-level ab initio/DFT calculations (Baer et al., Faraday Discuss 2013, 160, 89).

Ionic strength-modulated catalytic efficiency of a multienzyme cascade nanoconfined on charged hierarchical scaffolds

Understanding the effect of ionic strength on the efficiency of this enzyme cascade within charged hierarchical nanospace is not only fundamentally interesting, but also important for translating biochemical pathways to noncellular environments.

Computational study of elements of stability of a four-helix bundle protein biosurfactant

Biosurfactants are surface-active molecules produced principally by microorganisms. They are a sustainable alternative to chemically-synthesized surfactants, having the advantages of being non-toxic, highly functional, eco-friendly and biodegradable. However they are currently only used in a few industrial products due to costs associated with production and purification, which exceed those for commodity chemical surfactants. DAMP4, a member of a four-helix bundle biosurfactant protein family, can be produced in soluble form and at high yield in Escherichia coli, and can be recovered using a facile thermal phase-separation approach. As such, it encompasses an interesting synergy of biomolecular and chemical engineering with prospects for low-cost production even for industrial sectors. DAMP4 is highly functional, and due to its extraordinary thermal stability it can be purified in a simple two-step process, in which the combination of high temperature and salt leads to denaturation of all contaminants, whereas DAMP4 stays stable in solution and can be recovered by filtration. This study aimed to characterize and understand the fundamental drivers of DAMP4 stability to guide further process and surfactant design studies. The complementary use of experiments and molecular dynamics simulation revealed a broad pH and temperature tolerance for DAMP4, with a melting point of 122.4 °C, suggesting the hydrophobic core as the major contributor to thermal stability. Simulation of systematically created in silico variants of DAMP4 showed an influence of number and location of hydrophilic mutations in the hydrophobic core on stability, demonstrating a tolerance of up to three mutations before a strong loss in stability occurred. The results suggest a consideration of a balance of stability, functionality and kinetics for new designs according to their application, aiming for maximal functionality but at adequate stability to allow for cost-efficient production using thermal phase separation approaches.

Effects of the protein denaturant guanidinium chloride on aqueous hydrophobic contact-pair interactions

Guanidinium chloride (GdmCl) is one of the most common protein denaturants. Although GdmCl is well known in the field of protein folding, the mechanism by which it denatures proteins is not well understood. In fact, there are few studies looking at its effects on hydrophobic interactions. In this work the effect of GdmCl on hydrophobic interactions has been studied by observing how the denaturant influences model systems of phenyl and alkyl hydrophobic contact pairs. Contact pair formation is monitored through the use of fluorescence spectroscopy, i.e., measuring the intrinsic phenol fluorescence being quenched by carboxylate ions. Hydrophobic interactions are isolated from other interactions through a previously developed methodology. The results show that GdmCl does not significantly affect hydrophobic interactions between small moieties such as methyl groups and phenol; while on the other hand, the interaction of larger hydrophobes such as hexyl and heptyl groups with phenol is significantly destabilized.

Kinetic evidence for a two-stage mechanism of protein denaturation by guanidinium chloride

Dry molten globular (DMG) intermediates, an expanded form of the native protein with a dry core, have been observed during denaturant-induced unfolding of many proteins. These observations are counterintuitive because traditional models of chemical denaturation rely on changes in solvent-accessible surface area, and there is no notable change in solvent-accessible surface area during the formation of the DMG. Here we show, using multisite fluorescence resonance energy transfer, far-UV CD, and kinetic thiol-labeling experiments, that the guanidinium chloride (GdmCl)-induced unfolding of RNase H also begins with the formation of the DMG. Population of the DMG occurs within the 5-ms dead time of our measurements. We observe that the size and/or population of the DMG is linearly dependent on [GdmCl], although not as strongly as the second and major step of unfolding, which is accompanied by core solvation and global unfolding. This rapid GdmCl-dependent population of the DMG indicates that GdmCl can interact with the protein before disrupting the hydrophobic core. These results imply that the effect of chemical denaturants cannot be interpreted solely as a disruption of the hydrophobic effect and strongly support recent computational studies, which hypothesize that chemical denaturants first interact directly with the protein surface before completely unfolding the protein in the second step (direct interaction mechanism).

Concentration-dependent like-charge pairing of guanidinium ions and effect of guanidinium chloride on the structure and dynamics of water from all-atom molecular dynamics simulation

An all-atom molecular dynamics simulation shows concentration-dependent like-charge ion pairing of the guanidinium ion in an aqueous solution of guanidinium chloride. We have observed two types of like-charge ion pairing for guanidinium ions, namely, stacked ion pairs and solvent-separated ion pairs. Interestingly, both of these like-charge ion-pair formations are dependent on the concentration of guanidinium chloride in water. The probability of stacked like-charge ion-pair formation decreases, whereas, the probability of solvent-separated like-charge pairing increases as the concentration of guanidinium chloride increases, which is shown from radial distribution functions and is confirmed from the energy calculations. Besides like-charge ion-pair formation, we also investigated guanidinium chloride induced changes in water structure. Hydrogen-bond analysis indicates that guanidinium chloride does not alter the strict-hydrogen-bonding patterns of water, whereas, it breaks the bend-hydrogen bond and the non-hydrogen-bonding patterns. Tetrahedral order, nearest neighbor orientation, and distance distribution of water molecules around a probe water molecule show the extent of water structure distortion.

Molecular Insight into Different Denaturing Efficiency of Urea, Guanidinium, and Methanol: A Comparative Simulation Study

We have designed various nanoslit systems, whose opposing surfaces can be either hydrophobic, hydrophilic, or simply a water-vapor interface, for the molecular dynamics simulation of confined water with three different protein denaturants, i.e., urea, guanidinium chloride (GdmCl), and methanol, respectively. Particular attention is placed on the preferential adsorption of the denaturant molecules onto the opposing surfaces and associated resident time in the vicinal layer next to the surfaces, as well as their implication in the denaturing efficiency of different denaturant molecules. Our simulation results show that among the three denaturants, the occupancy of methanol in the vicinal layer is the highest while the residence time of Gdm is the longest. Although the occupancy and the residence time of urea in the vicinal layer is less than those of the other two denaturant molecules, urea entails "all-around" properties for being a highly effective denaturant. The distinct characteristics of three denaturants may suggest a different molecular mechanism for the protein denaturation. This comparative simulation by design allows us to gain additional insights, on the molecular level, into the denaturation effect and related hydrophobic effect.

Collapse of a hydrophobic polymer in a mixture of denaturants

The solvent quality of an aqueous mixture of two good solvents, urea and guanidinium chloride (GdmCl), for a hydrophobic polymer was investigated using atomistic molecular dynamics simulations. A counterintuitive collapse of the polymer was found, suggesting that mixing the two denaturants reduces the solvent quality. This cononsolvency of the polymer in the urea + GdmCl mixture is found to be caused by the preferential adsorption of urea on the polymer. The polymer collapses as a result of indirect long-range interactions between monomers resulting from the presence of urea clouds surrounding them. Surprisingly, urea behaves as the better solvent in the mixture not because there exists a stronger affinity of the polymer for urea. Instead, attractive interactions between two unlike denaturant molecules combined with the direct dispersion interactions of the polymer with both denaturants determine the solvent quality of the mixture.

Collapse of unfolded proteins in a mixture of denaturants

Both urea and guanidinium chloride (GdmCl) are frequently used as protein denaturants. Given that proteins generally adopt extended or unfolded conformations in either aqueous urea or GdmCl, one might expect that the unfolded protein chains will remain or become further extended due to the addition of another denaturant. However, a collapse of denatured proteins is revealed using atomistic molecular dynamics simulations when a mixture of denaturants is used. Both hen egg-white lysozyme and protein L are found to undergo collapse in the denaturant mixture. The collapse of the protein conformational ensembles is accompanied by a decreased solubility and increased non-native self-interactions of hydrophobic residues in the urea/GdmCl mixture. The increase of non-native interactions rather than the native contacts indicates that the proteins experience a simple collapse transition from the fully denatured states. During the protein collapse, the relatively stronger denaturant GdmCl displays a higher tendency to be absorbed onto the protein surface due to their stronger electrostatic interactions with proteins. At the same time, urea molecules also accumulate near the protein surface, resulting in an enhanced "local crowding" for the protein near its first solvation shell. This rearrangement of denaturants near the protein surface and crowded local environment induce the protein collapse, mainly by burying their hydrophobic residues. These findings from molecular simulations are then further explained by a simple analytical model based on statistical mechanics.

Free energetics of carbon nanotube association in pure and aqueous ionic solutions

Carbon nanotubes are a promising platform across a broad spectrum of applications ranging from separations technology, drug delivery, to bio(electronic) sensors. Proper dispersion of carbon nanotube materials is important to retaining the electronic properties of nanotubes. Experimentally it has been shown that salts can regulate the dispersing properties of CNTs in aqueous system with surfactants (Niyogi, S.; Densmore, C. G.; Doorn, S. K. J. Am. Chem. Soc.2009, 131, 1144-1153); details of the physicochemical mechanisms underlying such effects continue to be explored. We address the effects of inorganic monovalent salts (NaCl and NaI) on dispersion stability of carbon nanotubes.We perform all-atom molecular dynamics simulations using nonpolarizable interaction models to compute the potential of mean force between two (10,10) single-walled carbon nanotubes (SWNTs) in the presence of NaCl/NaI and compare to the potential of mean force between SWNTs in pure water. Addition of salts enhances stability of the contact state between two SWNT's on the order of 4 kcal/mol. The ion-specific spatial distribution of different halide anions gives rise to starkly different contributions to the free energy stability of nanotubes in the contact state. Iodide anion directly stabilizes the contact state to a much greater extent than chloride anion. The enhanced stability arises from the locally repulsive forces imposed on nanotubes by the surface-segregated iodide anion. Within the time scale of our simulations, both NaI and NaCl solutions stabilize the contact state by equivalent amounts. The marginally higher stability for contact state in salt solutions recapitulates results for small hydrophobic solutes in NaCl solutions (Athawale, M. V.; Sarupria, S.; Garde, S. J. Phys. Chem. B2008, 112, 5661-5670) as well as single-walled carbon nanotubes in NaCl and CaCl2 aqueous solutions.

Hydrophobicity of proteins and interfaces: insights from density fluctuations

Macroscopic characterizations of hydrophobicity (e.g., contact angle measurements) do not extend to the surfaces of proteins and nanoparticles. Molecular measures of hydrophobicity of such surfaces need to account for the behavior of hydration water. Theory and state-of-the-art simulations suggest that water density fluctuations provide such a measure; fluctuations are enhanced near hydrophobic surfaces and quenched with increasing surface hydrophilicity. Fluctuations affect conformational equilibria and dynamics of molecules at interfaces. Enhanced fluctuations are reflected in enhanced cavity formation, more favorable binding of hydrophobic solutes, increased compressibility of hydration water, and enhanced water-water correlations at hydrophobic surfaces. These density fluctuation-based measures can be used to develop practical methods to map the hydrophobicity/philicity of heterogeneous surfaces including those of proteins. They highlight that the hydrophobicity of a group is context dependent and is significantly affected by its environment (e.g., chemistry and topography) and especially by confinement. The ability to include information about hydration water in mapping hydrophobicity is expected to significantly impact our understanding of protein-protein interactions as well as improve drug design and discovery methods and bioseparation processes.

Effect of cosolvents on nano-confined water: a molecular dynamics study

We present results from atomistic molecular dynamics simulations to characterize the effects of cosolvents, such as urea and guanidinium (Gdm) salts, on the water confined in hydrophobic carbon nanotubes. We observed complete drying of the nanotube interiors of diameter ranging from 8 to 17 Å in urea. In contrast, the water population within nanotube cores smaller than 12 Å remains unaffected in GdmCl solution, whereas larger nanotube interiors become partially dehydrated with prevailing presence of stable Gdm(+)-Gdm(+) dimers. The molecular arrangement and the lifetime inside the nanotube were found to be characteristics of a particular cosolvent. In both urea and GdmCl solutions, preferential cosolvent intrusion resulting in nanotube dehydration is driven by the stronger dispersion interaction of cosolvent than water with the nanotube. The partial drying of the hydrophobic core is attributed to guanidinium's better hydration and weaker self-association propensity compared to urea, as well as to its moderate ion-pairing with strongly hydrated chloride ions. The Gdm(+) induced dehydration varies with the charge density of counter-ions, as the presence of high charge-density sulfate ions impedes penetration of guanidinium, and consequent dehydration of the nanotube. These findings provide important insights into the effect of cosolvents on the nano-confined water in a hydrophobic environment.

Effects of Hofmeister ions on the α-helical structure of proteins

The molecular conformation of proteins is sensitive to the nature of the aqueous environment. In particular, the presence of ions can stabilize or destabilize (denature) protein secondary structure. The underlying mechanisms of these actions are still not fully understood. Here, we combine circular dichroism (CD), single-molecule Förster resonance energy transfer, and atomistic computer simulations to elucidate salt-specific effects on the structure of three peptides with large α-helical propensity. CD indicates a complex ion-specific destabilization of the α-helix that can be rationalized by using a single salt-free computer simulation in combination with the recently introduced scheme of ion-partitioning between nonpolar and polar peptide surfaces. Simulations including salt provide a molecular underpinning of this partitioning concept. Furthermore, our single-molecule Förster resonance energy transfer measurements reveal highly compressed peptide conformations in molar concentrations of NaClO(4) in contrast to strong swelling in the presence of GdmCl. The compacted states observed in the presence of NaClO(4) originate from a tight ion-backbone network that leads to a highly heterogeneous secondary structure distribution and an overall lower α-helical content that would be estimated from CD. Thus, NaClO(4) denatures by inducing a molten globule-like structure that seems completely off-pathway between a fully folded helix and a coil state.