Comparison between Free and Immobilized Ion Effects on Hydrophobic Interactions: A Molecular Dynamics Study

Influence of Temperature and Salt Concentration on the Hydrophobic Interactions of Adamantane and Hexane

One of the definitions of hydrophobic interactions is the aggregation of nonpolar particles in a polar solvent, such as water. While this phenomenon appears to be very simple, it is crucial for many complex processes, such as protein folding, to take place. In this work, the hydrophobic association of adamantane and hexane at various temperatures and ionic strengths was studied using molecular dynamics simulations with the AMBER 16.0 program and the GAFF force field. The potentials of mean force of hydrophobic dimer formation, as well as the excess free energy, excess energy, excess entropy, and excess heat capacity corresponding to the formation of the contact minimum, were determined and analyzed. For both systems, the depth of the contact minimum in the potential of mean force was found to increase with both temperature and ionic strength. The excess heat capacity of the association at the contact minimum and T = 298 K was found to be negative and to decrease, while the excess entropy and energy were found to be positive and to increase for both systems, the changes being more pronounced for the hexane dimer. The excess heat capacity is also greater in absolute value for the hexane dimer.

Conservation and Identity Selection of Cationic Residues Flanking the Hydrophobic Regions in Intermediate Filament Superfamily

The interplay between the hydrophobic interactions generated by the nonpolar region and the proximal functional groups within nanometers of the nonpolar region offers a promising strategy to manipulate the intermolecular hydrophobic attractions in an artificial molecule system, but the outcomes of such modulations in the building of a native protein architecture remain unclear. Here we focus on the intermediate filament (IF) coiled-coil superfamily to assess the conservation of positively charged residue identity via a biostatistical approach. By screening the disease-correlated mutations throughout the IF superfamily, 10 distinct hotspots where a cation-to-cation substitution is associated with a pathogenic syndrome have been identified. The analysis of the local chemical context surrounding the hotspots revealed that the cationic diversity depends on their separation distance to the hydrophobic domain. The nearby cationic residues flanking the hydrophobic domain of a helix (separation <1 nm) are relatively conserved in evolution. In contrast, the cationic residues that are not adjacent to the hydrophobic domain (separation >1 nm) tolerate higher levels of variation and replaceability. We attribute this bias in the conservation degree of the cationic residue identity to reflect the interplay between the proximal cations and the hydrophobic interactions.

Opposite Regulatory Effects of Immobilized Cations on the Folding Vs. Assembly of Melittin

Ions are crucial in modulating the protein structure. For the free ions in bulk solution, ammonium is kosmotropic (structure forming) and guanidinium is chaotropic (structure breaking) to the protein structure within the Hofmeister series. However, the effect of immobilized ions on a protein surface is less explored. Herein, we explored the influence of two immobilized cations (ammonium in the side chain of lysine and guanidinium in the side chain of arginine) on the folding and assembly of melittin. Melittin adopts an α-helix structure and is driven by hydrophobic interactions to associate into a helical bundle. To test the influence of immobilized cations on the peptide structure, we designed the homozygous mutants exclusively containing ammonium (melittin-K) or guanidinium (melittin-R) and compared the differences of melittin-K vs. melittin-R in their folding, assembly, and molecular functions. The side chains of lysine and arginine differ in their influences on the folding and assembly of melittin. Specifically, the side chain of R increases the α-helical propensity of melittin relative to that of K, following an inverse Hofmeister series. In contrast, the side chain of K favors the assembly of melittin relative to the side chain of R in line with a direct Hofmeister series. The opposite regulatory effects of immobilized cations on the folding and assembly of melittin highlight the complexity of the noncovalent interactions that govern protein intermolecular architecture.

Cationic Side Chain Identity Directs the Hydrophobically Driven Self-Assembly of Amphiphilic β-Peptides in Aqueous Solution

Hydrophobic interactions mediated by nonpolar molecular fragments in water are influenced by local chemical and physical contexts in ways that are not yet fully understood. Here, we use globally amphiphilic (GA) β-peptides (GA-Lys and GA-Arg) with stable conformations to explore if replacement of β³-homolysine (βLys) with β³-homoarginine (βArg) influences the hydrophobically driven assembly of these peptides in bulk aqueous solution. The studies were conducted in 10 mM triethanolamine buffer at pH 7, where both βLys (ammonium) and βArg (guanidinium) side chains are substantially protonated. Comparisons of light scattering measurements and cryo-electron micrographs before and after the addition of 60 vol% MeOH indicate very different outcomes of the hydrophobically driven assembly of AcY-GA-Lys versus AcY-GA-Arg (AcY denotes an N-acetylated-β³-homotyrosine (βTyr) at each N-terminus). Nuclear magnetic resonance and analytical ultracentrifugation confirm that AcY-GA-Lys assembles into large aggregates in aqueous buffer, whereas AcY-GA-Arg at comparable concentrations forms only small oligomers. Titration of AcY-GA-Arg into aqueous solutions of AcY-GA-Lys reveals that the driving force for AcY-GA-Lys association is far stronger than that for AcY-GA-Arg association. We discuss these results in the light of past experimental observations involving single molecule force measurements with GA β-peptides and hydrophobically driven dimerization of conventional peptides that form a GA α-helix upon dimerization (but do not display the Lys versus Arg trend predicted by extrapolating from the earlier AFM studies with β-peptides). Overall, our results establish that the identity of proximal cationic groups, ammonium versus guanidinium, profoundly modulates the hydrophobically driven self-assembly of conformationally stable β-peptides in bulk aqueous solution.

A minireview on the perturbation effects of polar groups to direct nanoscale hydrophobic interaction and amphiphilic peptide assembly

Hydrophobic interaction provides the essential driving force for creating diverse native and artificial supramolecular architectures. Accumulating evidence leads to a hypothesis that the hydrophobicity of a nonpolar patch of a molecule is non-additive and susceptible to the chemical context of a judicious polar patch. However, the quantification of the hydrophobic interaction at the nanoscale remains a central challenge to validate the hypothesis. In this review, we aim to outline the recent efforts made to determine the hydrophobic interaction at a nanoscopic length scale. The advances achieved in the understanding of proximal polar groups perturbing the magnitude of hydrophobic interaction generated by the nonpolar patch are introduced. We will also discuss the influence of chemical heterogeneity on the modulation of amphiphilic peptide/protein assembly and molecular recognition.

Hydrophobic interaction provides the essential driving force for creating diverse native and artificial supramolecular architectures.

Aligned Assembly in a 2-D Gel of a Water-Soluble Peptide

We demonstrate the assembly of a compact, gel-like Langmuir-Blodgett film of rods formed by self-assembly of a β-sheet-forming water-soluble peptide, Ac-IKHLSVN-NH₂, at the surface of aqueous electrolytes. We characterize surface pressure hysteresis and demonstrate shear stiffening of the surface caused by area cycling, which we interpret as due to rearrangement and alignment of the rods. We show strong effects of the electrolyte on the assembly of the elementary rods, which can be related to the Hofmeister series and interpreted by effects on the interaction energies mediated by ions and water. Formation of β-sheet structures and assembly of these into surface-segregated semicrystalline gels was strongly promoted by ammonium sulfate electrolyte. With ammonium sulfate electrolyte as subphase for Langmuir-Blodgett film deposition, shear stiffening by surface area cycling resulted in very compact films on transfer to a substrate.

Modeling the nucleoporins that form the hairy pores

Sitting on the nuclear envelope, nuclear pore complexes (NPCs) control the molecular transport between the nucleus and the cytoplasm. Without definite open or close states, the NPC uses a family of intrinsically disordered nucleoporins called FG-Nups to construct a selective permeability barrier whose functional structure is unclear. Experimental advances have offered high-resolution molecular knowledge of the NPC scaffold and docking of the unfolded FG-Nups, however, the ‘hairy’ barrier structure still appears as blurred lobes even under the state-of-the-art microscopy. Without accurate experimental visualization, the molecular mechanism for the NPC-mediated transport remains a matter of debate. Modeling provides an alternative way to resolve this long-standing mystery. Here, we briefly review different methods employed in modeling the FG-Nups, arranging from all-atom molecular dynamics to mean-field theories. We discuss the advantage and limit of each modeling technique, and summarize the theoretical insights that, despite certain controversy, deepened our understanding of the hairy pore.

Hidden Isolated OH at the Charged Hydrophobic Interface Revealed by Two-Dimensional Heterodyne-Detected VSFG Spectroscopy

Water around hydrophobic groups mediates hydrophobic interactions that play key roles in many chemical and biological processes. Thus, the molecular-level elucidation of the properties of water in the vicinity of hydrophobic groups is important. We report on the structure and dynamics of water at two oppositely charged hydrophobic ion/water interfaces, that is, the tetraphenylborate-ion (TPB^- )/water and tetraphenylarsonium-ion (TPA⁺ )/water interfaces, which are clarified by two-dimensional heterodyne-detected vibrational sum-frequency generation (2D HD-VSFG) spectroscopy. The obtained 2D HD-VSFG spectra of the anionic TPB^- interface reveal the existence of distinct π-hydrogen bonded OH groups in addition to the usual hydrogen-bonded OH groups, which are hidden in the steady-state spectrum. In contrast, 2D HD-VSFG spectra of the cationic TPA⁺ interface only show the presence of usual hydrogen-bonded OH groups. The present study demonstrates that the sign of the interfacial charge governs the structure and dynamics of water molecules that face the hydrophobic region.

Long-range Regulation of Partially Folded Amyloidogenic Peptides

Neurodegeneration involves abnormal aggregation of intrinsically disordered amyloidogenic peptides (IDPs), usually mediated by hydrophobic protein-protein interactions. There is mounting evidence that formation of α-helical intermediates is an early event during self-assembly of amyloid-β42 (Aβ42) and α-synuclein (αS) IDPs in Alzheimer’s and Parkinson’s disease pathogenesis, respectively. However, the driving force behind on-pathway molecular assembly of partially folded helical monomers into helical oligomers assembly remains unknown. Here, we employ extensive molecular dynamics simulations to sample the helical conformational sub-spaces of monomeric peptides of both Aβ42 and αS. Our computed free energies, population shifts, and dynamic cross-correlation network analyses reveal a common feature of long-range intra-peptide modulation of partial helical folds of the amyloidogenic central hydrophobic domains via concerted coupling with their charged terminal tails (N-terminus of Aβ42 and C-terminus of αS). The absence of such inter-domain fluctuations in both fully helical and completely unfolded (disordered) states suggests that long-range coupling regulates the dynamicity of partially folded helices, in both Aβ42 and αS peptides. The inter-domain coupling suggests a form of intra-molecular allosteric regulation of the aggregation trigger in partially folded helical monomers. This approach could be applied to study the broad range of amyloidogenic peptides, which could provide a new path to curbing pathogenic aggregation of partially folded conformers into oligomers, by inhibition of sites far from the hydrophobic core.

Probing Peptide Sequences on Their Ability to Generate Affinity Sites in Molecularly Imprinted Polymers

An array of 4000 defined and addressable tripeptides on a polymer-coated glass slide is used to synthesize molecularly imprinted polymer (MIP) nanoparticles. This work is undertaken to systematically probe the impact of the peptide sequence on the ability to generate affinity MIPs. The polymer affinity is assessed by measuring the fluorescence of bound MIP nanoparticles at each peptide spot on the surface after washing the array to remove any low-affinity polymer. The generic composition commonly used in the preparation of MIPs against proteins seems to be equally suitable for imprinting hydrophobic and hydrophilic tripeptides. The amino acids frequently contributing to the formation of high-affinity MIPs include T, F, D, N, Y, W, and P. The amino acids that rarely contribute to the formation of high-affinity interactions with MIPs are G, V, A, L, I, and M. These observations are confirmed by computational modeling. The basic technique proposed here may be applicable in optimizing polymer compositions for the production of high-affinity MIPs or, more specifically, for the selection of appropriate amino acid sequences when peptide epitopes are used instead of whole protein imprinting.

Nanocompartmentalization of the Nuclear Pore Lumen

The nuclear pore complex (NPC) employs the intrinsically disordered regions (IDRs) from a family of phenylalanine-glycine-rich nucleoporins (FG-Nups) to control nucleocytoplasmic transport. It has been a long-standing mystery how the IDR-mediated mass exchange can be rapid yet selective. Here, we use a computational microscope to show that nanocompartmentalization of IDR subdomains leads to a remarkably elaborate gating structure as programmed by the amino acid sequences. In particular, we reveal a heterogeneous permeability barrier that combines an inner ring barrier with two vestibular condensates. Throughout the NPC, we find a polarized electrostatic potential and a diffuse thermoreversible FG network featuring mosaic FG territories with low FG-FG pairing fraction. Our theoretical anatomy of the central transporter sheds light into the sequence-structure-function relationship of the FG-Nups and provides a picture of nucleocytoplasmic mass exchange that allows a reconciliation of transport efficiency and specificity.

Molecular Order Affects Interfacial Water Structure and Temperature-Dependent Hydrophobic Interactions between Nonpolar Self-Assembled Monolayers

Understanding how material properties affect hydrophobic interactions-the water-mediated interactions that drive the association of nonpolar materials-is vital to the design of materials in contact with water. Conventionally, the magnitude of the hydrophobic interactions between extended interfaces is attributed to interfacial chemical properties, such as the amount of nonpolar solvent-exposed surface area. However, recent experiments have demonstrated that the hydrophobic interactions between uniformly nonpolar self-assembled monolayers (SAMs) also depend on molecular-level SAM order. In this work, we use atomistic molecular dynamics simulations to investigate the relationship between SAM order, water structure, and hydrophobic interactions to explain these experimental observations. The SAM-SAM hydrophobic interactions calculated from the simulations increase in magnitude as SAM order increases, matching experimental observations. We explain this trend by showing that the molecular-level order of the SAM impacts the nanoscale structure of interfacial water molecules, leading to an increase in water structure near disordered SAMs. These findings are consistent with a decrease in the solvation entropy of disordered SAMs, which is confirmed by measuring the temperature dependence of hydrophobic interactions using both simulations and experiments. This study elucidates how hydrophobic interactions can be influenced by an interfacial physical property, which may guide the design of synthetic materials with fine-tuned interfacial hydrophobicity.