Improving the accuracy and convergence of drug permeation simulations via machine-learned collective variables.

Understanding the permeation of biomolecules through cellular membranes is critical for many biotechnological applications, including targeted drug delivery, pathogen detection, and the development of new antibiotics. To this end, computer simulations are routinely used to probe the underlying mechanisms of membrane permeation. Despite great progress and continued development, permeation simulations of realistic systems (e.g., more complex drug molecules or biologics through heterogeneous membranes) remain extremely challenging if not intractable. In this work, we combine molecular dynamics simulations with transition-tempered metadynamics and techniques from the variational approach to conformational dynamics to study the permeation mechanism of a drug molecule, trimethoprim, through a multicomponent membrane. We show that collective variables (CVs) obtained from an unsupervised machine learning algorithm called time-structure based Independent Component Analysis (tICA) improve performance and substantially accelerate convergence of permeation potential of mean force (PMF) calculations. The addition of cholesterol to the lipid bilayer is shown to increase both the width and height of the free energy barrier due to a condensing effect (lower area per lipid) and increase bilayer thickness. Additionally, the tICA CVs reveal a subtle effect of cholesterol increasing the resistance to permeation in the lipid head group region, which is not observed when canonical CVs are used. We conclude that the use of tICA CVs can enable more efficient PMF calculations with additional insight into the permeation mechanism.

Stressed Lipid Droplets: How Neutral Lipids Relieve Surface Tension and Membrane Expansion Drives Protein Association.

Lipid droplets (LDs) are intracellular storage organelles composed of neutral lipids, such as triacylglycerol (TG), surrounded by a phospholipid (PL) monolayer decorated with specific proteins. Herein, we investigate the mechanism of protein association during LD and bilayer membrane expansion. We find that the neutral lipids play a dynamic role in LD expansion by further intercalating with the PL monolayer to create more surface-oriented TG molecules (SURF-TG). This interplay both reduces high surface tension incurred during LD budding or growth and also creates expansion-specific surface features for protein recognition. We then show that the autoinhibitory (AI) helix of CTP:phosphocholine cytidylyltransferase, a protein known to target expanding monolayers and bilayers, preferentially associates with large packing defects in a sequence-specific manner. Despite the presence of three phenylalanines, the initial binding with bilayers is predominantly mediated by the sole tryptophan due to its preference for membrane interfaces. Subsequent association is dependent on the availability of large, neighboring defects that can accommodate the phenylalanines, which are more probable in the stressed systems. Tryptophan, once fully associated, preferentially interacts with the glycerol moiety of SURF-TG in LDs. The calculation of AI binding free energy, hydrogen bonding and depth analysis, and in silico mutation experiments support the findings. Hence, SURF-TG can both reduce surface tension and mediate protein association, facilitating class II protein recruitment during LD expansion.

Small molecule therapeutics for COVID-19: repurposing of inhaled furosemide.

The novel coronavirus SARS-CoV-2 has become a global health concern. The morbidity and mortality of the potentially lethal infection caused by this virus arise from the initial viral infection and the subsequent host inflammatory response. The latter may lead to excessive release of pro-inflammatory cytokines, IL-6 and IL-8, as well as TNF-α ultimately culminating in hypercytokinemia ("cytokine storm"). To address this immuno-inflammatory pathogenesis, multiple clinical trials have been proposed to evaluate anti-inflammatory biologic therapies targeting specific cytokines. However, despite the obvious clinical utility of such biologics, their specific applicability to COVID-19 has multiple drawbacks, including they target only one of the multiple cytokines involved in COVID-19's immunopathy. Therefore, we set out to identify a small molecule with broad-spectrum anti-inflammatory mechanism of action targeting multiple cytokines of innate immunity. In this study, a library of small molecules endogenous to the human body was assembled, subjected to in silico molecular docking simulations and a focused in vitro screen to identify anti-pro-inflammatory activity via interleukin inhibition. This has enabled us to identify the loop diuretic furosemide as a candidate molecule. To pre-clinically evaluate furosemide as a putative COVID-19 therapeutic, we studied its anti-inflammatory activity on RAW264.7, THP-1 and SIM-A9 cell lines stimulated by lipopolysaccharide (LPS). Upon treatment with furosemide, LPS-induced production of pro-inflammatory cytokines was reduced, indicating that furosemide suppresses the M1 polarization, including IL-6 and TNF-α release. In addition, we found that furosemide promotes the production of anti-inflammatory cytokine products (IL-1RA, arginase), indicating M2 polarization. Accordingly, we conclude that furosemide is a reasonably potent inhibitor of IL-6 and TNF-α that is also safe, inexpensive and well-studied. Our pre-clinical data suggest that it may be a candidate for repurposing as an inhaled therapy against COVID-19.

Effect of Cholesterol on the Structure of Networked Water at the Surface of a Model Lipid Membrane.

Using all-atom molecular dynamics simulations and network analysis, we investigated the effect of membrane cholesterol level on the structure of organized water at the interface between bulk water and a model lipid membrane. Irrespective of membrane cholesterol content, interfacial water structure is largely perturbed by the presence of the membrane surface due to water-phospholipid interactions, which deplete the chance of hydrogen bonding among water molecules. In contrast, the addition of cholesterol suppresses the disturbing effect of the membrane on water-water hydrogen bonding as cholesterol provides a more bulklike environment for the interfacial water molecules, as evidenced by enhancement of local water density, a reduction in their orientational bias, and increases in both the number of hydrogen bonds and the topological complexity of the hydrogen bond network.

Understanding the effect of nanoconfinement on the structure of water hydrogen bond networks.

Using an integrated approach of network theory and atomistic molecular dynamics simulations, we performed a detailed topological analysis on hydrogen bond networks of water confined between either two graphene sheets or two lipid bilayers to explore the structural perturbation of these networks under nanoscale confinement. The hydrogen bond network structure can be perturbed to a considerable extent when water is confined by such surfaces, yet no small-world behaviour is observed. The presence of ions also reduces the network complexity but its effect may be small depending on the type of confining surfaces. We developed an information flow model to evaluate the fluctuating nature of hydrogen bond networks and to characterise the dynamic, long-distance hydrogen-bonded chains in water. We found that the length and directionality of the hydrogen bond "trails" are highly susceptible to the type of confining surfaces and the degree of confinement. In particular, the endpoints of the hydrogen bond trails are not completely random in confined water, in which inherent distributions are determined by the density of water and the density of hydrogen bonds. This work forms the basis for the study of the pure effect of hydrogen bond network topology on various transport processes, such as proton transfer, that occur along a sequence of hydrogen bonds in a biochemical system. Our results suggest that a combined effect of the structure and lifetime of the hydrogen bond network of interfacial water may contribute to high lateral proton diffusivity at the surface of a lipid membrane.

Understanding Water Structure in an Ion-Pair Solvation Shell in the Vicinity of a Water/Membrane Interface.

The anomalous properties of interfacial water at the surface of a lipid membrane and their implications on nearby chemical processes are well recognized. However, we have found that ion pairing thermodynamics may not be significantly affected by interfacial water in a classical, nonpolarizable force field. To trace the root cause of such a counterintuitive finding, we performed atomistic molecular dynamics simulations to explore the impact of polarizable interactions and characterize the hydration structure of a sodium chloride (NaCl) ion pair at the surface of a model lipid membrane and in a bulk phase. Our study reveals that the effect of the aqueous interface on the first solvation shell of the ion pair and thus on the ion pairing thermodynamics becomes more pronounced in the polarizable model, and that the free energy profile along the interionic distance cannot capture the difference in the degree of solvent participation in ion pairing at the water/membrane interface. This study also forms the basis for the future design of a reaction coordinate that takes the behavior of the interfacial water into account.

Macroion-Solvent Interactions in Charged Droplets.

Macroion-droplet interactions play a critical role in many settings such as ionization techniques of samples in mass spectrometry analysis and atmospheric aerosols. The droplets under investigation are composed of a polar solvent, primarily water, a charged macroion, and, possibly, buffer ions. We present highlights of our research on the relation between the charge state of a macroion and the droplet morphologies. We have determined that, depending on the charge on the macroion and certain macroscopic properties of the solvent, such as its dielectric constant and surface tension, a droplet may obtain striking conformations such as droplets with extruded tails, "pearl-necklace" conformations, and multipoint "star" shapes. The shapes of the droplet containing the macroion influence the charging mechanism of the macroion in a reciprocal manner. Understanding of the macroion-droplet interactions plays a central role in explaining the origin and the magnitude of the charge in spectra obtained in electrospray ionization mass spectrometry experiments and in controlling the stability of complexes of nucleic acids and proteins in droplets.

What factors determine the stability of a weak protein-protein interaction in a charged aqueous droplet?

Maintaining the interface of a weak transient protein complex transferred from bulk solution to the gaseous state via evaporating droplets is a critical question in the detection of the complex association (dissociation) constant by using electrospray ionization mass spectrometry (ESI-MS). Here we explore the factors that may affect the stability of a protein-protein interaction (PPI) using atomistic molecular dynamics (MD) modelling of a complex of ubiquitin (Ub) and the ubiquitin-associated domain (UbA) (RCSB PDB code ) and a non-covalent complex of diubiquitin (RCSB PDB code ) in aqueous droplets. A general method is presented to determine the protonation states of the complexes we investigate in particular, and that of a protein in general, under various pH conditions that an evaporating droplet acquires due to its change in size. We find that the combination of high temperature and high charge states of the protein complexes may destabilize the interface by creating new interfaces instead of a direct rupture of the initial stable interface. We provide evidence that highly charged protein complexes are found in droplets that form conical extrusions of the solvent on the surface due to charge-induced instability. This distinct droplet morphology leads to a higher solvent evaporation rate that assists in transferring the complex in the gaseous state without dissociation. The conical solvent protrusions expose on the droplet surface certain amino acids that otherwise would be solvated in a droplet with the protein complex of low charge states. The new vapor-protein interface does not have a direct effect on the stability of the PPI. A common way in experiments to stabilize the protein complexes in droplets is to reduce the protonation state of the proteins. Here we find that weakly bound protein complexes even at high protonation states can be stabilized by the presence of a small number of counterions, without affecting the protonation state of the protein. Our findings may provide guiding principles in ESI-MS experiments to stabilize weak transient PPIs.

When droplets become stars: charged dielectric droplets beyond the Rayleigh limit.

When a nano-drop comprising a single spherical central ion and dielectric solvent is charged above a well-defined threshold, it acquires a stable star morphology. In contrast, conducting droplets, will undergo fission. Here we report combined atomistic molecular dynamics and continuum modelling study of star formation of droplets that contain a highly charged ion. We assume that in the continuum model the dielectric response is linear. In this linear continuum model, which is an extension of Rayleigh model, the energy of the drop is comprised of terms analogous to those in Rayleigh model, which are surface energy and electrostatic energy of dielectric droplet charged by a central point charge. We present the stability analysis of the continuum model to determine the threshold of instability. Indeed we find that the model accounts well for the onset of the instabilities. Molecular dynamics show that the number of points of the star-shaped nano-drops depends only on the surface tension, dielectric constant and size of the droplet, and on the magnitude of the charge of the central ion, but not on its sign. Intuitively, it is expected that when a spherical dielectric drop becomes unstable it would transform into a non-spherical finite shape of the same volume as the initial spherical shape with the point charge located in the drop interior. To test whether the extended Rayleigh model can account for the observed droplet shapes, we performed numerical simulations of the linear continuum model. Contrary to the expectations, the simulations of the extended Rayleigh model does not reproduce the stable star shapes found in the atomistic simulations, not even when we account for the bending rigidity and spontaneous curvature of the surface. We argue that the assumption that the dielectric response is linear breaks down if the droplet surface approaches the central macro-ion, where the electric field strength is such that dielectric saturation sets in. We envisage that for certain solvents, these stars could be made permanent by cross-linking, opening the way to the production of a novel class of highly-non-convex colloids.

Charging and Release Mechanisms of Flexible Macromolecules in Droplets.

We study systematically the charging and release mechanisms of a flexible macromolecule, modeled by poly(ethylene glycol) (PEG), in a droplet by using molecular dynamics simulations. We compare how PEG is solvated and charged by sodium Na+ ions in a droplet of water (H2O), acetonitrile (MeCN), and their mixtures. Initially, we examine the location and the conformation of the macromolecule in a droplet bearing no net charge. It is revealed that the presence of charge carriers do not affect the location of PEG in aqueous and MeCN droplets compared with that in the neutral droplets, but the location of the macromolecule and the droplet size do affect the PEG conformation. PEG is charged on the surface of a sodiated aqueous droplet that is found close to the Rayleigh limit. Its charging is coupled to the extrusion mechanism, where PEG segments leave the droplet once they coordinate a Na+ ion or in a correlated motion with Na+ ions. In contrast, as PEG resides in the interior of a MeCN droplet, it is sodiated inside the droplet. The compact macro-ion transitions through partially unwound states to an extended conformation, a process occurring during the final stage of desolvation and in the presence of only a handful of MeCN molecules. For charged H2O/MeCN droplets, the sodiation of PEG is determined by the H2O component, reflecting its slower evaporation and preference over MeCN for solvating Na+ ions. We use the simulation data to construct an analytical model that suggests that the droplet surface electric field may play a role in the macro-ion-droplet interactions that lead to the extrusion of the macro-ion. This study provides the first evidence of the effect of the surface electric field by using atomistic simulations. Graphical Abstract ᅟ.

Advances in Modeling the Stability of Noncovalent Complexes in Charged Droplets with Applications in Electrospray Ionization-MS Experiments.

Electrospray ionization mass spectrometry is used extensively to measure the equilibrium constant of noncovalent complexes. In this Perspective, we attempt to present an accessible introduction to computational methodologies that can be applied to determine the stability of weak noncovalent complexes in their journey from bulk solution into the gaseous state. We demonstrate the usage of the methods on two typical examples of noncovalent complexes drawn from a broad class of nucleic acids and transient protein complexes found in aqueous droplets. We conclude that this new emerging direction in the use of simulations can lead to estimates of equilibrium constant corrections due to complex dissociation in the carrier droplet and finding of agents that may stabilize the protein interfaces.

Stability of a Transient Protein Complex in a Charged Aqueous Droplet with Variable pH.

Electrospray ionization mass spectrometry (ESI-MS) has the potential to become a high-throughput robust experimental method for the detection of protein-protein equilibrium constants. Poorly understood processes that affect the stability of weak noncovalent protein complexes in the intervening droplet environment are a significant factor that precludes the advancement of the method. We use molecular dynamics to study the stability of a ubiquitin and ubiquitin-associated domain complex (RCSB PDB code 2MRO ) in an aqueous droplet with changing size and charge concentration. We present evidence that a weak protein complex changes conformation and may dissociate in shrinking droplets. Then, the droplets containing these dissociated proteins divide. Our findings suggest that in some cases ESI-MS does not measure the correct association constants. The study intends to stimulate research for systematic development of experimental protocols that stabilize weakly bound protein interfaces in droplets.

Effect of solvent on the charging mechanisms of poly(ethylene glycol) in droplets.

We examine the effect of solvent on the charging mechanisms of a macromolecule in a droplet by using molecular dynamics simulations. The droplet contains excess charge that is carried by sodium ions. To investigate the principles of the charging mechanisms of a macromolecule in a droplet, we simulate aqueous and methanol droplets that contain a poly(ethylene glycol) (PEG) molecule. We find that the solvent plays a critical role in the charging mechanism and in the manner that the sodiated PEG emerges from a droplet. In the aqueous droplets, the sodiated PEG is released from the droplet while it is being charged at a droplet charge state below the Rayleigh limit. The charging of PEG occurs on the surface of the droplet. In contrast to the aqueous droplets, in the methanol droplet, the sodiated PEG resides in the interior of the droplet and it may become charged at any location in the droplet, interior or surface. The sodiated PEG emerges from the droplet by drying-out of the solvent. Even though these two mechanisms appear to be phenomenologically similar to the widely accepted ion-evaporation and charge-residue mechanisms, they have fundamental differences from those. An integral part of the mechanism that the macromolecular ions emerge from droplets is the droplet morphology. Droplet morphologies give rise to different solvation interactions between the solvent and the macromolecule. In the water-sodiated PEG system, we find the extrusion of the PEG morphology, while in methanol-sodiated droplet, we find the "pearl-on-the-necklace" morphology and the extrusion of the sodiated PEG in the last stage of the desolvation process. These findings provide insight into the mechanisms that macromolecules acquire their charge in droplets produced in electrospray ionization experiments.