Small-angle X-ray scattering unveils the internal structure of lipid nanoparticles.

Lipid nanoparticles own a remarkable potential in nanomedicine, only partially disclosed. While the clinical use of liposomes and cationic lipid-nucleic acid complexes is well-established, liquid lipid nanoparticles (nanoemulsions), solid lipid nanoparticles, and nanostructured lipid carriers have even greater possibilities. However, they face obstacles in being used in clinics due to a lack of understanding about the molecular mechanisms controlling their drug loading and release, interactions with the biological environment (such as the protein corona), and shelf-life stability. To create effective drug delivery carriers and successfully translate bench research to clinical settings, it is crucial to have a thorough understanding of the internal structure of lipid nanoparticles. Through synchrotron small-angle X-ray scattering experiments, we determined the spatial distribution and internal structure of the nanoparticles' lipid, surfactant, and the bound water in them. The nanoparticles themselves have a barrel-like shape that consists of coplanar lipid platelets (specifically cetyl palmitate) that are covered by loosely spaced polysorbate 80 surfactant molecules, whose polar heads retain a large amount of bound water. To reduce the interface cost of bound water with unbound water without stacking, the platelets collapse onto each other. This internal structure challenges the classical core-shell model typically used to describe solid lipid nanoparticles and could play a significant role in drug loading and release, biological fluid interaction, and nanoparticle stability, making our findings valuable for the rational design of lipid-based nanoparticles.

Hydrophobic Homopolymer's Coil-Globule Transition and Adsorption onto a Hydrophobic Surface under Different Conditions.

Unstructured proteins can modulate cellular responses to environmental conditions by undergoing coil-globule transitions and phase separation. However, the molecular mechanisms of these phenomena still need to be fully understood. Here, we use Monte Carlo calculations of a coarse-grained model incorporating water's effects on the system's free energy. Following previous studies, we modeled an unstructured protein as a polymer chain. Because we are interested in investigating how it responds to thermodynamic changes near a hydrophobic surface under different conditions, we chose an entirely hydrophobic sequence to maximize the interaction with the interface. We show that a slit pore confinement without top-down symmetry enhances the unfolding and adsorption of the chain in both random coil and globular states. Moreover, we demonstrate that the hydration water modulates this behavior depending on the thermodynamic parameters. Our findings provide insights into how homopolymers and possibly unstructured proteins can sense and adjust to external stimuli such as nanointerfaces or stresses.

Size-Pore-Dependent Methanol Sequestration from Water-Methanol Mixtures by an Embedded Graphene Slit.

The separation of liquid mixture components is relevant to many applications-ranging from water purification to biofuel production-and is a growing concern related to the UN Sustainable Development Goals (SDGs), such as "Clean water and Sanitation" and "Affordable and clean energy". One promising technique is using graphene slit-pores as filters, or sponges, because the confinement potentially affects the properties of the mixture components in different ways, favoring their separation. However, no systematic study has shown how the size of a pore changes the thermodynamics of the surrounding mixture. Here, we focus on water-methanol mixtures and explore, using Molecular Dynamics simulations, the effects of a graphene pore, with size ranging from 6.5 to 13 Å, for three compositions: pure water, 90%-10%, and 75%-25% water-methanol. We show that tuning the pore size can change the mixture pressure, density and composition in bulk due to the size-dependent methanol sequestration within the pore. Our results can help in optimizing the graphene pore size for filtering applications.

Correction: March et al. Protein Unfolding and Aggregation near a Hydrophobic Interface. Polymers 2021, 13, 156.

In the original publication [...].

Many-Body Contributions in Water Nanoclusters.

Many-body interactions in water are known to be important but difficult to treat in atomistic models and often are included only as a correction. Polarizable models treat them explicitly, with long-range many-body potentials, within their classical approximation. However, their calculation is computationally expensive. Here, we evaluate how relevant the contributions to the many-body interaction associated with different coordination shells are. We calculate the global energy minimum, and the corresponding configuration, for nanoclusters of up to 20 water molecules. We find that including the first coordination shell, i.e., the five-body term of the central molecule, is enough to approximate within 5% the global energy minimum and its structure. We show that this result is valid for three different polarizable models, the Dang-Chang, the MB-pol, and the Kozack-Jordan potentials. This result suggests a strategy to develop many-body potentials for water that are reliable and, at the same time, computationally efficient.

Structure and dynamics of nanoconfined water and aqueous solutions.

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Nanoconfined Fluids: Uniqueness of Water Compared to Other Liquids.

Nanoconfinement can drastically change the behavior of liquids, puzzling us with counterintuitive properties. It is relevant in applications, including decontamination and crystallization control. However, it still lacks a systematic analysis for fluids with different bulk properties. Here we address this gap. We compare, by molecular dynamics simulations, three different liquids in a graphene slit pore: (1) A simple fluid, such as argon, described by a Lennard-Jones potential; (2) an anomalous fluid, such as a liquid metal, modeled with an isotropic core-softened potential; and (3) water, the prototypical anomalous liquid, with directional HBs. We study how the slit-pore width affects the structure, thermodynamics, and dynamics of the fluids. All the fluids show similar oscillating properties by changing the pore size. However, their free-energy minima are quite different in nature: (i) are energy-driven for the simple liquid; (ii) are entropy-driven for the isotropic core-softened potential; and (iii) have a changing nature for water. Indeed, for water, the monolayer minimum is entropy driven, at variance with the simple liquid, while the bilayer minimum is energy driven, at variance with the other anomalous liquid. Also, water has a large increase in diffusion for subnm slit pores, becoming faster than bulk. Instead, the other two fluids have diffusion oscillations much smaller than water, slowing down for decreasing slit-pore width. Our results, clarifying that water confined at the subnm scale behaves differently from other (simple or anomalous) fluids under similar confinement, are possibly relevant in nanopores applications, for example, in water purification from contaminants.

Advances in the study of supercooled water.

In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.

Redefining the concept of hydration water near soft interfaces.

Water determines the properties of biological systems. Therefore, understanding the nature of the mutual interaction between water and biosystems is of primary importance for a proper assessment of any biological activity, e.g., the efficacy of new drugs or vaccines. A convenient way to characterize the interactions between biosystems and water is to analyze their impact on water density and dynamics in the proximity of the interfaces. It is commonly accepted that water bulk density and dynamical properties are recovered at distances of the order of 1 nm away from the surface of biological systems. This notion leads to the definition of hydration or biological water as the nanoscopic layer of water covering the surface of biosystems and to the expectation that all the effects of the water-interface interaction are limited to this thin region. Here, we review some of our latest contributions, showing that phospholipid membranes affect the water dynamics, structural properties, and hydrogen bond network at a distance that is more than twice as large as the commonly evoked ∼1nm thick layer and of the order of 2.4 nm. Furthermore, we unveil that at a shorter distance ∼0.5nm from the membrane, instead, there is an additional interface between lipid-bound and unbound water. Bound water has a structural role in the stability of the membrane. Our results imply that the concept of hydration water should be revised or extended and pave the way to a deeper understanding of the mutual interactions between water and biological systems.

Protein Unfolding and Aggregation near a Hydrophobic Interface.

The behavior of proteins near interfaces is relevant for biological and medical purposes. Previous results in bulk show that, when the protein concentration increases, the proteins unfold and, at higher concentrations, aggregate. Here, we study how the presence of a hydrophobic surface affects this course of events. To this goal, we use a coarse-grained model of proteins and study by simulations their folding and aggregation near an ideal hydrophobic surface in an aqueous environment by changing parameters such as temperature and hydrophobic strength, related, e.g., to ions concentration. We show that the hydrophobic surface, as well as the other parameters, affect both the protein unfolding and aggregation. We discuss the interpretation of these results and define future lines for further analysis, with their possible implications in neurodegenerative diseases.

Network Topology in Water Nanoconfined between Phospholipid Membranes.

Water provides the driving force for the assembly and stability of many cellular components. Despite its impact on biological functions, a nanoscale understanding of the relationship between its structure and dynamics under soft confinement has remained elusive. As expected, water in contact with biological membranes recovers its bulk density and dynamics at ∼1 nm from phospholipid headgroups but surprisingly enhances its intermediate range order (IRO) over a distance, at least, twice as large. Here, we explore how the IRO is related to the water's hydrogen-bond network (HBN) and its coordination defects. We characterize the increased IRO by an alteration of the HBN up to more than eight coordination shells of hydration water. The HBN analysis emphasizes the existence of a bound-unbound water interface at ∼0.8 nm from the membrane. The unbound water has a distribution of defects intermediate between bound and bulk water, but with density and dynamics similar to bulk, while bound water has reduced thermal energy and many more HBN defects than low-temperature water. This observation could be fundamental for developing nanoscale models of biological interactions and for understanding how alteration of the water structure and topology, for example, due to changes in extracellular ions concentration, could affect diseases and signaling. More generally, it gives us a different perspective to study nanoconfined water.

In Silico Evidence That Protein Unfolding is a Precursor of Protein Aggregation.

We present a computational study on the folding and aggregation of proteins in an aqueous environment, as a function of its concentration. We show how the increase of the concentration of individual protein species can induce a partial unfolding of the native conformation without the occurrence of aggregates. A further increment of the protein concentration results in the complete loss of the folded structures and induces the formation of protein aggregates. We discuss the effect of the protein interface on the water fluctuations in the protein hydration shell and their relevance in the protein-protein interaction.

Understanding the Kinetics of Protein-Nanoparticle Corona Formation.

When a pristine nanoparticle (NP) encounters a biological fluid, biomolecules spontaneously form adsorption layers around the NP, called "protein corona". The corona composition depends on the time-dependent environmental conditions and determines the NP's fate within living organisms. Understanding how the corona evolves is fundamental in nanotoxicology as well as medical applications. However, the process of corona formation is challenging due to the large number of molecules involved and to the large span of relevant time scales ranging from 100 µs, hard to probe in experiments, to hours, out of reach of all-atoms simulations. Here we combine experiments, simulations, and theory to study (i) the corona kinetics (over 10-3-103 s) and (ii) its final composition for silica NPs in a model plasma made of three blood proteins (human serum albumin, transferrin, and fibrinogen). When computer simulations are calibrated by experimental protein-NP binding affinities measured in single-protein solutions, the theoretical model correctly reproduces competitive protein replacement as proven by independent experiments. When we change the order of administration of the three proteins, we observe a memory effect in the final corona composition that we can explain within our model. Our combined experimental and computational approach is a step toward the development of systematic prediction and control of protein-NP corona composition based on a hierarchy of equilibrium protein binding constants.

Effects of confinement between attractive and repulsive walls on the thermodynamics of an anomalous fluid.

We study via molecular-dynamics simulations the thermodynamics of an anomalous fluid confined in a slit pore with one wall structured and attractive and another unstructured and repulsive. We find that the phase diagram of the homogeneous part of the confined fluid is shifted to higher temperatures, densities, and pressures with respect to the bulk, but it can be rescaled on the bulk case. We calculate a moderate increase of mobility of the homogeneous confined fluid that we interpret as a consequence of the layering due to confinement and the collective modes due to long-range correlations. We show that, as in bulk, the confined fluid has structural, diffusion, and density anomalies that order in the waterlike hierarchy, and a liquid-liquid critical point (LLCP). The overall anomalous region moves to higher temperatures, densities, and pressure, and the LLCP displaces to higher temperature compared to bulk. Motivated by experiments, we calculate also the phase diagram not just for the homogeneous part of the confined fluid but for the entire fluid in the pore, and we show that it is shifted toward higher pressures but preserves the thermodynamics, including the LLCP. Because our model has waterlike properties, we argue that in experiments with supercooled water confined in slit pores with a width of >3 nm if hydrophilic and of >1.5 nm if hydrophobic, the existence of the LLCP could be easier to test than in bulk, where it is not directly accessible.

Structural Interpretation of the Large Slowdown of Water Dynamics at Stacked Phospholipid Membranes for Decreasing Hydration Level: All-Atom Molecular Dynamics.

Hydration water determines the stability and function of phospholipid membranes as well as the interaction of membranes with other molecules. Experiments and simulations have shown that water dynamics slows down dramatically as the hydration decreases, suggesting that the interfacial water that dominates the average dynamics at low hydration is slower than water away from the membrane. Here, based on all-atom molecular dynamics simulations, we provide an interpretation of the slowdown of interfacial water in terms of the structure and dynamics of water-water and water-lipid hydrogen bonds (HBs). We calculate the rotational and translational slowdown of the dynamics of water confined in stacked phospholipid membranes at different levels of hydration, from completely hydrated to poorly hydrated membranes. For all hydrations, we analyze the distribution of HBs and find that water-lipids HBs last longer than water-water HBs and that at low hydration most of the water is in the interior of the membrane. We also show that water-water HBs become more persistent as the hydration is lowered. We attribute this effect (i) to HBs between water molecules that form, in turn, persistent HBs with lipids; (ii) to the hindering of the H-bonding switching between water molecules due to the lower water density at the interface; and (iii) to the higher probability of water-lipid HBs as the hydration decreases. Our interpretation of the large dynamic slowdown in water under dehydration is potentially relevant in understanding membrane biophysics at different hydration levels.

Contribution of Water to Pressure and Cold Denaturation of Proteins.

The mechanisms of cold and pressure denaturation of proteins are matter of debate and are commonly understood as due to water-mediated interactions. Here, we study several cases of proteins, with or without a unique native state, with or without hydrophilic residues, by means of a coarse-grain protein model in explicit solvent. We show, using Monte Carlo simulations, that taking into account how water at the protein interface changes its hydrogen bond properties and its density fluctuations is enough to predict protein stability regions with elliptic shapes in the temperature-pressure plane, consistent with previous theories. Our results clearly identify the different mechanisms with which water participates to denaturation and open the perspective to develop advanced computational design tools for protein engineering.

Optimization of crystal nucleation close to a metastable fluid-fluid phase transition.

The presence of a metastable fluid-fluid critical point is thought to dramatically influence the crystallization pathway, increasing the nucleation rate by many orders of magnitude over the predictions of classical nucleation theory. We use molecular dynamics simulations to study the kinetics of crystallization in the vicinity of this metastable critical point and throughout the metastable fluid-fluid phase diagram. To quantitatively understand how the fluid-fluid phase separation affects the crystal nucleation, we evaluate accurately the kinetics and reconstruct the thermodynamic free-energy landscape of crystal formation. Contrary to expectations, we find no special advantage of the proximity of the metastable critical point on the crystallization rates. However, we find that the ultrafast formation of a dense liquid phase causes the crystallization to accelerate both near the metastable critical point and almost everywhere below the fluid-fluid spinodal line. These results unveil three different scenarios for crystallization that could guide the optimization of the process in experiments.

Structural behavior and dynamics of an anomalous fluid between attractive and repulsive walls: templating, molding, and superdiffusion.

Confinement can modify the dynamics, the thermodynamics, and the structural properties of liquid water, the prototypical anomalous liquid. By considering a generic model for anomalous liquids, suitable for describing solutions of globular proteins, colloids, or liquid metals, we study by molecular dynamics simulations the effect that an attractive wall with structure and a repulsive wall without structure have on the phases, the crystal nucleation, and the dynamics of the fluid. We find that at low temperatures the large density of the attractive wall induces a high-density, high-energy structure in the first layer ("templating" effect). In turn, the first layer induces a "molding" effect on the second layer determining a structure with reduced energy and density, closer to the average density of the system. This low-density, low-energy structure propagates further through the layers by templating effect and can involve all the existing layers at the lowest temperatures investigated. Therefore, although the high-density, high-energy structure does not self-reproduce further than the first layer, the structured wall can have a long-range influence thanks to a sequence of templating, molding, and templating effects through the layers. We find that the walls also have an influence on the dynamics of the liquid, with a stronger effect near the attractive wall. In particular, we observe that the dynamics is largely heterogeneous (i) among the layers, as a consequence of the sequence of structures caused by the walls presence, and (ii) within the same layer, due to superdiffusive liquid veins within a frozen matrix of particles near the walls at low temperature and high density. Hence, the partial freezing of the first layer does not correspond necessarily to an effective reduction of the channel's section in terms of transport properties, as suggested by other authors.

Critical behavior of a water monolayer under hydrophobic confinement.

The properties of water can have a strong dependence on the confinement. Here, we consider a water monolayer nanoconfined between hydrophobic parallel walls under conditions that prevent its crystallization. We investigate, by simulations of a many-body coarse-grained water model, how the properties of the liquid are affected by the confinement. We show, by studying the response functions and the correlation length and by performing finite-size scaling of the appropriate order parameter, that at low temperature the monolayer undergoes a liquid-liquid phase transition ending in a critical point in the universality class of the two-dimensional (2D) Ising model. Surprisingly, by reducing the linear size L of the walls, keeping the walls separation h constant, we find a 2D-3D crossover for the universality class of the liquid-liquid critical point for L=h ≃ 50, i.e. for a monolayer thickness that is small compared to its extension. This result is drastically different from what is reported for simple liquids, where the crossover occurs for L=h ≃ 5, and is consistent with experimental results and atomistic simulations. We shed light on these findings showing that they are a consequence of the strong cooperativity and the low coordination number of the hydrogen bond network that characterizes water.

Effect of pressure on the anomalous response functions of a confined water monolayer at low temperature.

We study a coarse-grained model for a water monolayer that cannot crystallize due to the presence of confining interfaces, such as protein powders or inorganic surfaces. Using both Monte Carlo simulations and mean field calculations, we calculate three response functions: the isobaric specific heat C(P), the isothermal compressibility K(T), and the isobaric thermal expansivity α(P). At low temperature T, we find two distinct maxima in C(P), K(T), and ∣α(P)∣, all converging toward a liquid-liquid critical point (LLCP) with increasing pressure P. We show that the maximum in C(P) at higher T is due to the fluctuations of hydrogen (H) bond formation and that the second maximum at lower T is due to the cooperativity among the H bonds. We discuss a similar effect in K(T) and ∣α(P)∣. If this cooperativity were not taken into account, both the lower-T maximum and the LLCP would disappear. However, comparison with recent experiments on water hydrating protein powders provides evidence for the existence of the lower-T maximum, supporting the hypothesized LLCP at positive P and finite T. The model also predicts that when P moves closer to the critical P the C(P) maxima move closer in T until they merge at the LLCP. Considering that other scenarios for water are thermodynamically possible, we discuss how an experimental measurement of the changing separation in T between the two maxima of C(P) as P increases could determine the best scenario for describing water.