Interaction mechanism between the focused ultrasound and lipid membrane at the molecular level

Ultrasound pulse repetition frequency preferentially activates different neuron populations independent of cell type

Objective. Transcranial ultrasound (US) stimulation serves as an external input to a neuron, and thus the evoked response relies on neurons' intrinsic properties. Neural activity is limited to a couple hundred hertz and often exhibits preference to input frequencies. Accordingly, US pulsed at specific physiologic pulse repetition frequencies (PRFs) may selectively engage neurons with the corresponding input frequency preference. However, most US parametric studies examine the effects of supraphysiologic PRFs. It remains unclear whether pulsing US at different physiologic PRFs could activate distinct neurons in the awake mammalian brain.Approach. We recorded cellular calcium responses of individual motor cortex neurons to US pulsed at PRFs of 10, 40, and 140 Hz in awake mice. We compared the evoked responses across these PRFs in the same neurons. To further understand the cell-type dependent effects, we categorized the recorded neurons as parvalbumin positive fast spiking interneurons or putative excitatory neurons and analyzed single-cell mechanosensitive channel expression in mice and humans using the Allen Brain Institute's RNA-sequencing databases.Main results. We discovered that many neurons were preferentially activated by only one PRF and different PRFs selectively engaged distinct neuronal populations. US-evoked cellular calcium responses exhibited the same characteristics as those naturally occurring during spiking, suggesting that US increases intrinsic neuronal activity. Furthermore, evoked responses were similar between fast-spiking inhibitory neurons and putative excitatory neurons. Thus, variation in individual neuron's cellular properties dominates US-evoked response heterogeneity, consistent with our observed cell-type independent expression patterns of mechanosensitive channels across individual neurons in mice and humans. Finally, US transiently increased network synchrony without producing prolonged over-synchronization that could be detrimental to neural circuit functions.Significance. These results highlight the feasibility of activating distinct neuronal subgroups by varying PRF and the potential to improve neuromodulation effects by combining physiologic PRFs.

Ultrasound enhanced diffusion in hydrogels: An experimental and non-equilibrium molecular dynamics study

Focused ultrasound has experimentally been found to enhance the diffusion of nanoparticles; our aim with this work is to study this effect closer using both experiments and non-equilibrium molecular dynamics. Measurements from single particle tracking of 40 nm polystyrene nanoparticles in an agarose hydrogel with and without focused ultrasound are presented and compared with a previous experimental study using 100 nm polystyrene nanoparticles. In both cases, we observed an increase in the mean square displacement during focused ultrasound treatment. We developed a coarse-grained non-equilibrium molecular dynamics model with an implicit solvent to investigate the increase in the mean square displacement and its frequency and amplitude dependencies. This model consists of polymer fibers and two sizes of nanoparticles, and the effect of the focused ultrasound was modeled as an external oscillating force field. A comparison between the simulation and experimental results shows similar mean square displacement trends, suggesting that the particle velocity is a significant contributor to the observed ultrasound-enhanced mean square displacement. The resulting diffusion coefficients from the model are compared to the diffusion equation for a two-time continuous time random walk. The model is found to have the same frequency dependency. At lower particle velocity amplitude values, the model has a quadratic relation with the particle velocity amplitude as described by the two-time continuous time random walk derived diffusion equation, but at higher amplitudes, the model deviates, and its diffusion coefficient reaches the non-hindered diffusion coefficient. This observation suggests that at higher ultrasound intensities in hydrogels, the non-hindered diffusion coefficient can be used.

Ultrasound pulse repetition frequency preferentially activates different neuron populations independent of cell type

Transcranial ultrasound activates mechanosensitive cellular signaling and modulates neural dynamics. Given that intrinsic neuronal activity is limited to a couple hundred hertz and often exhibits frequency preference, we examined whether pulsing ultrasound at physiologic pulse repetition frequencies (PRFs) could selectively influence neuronal activity in the mammalian brain. We performed calcium imaging of individual motor cortex neurons, while delivering 0.35 MHz ultrasound at PRFs of 10, 40, and 140 Hz in awake mice. We found that most neurons were preferentially activated by only one of the three PRFs, highlighting unique cellular effects of physiologic PRFs. Further, ultrasound evoked responses were similar between excitatory neurons and parvalbumin positive interneurons regardless of PRFs, indicating that individual cell sensitivity dominates ultrasound-evoked effects, consistent with the heterogeneous mechanosensitive channel expression we found across single neurons in mice and humans. These results highlight the feasibility of tuning ultrasound neuromodulation effects through varying PRFs.

Molecular insights into the effects of focused ultrasound mechanotherapy on lipid bilayers: Unlocking the keys to design effective treatments

Administration of focused ultrasounds (US) represents an attractive complement to classical therapies for a wide range of maladies, from cancer to neurological pathologies, as they are non-invasive, easily targeted, their dosage is easy to control, and they involve low risks. Different mechanisms have been proposed for their activity but the direct effect of their interaction with cell membranes is not well understood at the molecular level. This is in part due to the difficulty of designing experiments able to probe the required spatio-temporal resolutions. Here we use Molecular Dynamics (MD) simulations at two resolution levels and machine learning (ML) classification tools to shed light on the effects that focused US mechanotherapy methods have over a range of lipid bilayers. Our results indicate that the dynamic-structural response of the membrane models to the mechanical perturbations caused by the sound waves strongly depends on the lipid composition. The analyses performed on the MD trajectories contribute to a better understanding of the behavior of lipid membranes, and to open up a path for the rational design of new therapies for the long list of diseases characterized by specific lipid profiles of pathological membrane cells.

Computational Insight of Phase Transformation and Drug Release Behaviour of Doxycycline-Loaded Ibuprofen-Based In-Situ Forming Gel

This research investigates the gel formation behaviour and drug-controlling performance of doxycycline-loaded ibuprofen-based in-situ forming gels (DH-loaded IBU-based ISGs) for potential applications in periodontal treatment. The investigation begins by exploring the physical properties and gel formation behaviour of the ISGs, with a particular focus on determining their sustained release capabilities. To gain a deeper understanding of the molecular interactions and dynamics within the ISGs, molecular dynamic (MD) simulations are employed. The effects of adding IBU and DH on reducing surface tension and water tolerance properties, thus affecting molecular properties. The phase transformation phenomenon is observed around the interface, where droplets of ISGs move out to the water phase, leading to the precipitation of IBU around the interface. The optimization of drug release profiles ensures sustained local drug release over seven days, with a burst release observed on the first day. Interestingly, different organic solvents show varying abilities to control DH release, with dimethyl sulfoxide (DMSO) demonstrating superior control compared to N-Methyl-2-pyrrolidone (NMP). MD simulations using AMBER20 software provide valuable insights into the movement of individual molecules, as evidenced by root-mean-square deviation (RMSD) values. The addition of IBU to the system results in the retardation of IBU molecule movement, particularly evident in the DMSO series, with the diffusion constant value of DH reducing from 1.2452 to 0.3372 and in the NMP series from 0.3703 to 0.2245 after adding IBU. The RMSD values indicate a reduction in molecule fluctuation of DH, especially in the DMSO system, where it decreases from over 140 to 40 Å. Moreover, their radius of gyration is influenced by IBU, with the DMSO system showing lower values, suggesting an increase in molecular compactness. Notably, the DH-IBU configuration exhibits stable pairing through H-bonding, with a higher amount of H-bonding observed in the DMSO system, which is correlated with the drug retardation efficacy. These significant findings pave the way for the development of phase transformation mechanistic studies and offer new avenues for future design and optimization formulation in the ISG drug delivery systems field.

Molecular dynamics simulation of cancer cell membrane perforated by shockwave induced bubble collapse

It has been widely accepted that cancer cells are softer than their normal counterparts. This motivates us to propose, as a proof-of-concept, a method for the efficient delivery of therapeutic agents into cancer cells, while normal cells are less affected. The basic idea of this method is to use a water jet generated by the collapse of the bubble under shockwaves to perforate pores in the cell membrane. Given a combination of shockwave and bubble parameters, the cancer membrane is more susceptible to bending, stretching, and perforating than the normal membrane because the bending modulus of the cancer cell membrane is smaller than that of the normal cell membrane. Therefore, the therapeutic agent delivery into cancer cells is easier than in normal cells. Adopting two well-studied models of the normal and cancer membranes, we perform shockwave induced bubble collapse molecular dynamics simulations to investigate the difference in the response of two membranes over a range of shockwave impulse 15-30 mPa s and bubble diameter 4-10 nm. The simulation shows that the presence of bubbles is essential for generating a water jet, which is required for perforation; otherwise, pores are not formed. Given a set of shockwave impulse and bubble parameters, the pore area in the cancer membrane is always larger than that in the normal membrane. However, a too strong shockwave and/or too large bubble results in too fast disruption of membranes, and pore areas are similar between two membrane types. The pore closure time in the cancer membrane is slower than that in the normal membrane. The implications of our results for applications in real cells are discussed in some details. Our simulation may be useful for encouraging future experimental work on novel approaches for cancer treatment.

Engineering of 2D nanomaterials to trap and kill SARS-CoV-2: a new insight from multi-microsecond atomistic simulations

In late 2019, coronavirus disease 2019 (COVID-19) was caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Spike protein is one of the surface proteins of SARS-CoV-2 that is essential for its infectious function. Therefore, it received lots of attention for the preparation of antiviral drugs, vaccines, and diagnostic tools. In the current study, we use computational methods of chemistry and biology to study the interaction between spike protein and its receptor in the body, angiotensin-I-converting enzyme-2 (ACE2). Additionally, the possible interaction of two-dimensional (2D) nanomaterials, including graphene, bismuthene, phosphorene, p-doped graphene, and functionalized p-doped graphene, with spike protein is investigated. The functionalized p-doped graphene nanomaterials were found to interfere with spike protein better than the other tested nanomaterials. In addition, the interaction of the proposed nanomaterials with the main protease (Mpro) of SARS-CoV-2 was studied. Functionalized p-doped graphene nanomaterials showed more capacity to prevent the activity of Mpro. These 2D nanomaterials efficiently reduce the transmissibility and infectivity of SARS-CoV-2 by both the deformation of the spike protein and inhibiting the Mpro. The results suggest the potential use of 2D nanomaterials in a variety of prophylactic approaches, such as masks or surface coatings, and would deserve further studies in the coming years.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Molecular Study of Ultrasound-Triggered Release of Fluorescein from Liposomes

Several investigations have suggested that ultrasound triggers the release of drugs encapsulated into liposomes at acoustic pressures low enough to avoid cavitation or high hyperthermia. However, the mechanism leading to this triggered release as well as the adequate composition of the liposome membrane remains unknown. Here, we investigate the ultrasound-triggered release of fluorescein disodium salt encapsulated into liposomes made of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) or 1,2-distearoylphosphatidyl-ethanolamine (DSPC) lipids with various concentrations of cholesterol (from 0 to 44 mol %). The passive release of encapsulated fluorescein was first characterized. It was observed to be higher when the membrane is in a fluid phase and increased with temperature but decreased upon addition of cholesterol. Next, the release of fluorescein was measured at different acoustic frequencies (0.8, 1.1, and 3.3 MHz) and peak-to-peak pressures (0, 2, 2.5, 5, and 8 MPa). Measurements were performed at temperatures where DOPC and DSPC liposomes were, respectively, in the fluid or gel phase. We found that the release rate did not depend on the ultrasound frequency. For DOPC liposomes, the ultrasound-triggered release of fluorescein decreased with increasing concentration of cholesterol in liposomes, while the behavior was more complex for DSPC liposomes. Overall, the triggered release from DSPC liposomes was up to ten times less than DOPC liposomes. Molecular dynamics simulations performed on a pure DOPC membrane showed that a membrane experiences, under a directional pressure of ±2.4 MPa, various changes in properties such as the area per lipid (APL). An increase in the APL was notably observed when the simulation box was laterally stretched or perpendicularly compressed, which was accompanied by an increase in the number of water molecules crossing the membrane. This suggests that ultrasound most probably enhances the diffusion of encapsulated molecules at small acoustic pressures by increasing the distance between lipids.

Molecular mechanism of ultrasound interaction with a blood brain barrier model

The brain is strictly protected by the blood brain barrier preventing the crossing of therapeutics to treat brain diseases. The high and low intensity focused ultrasound methods have been used to temporarily open the blood brain barrier, facilitating the transport of drugs. The methods are very promising because the opening is transient, localized, and noninvasive. However, the molecular mechanism of the opening is unknown, and this limits the development and application of these methods. With this in mind, we carry out a molecular dynamics simulation study to understand the interaction of ultrasound with the cell membrane and the tight junction. Our minimal blood brain barrier model is composed of two lipid bilayers, mimicking two portions of neighboring cells, connected together by a tight junction formed by a pair of two cis-dimers of the claudin-5 protein. Using an experimental ultrasound frequency of 50 MHz, simulations show that at low intensities, ultrasound does not impact the structure of the cell membranes and tight junction, implying that the direct interaction of ultrasound with the blood brain barrier is not responsible for the experimentally observed opening. At high intensities, the ultrasound pulls the monolayers of individual cell membrane lipid bilayers apart, creating air compartments inside the bilayers. This reduces the free energy barrier for the translocation of drugs across the lipid bilayer and enhances drug permeability. At very high intensities, the two monolayers are largely separated, resulting in cell damage and implying that the blood brain barrier is primarily opened at the experimentally observed damaged areas.