How a few help all: cooperative crossing of lipid membranes by COSAN anions.

Experimental results show that the presence of a concentration gradient of certain nano-ions (most notably cobaltabisdicarbollide ([o-COSAN]- anions), induce a current across intact artificial phospholipid bilayers in spite of the high Born free energy estimated for these ions. The mechanism underlying this observed translocation of nano-anions across membranes has yet to be determined. Here we show, using molecular dynamics simulations, that the permeation of [o-COSAN]- anions across a lipid bilayer proceeds in a cooperative manner. Single nano-ions can enter the bilayer but permeation is hampered by a free energy barrier of about 8kBT. The interaction between these nano-ions inside a leaflet induces a flip-flop translocation mechanism with the formation of transient, elongated structure inside the membrane. This cooperative flip-flop allows an efficient distribution of [o-COSAN]- anions in both leaflets of the bilayer. These results suggest the existence of a new mechanism for permeation of nano-ions across lipid membranes, relevant for those that have the appropriate self-assembly character.

Computer simulations of the interaction between SARS-CoV-2 spike glycoprotein and different surfaces.

A prominent feature of coronaviruses is the presence of a large glycoprotein spike protruding from a lipidic membrane. This glycoprotein spike determines the interaction of coronaviruses with the environment and the host. In this paper, we perform all atomic molecular dynamics simulations of the interaction between the SARS-CoV-2 trimeric glycoprotein spike and surfaces of materials. We considered a material with high hydrogen bonding capacity (cellulose) and a material capable of strong hydrophobic interactions (graphite). Initially, the spike adsorbs to both surfaces through essentially the same residues belonging to the receptor binding subunit of its three monomers. Adsorption onto cellulose stabilizes in this configuration, with the help of a large number of hydrogen bonds developed between cellulose and the three receptor-binding domains of the glycoprotein spike. In the case of adsorption onto graphite, the initial adsorption configuration is not stable and the surface induces a substantial deformation of the glycoprotein spike with a large number of adsorbed residues not pertaining to the binding subunits of the spike monomers.

Surface specificity and mechanistic pathway of de-fluorination of C60F48 on coinage metals.

We provide experimental and theoretical understanding on fundamental processes taking place at room temperature when a fluorinated fullerene dopant gets close to a metal surface. By employing scanning tunneling microscopy and photoelectron spectroscopies, we demonstrate that the on-surface integrity of C60F48 depends on the interaction with the particular metal it approaches. Whereas on Au(111) the molecule preserves its chemical structure, on more reactive surfaces such as Cu(111) and Ni(111), molecules interacting with the bare metal surface lose the halogen atoms and transform to C60. Though fluorine-metal bonding can be detected depending on the molecular surface density, no ordered fluorine structures are observed. We show the implications of the metal-dependent de-fluorination in the electronic structure of the molecules and the energy alignment at the molecule-metal interface. Molecular dynamics simulations with ReaxFF reactive force field corroborate the experimental facts and provide a detailed mechanistic picture of the surface-induced de-fluorination, which involves the rotation of the molecule on the surface. Outstandingly, a thermodynamic analysis indicates that the effect of the metal surface is lowering and diminishing the energy barrier for C-F cleave, demonstrating the catalytic role of the surface. The present study contributes to in-depth knowledge of the mechanisms that affect the degree of stability of chemical species on surfaces, which is essential to advance our understanding of the chemical reactivity of metals and their role in on-surface chemical reactions.

Atomistic Simulations of COSAN: Amphiphiles without a Head-and-Tail Design Display "Head and Tail" Surfactant Behavior.

Cobaltabisdicarbollide (COSAN) anions have an unexpectedly rich self-assembly behavior, which can lead to vesicles and micelles without having a classical surfactant molecular architecture. This was rationalized by the introduction of new terminology and novel driving forces. A key aspect in the interpretation of COSAN behavior is the assumption that the most stable form of these ions is the transoid rotamer, which lacks a "hydrophilic head" and a "hydrophobic tail". Using implicit solvent DFT calculations and MD simulations we show that in water, 1) the cisoid rotamer is the most stable form of COSAN and 2) this cisoid rotamer has a well-defined hydrophilic polar region ("head") and a hydrophobic apolar region ("tail"). In addition, our simulations show that the properties of this rotamer in water (interfacial affinity, micellization) match those expected for a classical surfactant. Therefore, we conclude that the experimental results for the COSAN ions can now be understood in terms of its amphiphilic molecular architecture.

Effect of collagenase-gelatinase ratio on the mechanical properties of a collagen fibril: a combined Monte Carlo-molecular dynamics study.

Loading in cartilage is supported primarily by fibrillar collagen, and damage will impair the function of the tissue, leading to pathologies such as osteoarthritis. Damage is initiated by two types of matrix metalloproteinases, collagenase and gelatinase, that cleave and denature the collagen fibrils in the tissue. Experimental and modeling studies have revealed insights into the individual contributions of these two types of MMPs, as well as the mechanical response of intact fibrils and fibrils that have experienced random surface degradation. However, no research has comprehensively examined the combined influences of collagenases and gelatinases on collagen degradation nor studied the mechanical consequences of biological degradation of collagen fibrils. Such preclinical examinations are required to gain insights into understanding, treating, and preventing degradation-related cartilage pathology. To develop these insights, we use sequential Monte Carlo and molecular dynamics simulations to probe the effect of enzymatic degradation on the structure and mechanics of a single collagen fibril. We find that the mechanical response depends on the ratio of collagenase to gelatinase-not just the amount of lost fibril mass-and we provide a possible mechanism underlying this phenomenon. Overall, by characterizing the combined influences of collagenases and gelatinases on fibril degradation and mechanics at the preclinical research stage, we gain insights that may facilitate the development of targeted interventions to prevent the damage and loss of mechanical integrity that can lead to cartilage pathology.

Molecular insight into the wetting behavior and amphiphilic character of cellulose nanocrystals.

The study of nanocellulose is a field of growing interest due to its many applications and its use in the development of biocompatible and eco-friendly materials. In spite of the vast number of studies in the field, many questions about the role of the molecular structure in the properties of cellulose are still subject of debate. One of these fundamental questions is the possible amphiphilic nature of cellulose and the relative role of hydrogen bonding and hydrophobic effect on the interactions of cellulose. In this work we present an extensive molecular dynamics simulation study of this question by analyzing the wetting of cellulose with water and organic solvent, its interaction with hydrophilic and hydrophobic ions and its interaction with a protein (human epidermal growth factor, hEGF). We consider two characteristic cellulose crystal planes of Iß cellulose with very different roughness, different hydrogen bonding capability and different exposure of cellulose hydrophobic groups (the (010) plane which has exposed -OH groups and the (100) plane with buried -OH groups). Our results show that both surfaces are simultaneously hydrophilic and lipophilic, with both surfaces having very similar contact angles. In spite of the global similarity of wetting of both surfaces, the molecular details of wetting are very different and substantial local wetting heterogeneities (which strongly depend on the surface) appear for both solvents. We also observe a weak interaction of both surfaces with hydrophobic and hydrophilic solutes. These weak interactions are attributed to the simultaneous lipophilic and hydrophilic character of both (100) and (010) cellulose surfaces. Interestingly, we found a substantial interaction of both cellulose planes with polar and apolar residues of the hEGF protein.

Insights into Preformed Human Serum Albumin Corona on Iron Oxide Nanoparticles: Structure, Effect of Particle Size, Impact on MRI Efficiency, and Metabolization.

In the past decade, profuse research efforts explored the uses of iron oxide particles in nanomedicine. To a great extent, the efficiency and fate of those magnetic nanoparticles depend on how their surfaces interface with the proteins in a physiological environment. It is well reported how an ungoverned protein corona can be detrimental to cellular uptake and targeting efficiency and how it can modify the nanoparticles biodistribution. Novel strategies are emerging to achieve enhanced and more reproducible performances of engineered nanoparticles with a custom-built protein corona. Here we report on a generalized protocol to preform a monolayer of human serum albumin (HSA) on superparamagnetic iron oxide nanoparticles (SPIONs) of different sizes. The resulting molecular structures are described by molecular dynamics simulations of the hybrid nanoconjugates. The simulations outcomes regarding the number of proteins in the corona and their monolayer arrangement on the particle surface are in agreement with the results obtained from dynamic light scattering and electronic microscopy analysis. Using tryptophan fluorescence quenching, we revealed the existence of a strong interaction between the SPIONs and the HSA which endorses the robustness of the protein-nanoparticle conjugates in this system. Moreover, we evaluated the effect of the HSA corona on the SPIONs efficiency as magnetic resonance imaging (MRI) contrast agents in water, human serum, and saline media. The protein corona did not affect the efficiency of the SPIONs as T2 contrast agents but reduce their T1 efficiency. In addition, we observed a greater stability for HSA-SPIONs nanoconjugates in saline and in acid media, preventing nanoparticle dissolution in extreme gastric conditions.

Coarse-grained molecular dynamics simulation of the interface behaviour and self-assembly of CTAB cationic surfactants.

In this work we study the behaviour at interfaces and the micelle self-assembly of a cationic surfactant (CTAB) by Molecular Dynamics (MD) simulations of coarse-grained models. We consider both the standard (with explicit water) Martini force field and the implicit solvent version of the Martini force field (Dry Martini). First, we study the behaviour of CTAB at a water/vacuum interface, at a water/organic solvent interface and in a pre-assembled CTAB micelle using both standard and Dry Martini and all-atomic simulations. Our results indicate that there are significant quantitative differences between the predictions of the two models. Interestingly, implicit solvent simulations with Dry Martini show good quantitative agreement with all-atomic MD simulations, better than explicit solvent Martini MD simulations. The computational efficiency of the Martini and Dry Martini models allowed us to study the self-assembly of CTAB in a large system with many micelles. We observe the self-assembly of CTAB into micelles and also the exchange of CTAB molecules between micelles by events such as micelle fusion and fission which are difficult to observe in all-atomic MD simulations due to the time and length scales involved. Under the studied conditions, both Martini models predict a rather different self-assembly behaviour. The standard Martini model predicts a final equilibrium state with spherical micelles with an average size of ≈70 CTAB molecules. In contrast, the Dry Martini model predicts the formation of large tubular micelles with ≈330 CTAB molecules. Compared with experiments, standard Martini and Dry Martini underestimate and overestimate, respectively, the micelle size.

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.

Behavior of ligand binding assays with crowded surfaces: Molecular model of antigen capture by antibody-conjugated nanoparticles.

Ligand-receptor binding is of utmost importance in several biologically related disciplines. Ligand binding assays (LBA) use the high specificity and high affinity of ligands to detect, target or measure a specific receptors. One particular example of ligand binding assays are Antibody conjugated Nanoparticles (AcNPs), edge-cutting technologies that are present in several novel biomedical approaches for imaging, detection and treatment of diseases. However, the nano-confinement in AcNPs and LBA nanostructures introduces extra complexity in the analysis of ligand-receptor equilibriums. Because antibodies are large voluminous ligands, the effective affinity in AcNPs is often determined by antibody orientation and surface coverage. Moreover, antibodies have two binding sites introducing an extra ligand-receptor binding equilibrium. As consequence of all this, experimental or theoretical studies providing a guidelines for the prediction of the binding behavior in AcNPs are scarce. In this work, we present a set of theoretical calculations to shed light into the complex binding behavior of AcNPs and its implications in biomedical applications. To investigate the ligand-receptor binding on AcNPs, we have used a molecular theory that predicts the probability of different molecular conformations of the system depending on the local environment. We have considered two different pathways for designing these devices: covalently conjugated antibodies and streptavidin-biotin conjugated antibodies. We also explore the effects of surface coverage, bulk concentrations, nanoparticle size and antibody-antigen affinity. Overall, this work offers a series of theoretical predictions that can be used as a guide in the design of antibody conjugated nanoparticles for different applications.

Mechanical properties of a collagen fibril under simulated degradation.

Collagen fibrils are a very important component in most of the connective tissue in humans. An important process associated with several physiological and pathological states is the degradation of collagen. Collagen degradation is usually mediated by enzymatic and non-enzymatic processes. In this work we use molecular dynamics simulations to study the influence of simulated degradation on the mechanical properties of the collagen fibril. We applied tensile stress to the collagen fiber at different stages of degradation. We compared the difference in the fibril mechanical priorities due the removal of enzymatic crosslink, surface degradation and volumetric degradation. As anticipated, our results indicated that, regardless of the degradation scenario, fibril mechanical properties is reduced. The type of degradation mechanism (crosslink, surface or volumetric) expressed differential effect on the change in the fibril stiffness. Our simulation results showed dramatic change in the fibril stiffness with a small amount of degradation. This suggests that the hierarchical structure of the fibril is a key component for the toughness and is very sensitive to changes in the organization of the fibril. The overall results are intended to provide a theoretical framework for the understanding the mechanical behavior of collagen fibrils under degradation.

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.

Controlling the hydration rate of a hydrophilic matrix in the core of an intravaginal ring determines antiretroviral release.

Intravaginal ring technology is generally limited to releasing low molecular weight species that can diffuse through the ring elastomer. To increase the diversity of drugs that can be delivered from intravaginal rings, we designed an IVR that contains a drug matrix encapsulated in the core of the IVR whereby the mechanism of drug release is uncoupled from the interaction of the drug with the ring elastomer. We call the device a flux controlled pump, and it is comprised of compressed pellets of a mixture of drug and hydroxypropyl cellulose within the hollow core of the ring. The pump orifice size and chemistry of the polymer pellets control the rate of hydration and diffusion of the drug-containing hydroxypropyl cellulose gel from the device. A mechanistic model describing the hydration and diffusion of the hydroxypropyl cellulose matrix is presented. Good agreement between the quantitative model predictions and the experimental studies of drug release was obtained. We achieved controlled release rates of multiple antiretrovirals ranging from µg/d to mg/d by altering the orifice design, drug loading, and mass of pellets loaded in the device. This device could provide an adaptable platform for the vaginal drug delivery of many molecules.

Cooperative dynamic and diffusion behavior above and below the dynamical crossover of supercooled water.

Using extensive molecular dynamics simulations combined with a novel approach to analyze the molecular displacements we analyzed the change in the dynamics above and below the crossover temperature T(x) for supercooled water. Our findings suggest that the crossover from fragile to strong glass former occurring at T(x) is related with a change in the diffusion mechanism evidencing the presence of jump-like diffusion at lower temperatures. Also we observe that fluctuations of the local environments are intimately connected with fluctuations in the size and the amount of cooperative cluster of mobile molecules, and in particular we find a highly cooperative nature of the motion at low temperatures.

The water supercooled regime as described by four common water models.

The temperature scale of simple water models in general does not coincide with the natural one. Therefore, in order to make a meaningful evaluation of different water models, a temperature rescaling is necessary. In this paper, we introduce a rescaling using the melting temperature and the temperature corresponding to the maximum of the heat capacity to evaluate four common water models (TIP4P-Ew, TIP4P-2005, TIP5P-Ew and Six-Sites) in the supercooled regime. Although all the models show the same general qualitative behavior, the TIP5P-Ew appears as the best representation of the supercooled regime when the rescaled temperature is used. We also analyze, using thermodynamic arguments, the critical nucleus size for ice growth. Finally, we speculate on the possible reasons why atomistic models do not usually crystalize while the coarse grained mW model do crystallize.