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     

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     

Binding of nucleobases with graphene and carbon nanotube: a review of computational studies

Functionalized carbon nanotubes (CNTs) constitute a new class of nanostructured materials that have vast applications in CNT purification and separation, biosensing, drug delivery, etc. Hybrids formed from the functionalization of CNT with biological molecules have shown interesting properties and have attracted great attention in recent years. Of particular interest is the hybridization of single- or double-stranded nucleic acid (NA) with CNT. Nucleobases, as the building blocks of NA, interact with CNT and contribute strongly to the stability of the NA–CNT hybrids and their properties. In this work, we present a thorough review of previous studies on the binding of nucleobases with graphene and CNT, with a focus on the simulation works that attempted to evaluate the structure and strength of binding. Discrepancies among these works are identified, and factors that might contribute to such discrepancies are discussed.     

Theoretical study of the interaction between X (H,F) and graphene

This study investigates the interaction between X (X = H and F) and graphene C₅₄H₁₈ (D_6 h), and the potential energy surface of the graphene radical. The calculations on the structures and energies are further discussed thermodynamically and kinetically using the density function theory method at the B3LYP/6-31G (d) level. Our findings show that there are four distinct isomers of C₅₄H₁₈–X. C₅₄H₁₈–H² and C₅₄H₁₈–F⁴ are the most stable isomers in their own systems. In addition, the transition states, as well as reaction pathways of H transferring between different key points on representative patch, are given to explore the possible reaction mechanism. Finally, the stability of C₅₄H₁₈–X₂ is discussed through the density functional theory.     

Molecular dynamics simulation of dispersion improvement of graphene sheets in nanofluids by steric hindrance resulting from functional groups

Nanofluids are candidate materials for thermal management of heat transfer equipment. Practical applications of thermally enhanced nanofluids contribute to the reduction of weight of systems, leading to improved energy efficiency. Microsize particles sink into the systems because of gravity, therefore rendering the addition meaningless in terms of improving thermal properties. However, nanoparticles can be buoyant, leading to Brownian motion in the fluid, when they do not aggregate with each other. The most important factor in nanofluids is long-term stability of the dispersion in the fluid. Numerous studies have reported the dispersion stability; functional groups attached to nanoparticles play a role in causing steric hindrance and have an affinity for the surrounding fluid, resulting in preserving the dispersion. We investigate the structural effects on dispersion by molecular dynamics simulations of nanofluid containing graphene sheets with functional groups of varying lengths at the surface. The results demonstrate that short functional groups were too short to cause significant steric hindrance, while relatively longer functional groups tended to stack onto the graphene sheets, leading to trapping due to strong van der Waals interactions. Additionally, we discuss the minimum number of functional groups necessary for maintaining dispersion through calculations of the area of a single functional group.     

Wrinkling and thermal conductivity of one graphene sheet under shear

Tersoff-potential - based molecular dynamics method is used to simulate wrinkling deformation of one graphene sheet under shear, and the obtained deformation is compared with analytical solutions of macro-membrane. Furthermore, thermal conductivity of the wrinkled graphene at different temperatures is calculated. It is found that (1) the wrinkling deformation of graphene sheet under shear is close to the analytical solutions of macro-membrane under shear, which implies that the solutions of macro-membrane are applicable to predict the wrinkling deformation of graphene sheets under shear; (2) the more serious the wrinkling of the graphene under shear is, the stronger the phonon scattering is and, therefore, the lower the thermal conductivity of the wrinkled graphene is; (3) within the temperature range of 400–700 K, the thermal conductivity of graphene sheet decreases with increase in temperature.     

Formation of binocular-like structure using graphene nanosheet and carbon nanotubes

Binocular-like carbon structure was successfully fabricated using forced-field-based molecular dynamics simulations. The mechanism for the self-scrolling of graphene sheet using carbon nanotubes as template was analysed. By analysing the self-scrolling process, it was found that the van der Waals interactions were responsible for the formation of the binocular-like carbon nanostructures. Furthermore, the position and number of carbon tubes influenced the self-scrolling process. Moreover, the mechanism was also discussed. Our theoretical results will provide researchers a powerful guide and helpful assistance in designing better targeted programmes in experiments.     

Molecular dynamics simulation of lactate dehydrogenase adsorption onto pristine and carboxylic-functionalized graphene

ABSTRACT

Lactate dehydrogenase (LDH) is a tetrameric enzyme which is composed of two subunits known as LDHA and LDHB, which are encoded by the LDHA and LDHB genes respectively. LDH catalyses the last step in anaerobic glycolysis through the reversible conversion of pyruvate to lactate via coupled oxidation of NADH cofactor. The LDHA plays an important regulatory role in anaerobic glycolysis, by catalysing the final step of the process. Therefore, it is likely that increases in the expression level of LDHA in cancer cells could facilitate the efficiency of anaerobic glycolysis. Measuring the level of serum LDHA is a key step in the diagnosis of many cancer types. In this study, the adsorption, stability, and dynamics of LDHA on the surface of pristine graphene (PG) and carboxylated graphene (COOH-Graphene) were investigated using its molecular dynamics simulation. Variations in root mean square deviation, root mean square fluctuation, solvent accessible surface area and adsorption energy of the LDHA during the simulation were calculated to analyse the effect of PG and COOH-Graphene on the overall conformation of LDHA. Results showed that the adsorption of LDHA on COOH-Graphene is mostly mediated by electrostatic interactions, whereas on the PG, both Van der Waals and π-π interactions are prominent.     

Identification of InhA inhibitors: A combination of virtual screening,molecular dynamics simulations and quantum chemical studies

In the present work, multiple pharmacophore-based virtual screening of the SPECS natural product database was carried out to identify novel inhibitors of the validated biological target, InhA. The pharmacophore models were built from the five different groups of the co-crystallized ligands present within the active site. The generated models with the same features from each group were pooled and subjected to the test set validation, receiver–operator characteristic analysis and Güner–Henry studies. A set of five hypotheses with sensitivity > 0.5, specificity > 0.5, area under curve (AUC) > 0.7, and goodness of hit score > 0.7 were retrieved and exploited for the virtual screening. The common hits (87 molecules) obtained from these hypotheses were processed via drug-likeness filters. The filtered molecules (27 molecules) were compared for the binding modes and the top scored molecules (12 molecules) along with the reference (triclosan (TCL), docking score = ?11.65 kcal/mol) were rescored and reprioritized via molecular mechanics-generalized Born surface area approach. Eventually, the stability of reprioritized (10 molecules) docked complexes was scrutinized via molecular dynamics simulations. Moreover, the quantum chemical studies of the dynamically stable compounds (9 molecules) were performed to understand structural features essential for the activity. Overall, the protocol resulted in the recognition of nine lead compounds that can be targeted against InhA.     

van der Waals potential barrier for cobaltocene encapsulation into single-walled carbon nanotubes: classical molecular dynamics and ab initio study

In this work, we carried out geometry optimisations and classical molecular dynamics for the problem of cobaltocene (CC) encapsulation into different carbon nanotubes (CNTs) ((7,7), (8,8), (13,0) and (14,0) tubes were used). CCs are molecules composed of two aromatic pentagonal rings (C5H5) sandwiching one cobalt atom. From our simulation results, we observed that CC was encapsulated into CNTs (8,8), (13,0) and (14,0). However, for CNT (7,7), the encapsulation could not occur, in disaggrement with some previous works in the literature. Our results show that the encapsulation process is mainly governed by van der Waals potential barriers.     

Molecular modeling,dynamics studies and density functional theory approaches to identify potential inhibitors of SIRT4 protein from Homo sapiens : a novel target for the treatment of type 2 diabetes

Type 2 diabetes is one of the biggest health challenges in the world and WHO projects it to be the 7th leading cause of death in 2030. It is a chronic condition affecting the way our body metabolizes sugar. Insulin resistance is high risk factor marked by expression of Lipoprotein Lipases and Peroxisome Proliferator-Activated Receptor that predisposes to type 2 diabetes. AMP-dependent protein kinase in AMPK signaling pathway is a central sensor of energy status. Deregulation of AMPK signaling leads to inflammation, oxidative stress, and deactivation of autophagy which are implicated in pathogenesis of insulin resistance. SIRT4 protein deactivates AMPK as well as directly inhibits insulin secretion. SIRT4 overexpression leads to dyslipidimeia, decreased fatty acid oxidation, and lipogenesis which are the characteristic features of insulin resistance promoting type 2 diabetes. This makes SIRT4 a novel therapeutic target to control type 2 diabetes. Virtual screening and molecular docking studies were performed to obtain potential ligands. To further optimize the geometry of protein–ligand complexes Quantum Polarized Ligand Docking was performed. Binding Free Energy was calculated for the top three ligand molecules. In view of exploring the stereoelectronic features of the ligand, density functional theory approach was implemented at B3LYP/6-31G* level. 30 ns MD simulation studies of the protein–ligand complexes were done. The present research work proposes ZINC12421989 as potential inhibitor of SIRT4 with docking score (?7.54 kcal/mol), docking energy (?51.34 kcal/mol), binding free energy (?70.21 kcal/mol), and comparatively low energy gap (?0.1786 eV) for HOMO and LUMO indicating reactivity of the lead molecule.     

Molecular simulation for separation of ethylene and ethane by functionalised graphene membrane

ABSTRACT

Graphene is an excellent adsorbent and a membrane material for separation which has attracted wide attention in recent years. Moreover, compared with typical polymer materials, porous graphene has exhibited superior performance. In this paper, molecular dynamics and quantum mechanics were used to explore the appropriate pore size and separation mechanism of graphene. The 2N-Pore-13 (modified by N and H atom) membrane can prevent the penetration of ethane while maintaining high ethylene flux. The permeation rate of ethylene reached 3.7×10⁶ GPU in 5N-Pore-13 membrane, while the one of ethane was only 227 GPU. The mechanism is based on the fact that molecular structure of ethylene is two-dimensional, so that ethylene can get closer to membrane surface when it is adsorbed. When passing through the pores, ethylene has lower enthalpy and entropy barrier.     

Molecular phenomena in colloidal nanostructure synthesis

We review our recent studies of molecular phenomena that can affect the shapes of colloidal nanocrystals grown in solution-phase syntheses. We first present an overview of our first-principles studies with density functional theory aimed at understanding the workings of polymeric structure-directing agents (SDAs). We demonstrate that polyvinylpyrrolidone (PVP), a successful SDA for producing {100}-faceted Ag nanostrutures, possesses selective binding to Ag(100) at the segment level and it is also sufficiently stiff to achieve the correlated binding of neighbouring chain segments. In contrast, polyethylene oxide, which is not so successful, has weaker segmental selectivity and is not as stiff as PVP. We also review our studies of how solvent can influence the growth of anisotropic nanocrystals by directing nanocrystal aggregation.     

Vibrational circular dichroism and IR spectral analysis as a test of theoretical conformational modeling for a cyclic hexapeptide

A model cyclohexapeptide, cyclo-(Phe-(D)Pro-Gly-Arg-Gly-Asp) was synthesized and its IR and VCD spectra were used as a test of density functional theory (DFT) level predictions of spectral intensities for a peptide with a nonrepeating but partially constricted conformation. Peptide structure and flexibility was estimated by molecular dynamics (MD) simulations and the spectra were simulated using full quantum mechanical (QM) approaches for the complete peptide and for simplified models with truncated side chains. After simulated annealing, the backbone conformation of the ring structure is relatively stable, consisting of a normal beta-turn and a tight loop (no H-bond) which does not vary over short trajectories. Only in quite long MD runs at high temperatures do other conformations appear. MD simulations were carried out for the cyclic peptide in water and in TFE, which match experimental solvents, as well as with and without protonation of the Asp carboxyl group. DFT spectral simulations were made using the annealed structure and were extended to include basis set variation, to determine an optimal computational approach, and solvent simulation with a polarized continuum model (PCM). Stepwise full DFT simulation of spectra was done for various sequences with the same backbone geometry but based on (1) solely Gly residues, (2) Ala substitution except Gly and Pro, and (3) complete sequences with side chains. Additionally, a selection of structures was used to compute IR and VCD spectra with the optimal method to determine structural variation effects. The side chains, especially the Asp-COOH and Arg-NH(2) transitions, had an impact on the computed amide frequencies, IR intensities and VCD pattern. Since experimentally these groups would have little chirality, due to conformational variation, they do not impact the observed VCD spectra. Correcting for frequency shifts, the Ala model for the cyclopeptide gives the clearest representation of the amide VCD. The experimental sign pattern for the amide I' band in D(2)O and also the sharper, more intense amide I VCD band in TFE was seen to some degree in one conformer with Type II' turns, but the data favor a mix of structures.     

Molecular dynamics simulations of the aggregation behaviour of overlapped graphene sheets in linear aliphatic hydrocarbons

Molecular dynamics simulations were used to investigate the aggregation of two partially overlapped graphene sheets in hexane, dodecane and eicosane. When partially overlapped graphene sheets are adjacent to one another, they will expel the adsorbed layers of the solvent molecules on the graphene surface, and the amount of overlap will increase. When the overlapped regions of the graphene sheets are separated by solvent molecules, they cannot expel the adsorption layers between them, and so the sheets remain separated. The driving force for aggregation is the van der Waals interaction between the two graphene sheets, while the van der Waals interaction between the graphene sheets and the solvent molecules inhibits graphene aggregation. The diffusion rate of the hydrocarbon molecules with shorter chain lengths is higher. Thus, they diffuse faster during graphene aggregation, which leads to a higher rate of graphene overlapping in the shorter hydrocarbons. This work provides useful insights into graphene aggregation in linear hydrocarbon solvents of varying lengths at the nanoscale.     

Design of novel PhMTNA inhibitors,targeting neurological disorder through homology modeling,molecular docking,and dynamics approaches

Vanishing white matter (VWM) is a hereditary human disease, mostly prevalent in childhood caused by the defects in the eukaryotic initiation factor beta subunits. It is the first disease involved in the translation initiation factor, eIF2B. There is no specific treatment for VWM which mainly affect the brain and ovaries. The gray matter remains normal in all characteristics while the white matter changes texture, coming to the pathophysiology, many initiation factors are involved in the initiation of translation of mRNAs into polypeptides. In this study, the three-dimensional structure of PhMTNA protein was modeled and the stability ascertained through Molecular dynamic simulation (MDS) for 100?ns. The active site residues are conserved with the reported BsMTNA structure which is also confirmed through sitemap prediction. Through virtual screening and induced fit docking, top five leads against PhMTNA protein was identified based on their binding mode and affinity. ADME properties and DFT (Density Functional Theory) studies of these compounds were studied. In addition to that, computational mutagenesis studies were performed to identify the hotspot residues involved in the protein-ligand interactions. Overall analysis showed that the compound NCI_941 has a highest binding energy of ?46.256?kcal mol^?1 in the Arg57Ala mutant. Thus, the results suggest that NCI_941 would act as a potent inhibitor against PhMTNA protein.     

Computational insights into the binding pattern of mitochondrial calcium uniporter inhibitor through homology modeling,molecular dynamics simulation,binding free energy prediction and density functional theory calculation

Abstract

The mitochondrial calcium uniporter (MCU) is the critical protein of the inner mitochondrial membrane that is the primary mediator for calcium uptake into the mitochondrial matrix. Herein we built the optimal homology model of human MCU which was refined through all-atom molecular dynamics simulation. Then, the binding mode of known inhibitor was predicted through molecular docking method, along with molecular dynamics simulation and binding free energy calculation to verify the docking result and stability of the protein-inhibitor complex. Finally, density functional theory (DFT) calculation enhanced our understanding of the molecular interaction of MCU inhibitor. Our research would provide a deeper insight into the interactions between human MCU and its inhibitor, which boosts to develop novel therapy against MCU related disease.

Communicated by Ramaswamy H. Sarma     

The influence of nicotine on pioglitazone encapsulation into carbon nanotube: the investigation of molecular dynamic and density functional theory

In this work, molecular dynamics simulations of the insertion of pioglitazone into the nanotube with chirality (10, 10) at 400 K and 1 bar in the presence and absence of nicotine molecules and in different drug concentrations have been studied. The main aim is consideration of the effect of nicotine in the drug encapsulation process. The results indicate that encapsulation of pioglitazone could be attributed to the water flow via van der Waals and hydrophilic interactions. Because of the existence of the partial π–π interactions between aromatic rings of pioglitazone and the conjugated aromatic rings of nanotube, pioglitazone molecule can enter inside the nanotube. Some physical properties such as hydrogen bonding, number of contacts, also, the diffusion coefficient of the pioglitazone and water molecules, and variation of the center of mass have been calculated during the simulation. Furthermore, computing the electronic structure has also been done on model systems for quantitative determination of the adsorption energy (E_ads). The B3LYP/6-31G* level calculations on four different configurations of pioglitazone/carbon nanotube (CNT) and nicotine/CNT show that the interaction of drug with the inside of the nanotube is stronger than the other forms.