Hydrophobicity of protein surfaces: Separating geometry from chemistry

Homologous polydopamine ameliorates haemolysis of melittin for enhancing its anticancer efficacy

Despite exhibiting potent anticancer activity, the strong hemolytic properties of melittin (MEL) significantly restrict its delivery efficiency and clinical applications. To address this issue, we have devised a strategy wherein homologous dopamine (DA), an essential component of bee venom, is harnessed as a vehicle for the synthesis of MEL-polydopamine (PDA) nanoparticles (MP NPs). The ingenious approach lies in the fact that MEL is a basic polypeptide, and the polymerization of DA is also conducted under alkaline conditions, indicating the distinctive advantages of PDA in MEL encapsulation. Furthermore, MP NPs are modified with folic acid to fabricate tumor-targeted nanomedicine (MPF NPs). MPF NPs can ameliorate the hemolysis of MEL in drug delivery and undergo degradation triggered by high levels of reactive oxygen species (ROS) within solid tumors, thereby facilitating MEL release and subsequent restoration of anticancer activity. After cellular uptake, MPF NPs induce cell apoptosis through the PI3K/Akt-mediated p53 signaling pathway. The tumor growth inhibitory rate of MPF NPs in FA receptor-positive 4T1 and CT26 xenograft mice reached 78.04% and 81.66%, which was significantly higher compared to that in FA receptor-negative HepG2 xenograft mice (45.79%). Homologous vehicles provide a new perspective for nanomedicine design.

Inverse design of pore wall chemistry and topology through active learning of surface group interactions

Design of next-generation membranes requires a nanoscopic understanding of the effect of biologically inspired heterogeneous surface chemistries and topologies (roughness) on local water and solute behavior. In particular, the rejection of small, neutral solutes, such as boric acid, poses a heretofore unsolved challenge. In prior work, a computational inverse design technique using an evolutionary optimization successfully uncovered new surface design strategies for optimized transport of water over solutes in smooth, model pores consisting of two surface chemistries. However, extending such an approach to more complex (and realistic) scenarios involving many surface chemistries as well as surface roughness is challenging due to the expanded design space. In this work, we develop a new approach that uses active learning to optimize in a reduced feature space of surface group interactions, finding parameters that lead to their assembly into ordered, optimal patterns. This approach rapidly identifies novel surface functionalizations that maximize the difference in water and boric acid transport through the nanopore. Moreover, we find that the roughness of the nanopore wall, independent of its chemistry, can be leveraged to enhance transport selectivity: oscillations in the pore wall diameter optimally inhibit boric acid transport by creating energetic wells from which the solute must escape to transport down the pore. This proof-of-concept demonstrates the potential for active learning strategies, in concert with molecular simulations, to rapidly navigate complex design spaces of aqueous interfaces and is promising as a tool for engineering water-mediated surface interactions for a broad range of applications.

An Affordable Topography-Based Protocol for Assigning a Residue's Character on a Hydropathy (PARCH) Scale

The hydropathy of proteins or quantitative assessment of protein-water interactions has been a topic of interest for decades. Most hydropathy scales use a residue-based or atom-based approach to assign fixed numerical values to the 20 amino acids and categorize them as hydrophilic, hydroneutral, or hydrophobic. These scales overlook the protein's nanoscale topography, such as bumps, crevices, cavities, clefts, pockets, and channels, in calculating the hydropathy of the residues. Some recent studies have included protein topography in determining hydrophobic patches on protein surfaces, but these methods do not provide a hydropathy scale. To overcome the limitations in the existing methods, we have developed a Protocol for Assigning a Residue's Character on the Hydropathy (PARCH) scale that adopts a holistic approach to assigning the hydropathy of a residue. The parch scale evaluates the collective response of the water molecules in the protein's first hydration shell to increasing temperatures. We performed the parch analysis of a set of well-studied proteins that include the following─enzymes, immune proteins, and integral membrane proteins, as well as fungal and virus capsid proteins. Since the parch scale evaluates every residue based on its location, a residue may have very different parch values inside a crevice versus a surface bump. Thus, a residue can have a range of parch values (or hydropathies) dictated by the local geometry. The parch scale calculations are computationally inexpensive and can compare hydropathies of different proteins. The parch analysis can affordably and reliably aid in designing nanostructured surfaces, identifying hydrophilic and hydrophobic patches, and drug discovery.

Relationships between Water's Structure and Solute Affinity at Polypeptoid Brush Surfaces

Excellent antifouling surfaces are generally thought to create a tightly bound layer of water that resists solute adsorption, and highly hydrophilic surfaces such as those with zwitterionic functionalities are of significant current interest as antifoulant strategies. However, despite significant proofs-of-concept, we still lack a fundamental understanding of how the nanoscopic structure of this hydration layer translates to reduced fouling, how surface chemistry can be tuned to achieve antifouling through hydration water, and why, in particular, zwitterionic surfaces seem so promising. Here, we use molecular dynamics simulations and free energy calculations to investigate the molecular relationships among surface chemistry, hydration water structure, and surface-solute affinity across a variety of surface-decorated chemistries. Specifically, we consider polypeptoid-decorated surfaces that display well-known experimental antifouling capabilities and that can be synthesized sequence specifically, with precise backbone positioning of, e.g., charged groups. Through simulations, we calculate the affinities of a range of small solutes to polypeptoid brush surfaces of varied side-chain chemistries. We then demonstrate that measures of the structure of surface hydration water in response to a particular surface chemistry signal solute-surface affinity; specifically, we find that zwitterionic chemistries produce solute-surface repulsion through highly coordinated hydration water while suppressing tetrahedral structuring around the solute, in contrast to uncharged surfaces that show solute-surface affinity. Based on the relationship of this structural perturbation to the affinity of small-molecule solutes, we propose a molecular mechanism by which zwitterionic surface chemistries enhance solute repulsion, with broader implications for the design of antifouling surfaces.

UNILIPID, a Methodology for Energetically Accurate Prediction of Protein Insertion into Implicit Membranes of Arbitrary Shape

The insertion of proteins into membranes is crucial for understanding their function in many biological processes. In this work, we present UNILIPID, a universal implicit lipid-protein description as a methodology for dealing with implicit membranes. UNILIPID is independent of the scale of representation and can be applied at the level of all atoms, coarse-grained particles down to the level of a single bead per amino acid. We provide example implementations for these scales and demonstrate the versatility of our approach by accurately reflecting the free energy of transfer for each amino acid. In addition to single membranes, we describe the analytical implementation of double membranes and show that UNILIPID is well suited for modeling at multiple scales. We generalize to membranes of arbitrary shape. With UNILIPID, we provide a methodological framework for a simple and general parameterization tuned to reproduce a selected reference hydrophobicity scale. The software we provide along with the methodological description is optimized for specific user features such as real-time response, live visual analysis, and virtual reality experiences.

Interactions of SARS-CoV-2 and MERS-CoV fusion peptides measured using single-molecule force methods

We address the challenge of understanding how hydrophobic interactions are encoded by fusion peptide (FP) sequences within coronavirus (CoV) spike proteins. Within the FPs of severe acute respiratory syndrome CoV 2 and Middle East respiratory syndrome CoV (MERS-CoV), a largely conserved peptide sequence called FP1 (SFIEDLLFNK and SAIEDLLFDK in SARS-2 and MERS, respectively) has been proposed to play a key role in encoding hydrophobic interactions that drive viral-host cell membrane fusion. Although a non-polar triad (Leu-Leu-Phe (LLF)) is common to both FP1 sequences, and thought to dominate the encoding of hydrophobic interactions, FP1 from SARS-2 and MERS differ in two residues (Phe 2 versus Ala 2 and Asn 9 versus Asp 9, respectively). Here we explore whether single-molecule force measurements can quantify hydrophobic interactions encoded by FP1 sequences, and then ask whether sequence variations between FP1 from SARS-2 and MERS lead to significant differences in hydrophobic interactions. We find that both SARS-2 and MERS wild-type FP1 generate measurable hydrophobic interactions at the single-molecule level, but that SARS-2 FP1 encodes a substantially stronger hydrophobic interaction than its MERS counterpart (1.91 ± 0.03 nN versus 0.68 ± 0.03 nN, respectively). By performing force measurements with FP1 sequences with single amino acid substitutions, we determine that a single-residue mutation (Phe 2 versus Ala 2) causes the almost threefold difference in the hydrophobic interaction strength generated by the FP1 of SARS-2 versus MERS, despite the presence of LLF in both sequences. Infrared spectroscopy and circular dichroism measurements support the proposal that the outsized influence of Phe 2 versus Ala 2 on the hydrophobic interaction arises from variation in the secondary structure adopted by FP1. Overall, these insights reveal how single-residue diversity in viral FPs, including FP1 of SARS-CoV-2 and MERS-CoV, can lead to substantial changes in intermolecular interactions proposed to play a key role in viral fusion, and hint at strategies for regulating hydrophobic interactions of peptides in a range of contexts.

Changes in cell surface properties of Pseudomonas fluorescens by adaptation to NaCl induced hypertonic stress

Determination of the effect of water stress on the surface properties of bacteria is crucial to study bacterial induced soil water repellency. Changes in the environmental conditions may affect several properties of bacteria such as the cell hydrophobicity and morphology. Here, we study the influence of adaptation to hypertonic stress on cell wettability, shape, adhesion, and surface chemical composition of Pseudomonas fluorescens. From this we aim to discover possible relations between the changes in wettability of bacterial films studied by contact angle and single cells studied by atomic and chemical force microscopy (AFM, CFM), which is still lacking. We show that by stress the adhesion forces of the cell surfaces towards hydrophobic functionalized probes increase while they decrease towards hydrophilic functionalized tips. This is consistent with the contact angle results. Further, cell size shrunk and protein content increased upon stress. The results suggest two possible mechanisms: Cell shrinkage is accompanied by the release of outer membrane vesicles by which the protein to lipid ratio increases. The higher protein content increases the rigidity and the number of hydrophobic nano-domains per surface area.

Learning the relationship between nanoscale chemical patterning and hydrophobicity

The hydrophobicity of proteins and similar surfaces, which display chemical heterogeneity at the nanoscale, drives countless aqueous interactions and assemblies. However, predicting how surface chemical patterning influences hydrophobicity remains a challenge. Here, we address this challenge by using molecular simulations and machine learning to characterize and model the hydrophobicity of a diverse library of patterned surfaces, spanning a wide range of sizes, shapes, and chemical compositions. We find that simple models, based only on polar content, are inaccurate, whereas complex neural network models are accurate but challenging to interpret. However, by systematically incorporating chemical correlations between surface groups into our models, we are able to construct a series of minimal models of hydrophobicity, which are both accurate and interpretable. Our models highlight that the number of proximal polar groups is a key determinant of hydrophobicity and that polar neighbors enhance hydrophobicity. Although our minimal models are trained on particular patch size and shape, their interpretability enables us to generalize them to rectangular patches of all shapes and sizes. We also demonstrate how our models can be used to predict hot-spot locations with the largest marginal contributions to hydrophobicity and to design chemical patterns that have a fixed polar content but vary widely in their hydrophobicity. Our data-driven models and the principles they furnish for modulating hydrophobicity could facilitate the design of novel materials and engineered proteins with stronger interactions or enhanced solubilities.

Engineered Design of the E-Helix Structure on Ferritin Nanoparticles

Insertion of an immunogenic epitope at the C-terminus of ferritin has shown the potential to produce a stable and efficacious vaccine. There is however limited understanding of how C-terminus insertion affects ferritin protein stability. The E-helix at the C-terminus has attracted interest because there are contradictory reports as to whether it has a role in protein stabilization. Here, we report, for the first time, combining molecular dynamics simulation (MDS) with experiment to engineer the design of the E-helix at the C-terminus of engineered human ferritin heavy chain (F₁) inserted with Epstein-Barr nuclear antigen 1 (EBNA1, E₁) and flexible linker (L₃) residues (to afford F₁L₃E₁). Hot spots on the E-helix of the C-terminus were predicted by MDS at aa 167 (Glu) and aa 171 (Asp). Five (5) variants of F₁L₃E₁ were constructed by considering hot spots and alteration of electrostatic or hydrophobic interfaces, namely, (1) C1, hot spots substituted with noncharged residue Gln; (2) C2, hot spots substituted with positively charged residue Arg; (3) C3, hydrophobic residues substituted with the most hydrophobic residues Val and Ile; (4) C4, hydrophobic residues substituted with the most hydrophilic residues Gln and Asn; and (5) C5, a heptad repeat structure in the E-helix disrupted by substituting "a" and "d" heptad residues with noncharged polar residue Gln. It was found that the E-helix is essential to maintain integrated protein stability and that changing the hydrophobic interface (C3 and C4) had more significant effects on protein folding and stability than changing the electrostatic interface (C1 and C2). It was confirmed by both MDS and experiment that variants C1, C2, and C5 were able to fold to form stable conformational structures with protein surface hydrophobicity similar to that of F₁L₃E₁. However, they are less thermally stable than F₁L₃E₁. Significant changes in hydrophobicity drove significant protein aggregation for variants C3 and C4. It is concluded that the molecular design of the C-terminus in engineered ferritin, especially the E-helix, is important to ensure the epitope-based chimeric vaccine is safe (aggregate free) and efficacious.

Sequence Modulates Polypeptoid Hydration Water Structure and Dynamics

We use molecular dynamics simulations to investigate the effect of polypeptoid sequence on the structure and dynamics of its hydration waters. Polypeptoids provide an excellent platform to study small-molecule hydration in disordered polymers, as they can be precisely synthesized with a variety of sidechain chemistries. We examine water behavior near a set of peptoid oligomers in which the number and placement of nonpolar versus polar sidechains are systematically varied. To do this, we leverage a new computational workflow enabling accurate sampling of polypeptoid conformations. We find that the hydration waters are less dense, are more tetrahedral, and have slower dynamics compared to bulk water. The magnitude of these shifts increases with the number of nonpolar groups. We also find that shifts in the water structure and dynamics are strongly correlated, suggesting that experimental insight into the dynamics of hydration water obtained by Overhauser dynamic nuclear polarization (ODNP) also contains information about water structural properties. We then demonstrate the ability of ODNP to probe site-specific dynamics of hydration water near these model peptoid systems.

Identifying hydrophobic protein patches to inform protein interaction interfaces

Interactions between proteins lie at the heart of numerous biological processes and are essential for the proper functioning of the cell. Although the importance of hydrophobic residues in driving protein interactions is universally accepted, a characterization of protein hydrophobicity, which informs its interactions, has remained elusive. The challenge lies in capturing the collective response of the protein hydration waters to the nanoscale chemical and topographical protein patterns, which determine protein hydrophobicity. To address this challenge, here, we employ specialized molecular simulations wherein water molecules are systematically displaced from the protein hydration shell; by identifying protein regions that relinquish their waters more readily than others, we are then able to uncover the most hydrophobic protein patches. Surprisingly, such patches contain a large fraction of polar/charged atoms and have chemical compositions that are similar to the more hydrophilic protein patches. Importantly, we also find a striking correspondence between the most hydrophobic protein patches and regions that mediate protein interactions. Our work thus establishes a computational framework for characterizing the emergent hydrophobicity of amphiphilic solutes, such as proteins, which display nanoscale heterogeneity, and for uncovering their interaction interfaces.

Behavior of Water Near Multimodal Chromatography Ligands and Its Consequences for Modulating Protein-Ligand Interactions

Multimodal chromatography is a powerful approach for purifying proteins that uses ligands containing multiple modes of interaction. Recent studies have shown that selectivity in multimodal chromatographic separations is a function of the ligand structure and geometry. Here, we performed molecular dynamics simulations to explore how the ligand structure and geometry affect ligand-water interactions and how these differences in solution affect the nature of protein-ligand interactions. Our investigation focused on three chromatography ligands: Capto MMC, Nuvia cPrime, and Prototype 4, a structural variant of Nuvia cPrime. First, the solvation characteristics of each ligand were quantified via three metrics: average water density, fluctuations, and residence time. We then explored how solvation was perturbed when the ligand was bound to the protein surface and found that the probability of the phenyl ring dewetting followed the order: Capto MMC > Prototype 4 > Nuvia cPrime. To explore how these differences in dewetting affect protein-ligand interactions, we calculated the probability of each ligand binding to different types of residues on the protein surface and found that the probability of binding to a hydrophobic residue followed the same order as the dewetting behavior. This study illustrates the role that wetting and dewetting play in modulating protein-ligand interactions.

Size dependence of hydrophobic hydration at electrified gold/water interfaces

Hydrophobic hydration at metal/water interfaces actively contributes to the energetics of electrochemical reactions, e.g. [Formula: see text] and [Formula: see text] reduction, where small hydrophobic molecules are involved. In this work, constant applied potential molecular dynamics is employed to study hydrophobic hydration at a gold/water interface. We propose an adaptation of the Lum-Chandler-Weeks (LCW) theory to describe the free energy of hydrophobic hydration at the interface as a function of solute size and applied voltage. Based on this model we are able to predict the free energy cost of cavity formation at the interface directly from the free energy cost in the bulk plus an interface-dependent correction term. The interfacial water network contributes significantly to the free energy, yielding a preference for outer-sphere adsorption at the gold surface for ideal hydrophobes. We predict an accumulation of small hydrophobic solutes of sizes comparable to CO or [Formula: see text], while the free energy cost to hydrate larger hydrophobes, above 2.5-Å radius, is shown to be greater at the interface than in the bulk. Interestingly, the transition from the volume dominated to the surface dominated regimes predicted by the LCW theory in the bulk is also found to take place for hydrophobes at the Au/water interface but occurs at smaller cavity radii. By applying the adapted LCW theory to a simple model addition reaction, we illustrate some implications of our findings for electrochemical reactions.

Structural aspects of an energy-based water classification index and the structure-dynamics link in glassy relaxation

An energy-based structural indicator for water, [Formula: see text], has been recently introduced by our group. In turn, in this work we aim at: (1) demonstrating that [Formula: see text] is indeed able to correctly classify water molecules between locally structured tetrahedral (T) and locally distorted (D) ones, circumventing the usual problem of certain previous indicators of overestimating the distorted state; (2) correlating [Formula: see text] with dynamic propensity, a measure of the molecular mobility tendency, in order to seek for the existence of a connection between structure and dynamics within the supercooled regime. More specifically, in the first part of this work we will show that [Formula: see text] accurately discriminates between merely thermally deformed local molecular arrangements and truly distorted molecules (defects). This fact will be made evident not only from radial distribution function results but also from the dynamic propensity distributions of the different kinds of molecules. In turn, we shall devote the second part of this work to finding correlations between T and D molecules with low- and high-dynamic-propensity molecules, respectively, thus revealing the existence of a link between local structure and dynamics, while also making evident the dominant role of the D molecules (defects) in the structural relaxation. Moreover, the availability of a proper molecular classification technique will enable us to study the timescale of such influence of structure on dynamics by defining a modified dynamic propensity measure and by applying it to the structured and unstructured water molecular states.

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.

Affinity of small-molecule solutes to hydrophobic, hydrophilic, and chemically patterned interfaces in aqueous solution

Performance of membranes for water purification is highly influenced by the interactions of solvated species with membrane surfaces, including surface adsorption of solutes upon fouling. Current efforts toward fouling-resistant membranes often pursue surface hydrophilization, frequently motivated by macroscopic measures of hydrophilicity, because hydrophobicity is thought to increase solute-surface affinity. While this heuristic has driven diverse membrane functionalization strategies, here we build on advances in the theory of hydrophobicity to critically examine the relevance of macroscopic characterizations of solute-surface affinity. Specifically, we use molecular simulations to quantify the affinities to model hydroxyl- and methyl-functionalized surfaces of small, chemically diverse, charge-neutral solutes represented in produced water. We show that surface affinities correlate poorly with two conventional measures of solute hydrophobicity, gas-phase water solubility and oil-water partitioning. Moreover, we find that all solutes show attraction to the hydrophobic surface and most to the hydrophilic one, in contrast to macroscopically based hydrophobicity heuristics. We explain these results by decomposing affinities into direct solute interaction energies (which dominate on hydroxyl surfaces) and water restructuring penalties (which dominate on methyl surfaces). Finally, we use an inverse design algorithm to show how heterogeneous surfaces, with multiple functional groups, can be patterned to manipulate solute affinity and selectivity. These findings, importantly based on a range of solute and surface chemistries, illustrate that conventional macroscopic hydrophobicity metrics can fail to predict solute-surface affinity, and that molecular-scale surface chemical patterning significantly influences affinity-suggesting design opportunities for water purification membranes and other engineered interfaces involving aqueous solute-surface interactions.

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.

Dynamic behaviors of interfacial water on the self-assembly monolayer (SAM) heterogeneous surface

Dynamic behaviors of water molecules near the surface with mixed hydrophobic and hydrophilic areas are studied by molecular dynamics simulation. More specifically, the diffusion coefficient and hydrogen bond lifetime of interfacial water on the self-assembly monolayer composed of hydrophobic and hydrophilic groups and their dependence on the mixing ratio are studied. The diffusion dramatically slows down, and the hydrogen bond lifetime considerably increases when a few hydrophilic groups are added to the hydrophobic surface. When the percentage of hydrophilic groups increases to 25%, the behavior of interfacial water is similar to the case of the pure hydrophilic surface. The sensitivity to the hydrophilic group can be attributed to the fact that the grafted hydrophilic groups can not only retard the directly bound water molecules but also affect indirectly bound water by stabilizing hydrogen bonds among interfacial water molecules.

Exploring the Free-Energy Landscape and Thermodynamics of Protein-Protein Association

We report the free-energy landscape and thermodynamics of the protein-protein association responsible for the drug-induced multimerization of HIV-1 integrase (IN). Allosteric HIV-1 integrase inhibitors promote aberrant IN multimerization by bridging IN-IN intermolecular interactions. However, the thermodynamic driving forces and kinetics of the multimerization remain largely unknown. Here, we explore the early steps in the IN multimerization by using umbrella sampling and unbiased molecular dynamics simulations in explicit solvent. In direct simulations, the two initially separated dimers spontaneously associate to form near-native complexes that resemble the crystal structure of the aberrant tetramer. Most strikingly, the effective interaction of the protein-protein association is very short-ranged: the two dimers associate rapidly within tens of nanoseconds when their binding surfaces are separated by d ≤ 4.3 Å (less than two water diameters). Beyond this distance, the oligomerization kinetics appears to be diffusion controlled with a much longer association time. The free-energy profile also captured the crucial role of allosteric IN inhibitors in promoting multimerization and explained why several C-terminal domain mutations are remarkably resistant to the drug-induced multimerization. The results also show that at small separation, the protein-protein binding process contains two consecutive phases with distinct thermodynamic signatures. First, interprotein water molecules are expelled to the bulk, resulting in a small increase in entropy, as the solvent entropy gain from the water release is nearly cancelled by the loss of side-chain entropies as the two proteins approach each other. At shorter distances, the two dry binding surfaces adapt to each other to optimize their interaction energy at the expense of further protein configurational entropy loss. Although the binding interfaces feature clusters of hydrophobic residues, overall, the protein-protein association in this system is driven by enthalpy and opposed by entropy.

Dynamical control by water at a molecular level in protein dimer association and dissociation

Water, often termed as the "lubricant of life," is expected to play an active role in navigating protein dissociation-association reactions. In order to unearth the molecular details, we first compute the free-energy surface (FES) of insulin dimer dissociation employing metadynamics simulation, and then carry out analyses of insulin dimerization and dissociation using atomistic molecular-dynamics simulation in explicit water. We select two sets of initial configurations from 1) the dissociated state and 2) the transition state, and follow time evolution using several long trajectories (∼1-2 μs). During the process we not only monitor configuration of protein monomers, but also the properties of water. Although the equilibrium structural properties of water between the two monomers approach bulklike characteristics at a separation distance of ∼5 nm, the dynamics differ considerably. The complex association process is observed to be accompanied by several structural and dynamical changes of the system, such as large-scale correlated water density fluctuations, coupled conformational fluctuation of protein monomers, a dewettinglike transition with the change of intermonomeric distance R_MM from ∼4 to ∼2 nm, orientation of monomers and hydrophobic hydration in the monomers. A quasistable, solvent-shared, protein monomer pair (SSPMP) forms at around 2 nm during association process which is a local free-energy minimum having ∼50-60% of native contacts. Simulations starting with arrangements sampled from the transition state (TS) of the dimer dissociation reveal that the final outcome depends on relative orientation of the backbone in the "hotspot" region.

Determination and evaluation of the nonadditivity in wetting of molecularly heterogeneous surfaces

Every folded protein presents an interface with water that is composed of domains of varying hydrophilicity/-phobicity. Many simulation studies have highlighted the nonadditivity in the wetting of such nanostructured surfaces in contrast with the accepted theoretical formula that is additive. We present here an experimental study on surfaces of identical composition but different organization of hydrophobic and hydrophilic domains. We prove that the interfacial energy of such surfaces differs by ∼20% and that a significant difference in the interfacial water H-bonding structure can be measured. As a result, in combination with molecular-dynamics simulations, we propose a model that captures the wetting of molecularly heterogeneous surfaces, showing the importance of local structure (first-nearest neighbors) in determining the wetting properties.

The interface between water and folded proteins is very complex. Proteins have “patchy” solvent-accessible areas composed of domains of varying hydrophobicity. The textbook understanding is that these domains contribute additively to interfacial properties (Cassie’s equation, CE). An ever-growing number of modeling papers question the validity of CE at molecular length scales, but there is no conclusive experiment to support this and no proposed new theoretical framework. Here, we study the wetting of model compounds with patchy surfaces differing solely in patchiness but not in composition. Were CE to be correct, these materials would have had the same solid–liquid work of adhesion (W_SL) and time-averaged structure of interfacial water. We find considerable differences in W_SL, and sum-frequency generation measurements of the interfacial water structure show distinctively different spectral features. Molecular-dynamics simulations of water on patchy surfaces capture the observed behaviors and point toward significant nonadditivity in water density and average orientation. They show that a description of the molecular arrangement on the surface is needed to predict its wetting properties. We propose a predictive model that considers, for every molecule, the contributions of its first-nearest neighbors as a descriptor to determine the wetting properties of the surface. The model is validated by measurements of W_SL in multiple solvents, where large differences are observed for solvents whose effective diameter is smaller than ∼6 Å. The experiments and theoretical model proposed here provide a starting point to develop a comprehensive understanding of complex biological interfaces as well as for the engineering of synthetic ones.

Recovering superhydrophobicity in nanoscale and macroscale surface textures

Here, we investigate the complete drying of hydrophobic cavities in order to elucidate the dependence of drying on the size, the geometry, and the degree of hydrophobicity of the confinement. Two complementary theoretical approaches are adopted: a macroscopic one based on classical capillarity and a microscopic classical density functional theory. This combination allows us to pinpoint unique drying mechanisms at the nanoscale and to clearly differentiate them from the mechanisms operational at the macroscale. Nanoscale hydrophobic cavities allow the thermodynamic destabilization of the confined liquid phase over an unexpectedly broad range of conditions, including pressures as large as 10 MPa and contact angles close to 90°. On the other hand, for cavities on the micron scale, such destabilization occurs only for much larger contact angles and close to liquid-vapor coexistence. These scale-dependent drying mechanisms are used to propose design criteria for hierarchical superhydrophobic surfaces capable of spontaneous self-recovery over a broad range of operating conditions. In particular, we detail the requirements under which it is possible to realize perpetual superhydrophobicity at positive pressures on surfaces with micron-sized textures by exploiting drying, facilitated by nanoscale coatings. Concerning the issue of superhydrophobicity, these findings indicate a promising direction both for surface fabrication and for the experimental characterization of perpetual surperhydrophobicity. From a more basic perspective, the present results have an echo on a wealth of biological problems in which hydrophobic confinement induces drying, such as in protein folding, molecular recognition, and hydrophobic gating.

Structural features of high-local-density water molecules: Insights from structure indicators based on the translational order between the first two molecular shells

The two-liquids scenario for liquid water assumes the existence of two competing preferential local molecular structural states characterized by either low or high local density. While the former is expected to present good local order thus involving privileged structures, the latter is usually regarded as conforming a high-entropy unstructured state. A main difference in the local arrangement of such "classes" of water molecules can be inferred from the degree of translational order between the first and second molecular shells. This is so, since the low-local-density molecules present a clear gap between the first two shells while in the case of the high-local-density ones, one or more molecules from the second shell have collapsed toward the first one, thus populating the intershell region. Some structural indicators, like the widely employed local structure index and the recently introduced ζ index, have been devised precisely on the basis of this observation, being successful in detecting well-structured low-local-density molecules. However, the nature of the high-local-density state has been mainly disregarded over the years. In this work we employ molecular dynamics simulations for two water models (the extended simple point charge model and the five-site model) at the liquid and supercooled regimes combined with the inherent dynamics approach (energy minimizations of the instantaneous configurations) in order to both rationalize the detailed structural and topological information that these indicators provide and to advance in our understanding of the high-density state.

Spatially Heterogeneous Water Properties at Disordered Surfaces Decrease the Hydrophobicity of Nonpolar Self-Assembled Monolayers

Understanding the relationship between hydrophobicity and the properties of functionalized surfaces is vital to the design of materials that interact in aqueous environments. In this Letter, we use atomistic molecular dynamics simulations to investigate the effects of surface order on the hydrophobicity of self-assembled monolayers (SAMs) containing nonpolar ligands. We find that the interfacial hydrophobicity is highly correlated with SAM order and, strikingly, poorly correlated with the solvent-accessible surface area, which typically has been related to interfacial hydrophobicity. Analysis of spatial variations in both SAM and water properties reveals that the SAM-water interface is pinned near regions of disordered SAM surfaces with increased free volume, decreasing the overall interfacial hydrophobicity. Spatial variations in ligand end group positions at disordered SAM surfaces thus translate to spatial variations in hydrophobicity, yielding heterogeneous surface properties. These findings provide new insights into how surface order can alter the hydrophobicity of chemically uniform surfaces.

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.

Characterizing surface wetting and interfacial properties using enhanced sampling (SWIPES)

We introduce an accurate and efficient method for characterizing surface wetting and interfacial properties, such as the contact angle made by a liquid droplet on a solid surface, and the vapor-liquid surface tension of a fluid. The method makes use of molecular simulations in conjunction with the indirect umbrella sampling technique to systematically wet the surface and estimate the corresponding free energy. To illustrate the method, we study the wetting of a family of Lennard-Jones surfaces by water. For surfaces with a wide range of attractions for water, we estimate contact angles using our method, and compare them with contact angles obtained using droplet shapes. Notably, our method is able to capture the transition from partial to complete wetting as surface-water attractions are increased. Moreover, the method is straightforward to implement and is computationally efficient, providing accurate contact angle estimates in roughly 5 nanoseconds of simulation time.

Protein G selects two binding sites for carbon nanotube with dissimilar behavior; a molecular dynamics study

BACKGROUND

Study of nanostructure-protein interaction for development of various types of nano-devices is very essential. Among carbon nanostructures, carbon nanotube (CNT) provides a suitable platform for functionalization by proteins. Previous studies have confirmed that the CNT induces changes in the protein structure.

METHODS

Molecular dynamics (MD) simulation study was employed to illustrate the changes occurring in the protein G (PGB) in the presence of a CNT. In order to predict the PGB surface patches for the CNT, Autodock tools were utilized.

RESULTS

Docking results indicate the presence of two different surface patches with diverse amino acids: the dominant polar residues in the first (PGB-CNT1) and the aromatic residues in the second (PGB-CNT2) surface patch. Displacement of amino acids in the PGB-CNT2 complex occurred during the simulation and it caused an increase in its stability at the end of simulation. The amino acids' displacements diminished the PGB α-helix structure by breakage of hydrogen bonds and generated more transient structures. Principal component analysis determined that the interaction of the CNT with the second surface patch of the PGB raised the extent and modes of the PGB motions. In contrast, insignificant structural changes induced in the PGB while the CNT bonded through the first surface patch.

CONCLUSION

Even though neither of the PGB-CNT complexes could prevent structural changes in the PGB, development of the PGB-CNT1 complex induce slight structural changes in its fragment of crystallizable receptor (FCR). Dissimilar structural changes induced in the PGB-CNT complexes are possibly related to various characteristics of the PGB binding sites.

Interactions between Metal Oxides and Biomolecules: from Fundamental Understanding to Applications

Metallo-oxide (MO)-based bioinorganic nanocomposites promise unique structures, physicochemical properties, and novel biochemical functionalities, and within the past decade, investment in research on materials such as ZnO, TiO₂, SiO₂, and GeO₂ has significantly increased. Besides traditional approaches, the synthesis, shaping, structural patterning, and postprocessing chemical functionalization of the materials surface is inspired by strategies which mimic processes in nature. Would such materials deliver new technologies? Answering this question requires the merging of historical knowledge and current research from different fields of science. Practically, we need an effective defragmentation of the research area. From our perspective, the superficial accounting of material properties, chemistry of the surfaces, and the behavior of biomolecules next to such surfaces is a problem. This is particularly of concern when we wish to bridge between technologies in vitro and biotechnologies in vivo. Further, besides the potential practical technological efficiency and advantages such materials might exhibit, we have to consider the wider long-term implications of material stability and toxicity. In this contribution, we present a critical review of recent advances in the chemistry and engineering of MO-based biocomposites, highlighting the role of interactions at the interface and the techniques by which these can be studied. At the end of the article, we outline the challenges which hamper progress in research and extrapolate to developing and promising directions including additive manufacturing and synthetic biology that could benefit from molecular level understanding of interactions occurring between inanimate (abiotic) and living (biotic) materials.

The Excess Chemical Potential of Water at the Interface with a Protein from End Point Simulations

We use end point simulations to estimate the excess chemical potential of water in the homogeneous liquid and at the interface with a protein in solution. When the pure liquid is taken as the reference, the excess chemical potential of interfacial water is the difference between the solvation free energy of a water molecule at the interface and in the bulk. Using the homogeneous liquid as an example, we show that the solvation free energy for growing a water molecule can be estimated by applying UWHAM to the simulation data generated from the initial and final states (i.e., "the end points") instead of multistate free energy perturbation simulations because of the possible overlaps of the configurations sampled at the end points. Then end point simulations are used to estimate the solvation free energy of water at the interface with a protein in solution. The estimate of the solvation free energy at the interface from two simulations at the end points agrees with the benchmark using 32 states within a 95% confidence interval for most interfacial locations. The ability to accurately estimate the excess chemical potential of water from end point simulations facilitates the statistical thermodynamic analysis of diverse interfacial phenomena. Our focus is on analyzing the excess chemical potential of water at protein receptor binding sites with the goal of using this information to assist in the design of tight binding ligands.

Density Fluctuation in Aqueous Solutions and Molecular Origin of Salting-Out Effect for CO2

Using molecular dynamics simulation, we studied the density fluctuations and cavity formation probabilities in aqueous solutions and their effect on the hydration of CO₂. With increasing salt concentration, we report an increased probability of observing a larger than the average number of species in the probe volume. Our energetic analyses indicate that the van der Waals and electrostatic interactions between CO₂ and aqueous solutions become more favorable with increasing salt concentration, favoring the solubility of CO₂ (salting in). However, due to the decreasing number of cavities forming when salt concentration is increased, the solubility of CO₂ decreases. The formation of cavities was found to be the primary control on the dissolution of gas, and is responsible for the observed CO₂ salting-out effect. Our results provide the fundamental understanding of the density fluctuation in aqueous solutions and the molecular origin of the salting-out effect for real gas.

Molecular Shape and the Hydrophobic Effect

This review focuses on papers published since 2000 on the topic of the properties of solutes in water. More specifically, it evaluates the state of the art of our understanding of the complex relationship between the shape of a hydrophobe and the hydrophobic effect. To highlight this, we present a selection of references covering both empirical and molecular dynamics studies of small (molecular-scale) solutes. These include empirical studies of small molecules, synthetic hosts, crystalline monolayers, and proteins, as well as in silico investigations of entities such as idealized hard and soft spheres, small solutes, hydrophobic plates, artificial concavity, molecular hosts, carbon nanotubes and spheres, and proteins.

Hydrophobicity of proteins and nanostructured solutes is governed by topographical and chemical context

Hydrophobic interactions drive many important biomolecular self-assembly phenomena. However, characterizing hydrophobicity at the nanoscale has remained a challenge due to its nontrivial dependence on the chemistry and topography of biomolecular surfaces. Here we use molecular simulations coupled with enhanced sampling methods to systematically displace water molecules from the hydration shells of nanostructured solutes and calculate the free energetics of interfacial water density fluctuations, which quantify the extent of solute-water adhesion, and therefore solute hydrophobicity. In particular, we characterize the hydrophobicity of curved graphene sheets, self-assembled monolayers (SAMs) with chemical patterns, and mutants of the protein hydrophobin-II. We find that water density fluctuations are enhanced near concave nonpolar surfaces compared with those near flat or convex ones, suggesting that concave surfaces are more hydrophobic. We also find that patterned SAMs and protein mutants, having the same number of nonpolar and polar sites but different geometrical arrangements, can display significantly different strengths of adhesion with water. Specifically, hydroxyl groups reduce the hydrophobicity of methyl-terminated SAMs most effectively not when they are clustered together but when they are separated by one methyl group. Hydrophobin-II mutants show that a charged amino acid reduces the hydrophobicity of a large nonpolar patch when placed at its center, rather than at its edge. Our results highlight the power of water density fluctuations-based measures to characterize the hydrophobicity of nanoscale surfaces and caution against the use of additive approximations, such as the commonly used surface area models or hydropathy scales for characterizing biomolecular hydrophobicity and the associated driving forces of assembly.

Dynamics of Water Monolayers Confined by Chemically Heterogeneous Surfaces: Observation of Surface-Induced Anisotropic Diffusion

Water present in confining geometries plays key roles in many systems of scientific and technological relevance. Prominent examples are living cells and nanofluidic devices. Despite its importance, a complete understanding of the dynamics of water in nanoscale confinement remains elusive. In this work, we use molecular dynamics (MD) simulation to investigate the diffusive dynamics of water monolayers confined in chemically heterogeneous silica slit pores. The effect of chemical heterogeneity is systematically investigated through the fraction f_SiOH of randomly distributed surface sites that possess hydroxyl functional groups. Partial hydroxylation results in heterogeneous surfaces comprising nanoscale hydrophobic and hydrophilic regions. We find that the in-plane diffusivity of water increases monotonically with f_SiOH; at low surface hydroxylation (f_SiOH ≤ 50%), slow water dynamics arise due to the formation of icelike structures in the hydrophobic regions, while at f_SiOH ≥ 75%, surface-water H-bonds in the hydrophilic regions result in faster dynamics. We show that surface patterning with ordered hydrophobic and hydrophilic "stripes" can be used to induce one-dimensional diffusion, with water diffusing through the slit pore preferentially along the direction of the hydrophilic surface patterns.

Studies on electrostatic interactions within model nano-confined aqueous environments of different chemical nature

We study the potential of mean force for pairs of parallel flat surfaces with attractive electrostatic interactions by employing model systems functionalized with different charged, hydrophobic and hydrophilic groups. We study the way in which the local environment (hydrophobic or hydrophilic moieties) modulates the interaction between the attractive charged groups on the plates by removing or attracting nearby water and thus screening or not the electrostatic interaction. To explicitly account for the role of the solvent and the local hydrophobicity, we also perform studies in vacuo. Additionally, the results are compared to that for non-charged plates in order to single out and rationalize the non-additivity of the different non-covalent interactions. Our simulations demonstrate that the presence of neighboring hydrophobic groups promote water removal in the vicinity of the charged groups, thus enhancing charge attraction upon self-assembly. This role of the local hydrophobicity modulating electrostatic interactions is consistent with recent qualitative descriptions in the protein binding context.

Molecular Determinants of Temperature Dependence of Protein Volume Change upon Unfolding

Pressure is a well-known environmental stressor that can either stabilize or destabilize proteins. The volumetric change upon protein unfolding determines the effect of pressure on protein stability, where negative volume changes destabilized proteins at high pressures. High temperature often accompanies high pressure, for example, in the ocean depths near hydrothermal vents or near faults, so it is important to study the effect of temperature on the volumetric change upon unfolding. We previously detailed the magnitude and sign of the molecular determinants of volumetric change, allowing us to quantitatively predict the volumetric change upon protein unfolding. Here, we present a comprehensive analysis of the temperature dependence of the volumetric components of proteins, showing that hydration volume is the primary component that defines expansivities of the native and unfolded states and void volume only contributes slightly to the folded state expansivity.

Water is an active matrix of life for cell and molecular biology

Szent-Győrgi called water the "matrix of life" and claimed that there was no life without it. This statement is true, as far as we know, on our planet, but it is not clear whether it must hold throughout the cosmos. To evaluate that question requires a close consideration of the many varied and subtle roles that water plays in living cells-a consideration that must be free of both an assumed essentialism that gives water an almost mystical life-giving agency and a traditional tendency to see it as a merely passive solvent. Water is a participant in the "life of the cell," and here I describe some of the features of that active agency. Water's value for molecular biology comes from both the structural and dynamic characteristics of its status as a complex, structured liquid as well as its nature as a polar, protic, and amphoteric reagent. Any discussion of water as life's matrix must, however, begin with an acknowledgment that our understanding of it as both a liquid and a solvent is still incomplete.