Effect of pristine and functionalized single- and multi-walled carbon nanotubes on CO<Subscript>2</Subscript> 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 CO₂ separation from other gases. According to the calculated results, the highest permeability and diffusivity for CO₂ 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 CO₂ separation performance of PIM-1 more than other CNTs studied in this work.     

Predicting doxorubicin drug delivery by single-walled carbon nanotube through cell membrane in the absence and presence of nicotine molecules: a molecular dynamics simulation study

Abstract

In order to study the interaction of the anticancer agent Doxorubicin with the single-walled carbon nanotubes with different diameters as drug delivery systems, the molecular dynamics (MD) simulations have been used. Also, for design and development of intracellular Doxorubicin drug delivery systems, a series of steered MD simulations are applied to explore the possibility of encapsulated Doxorubicin–carbon nanotube penetration through a lipid bilayer in presence and absence of Nicotine molecules at different pulling rates. Our simulation results showed that in spite of the adsorption of drug molecules on the outer sidewall of the nanotubes, the spontaneous localization of one Doxorubicin molecule into the cavity of the nanovectors with larger diameters is observed. It is found that the presence of Nicotine molecules in extracellular medium increases the required force for pulling nanotube-encapsulated drug as well as the required time for penetration process, especially at higher velocity. Also, the entering process of the Nicotine molecules into the carbon nanotube causes that the encapsulated drug molecule is fully released in the hydrophobic phase of the lipid bilayer.

Communicated by Ramaswamy H. Sarma     

Assessment of the chitosan-functionalized graphene oxide as a carrier for loading thioguanine,an antitumor drug and effect of urea on adsorption process: Combination of DFT computational and molecular dynamics simulation studies

In this study, the interaction thioguanine (TG) anticancer drug with the functionalized graphene oxide (GO) nanosheet surface is theoretically studied in both gas phase and separately in physiological media using the density functional theory (DFT) calculations. DFT calculations indicated the adsorption and solvation energies are negative for f-GONS/TG complexes which propose the adsorption process of TG molecule onto the f-GONS surface is possible from the energetic viewpoint. QTAIM calculations confirm the nature of partially covalent-partially electrostatic between drug and nanosheet. These results are sorely relevant that an approach for loading of TG molecule is the chemical modification of GO using covalent functionalization which can serve as a nanocarrier to load drug molecules. Moreover, to understand the effect of urea on the nature of the interaction between TG and f-GONS, molecular dynamics (MD) simulation was employed. The results indicated that in the presence of urea the adsorption process gets affected and leads to instability of system, while the affinity of the TG for adsorption onto GO surface is increased in pure water.

Communicated by Ramaswamy H. Sarma     

Comparative prediction of binding affinity of Hydroxyurea anti-cancer to boron nitride and carbon nanotubes as smart targeted drug delivery vehicles

Abstract

In this study, the adsorption of Hydroxyurea (HU) onto the inner and outer surfaces of boron nitride and carbon nanotubes (CNTs) was investigated using the density functional theory calculations and molecular dynamics (MDs) simulations in aqueous solution. The values of the adsorption energy show that HU molecule is preferentially adsorbed inside of boron nitride and CNTs with the molecular axis parallel to the tubes axis, which means that the cavity of nanotubes is favorable for encapsulation of this drug. Also, it was found that the HU/boron nitride nanotube (BNNT) system is more stable than the HU/CNT system. The stability of the complexes of HU/ BNNT attributed to the formation of the intermolecular hydrogen bonds between the H atoms of HU molecule and the N atoms of BNNT, which is confirmed by Bader’s quantum theory of atoms in molecules. The natural bond orbital analysis shows the charge transfers occur from HU molecule to nanotubes in all complexes. Moreover, the adsorption of HU molecule on the surfaces of the nanotubes was investigated by explicit water models. Also, the adsorption behavior of HU on the functionalized boron nitride and CNTs is investigated to design and develop new nanocarriers for biomedical applications. Furthermore, MDs simulations are examined in the presence of one and two drug molecules. The obtained results illustrate that the lowest value of Lennard–Jones (L–J) energy between drug and nanotubes exist in the simulation system with two drug molecules.     

Molecular insight into adsorption affinities of Carmustine drug on boron and nitrogen doped functionalized single-walled carbon nanotubes using density functional theory including dispersion correction calculations and molecular dynamics simulation

Abstract

We report a quantum mechanics calculation and molecular dynamics simulation study of Carmustine drug (BNU) adsorption on the surface of nitrogen (N) and boron (B) doped-functionalized single-walled carbon nanotubes. The stability of the optimized complexes is determined on the basis of relative adsorption energy (ΔE_ads). The ΔE_ads results claim that drug molecule tends to adsorb on the nitrogen and boron doped functionalized tubes with the energy values in the range of ?61.177 to ?95.806?kJ/mol. Based on the obtained results, it is observed that N-doping compared with B-doping has improved more effectively drug absorption on the surface of functionalized nanotube. The results of Atoms in Molecule calculations indicate that drug adsorbs molecularly via hydrogen bonds interactions on the surface doped-functionalized carbon nanotubes. Moreover, molecular dynamics simulation is performed to investigate the dynamics behavior of the drug molecules on the nitrogen-doped functionalized carbon nanotube (f-NNT) and functionalized carbon nanotube (f-CNT). The higher average calculated electrostatic and van der Waals energies as well as higher number of intermolecular hydrogen bonds in BNU-f-NNT in comparison with BNU-f-CNT model suggest the more effectual interaction between drug molecules and nitrogen-doped functionalized carbon nanotube.

Communicated by Ramaswamy H. Sarma     

Carbon nanotube prevents the secondary structure formation of amyloid-β trimers: an all-atom molecular dynamics study

Abstract

Increasing evidence shows that the formation of misfolded aggregates amyloid-β (Aβ) peptide is associated with the Alzheimer’s disease (AD). Recent experiments reveal a significant correlation of the amount of trimer species bound to neurons with increasing neuro-toxicity. Our previous computational study demonstrated that carbon nanotubes (CNT) can inhibit effectively the formation of β-sheet-rich oligomers of Aβ(16-22) – a hydrophobic peptide essential for Aβ fibrillization. However, the influence of CNT on the oligomers formed by full-length Aβ remains elusive. In this study, we have investigated the conformational dynamics of Aβ(1-42) trimer, built from an NMR structure of α-helical monomer, in the absence and presence of a single-walled carbon nanotube (SWCNT). Our simulations demonstrate that SWCNT can significantly hinder the trimerisation and prevents the secondary structure formation of Aβ(1-42) peptide. Hydrophobic and aromatic stacking interactions between SWCNT and Aβ play important roles in the secondary structure formation of the Aβ trimer. This study reveals a complete picture of the detailed preventable mechanism of full-length Aβ(1-42) by SWCNT, providing theoretical insights into the development of drug candidates of AD.     

On the vibrational behavior of single- and double-walled carbon nanotubes under the physical adsorption of biomolecules in the aqueous environment: a molecular dynamics study

The adsorption of biomolecules on the walls of carbon nanotubes (CNTs) in an aqueous environment is of great importance in the field of nanobiotechnology. In this study, molecular dynamics (MD) simulations were performed to understand the mechanical vibrational behavior of single- and double-walled carbon nanotubes (SWCNTs and DWCNTs) under the physical adsorption of four important biomolecules (L-alanine, guanine, thymine, and uracil) in vacuum and an aqueous environment. It was observed that the natural frequencies of these CNTs in vacuum reduce under the physical adsorption of biomolecules. In the aqueous environment, the natural frequency of each pure CNT decreased as compared to its natural frequency in vacuum. It was also found that the frequency shift for functionalized CNTs as compared to pure CNTs in the aqueous environment was dependent on the radius and the number of walls of the CNT, and could be positive or negative.     

Simple benzene derivatives adsorption on defective single-walled carbon nanotubes: a first-principles van der Waals density functional study

We have investigated the interaction between open-ended zig-zag single-walled carbon nanotube (SWCNT) and a few benzene derivatives using the first-principles van der Waals density functional (vdW-DF) method, involving full geometry optimization. Such sp ²-like materials are typically investigated using conventional DFT methods, which significantly underestimate non-local dispersion forces (vdW interactions), therefore affecting interactions between respected molecules. Here, we considered the vdW forces for the interacting molecules that originate from the interacting π electrons of the two systems. The ?0.54 eV adsorption energy reveals that the interaction of benzene with the side wall of the SWCNT is typical of the strong physisorption and comparable with the experimental value for benzene adsorption onto the graphene sheet. It was found that aromatics are physisorbed on the sidewall of perfect SWCNTs, as well as at the edge site of the defective nanotube. Analysis of the electronic structures shows that no orbital hybridization between aromatics and nanotubes occurs in the adsorption process. The results are relevant in order to identify the potential applications of noncovalent functionalized systems.

Molecular dynamics simulation of doxorubicin adsorption on a bundle of functionalized CNT

In this study, molecular dynamics simulation is used to investigate the adsorption of an anticancer drug, doxorubicin, on bundles of functionalized single-walled carbon nanotubes (SWNTs) in an aqueous solution. Carboxylic group has been selected as the functional group. Molecular dynamics (MD) simulations are performed for both separated systems containing a SWNT bundle and a functionalized carbon nanotube bundle, and results are compared with existing experimental data. MD results show that doxorubicin can be adsorbed on CNTs using different methods such as entrapment within CNT bundle, attachment to the side wall of the CNT, and adsorption on the CNT inner cavity. For functionalized CNT, the adsorption of drugs on the functional groups is essential for predicting the enhancement of drug loading on the functionalized nanotubes. Furthermore, the adsorption behavior of doxorubicin on CNTs is fitted with Langmuir and Freundlich isotherm models. The results show that Langmuir model can predict the adsorption behavior of doxorubicin on CNTs more accurately than Freundlich model does. As predicted by this isotherm model, the adsorption process of doxorubicin on CNTs is relatively difficult, but it can be improved by increasing the functional groups on the CNTs surface.     

Mitogen activated protein kinase-1 and cell division control protein-42 are putative targets for the binding of novel natural lead molecules: a therapeutic intervention against Candida albicans

Abstract

Candida albicans, fungal yeast causes several lethal infections in immune-suppressed patients and recently emerged as drug-resistant pathogens worldwide. The present study aimed to screen putative drug targets of Candia albicans and to study the binding potential of novel natural lead compounds towards these targets by computational virtual screening and molecular dynamic (MD) simulation. Through extensive analysis of mitogen-activated protein kinase (MAPK) signalling pathways, mitogen-activated protein kinase-1 (HOG1) and cell division control protein-42 (CDC42) genes were prioritized as putative targets based on their virulent functions. The three-dimensional structures of these genes, not available in their native forms, were computationally modeled and validated. 76 lead molecules from various natural sources were screened and their drug likeliness and pharmacokinetic features were predicted. Among these ligands, two lead molecules that demonstrated ideal drug-likeliness and pharmacokinetic features were docked against HOG1 and CDC42 and their binding potential was compared with the binding of conventional drug Fluconazole with their usual target. The prediction was computationally validated by MD simulation. The current study revealed that Cudraxanthone-S present in Cudrania cochinchinensis and Scutifoliamide-B present in Piper scutifolium exhibited ideal drug likeliness, pharmacokinetics and binding potential to the prioritized targets in comparison with the binding of Fluconazole and their usual target. MD simulation showed that CDC42-Cudraxanthone-S and HOG1-Scutifoliamide-B complexes were exhibited stability throughout MD simulation. Thus, the study provides significant insight into employing HOG1 and CDC42 of MAPK as putative drug targets of C. albicans and Cudraxanthone-S and Scutifoliamide-B as potential inhibitors for drug discovery.

Communicated by Ramaswamy H. Sarma     

Probing the Diameter Limit of Single Walled Carbon Nanotubes in SWCNT: Fullerene Solar Cells

In this work, for the first time, the diameter limit of surfactant wrapped single walled carbon nanotubes (SWCNTs) in SWCNT:C₆₀ solar cells is determined through preparation of monochiral small and large diameter nanotube devices as well as those from polychiral mixtures. Through assignment of the different nanotube chiralities by photoluminescence and optical density measurements a diameter limit yielding 0% internal quantum efficiency (IQE) is determined. This work provides insights into the required net driving energy for SWCNT exciton dissociation onto C₆₀ and establishes a family of (n,m) species which can efficiently be utilized in polymer‐free SWCNT:C₆₀ solar cells. Using this approach the largest diameter nanotube with an IQE > 0% is found to be (8,6) with a diameter of 0.95 nm. Possible strategies to extend this diameter limit are then discussed.     

Phosphatidylserine Targets Single-Walled Carbon Nanotubes to Professional Phagocytes In Vitro and In Vivo

Broad applications of single-walled carbon nanotubes (SWCNT) dictate the necessity to better understand their health effects. Poor recognition of non-functionalized SWCNT by phagocytes is prohibitive towards controlling their biological action. We report that SWCNT coating with a phospholipid “eat-me” signal, phosphatidylserine (PS), makes them recognizable in vitro by different phagocytic cells - murine RAW264.7 macrophages, primary monocyte-derived human macrophages, dendritic cells, and rat brain microglia. Macrophage uptake of PS-coated nanotubes was suppressed by the PS-binding protein, Annexin V, and endocytosis inhibitors, and changed the pattern of pro- and anti-inflammatory cytokine secretion. Loading of PS-coated SWCNT with pro-apoptotic cargo (cytochrome c) allowed for the targeted killing of RAW264.7 macrophages. In vivo aspiration of PS-coated SWCNT stimulated their uptake by lung alveolar macrophages in mice. Thus, PS-coating can be utilized for targeted delivery of SWCNT with specified cargoes into professional phagocytes, hence for therapeutic regulation of specific populations of immune-competent cells.     

Functionalized Carbon Nanotubes in the Brain: Cellular Internalization and Neuroinflammatory Responses

The potential use of functionalized carbon nanotubes (f-CNTs) for drug and gene delivery to the central nervous system (CNS) and as neural substrates makes the understanding of their in vivo interactions with the neural tissue essential. The aim of this study was to investigate the interactions between chemically functionalized multi-walled carbon nanotubes (f-MWNTs) and the neural tissue following cortical stereotactic administration. Two different f-MWNT constructs were used in these studies: shortened (by oxidation) amino-functionalized MWNT (oxMWNT-NH₃ ⁺) and amino-functionalized MWNT (MWNT-NH₃ ⁺). Parenchymal distribution of the stereotactically injected f-MWNTs was assessed by histological examination. Both f-MWNT were uptaken by different types of neural tissue cells (microglia, astrocytes and neurons), however different patterns of cellular internalization were observed between the nanotubes. Furthermore, immunohistochemical staining for specific markers of glial cell activation (GFAP and CD11b) was performed and secretion of inflammatory cytokines was investigated using real-time PCR (qRT-PCR). Injections of both f-MWNT constructs led to a local and transient induction of inflammatory cytokines at early time points. Oxidation of nanotubes seemed to induce significant levels of GFAP and CD11b over-expression in areas peripheral to the f-MWNT injection site. These results highlight the importance of nanotube functionalization on their interaction with brain tissue that is deemed critical for the development nanotube-based vector systems for CNS applications.     

Inhibitor binding studies of Mycobacterium tuberculosis MraY (Rv2156c): Insights from molecular modeling,docking, and simulation studies

Abstract

Tuberculosis (TB) is a contagious disease caused by Mycobacterium tuberculosis (M.tb) or tubercule bacillus, and H37Rv is the most studied clinical strain. The recent development of resistance to existing drugs is a global health-care challenge to control and cure TB. Hence, there is a critical need to discover new drug targets in M.tb. The members of peptidoglycan biosynthesis pathway are attractive target proteins for antibacterial drug development. We have performed in silico analysis of M.tb MraY (Rv2156c) integral membrane protein and constructed the three-dimensional (3D) structure model of M.tb MraY based on homology modeling method. The validated model was complexed with antibiotic muraymycin D2 (MD2) and was used to generate structure-based pharmacophore model (e-pharmacophore). High-throughput virtual screening (HTVS) of Asinex database and molecular docking of hits was performed to identify the potential inhibitors based on their mode of interactions with the key residues involved in M.tb MraY–MD2 binding. The validation of these molecules was performed using molecular dynamics (MD) simulations for two best identified hit molecules complexed with M.tb MraY in the lipid bilayer, dipalmitoylphosphatidyl-choline (DPPC) membrane. The results indicated the stability of the complexes formed and retained non-bonding interactions similar to MD2. These findings may help in the design of new inhibitors to M.tb MraY involved in peptidoglycan biosynthesis.     

Structural dynamics behind variants in pyrazinamidase and pyrazinamide resistance

Abstract

Pyrazinamide (PZA) is an important component of first-line anti-tuberculosis (anti-TB) drugs. The anti-TB agent is activated into an active form, pyrazinoic acid (POA), by Mycobacterium tuberculosis (MTB) pncA gene encoding pyrazinamidase (PZase). The major cause of PZA-resistance has been associated with mutations in the pncA gene. We have detected several novel mutations including V131F, Q141P, R154T, A170P, and V180F (GeneBank Accession No. MH461111) in the pncA gene of PZA-resistant isolates during PZA drug susceptibility testing followed by pncA gene sequencing. Here, we investigated molecular mechanism of PZA-resistance by comparing the results of experimental and molecular dynamics. The mutants (MTs) and wild type (WT) PZase structures in apo and complex with PZA were subjected to molecular dynamic simulations (MD) at the 40?ns. Multiple factors, including root mean square deviations (RMSD), binding pocket, total energy, dynamic cross correlation, and root mean square fluctuations (RMSF) of MTs and WT were compared. The MTs attained a high deviation and fluctuation compared to WT. Binding pocket volumes of the MTs, were found, lower than the WT, and the docking scores were high than WT while shape complementarity scores were lower than that of the WT. Residual motion in MTs are seemed to be dominant in anti-correlated motion. Mutations at locations, V131F, Q141P, R154T, A170P, and V180F, might be involved in the structural changes, possibly affecting the catalytic property of PZase to convert PZA into POA. Our study provides useful information that will enhance the understanding for better management of TB. Abbreviations DST drug susceptibility testing

Δelec electrostatic energy

LJ Lowenstein–Jensen medium

MGIT mycobacterium growth indicator tubes

MTs mutants

MD molecular dynamic simulations

MTB Mycobacterium tuberculosis

NALC–NaOH N-acetyl-l-cysteine–sodium hydroxide

NIH National Institutes of Health

NPT amount of substance (N), pressure (P) temperature (T)

NVT moles (N), volume (V) temperature (T)

PZase pyrazinamidase

Δps polar solvation energy

PTRL Provincial Tuberculosis Reference Laboratory

RMSD root mean square deviations

RMSF root mean square fluctuations

ΔSASA solvent accessible surface area energy

TB tuberculosis

GTotal total binding free energy

ΔvdW Van der Waals energy

WT wild type

Communicated by Ramaswamy H. Sarma     

Prediction of chlortetracycline adsorption on the Fe3O4 nanoparticle using molecular dynamics simulation

Abstract

Molecular dynamics (MD) simulation was applied to investigate the adsorption mechanism of chlortetracycline (CTC) antibiotic molecule as the aqueous pollutant on the Fe₃O₄ nanoparticle (NP). Two different NP sizes with a diameter of about 1.4?nm and 3.5?nm were selected. Initially, the stability of both NPs in water was investigated by calculating radial distribution function curves of NP atoms. Simulation results confirmed the stable crystallographic structures of both NPs. However, small NP induce greater structural stabilization. Then, CTC molecules were adsorbed on NPs surface in various pollutant concentrations. Electrostatic and hydrogen bond were the major types of interactions between CTC molecules and the adsorbent surface. CTC molecules formed a complex with NP surface from their amine side chains; while they were parallel to each other in their aromatic rings and π-π bond between two CTC molecules was formed. Diffusion rate of CTC molecules could predict the adsorption mechanism. At lower concentration of CTC, CTC molecules tend to adsorb on the NP surface. At these concentrations, the diffusion rate of CTC was high. By increasing the CTC concentration, the pollutant agglomeration was enhanced which decreased the diffusion rate. At this time, the surface of NP was saturated. In addition, the results of isotherm curves showed that CTC adsorption on small NPs could be defined with both Langmuir and Freundlich isotherm models, while Freundlich isotherm model was more appropriate for larger NPs. In conclusion, observations confirmed that MD simulation could successfully predict the behavior of CTC adsorption on the Fe₃O₄ NP surface.

Communicated by Ramaswamy H. Sarma     

Virtual identification of novel PPARα/γ dual agonists by 3D-QSAR,molecule docking and molecular dynamics studies

Abstract

Peroxisome proliferator-activated receptors (PPARs) are considered important targets for the treatment of Type 2 diabetes (T2DM). To accelerate the discovery of PPAR α/γ dual agonists, the comparative molecular field analysis (CoMFA) were performed for PPARα and PPARγ, respectively. Based on the molecular alignment, highly predictive CoMFA model for PPARα was obtained with a cross-validated q² value of 0.741 and a conventional r² of 0.975 in the non-cross-validated partial least-squares (PLS) analysis, while the CoMFA model for PPARγ with a better predictive ability was shown with q² and r² values of 0.557 and 0.996, respectively. Contour maps derived from the 3D-QSAR models provided information on main factors towards the activity. Then, we carried out structural optimization and designed several new compounds to improve the predicted biological activity. To investigate the binding modes of the predicted compounds in the active site of PPARα/γ, a molecular docking simulation was carried out. Molecular dynamic (MD) simulations indicated that the predicted ligands were stable in the active site of PPARα/γ. Therefore, combination of the CoMFA and structure-based drug design results could be used for further structural alteration and synthesis and development of novel and potent dual agonists. Abbreviations DM diabetes mellitus

T2DM type 2 diabetes

PPARs peroxisome proliferator-activated receptors

LBDD ligand based drug design

3D-QSAR three-dimensional quantitative structure activity relationship

CoMFA comparative molecular field analysis

PLS partial least square

LOO leave-one-out

q² cross-validated correlation coefficient

ONC optimal number of principal components

r² non-cross-validated correlation coefficient

SEE standard error of estimate

F the Fischer ratio

r²_pred predictive correlation coefficient

DBD DNA binding domain

MD molecular dynamics

RMSD root-mean-square deviation

RMSF root mean square fluctuations

Communicated by Ramaswamy H. Sarma     

Performance Enhancement of Polymer‐Free Carbon Nanotube Solar Cells via Transfer Matrix Modeling

Polymer‐free (6,5) single‐walled carbon nanotubes (SWCNTs) prepared using the gel permeation approach are integrated into SWCNT:C₆₀ solar cells. Evaporation‐driven self‐assembly is used to form large‐area SWCNT thin films from the surfactant‐stabilized aqueous suspensions. The thicknesses of various layers within the solar cell are optimized by theoretical modeling using transfer matrix calculations, where the distribution of the electric field within the stack is matched to light absorption by the SWCNTs through either their primary (S₁₁) or secondary (S₂₂) absorption peaks, or a combination thereof. The validity of the model is verified experimentally through a detailed parameter study and then used to develop SWCNT:C₆₀ solar cells with high open‐circuit voltage (0.44 V) as well as a cutting‐edge internal quantum efficiency of up to 86% through the nanotube S₁₁ transition, over an active area of 0.105 cm².     

Understanding co-loading of doxorubicin and camptothecin on graphene and folic acid-conjugated graphene for targeting drug delivery: classical MD simulation and DFT calculation

Abstract

The surface modification ability is one of the remarkable characters of graphene (G) nanosheet. Based on this strategy, G surface is modified with folic acid (FA) to improve the targeting delivery of chemotherapy agents. The dual delivery strategy for the transport of doxorubicin (DOX) and camptothecin (CPT) by using G and folic acid functionalized G nanocarriers is examined. The density functional theory (DFT) and molecular dynamics (MD) simulation are employed to gain a deep insight into the nature of the drug and the carrier interactions. The obtained results indicate that the drug molecules spontaneously move toward the carriers and form stable complexes. In the graphene-based systems, the drug molecules form strong π-π interactions with the carrier surface. It is found that the FA functionalization of G (FA-G) not only improves targeting effect but also reinforces drug-carrier interaction. Furthermore, the MD and DFT results show that interaction of DOX molecules with G and FA-G is stronger than CPT. We believe that the results obtained from this study can be helpful to improve the drug effectiveness in cancer treatment.

Communicated by Ramaswamy H. Sarma     

Modelling Diffusion in Zeolites by Molecular Dynamics Simulations

Abstract

The diffusion of molecules sorbed in zeolites is of growing interest for understanding the mechanisms of chemical processes with regard to selectivity and reactivity [1].

MD simulations give insight into physical systems on the molecular level allowing to study and visualize the motion of molecules even beyond the possibilities of experiments [2,3]. Single system parameters can easily be varied to study their influence, also those parameters that are fixed in reality (e.g., the size of particles). We present a cross section of our recent work to illustrate the capabilities of MD: The self diffusion coefficients (D) of a mixture of methane and xenon in silicalite show remarkable deviations from those of the pure species. This is shown and confirmed by PFG NMR experiments [4].

Simulating ethane in zeolite A the mechanism of diffusion has been studied. The effects of rotation on the diffusion lead to cases where D decreases with growing temperature [5].

The independence of self diffusion on lattice vibrations is proven even for zeolites with windows of guest particle size comparing simulations with rigid and vibrating zeolite lattice [6].