Water at the nanoscale: From filling or dewetting hydrophobic pores and carbon nanotubes to "sliding" on graphene.

In this work, we study the effect of nanoconfinement on the hydration properties of model hydrophobic pores and carbon nanotubes, determining their wetting propensity and the conditions for geometrically induced dehydration. By employing a recently introduced water structural index, we aim at two main goals: (1) to accurately quantify the local hydrophobicity and predict the drying transitions in such systems, and (2) to provide a molecular rationalization of the wetting process. In this sense, we will further discuss the number and strength of the interactions required by the water molecules to promote wetting. In the case of graphene-like surfaces, an explanation for their unexpectedly significant hydrophilicity will also be provided. On the one hand, the structural index will show that the net attraction to the dense carbon network that a water molecule experiences through several simultaneous weak interactions is sufficient to give rise to hydrophilic behavior. On the other hand, we will show that an additional effect is also at play: the hydrating water molecule is retained on the surface by a smooth exchange of such simultaneous weak interactions, as if "sliding" on graphene.

Persistent Local Structural Defectiveness as an Early Time Predictor of Intermittent Glassy Relaxation Events in Supercooled Water.

To become a glass from the metastable supercooled state, a liquid experiences a dramatic dynamical slowing down within a narrow temperature window. However, the attainment of solid rigidity is not the result of breaking translational symmetry as in a crystal: the structure of the resulting amorphous solid strikingly resembles that of the liquid state. Moreover, the supercooled liquid is dynamically heterogeneous; that is, the dynamics varies by orders of magnitude from one region of the sample to another, but the establishment of the existence of strong structural differences between such regions has demanded hard efforts along the years. In this work, we focus precisely on such a structure-dynamics link for supercooled water showing that local regions with structural defectiveness are persistent during the structural relaxation of the system, hence acting as early time predictors of later intermittent glassy relaxation events.

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.

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.

Comparing the performance of two structural indicators for different water models while seeking for connections between structure and dynamics in the glassy regime.

In this work, we compare the performance of two structural indicators based on the degree of translational order up to the second coordination shell in three water models: SPC/E, TIP4P/2005, and TIP5P. Beyond directly contrasting their distributions for different temperatures to evidence their usefulness in estimating the fraction of structured and unstructured molecules and, when possible, their classification capability, we also correlate them with an indirect measure of structural constraint: the dynamic propensity. Furthermore, this procedure enables us to show the existence of evident correlations between structural and dynamical information. More specifically, we find that locally structured molecules display a preference for low dynamic propensity values and, more conspicuously, that locally unstructured molecules are extremely subject to high dynamic propensity. This result is particularly relevant for the supercooled regime where the establishment of firm links between the structure and dynamics has remained rather elusive since the occurrence of dynamics that vary in orders of magnitude upon supercooling usually contrast with barely noticeable overall structural changes.

Size dependence of dynamic fluctuations in liquid and supercooled water.

We study the evolution of dynamic fluctuations averaged over different space lengths and time scales to characterize spatially and temporally heterogeneous behavior of TIP4P/2005 water in liquid and supercooled states. Analyzing a 250 000 molecules simulated system, we provide evidence of the existence, upon supercooling, of a significant enhancement of spatially localized dynamic fluctuations stemming from regions of correlated mobile molecules. We show that both the magnitude of the departure from the value expected for the system-size dependence of an uncorrelated system and the system size at which such a trivial regime is finally recovered clearly increase upon supercooling. This provides a means to estimate an upper limit to the maximum length scale of influence of the regions of correlated mobile molecules. Notably, such an upper limit grows two orders of magnitude on cooling, reaching a value corresponding to a few thousand molecules at the lowest investigated temperature.

Design, synthesis and biological evaluation of 1,3-dihydroxyxanthone derivatives: Effective agents against acetylcholinesterase.

The present work concerns the rational design and development of new inhibitors of acetylcholinesterase (AChE) based on the privileged xanthone scaffold. In order to understand and rationalize the mode of action of these target structures a theoretical study was initially conducted. From the results of rational design, a new variety of amphiphilic xanthone derivatives were synthesized, structurally characterized and evaluated as potential anti-Alzheimer agents. The results showed that most of the synthesized compounds exhibited high AChE inhibitory activity at the micromolar range (IC50, 0.46-12.09µM). The synthetic xanthone 11 showed the best inhibitory effect on AChE and a molecular modeling study revealed that 11 targeted both the catalytic active site (CAS) and the peripheral anionic site (PAS) of AChE. Therefore, this compound could be considered asa potential lead compound towards new drugs for the treatment of Alzheimer's disease.

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.

Structure and dynamics of high- and low-density water molecules in the liquid and supercooled regimes.

By combining the local structure index with potential energy minimisations we study the local environment of the water molecules for a couple of water models, TIP5P-Ew and SPC/E, in order to characterise low- and high-density "species". Both models show a similar behaviour within the supercooled regime, with two clearly distinguishable populations of unstructured and structured molecules, the fraction of the latter increasing with supercooling. Additionally, for TIP5P-Ew, we find that the structured component vanishes quickly at the normal liquid regime (above the melting temperature). Thus, while SPC/E provides a fraction of structured molecules similar to that found in X-ray experiments, we show that TIP5P-Ew underestimates such value. Moreover, unlike SPC/E, we demonstrate that TIP5P-Ew does not follow the linear dependence of the logarithm of the structured fraction with inverse temperature, as predicted by the two-order parameter model. Finally, we link structure to dynamics by showing that there exists a strong correlation between structural fluctuation and dynamics in the supercooled state with spatial correlations in both static and dynamic quantities.

Hydrogen Bond Dynamic Propensity Studies for Protein Binding and Drug Design.

We study the dynamic propensity of the backbone hydrogen bonds of the protein MDM2 (the natural regulator of the tumor suppressor p53) in order to determine its binding properties. This approach is fostered by the observation that certain backbone hydrogen bonds at the p53-binding site exhibit a dynamical propensity in simulations that differs markedly form their state-value (that is, formed/not formed) in the PDB structure of the apo protein. To this end, we conduct a series of hydrogen bond propensity calculations in different contexts: 1) computational alanine-scanning studies of the MDM2-p53 interface; 2) the formation of the complex of MDM2 with the disruptive small molecule Nutlin-3a (dissecting the contribution of the different molecular fragments) and 3) the binding of a series of small molecules (drugs) with different affinities for MDM2. Thus, the relevance of the hydrogen bond propensity analysis for protein binding studies and as a useful tool to complement existing methods for drug design and optimization will be made evident.

Hydrophilic behavior of graphene and graphene-based materials.

Graphene and the graphene-based materials like graphite, carbon nanotubes, and fullerenes are not only usually regarded as hydrophobic but also have been widely employed as paradigms for the investigation of the behavior of water under nonpolar confinement, a question of major concern for fields ranging from biology to materials design. However, some experimental and theoretical insights seem to contradict, at least partially, such a picture. In this work, we will provide firm evidence for a neat hydrophilic nature of graphene surfaces. Our molecular dynamics studies will demonstrate that parallel graphene sheets present a strong tendency to remain fully hydrated for moderately long times (even when the equilibrium state is indeed the collapse of the plates), and thus, they are less prone to self-assembly than model hydrophobic surfaces we shall employ as control which readily undergo a hydrophobic collapse. Potential of mean force calculations will indeed make evident that the solvent exerts a repulsive contribution on the self-assembly of graphene surfaces. Moreover, we shall also quantify graphene hydrophilicity by means of the calculation of water density at two pressures and water density fluctuations. This latter study has never been performed on graphene and represents a means both to confirm and to quantify its neat hydrophilic behavior. We shall also make evident the relevance of the mildly attractive water-carbon interactions, since their artificial weakening will be shown to revert from typically hydrophilic to typically hydrophobic behavior.

Protein packing defects "heat up" interfacial water.

Ligands must displace water molecules from their corresponding protein surface binding site during association. Thus, protein binding sites are expected to be surrounded by non-tightly-bound, easily removable water molecules. In turn, the existence of packing defects at protein binding sites has been also established. At such structural motifs, named dehydrons, the protein backbone is exposed to the solvent since the intramolecular interactions are incompletely wrapped by non-polar groups. Hence, dehydrons are sticky since they depend on additional intermolecular wrapping in order to properly protect the structure from water attack. Thus, a picture of protein binding is emerging wherein binding sites should be both dehydrons rich and surrounded by easily removable water. In this work we shall indeed confirm such a link between structure and dynamics by showing the existence of a firm correlation between the degree of underwrapping of the protein chain and the mobility of the corresponding hydration water molecules. In other words, we shall show that protein packing defects promote their local dehydration, thus producing a region of "hot" interfacial water which might be easily removed by a ligand upon association.

Wrapping effects within a proposed function-rescue strategy for the Y220C oncogenic mutation of protein p53.

Soluble proteins must protect their structural integrity from water attack by wrapping interactions which imply the clustering of nonpolar residues around the backbone hydrogen bonds. Thus, poorly wrapped hydrogen bonds constitute defects which have been identified as promoters of protein associations since they favor the removal of hydrating molecules. More specifically, a recent study of our group has shown that wrapping interactions allow the successful identification of protein binding hot spots. Additionally, we have also shown that drugs disruptive of protein-protein interfaces tend to mimic the wrapping behavior of the protein they replace. Within this context, in this work we study wrapping three body interactions related to the oncogenic Y220C mutation of the tumor suppressor protein p53. Our computational results rationalize the oncogenic nature of the Y220C mutation, explain the binding of a drug-like molecule already designed to restore the function of p53 and provide clues to help improve this function-rescue strategy and to apply in other drug design or re-engineering techniques.

Wrapping mimicking in drug-like small molecules disruptive of protein-protein interfaces.

The discovery of small-molecule drugs aimed at disrupting protein-protein associations is expected to lead to promising therapeutic strategies. The small molecule binds to the target protein thus replacing its natural protein partner. Noteworthy, structural analysis of complexes between successful disruptive small molecules and their target proteins has suggested the possibility that such ligands might somehow mimic the binding behavior of the protein they replace. In these cases, the molecules show a spatial and "chemical" (i.e., hydrophobicity) similarity with the residues of the partner protein involved in the protein-protein complex interface. However, other disruptive small molecules do not seem to show such spatial and chemical correspondence with the replaced protein. In turn, recent progress in the understanding of protein-protein interactions and binding hot spots has revealed the main role of intermolecular wrapping interactions: three-body cooperative correlations in which nonpolar groups in the partner protein promote dehydration of a two-body electrostatic interaction of the other protein. Hence, in the present work, we study some successful complexes between already discovered small disruptive drug-like molecules and their target proteins already reported in the literature and we compare them with the complexes between such proteins and their natural protein partners. Our results show that the small molecules do in fact mimic to a great extent the wrapping behavior of the protein they replace. Thus, by revealing the replacement the small molecule performs of relevant wrapping interactions, we convey precise physical meaning to the mimicking concept, a knowledge that might be exploited in future drug-design endeavors.