On the engulfment of antifreeze proteins by ice.

Antifreeze proteins (AFPs) are remarkable biomolecules that suppress ice formation at trace concentrations. To inhibit ice growth, AFPs must not only bind to ice crystals, but also resist engulfment by ice. The highest supercooling, [Formula: see text], for which AFPs are able to resist engulfment is widely believed to scale as the inverse of the separation, [Formula: see text], between bound AFPs, whereas its dependence on the molecular characteristics of the AFP remains poorly understood. By using specialized molecular simulations and interfacial thermodynamics, here, we show that in contrast with conventional wisdom, [Formula: see text] scales as [Formula: see text] and not as [Formula: see text]. We further show that [Formula: see text] is proportional to AFP size and that diverse naturally occurring AFPs are optimal at resisting engulfment by ice. By facilitating the development of AFP structure-function relationships, we hope that our findings will pave the way for the rational design of AFPs.

Accumulation of Antifreeze Proteins on Ice Is Determined by Adsorption.

Antifreeze proteins (AFPs) facilitate the survival of diverse organisms in frigid environments by adsorbing to ice crystals and suppressing their growth. The rate of AFP accumulation on ice is determined by an interplay between AFP diffusion from the bulk solution to the ice-water interface and the subsequent adsorption of AFPs to the interface. To interrogate the relative importance of these two processes, here, we combine nonequilibrium fluorescence experiments with a reaction-diffusion model. We find that as diverse AFPs accumulate on ice, their concentration in the aqueous solution does not develop a gradient but remains equal to its bulk concentration throughout our experiments. These findings lead us to conclude that AFP accumulation on ice crystals, which are smaller than 100 µm in radius, is not limited by the diffusion of AFPs, but by the kinetics of AFP adsorption. Our results imply that mass transport limitations do not hinder AFPs from performing their biological function.

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling.

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer icephilicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface-ice contact at the expense of surface-water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

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.

PRC2-Inactivating Mutations in Cancer Enhance Cytotoxic Response to DNMT1-Targeted Therapy via Enhanced Viral Mimicry.

Polycomb repressive complex 2 (PRC2) has oncogenic and tumor-suppressive roles in cancer. There is clinical success of targeting this complex in PRC2-dependent cancers, but an unmet therapeutic need exists in PRC2-loss cancer. PRC2-inactivating mutations are a hallmark feature of high-grade malignant peripheral nerve sheath tumor (MPNST), an aggressive sarcoma with poor prognosis and no effective targeted therapy. Through RNAi screening in MPNST, we found that PRC2 inactivation increases sensitivity to genetic or small-molecule inhibition of DNA methyltransferase 1 (DNMT1), which results in enhanced cytotoxicity and antitumor response. Mechanistically, PRC2 inactivation amplifies DNMT inhibitor-mediated expression of retrotransposons, subsequent viral mimicry response, and robust cell death in part through a protein kinase R (PKR)-dependent double-stranded RNA sensor. Collectively, our observations posit DNA methylation as a safeguard against antitumorigenic cell-fate decisions in PRC2-loss cancer to promote cancer pathogenesis, which can be therapeutically exploited by DNMT1-targeted therapy. SIGNIFICANCE: PRC2 inactivation drives oncogenesis in various cancers, but therapeutically targeting PRC2 loss has remained challenging. Here we show that PRC2-inactivating mutations set up a tumor context-specific liability for therapeutic intervention via DNMT1 inhibitors, which leads to innate immune signaling mediated by sensing of derepressed retrotransposons and accompanied by enhanced cytotoxicity. See related commentary by Guil and Esteller, p. 2020. This article is highlighted in the In This Issue feature, p. 2007.

Tumor-intrinsic PRC2 inactivation drives a context-dependent immune-desert microenvironment and is sensitized by immunogenic viruses.

Immune checkpoint blockade (ICB) has demonstrated clinical success in "inflamed" tumors with substantial T cell infiltrates, but tumors with an immune-desert tumor microenvironment (TME) fail to benefit. The tumor cell-intrinsic molecular mechanisms of the immune-desert phenotype remain poorly understood. Here, we demonstrated that inactivation of the polycomb-repressive complex 2 (PRC2) core components embryonic ectoderm development (EED) or suppressor of zeste 12 homolog (SUZ12), a prevalent genetic event in malignant peripheral nerve sheath tumors (MPNSTs) and sporadically in other cancers, drove a context-dependent immune-desert TME. PRC2 inactivation reprogramed the chromatin landscape that led to a cell-autonomous shift from primed baseline signaling-dependent cellular responses (e.g., IFN-γ signaling) to PRC2-regulated developmental and cellular differentiation transcriptional programs. Further, PRC2 inactivation led to diminished tumor immune infiltrates through reduced chemokine production and impaired antigen presentation and T cell priming, resulting in primary resistance to ICB. Intratumoral delivery of inactivated modified vaccinia virus Ankara (MVA) enhanced tumor immune infiltrates and sensitized PRC2-loss tumors to ICB. Our results identify molecular mechanisms of PRC2 inactivation-mediated, context-dependent epigenetic reprogramming that underline the immune-desert phenotype in cancer. Our studies also point to intratumoral delivery of immunogenic viruses as an initial therapeutic strategy to modulate the immune-desert TME and capitalize on the clinical benefit of ICB.

Characterizing the Interplay between Polymer Solvation and Conformation.

Conformational transitions of flexible molecules, especially those driven by hydrophobic effects, tend to be hindered by desolvation barriers. For such transitions, it is thus important to characterize and understand the interplay between solvation and conformation. Using specialized molecular simulations, here we perform such a characterization for a hydrophobic polymer solvated in water. We find that an external potential, which unfavorably perturbs the polymer hydration waters, can trigger a coil-to-globule or collapse transition, and that the relative stabilities of the collapsed and extended states can be quantified by the strength of the requisite potential. Our results also provide mechanistic insights into the collapse transition, highlighting that the bottleneck to polymer collapse is the formation of a sufficiently large cluster, and the collective dewetting of such a cluster. We also study the collapse of the hydrophobic polymer in octane, a nonpolar solvent, and interestingly, we find that the mechanistic details of the transition are qualitatively similar to that in water.

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

Correction for 'Characterizing surface wetting and interfacial properties using enhanced sampling (SWIPES)' by Hao Jiang et al., Soft Matter, 2019, 15, 860-869, DOI: .

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.

Overcoming BET Inhibitor Resistance in Malignant Peripheral Nerve Sheath Tumors.

PURPOSE: BET bromodomain inhibitors have emerged as a promising therapy for numerous cancer types in preclinical studies, including neurofibromatosis type 1 (NF1)-associated malignant peripheral nerve sheath tumor (MPNST). However, potential mechanisms underlying resistance to these inhibitors in different cancers are not completely understood. In this study, we explore new strategy to overcome BET inhibitor resistance in MPNST.Experimental Design: Through modeling tumor evolution by studying genetic changes underlying the development of MPNST, a lethal sarcoma with no effective medical treatment, we identified a targetable addiction to BET bromodomain family member BRD4 in MPNST. This served as a controlled model system to delineate mechanisms of sensitivity and resistance to BET bromodomain inhibitors in this disease. RESULTS: Here, we show that a malignant progression-associated increase in BRD4 protein levels corresponds to partial sensitivity to BET inhibition in MPNST. Strikingly, genetic depletion of BRD4 protein levels synergistically sensitized MPNST cells to diverse BET inhibitors in culture and in vivo. CONCLUSIONS: Collectively, MPNST sensitivity to combination genetic and pharmacologic inhibition of BRD4 revealed the presence of a unique addiction to BRD4 in MPNST. Our discovery that a synthetic lethality exists between BET inhibition and reduced BRD4 protein levels nominates MPNST for the investigation of emerging therapeutic interventions such as proteolysis-targeting chimeras (PROTACs) that simultaneously target bromodomain activity and BET protein abundance.

Characterizing Solvent Density Fluctuations in Dynamical Observation Volumes.

Hydrophobic effects drive diverse aqueous assemblies, such as micelle formation or protein folding, wherein the solvent plays an important role. Consequently, characterizing the free energetics of solvent density fluctuations can lead to important insights into these processes. Although techniques such as the indirect umbrella sampling (INDUS) method can be used to characterize solvent fluctuations in static observation volumes of various sizes and shapes, characterizing how the solvent mediates inherently dynamic processes, such as self-assembly or conformational change, remains a challenge. In this work, we generalize the INDUS method to facilitate the enhanced sampling of solvent fluctuations in dynamical observation volumes, whose positions and shapes can evolve. We illustrate the usefulness of this generalization by characterizing water density fluctuations in dynamical volumes pertaining to the hydration of flexible solutes, the assembly of small hydrophobes, and conformational transitions in a model peptide. We also use the method to probe the dynamics of hard spheres.

Protein Hydration Waters Are Susceptible to Unfavorable Perturbations.

The interactions of a protein, its phase behavior, and, ultimately, its ability to function are all influenced by the interactions between the protein and its hydration waters. Here, we study proteins with a variety of sizes, shapes, chemistries, and biological functions and characterize their interactions with their hydration waters using molecular simulations and enhanced sampling techniques. We find that, akin to extended hydrophobic surfaces, proteins situate their hydration waters at the edge of a dewetting transition, making them susceptible to unfavorable perturbations. We also find that the strength of the unfavorable potential needed to trigger dewetting is roughly the same for all proteins studied here and depends primarily on the width of the hydration shell being perturbed. Our findings establish a framework for systematically classifying protein patches according to how favorably they interact with water.

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 Hydration Thermodynamics: The Influence of Flexibility and Salt on Hydrophobin II Hydration.

The solubility of proteins and other macromolecular solutes plays an important role in numerous biological, chemical, and medicinal processes. An important determinant of protein solubility is the solvation free energy of the protein, which quantifies the overall strength of the interactions between the protein and the aqueous solution that surrounds it. Here we present an all-atom explicit-solvent computational framework for the rapid estimation of protein solvation free energies. Using this framework, we estimate the hydration free energy of hydrophobin II, an amphiphilic fungal protein, in a computationally efficient manner. We further explore how the protein hydration free energy is influenced by enhancing flexibility and by the addition of sodium chloride, and find that it increases in both cases, making protein hydration less favorable.

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.

The role of nerve microenvironment for neurofibroma development.

Deregulation of RAS signaling in Neurofibromatosis type 1 (NF1) results in the development of multiple neurofibromas, complex tumor of the peripheral nerves with no effective medical treatment. There is increasing evidences that neurofibroma initiates through loss of NF1 function in the Schwann cell lineage, followed by a cascade of interactions with other cell types in the surrounding tumor microenvironment. In NF1 patients, neurofibromas always develop along peripheral nerves and do not migrate to distant organs, including the central nervous system. In this study, we sought to identify the contributions of these peripheral nerves in neurofibroma formation. Using in vivo and in vitro three-dimensional (3D) culturing system, we show that peripheral nerves are absolutely required for neurofibroma tumorigenesis and report a novel 3D skin raft culture system for neurofibroma formation in vitro to decipher tumor pathogenesis. This interaction between neoplastic Schwann cells and their surrounding neural microenvironment has important implications for understanding early cellular events that dictate tumorigenesis. It also provides fertile ground for the elucidation of intrinsic and extrinsic factors within the nerve microenvironment that likely play essential roles in neurofibroma development and, therefore, viable therapeutic targets in neurofibroma therapy.

Spontaneous recovery of superhydrophobicity on nanotextured surfaces.

Rough or textured hydrophobic surfaces are dubbed "superhydrophobic" due to their numerous desirable properties, such as water repellency and interfacial slip. Superhydrophobicity stems from an aversion of water for the hydrophobic surface texture, so that a water droplet in the superhydrophobic "Cassie state" contacts only the tips of the rough surface. However, superhydrophobicity is remarkably fragile and can break down due to the wetting of the surface texture to yield the "Wenzel state" under various conditions, such as elevated pressures or droplet impact. Moreover, due to large energetic barriers that impede the reverse transition (dewetting), this breakdown in superhydrophobicity is widely believed to be irreversible. Using molecular simulations in conjunction with enhanced sampling techniques, here we show that on surfaces with nanoscale texture, water density fluctuations can lead to a reduction in the free energetic barriers to dewetting by circumventing the classical dewetting pathways. In particular, the fluctuation-mediated dewetting pathway involves a number of transitions between distinct dewetted morphologies, with each transition lowering the resistance to dewetting. Importantly, an understanding of the mechanistic pathways to dewetting and their dependence on pressure allows us to augment the surface texture design, so that the barriers to dewetting are eliminated altogether and the Wenzel state becomes unstable at ambient conditions. Such robust surfaces, which defy classical expectations and can spontaneously recover their superhydrophobicity, could have widespread importance, from underwater operation to phase-change heat transfer applications.