Molecular dynamics insight of interaction between the functionalized-carbon nanotube and cancerous cell membrane in doxorubicin delivery.

BACKGROUND AND OBJECTIVE: Doxorubicin (DOX) is a known anticancer drug which is widely used in cancer therapy. Carbon nanotubes (CNTs) are among the most promising platforms for smart drug delivery applications. However, due to the toxicity and their low sulubility their application is limited and their functionalization with wide range of biomolecules are suggested. Therefore, the functionalized carbon nanotubes (f-CNT) with carboxyl (CNT-COO) and folic acid (CNT-COO-FA) were investigated as drug-carrier. METHODS: Molecular dynamics (MD) simulation along with the Density Functional Theory (DFT) methods are being used to study the drug loading process on functionalized carbon nanotubes. RESULTS: The results indicate that doxorubicin molecules interact more with CNT-COO-FA than CNT-COO. The embedded dipalmitoylphosphatidylcholine (DPPC) lipid bilayer with a folate receptor was considered a cancerous cell's representative model. Then the drug release from the f-CNTs near the lipid bilayer was simulated. The results showed that CNT-COO-FA with a pH and ligand-sensitive mechanism strongly interacts with cancerous cells, which led to higher drug release, in agreement with the experimental results. The conformational changes of the lipid bilayer and folate receptor during drug release were evaluated. The analysis showed that drug release from CNT-COO-FA has significantly changed lipid bilayer and receptor conformations. The obtained results were interpreted and justified by considering the molecular mechanisms which control the drug delivery in the studied systems. CONCLUSIONS: Based on the obtained results, CNT-COO-FA has a better performance during the drug release compared to CNT-COO in delivering doxorubicin. Both pH and ligand sensitive mechanisms are found to be responsible for higher drug delivery efficiency of CNT-COO-FA. In contrast, CNT-COO can only enhance drug delivery efficiently with a pH-sensitive mechanism.

CellSys: An open-source tool for building initial structures for bio-membranes and drug-delivery systems.

Since phospholipids are the most important components in the structure of biomembranes, they deserve to be considered with a lot of attention in both experimental and computational theoretical studies using molecular simulation methods related to the research in the fields of drug design and drug delivery where they involve knowledge about the interactions of drug molecules with cell membranes. To employ the molecular simulation approach for this purpose the essential requirement is having information about the initial structure of phospholipids and how they interact with the drugs. Therefore in this article, we introduce an open-source software package in Python programming language for utilizing data manipulation for generation and developing the initial structure of biomolecular cells to provide the needed information for investigation in drug delivery systems. In addition, the proposed software package can be used for the efficient storage of membrane structural data to be exploited in designing new drug delivery systems. To verify the performance of the code and the results of the simulations, several analyses have been done, such as the calculation of area per lipid and self-diffusion coefficient, in addition to lipid order parameter. The results were in complete agreement with the references.

Molecular dynamics simulation study of doxorubicin adsorption on functionalized carbon nanotubes with folic acid and tryptophan.

In this work, molecular dynamics (MD) simulation is used to study the adsorption of the anticancer drug, doxorubicin (DOX), on the wall or surface of pristine and functionalized carbon nanotubes (FCNTs) in an aqueous solution. Initially, the CNTs were functionalized by tryptophan (Trp) and folic acid (FA), and then the DOX molecules were added to the system. The simulation results showed that the drug molecules can intensely interact with the FCNTs at physiological pH. Furthermore, it was found that as a result of functionalization, the solubility of FCNTs in an aqueous solution increases significantly. The effect of pH variation on drug release from both pristine and FCNTs was also investigated. The obtained results indicated that in acidic environments due to protonation of functional groups (Trp) and as a result of repulsive interaction between the DOX molecule and functional groups, the release of DOX molecules from FCNT's surface is facilitated. The drug release is also strongly dependent on the pH and protonated state of DOX and FCNT.

Molecular Dynamics Assessment of Doxorubicin Adsorption on Surface-Modified Boron Nitride Nanotubes (BNNTs).

Loading therapeutic agents on nanocarriers in order to protect them during drug delivery and exclusively targeting damaged tissues has gained substantial significance in biology realms in the past decade. Boron nitride nanotubes have given a new lease on designing nano delivery systems by virtue of their unique properties. The studies are still ongoing to thoroughly identify their chemical characteristics. In this study, we probed into the efficacy of boron nitride nanotubes and the impact of their surface modification by hydroxyl and amine functional groups in interaction with an anticancer drug model, i.e., doxorubicin. Defining the altered electronic properties of the nanotubes as well as the distribution of partial charges were carried out through density functional theory calculations, following the simulation of the drug loading process via molecular dynamics algorithms. The primary outcomes are inferred from systematical energies, van der Waals and electrostatic interactions, radial distribution functions, the number of hydrogen bonds, mean square displacement, diffusion coefficients, and binding free energies. Negative values of van der Waals energies imply a rapid, exothermic adsorption process whereby the contribution of these driving forces is more dominant than electrostatic ones. Ultimately, the values of overall diffusion coefficients of drugs and binding free energies, performed by the MM/PBSA approach, corroborate that the hydroxyl and amine-functionalized nanotubes reinforce the binding strength of the complexes to an approximate extent.

Synthesis, characterization, and CO2 adsorption properties of metal organic framework Fe-BDC.

The iron-containing Metal-Organic Frameworks (MOFs) have attracted a great deal of attention in the areas of gas separation, catalytic conversion, and drug delivery, due to their high surface area and activity, as well as the non-toxicity of iron. In this study, Fe-based MOFs using BDC ligands, MIL-101(Fe), MIL-53(Fe) and Amino-MIL-101(Fe) are synthesized by a solvothermal method and characterized by conventional methods such as BET, SEM, and TGA. Afterwards, the synthesized MOFs are investigated from the point of view of the adsorbing capability of carbon dioxide at different pressures and temperatures, and also their resistance to water and solvent. The results showed that Amino-MIL-101(Fe) achieved more CO2 adsorption than MIL-101(Fe) and MIL-53(Fe), equal to 13 mmol g-1 at 4 MP. Although MIL-53(Fe) has the best temperature resistance, around 350 °C, Amino-MIL-101(Fe) is more stable against water and ethanol and its surface area is increased from 670 to 915 m2 g-1 after washing with ethanol. The adsorption study reveals that CO2 is adsorbed not only by a physical adsorption mechanism, but also by chemisorption of acidic carbon dioxide by basic NH2 agent in the structure of Amino-MIL-101(Fe).

Effect of salt concentration on properties of mixed carbonate-based electrolyte for Li-ion batteries: a molecular dynamics simulation study.

In this work, a computational framework is proposed by utilizing molecular dynamics simulation to explore the existing relation between molecular structure and ionic conductivity of the electrolyte system [LiPF6+(EC+DMC 1:1)] consisting of a mixture of cyclic ethylene carbonate (EC) and acyclic dimethyl carbonate (DMC) solvents and lithium hexafluorophosphate (LiPF6) salt to propose as a novel mixed organic solvent-based electrolytes to promote the performance of lithium-ion batteries (LIBs). To acquire a clear understanding of the structural and transport properties of the designed electrolytes, quantum chemistry (QC) calculations and molecular dynamics (MD) simulation are used. In the first step, the accurate molecular structures of the studied electrolytes in addition to their corresponding atomic partial charges are evaluated. The MD simulations are performed at 330 K varying the LiPF6 concentration (0.5 M to 2.2 M). Analysis of the obtained results indicated that ionic diffusivity and conductivity of the electrolytes are dependent on the structure of solvated ions and lithium salt (LiPF6) concentration. It is found that the obtained MD simulation results are in reasonable agreement with experimental results. Graphical abstract A representation of dependence of transport properties of electrolyte system [LiPF6 +(EC+DMC 1:1)] as function of salt concentration to be used in Lithium-ion batteries (LIBs).

Adsorption and encapsulation of the drug doxorubicin on covalent functionalized carbon nanotubes: A scrutinized study by using molecular dynamics simulation and quantum mechanics calculation.

Adsorption of the drug doxorubicin (DOX) onto covalent functionalized carbon nanotubes (CNTs) as drug carriers was studied by employing molecular dynamics (MD) simulation. CNT was covalently functionalized by the chemical groups: amine, carboxyl and hydroxyl and the change in the electrostatic charge of CNT as a result of functionalization was investigated by quantum mechanics calculations. The drug adsorption onto the functionalized CNTs (f-CNT) was examined by analyzing the evaluated radial probability of the drug by MD simulation. Overall consideration of the results demonstrated that surface functionalization enhances the loading capacity of CNT for the drug encapsulation, also agglomeration of unprotonated drug molecules has increased encapsulation capacity. Analysis of the obtained results indicated that carboxyl and amine f-CNTs can act as a pH sensitive drug carrier where their protonation in acidic condition can decrease the electrostatic interactions of the loaded drug with the f-CNT and as a result can promote the drug release.

Molecular Perspective Mechanism for Drug Loading on Carbon Nanotube-Dendrimer: A Coarse-Grained Molecular Dynamics Study.

The loading mechanism of the protein ubiquitin and the drug pyrene, as a representatives of large and small molecules, onto the drug carrier carbon nanotube-polyamidoamine (PAMAM) was studied by using coarse-grained molecular dynamics simulation. The results indicated that the optimum and stable drug delivery system for protein loading can be obtained by inserting the molecules in the sequence of: (i) PAMAM, (ii) protein, and (iii) PAMAM. Also, it was found that properly adjusting the weight ratio of PAMAM to the protein, defined as MwPAMAM/ Mwprotein (where Mw is the molecular weight) can lead to achieve a stable system for loading the protein. However, for pyrene loading, it was found that the insertion sequence has no significant effect and only encapsulation of the pyrene molecules into PAMAM and adjustment of the weight ratio of PAMAM to pyrene ( MwPAMAM/ Mwpyrene) can affect the stability of the drug delivery system.

Effect of drug amlodipine on the charged lipid bilayer cell membranes DMPS and DMPS + DMPC: a molecular dynamics simulation study.

In this work, the effects of the anti-hypertensive drug amlodipine in native and PEGylated forms on the malfunctioning of negatively charged lipid bilayer cell membranes constructed from DMPS or DMPS + DMPC were studied by molecular dynamics simulation. The obtained results indicate that amlodipine alone aggregates and as a result its diffusion into the membrane is retarded. In addition, due to their large size aggregates of the drug can damage the cell, rupturing the cell membrane. It is shown that PEGylation of amlodipine prevents this aggregation and facilitates its diffusion into the lipid membrane. The interaction of the drug with negatively charged membranes in the presence of an aqueous solution of NaCl, as the medium, is investigated and its effects on the membrane are considered by evaluating the structural properties of the membrane such as area per lipid, thickness, lipid chain order and electrostatic potential difference between bulk solution and lipid bilayer surface. The effect of these parameters on the diffusion of the drug into the cell is critically examined and discussed.

A coarse grained molecular dynamics simulation study on the structural properties of carbon nanotube-dendrimer composites.

By employing coarse grained (CG) molecular dynamics (MD) simulation, the effect of the size and hydrophilic/hydrophobic properties of the interior/exterior structures of the dendrimers in carbon nanotube (CNT)-dendrimer composites has been studied, to find a stable composite with high solubility in water and the capability to be used in drug delivery applications. For this purpose, composites consisting of core-shell dendrimer complexes including: [PPI{core}-PAMAM{shell}], [PAMAM{core}-polyethyleneglycol (PEG){shell}] and [PAMAM{core}-fattyacid (FTA){shell}] were constructed. A new CG model for the fatty acid (FTA) molecules as functionalized to the dendrimer was developed, which, unlike the previous models, could generate the structural conformations of the FTA properly. The obtained results indicated that the dendrimer complexes with short FTA chains can form stable composites with the CNT. Also, it was found that the pristine PAMAM and PPI-PAMAM with small PPI, and PAMAM-PEG dendrimers with short PEG chains, can distribute their chains into the water medium and interact with the CNT efficiently, to form a stable water-soluble CNT-dendrimer composite. The results demonstrated that the structural difference between the interior and exterior of a core-shell dendrimer complex can prevent the core and the interior layers of the dendrimer complex from interacting with the CNT. An overall analysis of the results manifested that the CNT-PAMAM:4-PEG:4 is the most stable composite, due to strong binding of the dendrimer with the CNT while also having high solubility in water, and its core retains its structure properly and unchanged, suitable for encapsulating drugs in the targeted delivery applications.

Interaction of drugs amlodipine and paroxetine with the metabolizing enzyme CYP2B4: a molecular dynamics simulation study.

The interactions of the drugs amlodipine and paroxetine, which are prescribed respectively for treatment of hypertension and depression, with the metabolizing enzyme cytochrome CYP2B4 as the drug target, have been studied by molecular dynamics (MD) simulation. Poly ethylene glycol was used to control the drugs' interactions with each other and with the target CYP2B4. Thirteen simulation systems were carefully designed, and the results obtained from MD simulations indicated that amlodipine in the PEGylated form prescribed with paroxetine in the nonPEGylated form promotes higher cytochrome stability and causes fewer fluctuations as the drugs approach the target CYP2B4 and interact with it. The simulation results led us to hypothesize that the combination of the drugs with a specific drug ratio, as proposed in this work, manifests more effective diffusivity and less instability while metabolizing with enzyme CYP2B4. Also, the active residues in the CYP2B4 enzyme that interact with the drugs were determined by MD simulation, which were consistent with the reported experimental results. Graphical Abstract Efficient drug-enzyme interactions, as a result of PEGylation.

Effect of pristine and functionalized single- and multi-walled carbon nanotubes on CO2 separation of mixed matrix membranes based on polymers of intrinsic microporosity (PIM-1): a molecular dynamics simulation study.

Molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) simulations were conducted to investigate the transport properties of carbon dioxide, methane, nitrogen, and oxygen through pure and mixed matrix membranes (MMMs) based on polymers of intrinsic microporosity (PIM-1). For this purpose, first, 0.5 to 3 wt% of pristine single-walled carbon nanotube (p-SWCNT) and multi-walled carbon nanotube (p-MWCNT) were embedded into the pure PIM-1, and then for better dispersion of CNT particles into the polymer matrix and to improve the performance of the resulting MMMs, polyethylene glycol (PEG) functionalized SWCNT and MWCNT (f-SWCNT and f-MWCNT, respectively) were loaded. The characterization of the obtained MMMs was carried out by using density, glass transition temperature, X-ray pattern, and fractional free volume calculations. Comparing the obtained results with the available reported experimental data, indicate the authenticity of the applied simulation approach. The simulation results exhibit that the pristine and PEG-functionalized CNT particles improve the transport properties such as diffusivity, solubility, and permeability of the PIM-1 membranes, without sacrificing their selectivity. Also, the MMMs incorporated with 2 wt% of the functionalized CNT particles indicate better performance for the CO2 separation from other gases. According to the calculated results, the highest permeability and diffusivity for CO2 are observed in the [PIM-1/f-SWCNT] MMM among the other membranes which represent that the loading of the f-SWCNTs can enhance the CO2 separation performance of PIM-1 more than other CNTs studied in this work.

Combination of anti-hypertensive drugs: a molecular dynamics simulation study.

The anti-hypertensive drugs amlodipine, atenolol and lisinopril, in ordinary and PEGylated forms, with different combined-ratios, were studied by molecular dynamics simulations using GROMACS software. Twenty simulation systems were designed to evaluate the interactions of drug mixtures with a dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane, in the presence of water molecules. In the course of simulations, various properties of the systems were investigated, including drug location, diffusion and mass distribution in the membrane; drug orientation; the lipid chain disorder as a result of drug penetration into the DMPC membrane; the number of hydrogen bonds; and drug surface area. According to the results obtained, combined drugs penetrate deeper into the DMPC lipid bilayer membrane, and the lipid chains remain ordered. Also, the combined PEGylated drugs, at a combination ratio of 1:1:1, enhance drug penetration into the DMPC membrane, reduce drug agglomeration, orient the drug in a proper angle for easy penetration into the membrane, and decrease undesirable lipotoxicity due to distorted membrane self-assembly and thickness. Graphical abstract ᅟ.

Hybrid Dendrimers of PPI(core)-PAMAM(shell): A Molecular Dynamics Simulation Study.

The structural properties of hybrid dendrimers PPI(core)-PAMAM(shell) for application in drug delivery are studied by coarse-grained molecular dynamics simulation, and their capacity to encapsulate drug guest molecules such as pyrene is investigated by changing the core (PPI) in the PPI-PAMAM hybrids. For this purpose, a coarse-grained model for PPI dendrimer is developed and is used to predict the structural properties as a function of PPI core size, such as the size of hybrid dendrimers, the depth of water penetration, the extent of back-folding of their chain terminals, the size and distribution of created cavities, and asphericity. The results show that the location of pyrene in the interior structure of the hybrids is independent of PPI core size and the branching chains create a barrier against the penetrating molecules in the shell of PPI. Then, by adding the PAMAM to the surface of PPI, this barrier is removed, and this will enhance the encapsulation capacity of the hybrid.

Molecular dynamics simulation of coarse-grained poly(L-lysine) dendrimers.

Poly(L-lysine) (PLL) dendrimer are amino acid based macromolecules and can be used as drug delivery agents. Their branched structure allows them to be functionalized by various groups to encapsulate drug agents into their structure. In this work, at first, an attempt was made on all-atom simulation of PLL dendrimer of different generations. Based on all-atom results, a course-grained model of this dendrimer was designed and its parameters were determined, to be used for simulation of three generations of PLL dendrimer, at two pHs. Similar to the all-atom, the coarse-grained results indicated that by increasing the generation, the dendrimer becomes more spherical. At pH 7, the dendrimer had larger size, whereas at pH 12, due to back folding of branching chains, they had the tendency to penetrate into the inner layers. The calculated radial probability and radial distribution functions confirm that at pH 7, the PLL dendrimer has more cavities and as a result it can encapsulate more water molecules into its inner structure. By calculating the moment of inertia and the aspect ratio, the formation of spherical structure for PLL dendrimer was confirmed.

Molecular dynamics simulation study of chitosan and gemcitabine as a drug delivery system.

By using molecular dynamics (MD) simulation, biodegradable biopolymer chitosan as a carrier for the drug gemcitabine was investigated and the effect of three initial drug concentrations (10, 40, and 80%) on its loading efficiency was studied. Then water was added to the systems of drug and biopolymer and the effects of water on the interactions of drug and chitosan and on the drug loading efficiency were examined. From the results it was found that the maximum loading of the drug occurred at 40% of the drug concentration. The radial distribution function calculations indicated that in the absence of water molecules, the drug molecules were located at shorter distance from chitosan and the loading efficiency of the drug in these systems was higher.

Interaction of PEGylated anti-hypertensive drugs, amlodipine, atenolol and lisinopril with lipid bilayer membrane: A molecular dynamics simulation study.

The interaction of PEGylated anti-hypertensive drugs, amlodipine, atenolol and lisinopril with lipid bilayer membrane dimyristoylphosphatidylcholine (DMPC) has been studied in nine different simulation systems consisting of 128 lipid molecules and appropriate number of water molecules by molecular dynamics method and by utilizing GROMACS software. The influences of PEGylation on the mentioned drugs and the differences in application of two types of spacer molecules on the performance of drugs and DMPC membrane have been evaluated and mass density of the components in the simulation box, mean square displacement (MSD), electrostatic potential, hydrogen bonding, radial distribution function (RDF), area per lipid, order parameter, and angle distribution of the component molecules including drug, DMPC and PEG has been investigated. Furthermore, umbrella sampling analysis indicated that, PEGylation of the drugs made amlodipine to behave more hydrophilic, whereas in case of lisinopril and atenolol, PEGylation made these drugs to behave more hydrophobic. In almost all of the simulated systems, PEGylation increased the diffusion coefficient of the drugs.

Molecular dynamic simulation of Ca2+-ATPase interacting with lipid bilayer membrane.

In biomedical and drug delivery treatments, protein Ca2+-ATPase in the lipid bilayer (plasma) membrane plays a key role by reducing multidrug resistance of the cancerous cells. The lipid bilayer membrane and the protein Ca2+-ATPase were simulated by utilising the Gromacs software and by applying the all-atom/united atom and coarse-grained models. The initial structure of Ca2+-ATPase was derived from X-ray diffraction and electron microscopy patterns and was placed in a simulated bilayer membrane of dipalmitoylphosphatidylcholine. The conformational changes were investigated by evaluating the root mean square deviation, root mean square fluctuation, order parameter, diffusion coefficients, partial density, thickness and area per lipid.

Aqueous poly(amidoamine) dendrimer G3 and G4 generations with several interior cores at pHs 5 and 7: a molecular dynamics simulation study.

Poly(amidoamine) (PAMAM) dendrimers play an important role in drug delivery systems, because the dendrimers are susceptible to gain unique features with modification of their structure such as changing their terminals or improving their interior core. To investigate the core improvement and the effect of core nature on PAMAM dendrimers, we studied two generations G3 and G4 PAMAM dendrimers with the interior cores of commonly used ethylendiamine (EDA), 1,5-diaminohexane (DAH), and bis(3-aminopropyl) ether (BAPE) solvated in water, as an aqueous dendrimer system, by using molecular dynamics simulation and applying a coarse-grained (CG) dendrimer force field. To consider the electrostatic interactions, the simulations were performed at two protonation states, pHs 5 and 7. The results indicated that the core improvement of PAMAM dendrimers with DAH produces the largest size for G3 and G4 dendrimers at both pHs 5 and 7. The increase in the size was also observed for BAPE core but it was not so significant as that for DAH core. By considering the internal structure of dendrimers, it was found that PAMAM dendrimer shell with DAH core had more cavities than with BAPE core at both pHs 5 and 7. Also the moment of inertia calculations showed that the generation G3 is more open-shaped and has higher structural asymmetry than the generation G4. Possessing these properties by G3, specially due to its structural asymmetry, make penetration of water beads into the dendrimer feasible. But for higher generation G4 with its relatively structural symmetry, the encapsulation efficiency for water molecules can be enhanced by changing its core to DAH or BAPE. It is also observed that for the higher generation G4 the effect of core modification is more profound than G3 because the core modification promotes the structural asymmetry development of G4 more significantly. Comparing the number of water beads that penetrate into the PAMAM dendrimers for EDA, DAH, and BAPE cores indicates a significant increase when their cores have been modified with DAH or BAPE and substantiates the effective influence of the core nature in the dendrimer encapsulation efficiency.

Molecular simulation study of PAMAM dendrimer composite membranes.

Pure polysulfone (PSF) and its composites with chitosan (CST), hyaluronic acid (HA), conventional poly(amidoamine), and hydroxyl poly(amidoamine) dendrimers as the membranes for separation of the gases, methane, carbon dioxide, hydrogen sulfide, nitrogen, and oxygen have been studied by molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) simulations. The transport properties (solubility, diffusivity, and permeability) of pure and gas mixtures in the membranes were calculated and the results of the simulations were compared with the available experimental data. The simulated structural properties of the pure and composite PSF membranes including occupied volume, free volume, surface area, fractional free volume (FFV), and radius of gyration (R g ) were evaluated and their effects on the separability of the gases by the membranes were analyzed and interpreted by the obtained results.