An interplay of excluded-volume and polymer-(co)solvent attractive interactions regulates polymer collapse in mixed solvents.

Cosolvent effects on the coil-globule transitions in aqueous polymer solutions are not well understood, especially in the case of amphiphilic cosolvents that preferentially adsorb on the polymer and lead to both polymer swelling and collapse. Although a predominant focus in the literature has been placed on the role of polymer-cosolvent attractive interactions, our recent work has shown that excluded-volume interactions (repulsive interactions) can drive both preferential adsorption of the cosolvent and polymer collapse via a surfactant-like mechanism. Here, we further study the role of polymer-(co)solvent attractive interactions in two kinds of polymer solutions, namely, good solvent (water)-good cosolvent (alcohol) (GSGC) and poor solvent-good cosolvent (PSGC) solutions, both of which exhibit preferential adsorption of the cosolvent and a non-monotonic change in the polymer radius of gyration with the addition of the cosolvent. Interestingly, at low concentrations, the polymer-(co)solvent energetic interactions oppose polymer collapse in the GSGC solutions and contrarily support polymer collapse in the PSGC solutions, indicating the importance of the underlying polymer chemistry. Even though the alcohol molecules are preferentially adsorbed on the polymer, the trends of the energetic interactions at low cosolvent concentrations are dominated by the polymer-water energetic interactions in both the cases. Therefore, polymer-(co)solvent energetic interactions can either reinforce or compensate the surfactant-like mechanism, and it is this interplay that drives coil-to-globule transitions in polymer solutions. These results have implications for rationalizing the cononsolvency transitions in real systems such as polyacrylamides in aqueous alcohol solutions where the understanding of microscopic driving forces is still debatable.

A cosolvent surfactant mechanism affects polymer collapse in miscible good solvents.

The coil-globule transition of aqueous polymers is of profound significance in understanding the structure and function of responsive soft matter. In particular, the remarkable effect of amphiphilic cosolvents (e.g., alcohols) that leads to both swelling and collapse of stimuli-responsive polymers has been hotly debated in the literature, often with contradictory mechanisms proposed. Using molecular dynamics simulations, we herein demonstrate that alcohols reduce the free energy cost of creating a repulsive polymer-solvent interface via a surfactant-like mechanism which surprisingly drives polymer collapse at low alcohol concentrations. This hitherto neglected role of interfacial solvation thermodynamics is common to all coil-globule transitions, and rationalizes the experimentally observed effects of higher alcohols and polymer molecular weight on the coil-to-globule transition of thermoresponsive polymers. Polymer-(co)solvent attractive interactions reinforce or compensate this mechanism and it is this interplay which drives polymer swelling or collapse.

Improved Temperature Behavior of PNIPAM in Water with a Modified OPLS Model.

We test the OPLS/AA force field for a single PNIPAM 40-mer in aqueous solution using replica exchange molecular dynamics simulations and find that the force field fails to reproduce the experimental temperature behavior. To resolve this issue, we apply a modification on the partial charges previously suggested to reproduce the liquid-liquid phase separation of NIPAM aqueous solutions. The modified force field features stronger amide-water electrostatic interactions than the original OPLS model, predicts a weaker water-mediated monomer-monomer attraction, and reproduces the experimental coil-globule collapse enthalpy of PNIPAM in water. We revisit the cononsolvency problem of PNIPAM in methanol/water mixtures with the modified model and show that the dependence of the coil-globule collapse enthalpy on methanol concentration follows the experimental trend of the lower critical solution temperature. The calculations with the modified force field confirm that polymer dehydration is the determining factor for chain collapse in the cononsolvency regime.

A multi-state coarse grained modeling approach for an intrinsically disordered peptide.

Many proteins display a marginally stable tertiary structure, which can be altered via external stimuli. Since a majority of coarse grained (CG) models are aimed at structure prediction, their success for an intrinsically disordered peptide's conformational space with marginal stability and sensitivity to external stimuli cannot be taken for granted. In this study, by using the LKα14 peptide as a test system, we demonstrate a bottom-up approach for constructing a multi-state CG model, which can capture the conformational behavior of this peptide in three distinct environments with a unique set of interaction parameters. LKα14 is disordered in dilute solutions; however, it strictly adopts the α-helix conformation upon aggregation or when in contact with a hydrophobic/hydrophilic interface. Our bottom-up approach combines a generic base model, that is unbiased for any particular secondary structure, with nonbonded interactions which represent hydrogen bonds, electrostatics, and hydrophobic forces. We demonstrate that by using carefully designed all atom potential of mean force calculations from all three states of interest, one can get a balanced representation of the nonbonded interactions. Our CG model behaves intrinsically disordered in bulk water, folds into an α-helix in the presence of an interface or a neighboring peptide, and is stable as a tetrameric unit, successfully reproducing the all atom molecular dynamics simulations and experimental results.

Computational Calorimetry of PNIPAM Cononsolvency in Water/Methanol Mixtures.

We revisit the mechanism for cononsolvency of PNIPAM in water/methanol mixtures. Using extensive molecular dynamics simulations, we calculate the calorimetric enthalpy of the PNIPAM collapse transition and observe a unique fingerprint of PNIPAM cononsolvency which is analyzed in terms of microscopic interactions. We find that polymer hydration is the determining factor for PNIPAM collapse in the cononsolvency regime. In particular, it is shown that methanol frustrates the ability of water to form hydrogen bonds with the amide proton and therefore causes polymer collapse.

Conformation and Aggregation of LKα14 Peptide in Bulk Water and at the Air/Water Interface.

Historically, the protein folding problem has mainly been associated with understanding the relationship between amino acid sequence and structure. However, it is known that both the conformation of individual molecules and their aggregation strongly depend on the environmental conditions. Here, we study the aggregation behavior of the model peptide LKα14 (with amino acid sequence LKKLLKLLKKLLKL) in bulk water and at the air/water interface. We start by a quantitative analysis of the conformational space of a single LKα14 in bulk water. Next, in order to analyze the aggregation tendency of LKα14, by using the umbrella sampling technique we calculate the potential of mean force for pulling a single peptide from an n-molecule aggregate. In agreement with the experimental results, our calculations yield the optimal aggregate size as four. This equilibrium state is achieved by two opposing forces: Coulomb repulsion between the lysine side chains and the reduction of solvent accessible hydrophobic surface area upon aggregation. At the vacuum/water interface, however, even dimers of LKα14 become marginally stable, and any larger aggregate falls apart instantaneously. Our results indicate that even though the interface is highly influential in stabilizing the α-helix conformation for a single molecule, it significantly reduces the attraction between two LKα14 peptides, along with their aggregation tendency.

Tipping the Scale from Disorder to Alpha-helix: Folding of Amphiphilic Peptides in the Presence of Macroscopic and Molecular Interfaces.

Secondary amphiphilicity is inherent to the secondary structural elements of proteins. By forming energetically favorable contacts with each other these amphiphilic building blocks give rise to the formation of a tertiary structure. Small proteins and peptides, on the other hand, are usually too short to form multiple structural elements and cannot stabilize them internally. Therefore, these molecules are often found to be structurally ambiguous up to the point of a large degree of intrinsic disorder in solution. Consequently, their conformational preference is particularly susceptible to environmental conditions such as pH, salts, or presence of interfaces. In this study we use molecular dynamics simulations to analyze the conformational behavior of two synthetic peptides, LKKLLKLLKKLLKL (LK) and EAALAEALAEALAE (EALA), with built-in secondary amphiphilicity upon forming an alpha-helix. We use these model peptides to systematically study their aggregation and the influence of macroscopic and molecular interfaces on their conformational preferences. We show that the peptides are neither random coils in bulk water nor fully formed alpha helices, but adopt multiple conformations and secondary structure elements with short lifetimes. These provide a basis for conformation-selection and population-shift upon environmental changes. Differences in these peptides' response to macroscopic and molecular interfaces (presented by an aggregation partner) can be linked to their inherent alpha-helical tendencies in bulk water. We find that the peptides' aggregation behavior is also strongly affected by presence or absence of an interface, and rather subtly depends on their surface charge and hydrophobicity.

A transferable coarse-grained model for diphenylalanine: how to represent an environment driven conformational transition.

One of the major challenges in the development of coarse grained (CG) simulation models that aim at biomolecular structure formation processes is the correct representation of an environment-driven conformational change, for example, a folding/unfolding event upon interaction with an interface or upon aggregation. In the present study, we investigate this transferability challenge for a CG model using the example of diphenylalanine. This dipeptide displays a transition from a trans-like to a cis-like conformation upon aggregation as well as upon transfer from bulk water to the cyclohexane/water interface. Here, we show that one can construct a single CG model that can reproduce both the bulk and interface conformational behavior and the segregation between hydrophobic/hydrophilic medium. While the general strategy to obtain nonbonded interactions in the present CG model is to reproduce solvation free energies of small molecules representing the CG beads in the respective solvents, the success of the model strongly depends on nontrivial decisions one has to make to capture the delicate balance between the bonded and nonbonded interactions. In particular, we found that the peptide's conformational behavior is qualitatively affected by the cyclohexane/water interaction potential, an interaction that does not directly involve the peptide at all but merely influences the properties of the hydrophobic/hydrophilic interface. Furthermore, we show that a small modification to improve the structural/conformational properties of the CG model could dramatically alter the thermodynamic properties.