Molecular dynamics simulations of doxorubicin in sphingomyelin-based lipid membranes

Molecular Dynamics Simulations Reveal Octanoylated Hyaluronic Acid Enhances Liposome Stability, Stealth and Targeting

Liposome-based drug delivery systems have been widely used in drug and gene delivery. However, issues such as instability, immune clearance, and poor targeting have significantly limited their clinical utility. Consequently, there is an urgent need for innovative strategies to improve liposome performance. In this study, we explore the interaction mechanisms of hyaluronic acid (HA), a linear anionic polysaccharide composed of repeating disaccharide units of d-glucuronic acid and N-acetyl-d-glucosamine connected by alternating β-1,3 and β-1,4 glycosidic linkages, and its octanoylated derivates (OHA) with liposomes using extensive coarse-grained molecular dynamics simulations. The octyl moieties of OHA spontaneously inserted into the phospholipid bilayer of liposomes, leading to their effective coating onto the surface of liposome and enhancing their structural stability. Furthermore, encapsulating liposome with OHA neutralized their surface potential, interfering with the formation of a protein corona known to contribute to liposomal immune clearance. Importantly, the encapsulated OHA maintained its selectivity and therefore targeting ability for CD44, which is often overexpressed in tumor cells. These molecular-scale findings shed light on the interaction mechanisms between HA and liposomes and will be useful for the development of next-generation liposome-based drug delivery systems.

Investigating the Influence of Photoswitchable Lipids on the Structure and Dynamics of Lipid Membranes: Fundamentals and Potential Applications

In this work, we delve into the impact of photoisomerization of photoswitchable lipids (PSLs) on the membrane structure and dynamics at a molecular level. Through all-atom molecular dynamics simulations, we explore how UV irradiation-induced trans-to-cis isomerization of these lipids, particularly the azobenzene-derivatized phosphatidylcholine (AzoPC) lipid, influences the structure and dynamics of a simplified lipid membrane, mimicking those of E. coli bacteria across different temperatures. Our findings align with previous experimental observations regarding membrane properties and offer insights into localized effects and microscopic heterogeneity. Additionally, we estimate the relaxation time scale of the lipid membrane following AzoPC photoisomerization. Moreover, we demonstrate the feasibility of photoactivated drug release, exemplified by the controlled liberation of doxorubicin, an anticancer agent, through the membrane, suggesting the potential of PSLs in engineering photoactivated liposomes, coined as photoazosomes, for precise targeted drug delivery applications.

Unravelling the interactions between small molecules and liposomal bilayers via molecular dynamics and thermodynamic modelling

Lipid-based drug delivery systems hold immense promise in addressing critical medical needs, from cancer and neurodegenerative diseases to infectious diseases. By encapsulating active pharmaceutical ingredients - ranging from small molecule drugs to proteins and nucleic acids - these nanocarriers enhance treatment efficacy and safety. However, their commercial success faces hurdles, such as the lack of a systematic design approach and the issues related to scalability and reproducibility. This work aims to provide insights into the drug-phospholipid interaction by combining molecular dynamic simulations and thermodynamic modelling techniques. In particular, we have made a connection between the structural properties of the drug-phospholipid system and the physicochemical performance of the drug-loaded liposomal nanoformulations. We have considered two prototypical drugs, felodipine (FEL) and naproxen (NPX), and one model hydrogenated soy phosphatidylcholine (HSPC) bilayer membrane. Molecular dynamic simulations revealed which regions within the phospholipid bilayers are most and least favoured by the drug molecules. NPX tends to reside at the water-phospholipid interface and is characterized by a lower free energy barrier for bilayer membrane permeation. Meanwhile, FEL prefers to sit within the hydrophobic tails of the phospholipids and is characterized by a higher free energy barrier for membrane permeation. Flory-Huggins thermodynamic modelling, small angle X-ray scattering, dynamic light scattering, TEM, and drug release studies of these liposomal nanoformulations confirmed this drug-phospholipid structural difference. The naproxen-phospholipid system has a lower free energy barrier for permeation, higher drug miscibility with the bilayer, larger liposomal nanoparticle size, and faster drug release in the aqueous medium than felodipine. We suggest that this combination of molecular dynamics and thermodynamics approach may offer a new tool for designing and developing lipid-based nanocarriers for unmet medical applications.

Insights into Transfer of Supramolecular Doxorubicin/Congo Red Aggregates through Phospholipid Membranes

Doxorubicin (DOX) is a commonly used chemotherapeutic drug, from the anthracycline class, which is genotoxic to neoplastic cells via a DNA intercalation mechanism. It is effective and universal; however, it also causes numerous side effects. The most serious of them are cardiotoxicity and a decrease in the number of myeloid cells. For this reason, targeted DOX delivery systems are desirable, since they would allow lowering the drug dose and therefore limiting systemic side effects. Recently, synthetic dyes, in particular Congo red (CR), have been proposed as possible DOX carriers. CR is a planar molecule, built of a central biphenyl moiety and two substituted naphthalene rings, connected with diazo bonds. In water, it forms elongated ribbon-shaped supramolecular structures, which are able to selectively interact with immune complexes. In our previous studies, we have shown that CR aggregates can intercalate DOX molecules. In this way, they preclude DOX precipitation in water solutions and increase its uptake by MCF7 breast cancer cells. In the present work, we further explore the interactions between DOX, CR, and their aggregates (CR/DOX) with phospholipid membranes. In addition to neutral molecules, the protonated doxorubicin form, DXP, is also studied. Molecular dynamics simulations are employed to study the transfer of CR, DOX, DXP, and their aggregates through POPC bilayers. Interactions of CR, DOX, and CR/DOX with model monolayers are studied with Langmuir trough measurements. This study shows that CR may support the transfer of doxorubicin molecules into the bilayer. Both electrostatic and van der Waals interactions with lipids are important in this respect. The former promote the initial stages of the insertion process, the latter keep guest molecules inside the bilayer.

Dynamic Simulations of Interaction of the PEG-DPPE Micelle-Encapsulated Short-Chain Ceramides with the Raft-Included Membrane

The lipid raft subdomains in cancer cell membranes play a key role in signal transduction, biomolecule recruitment, and drug transmembrane transport. Augmented membrane rigidity due to the formation of a lipid raft is unfavorable for the entry of drugs, a limiting factor in clinical oncology. The short-chain ceramide (CER) has been reported to promote drug entry into membranes and disrupt lipid raft formation, but the underlying mechanism is not well understood. We recently explored the carrier-membrane fusion dynamics of PEG-DPPE micelles in delivering doxorubicin (DOX). Based on the phase-segregated membrane model composed of DPPC/DIPC/CHOL/GM1/PIP2, we aim to explore the dynamic mechanism of the PEG-DPPE micelle-encapsulating DOXs in association with the raft-included cell membrane modulated by C8 acyl tail CERs. The results show that the lipid raft remains integrated and DOX-resistant subjected to free DOXs and the micelle-encapsulating ones. Addition of CERs disorganizes the lipid raft by pushing CHOL aside from DPPC. It subsequently allows for a good permeability for PEG-DPPE micelle-encapsulated DOXs, which penetrate deeper as CER concentration increases. GM1 is significant in guiding drugs' redistributing between bilayer phases, and the anionic PIP2 further helps DOXs attain the inner bilayer surface. These results elaborate on the perturbing effect of CERs on lipid raft stability, which provides a new comprehensive approach for further design of drug delivery systems.

The effect of polymer coating on nanoparticles' interaction with lipid membranes studied by coarse-grained molecular dynamics simulations

Nanoparticles' (NPs) permeation through cell membranes, whether it happens via passive or active transport, is an essential initial step for their cellular internalization. The NPs' surface coating impacts the way they translocate through the lipid bilayer and the spontaneity of the process. Understanding the molecular details of NPs' interaction with cell membranes allows the design of nanosystems with optimal characteristics for crossing the lipid bilayer: computer simulations are a powerful tool for this purpose. In this work, we have performed coarse-grained molecular dynamics simulations and free energy calculations on spherical titanium dioxide NPs conjugated with polymer chains of different chemical compositions. We have demonstrated that the hydrophobic/hydrophilic character of the chains, more than the nature of their terminal group, plays a crucial role in determining the NPs' interaction with the lipid bilayer and the thermodynamic spontaneity of NPs' translocation from water to the membrane. We envision that this computational work will be helpful to the experimental community in terms of the rational design of NPs for efficient cell membrane permeation.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

A comparative study of 5- fluorouracil, doxorubicin, methotrexate, paclitaxel for their inhibition ability for Mpro of nCoV: Molecular docking and molecular dynamics simulations

A new corona virus (nCoV) is aetiological agent responsible for the viral pneumonia epidemic. Three is no specific therapeutic medicines available for the treatment of this condition and also effective treatment choices are few. In this work author tried to investigate some repurposing drug such as 5- fluorouracil, doxorubicin, methotrexate and paclitaxel against the main protease (Mpro) of nCoV by the computational model. Molecular docking was performed to screen out the best compound and doxorubicin was found to have minimum binding energy −121.89 kcal/mol. To further study, MD simulations were performed at 300 K and the result successfully corroborate the energy obtained by molecular docking. Temperature dependent MD simulation of the best molecule that is doxorubicin obtained from docking result was performed to check the variation in structural changes in Mpro of nCoV at 290 K, 310 K, 320 K and 325 K. It is sound that doxorubicin binds effectively with Mpro of nCoV at 290 K. Further ADME properties of the 5- fluorouracil, doxorubicin, methotrexate and paclitaxel were also evaluated to understand the bioavailability.

Multi-scale modeling of folic acid-functionalized TiO2 nanoparticles for active targeting of tumor cells

Strategies based on the active targeting of tumor cells are emerging as smart and efficient nanomedical procedures. Folic acid (FA) is a vitamin and a well-established tumor targeting agent because of its strong affinity for the folate receptor (FR), which is an overexpressed protein on the cell membranes of the tumor cells. FA can be successfully anchored to several nanocarriers, including inorganic nanoparticles (NPs) based on transition metal oxides. Among them, TiO₂ is extremely interesting because of its excellent photoabsorption and photocatalytic properties, which can be exploited in photodynamic therapy. However, it is not yet clear in which respects direct anchoring of FA to the NP or the use of spacers, based on polyethylene glycol (PEG) chains, are different and whether one approach is better than the other. In this work, we combine Quantum Mechanics (QM) and classical Molecular Dynamics (MD) to design and optimize the FA functionalization on bare and PEGylated TiO₂ models and to study the dynamical behavior of the resulting nanoconjugates in a pure water environment and in physiological conditions. We observe that they are chemically stable, even under the effect of increasing temperature (up to 500 K). Using the results from long MD simulations (100 ns) and from free energy calculations, we determine how the density of FA molecules on the TiO₂ NP and the presence of PEG spacers impact on the actual exposure of the ligands, especially by affecting the extent of FA-FA intermolecular interactions, which are detrimental for the targeting ability of FA towards the folate receptor. This analysis provides a solid and rational basis for experimentalists to define the optimal FA density and the more appropriate mode of anchoring to the carrier, according to the final purpose of the nanoconjugate.

An in Silico Approach to Reveal the Nanodisc Formulation of Doxorubicin

Molecular dynamic behaviors of nanodisc (ND) formulations of free doxorubicin (DOX) and DOX conjugated lipid prodrug molecules were investigated by molecular dynamics (MD) simulations. We have unveiled how formulation design affects the drug release profile and conformational stability of ND assemblies. Our simulation results indicate that free DOX molecules loaded in the ND system experienced rapid dissociation due to the unfavorable orientation of DOX attached to the lipid surface. It is found that DOX tends to form aggregates with higher drug quantities. In contrast, lipidated DOX-prodrugs incorporated in ND formulations exhibited sufficient ND conformational stability. The drug loading capacity is dependent on the type of lipid molecules grafted on the DOX-prodrug, and the drug loading quantities in a fixed area of NDs follow the order: DOX-BMPH-MP > DOX-BMPH-TC > DOX-BMPH-PTE. To gain further insight into the dynamic characteristics of ND formulations governed by different kinds of lipidation, we investigated the conformational variation of ND components, intermolecular interactions, the solvent accessible surface area, and individual MSP1 residue flexibility. We found that the global conformational stability of DOX-prodrug-loaded ND assemblies is influenced by the molecular flexibility and lipidated forms of DOX-prodrug. We also found that the spontaneous self-aggregation of DOX-prodrugs with increasing quantities on ND could reduce the membrane fluidity and enhance the conformational stability of ND formulations.