Search results for: molecular docking simulations

Authors: Ling Dai

Abstract:

Surface-sensing devices such as atomic force microscope have been widely used to characterize the surface structure and properties of nanoscale polymer films. However, using molecular dynamics simulations, we show that there is intrinsic and unavoidable inelastic deformation at polymer surfaces induced by the sensing tip. For linear chain polymers like perfluoropolyether, such tip-induced deformation derives from the differences in the atomic interactions which are atomic specie-based Van der Waals interactions, and resulting in atomic shuffling and causing inelastic alternation in both molecular structures and mechanical properties at the regions of the polymer surface. For those aromatic chain polymers like epoxy, the intrinsic deformation is depicted as the intra-chain rotation of aromatic rings and kinking of linear atomic connections. The present work highlights the need to reinterpret the data obtained from surface-sensing tests by considering this intrinsic inelastic deformation occurring at polymer surfaces.

Keywords: polymer, surface, nano, molecular dynamics

Procedia PDF Downloads 339

Authors: Gagan Dhawan

Abstract:

Mycobacterium tuberculosis was described as ‘captain of death’ with an inherent property of multiple drug resistance majorly caused by the competent mechanism of efflux pumps. In this study, various open source tools combining chemo-informatics with bioinformatics were used for efficient in-silico drug designing. The efflux pump, Rv1218c, belonging to the ABC transporter superfamily, which is predicted to be a tetronasin-transporter in M. tuberculosis was targeted. Recent studies have shown that Rv1218c forms a complex with two more efflux pumps (Rv1219c and Rv1217c) to provide multidrug resistance to the bacterium. The 3D structure of the protein was modeled (as the structure was unavailable in the previously collected databases on this gene). The TMHMM analysis of this protein in TubercuList has shown that this protein is present in the outer membrane of the bacterium. Virtual screening of compounds from various publically available chemical libraries was performed on the M. tuberculosis protein using various open source tools. These ligands were further assessed where various physicochemical properties were evaluated and analyzed. On comparison of different physicochemical properties, toxicity and docking, the ligand 2-(hydroxymethyl)-6-[4, 5, 6-trihydroxy-2-(hydroxymethyl) tetrahydropyran-3-yl] oxy-tetrahydropyran-3, 4, 5-triol was found to be best suited for further studies.

Keywords: drug resistance, efflux pump, molecular docking, virtual screening

Procedia PDF Downloads 356

Authors: Misbah Khan, Jian Wen, Muhammad Asif Shakoori

Abstract:

The Organic Rankine Cycle (ORC) plays an important role in converting low-temperature heat sources into electrical power by using refrigerants as working fluids. The thermophysical properties of working fluids are essential for designing ORC. HFO-1336mzz(Z) (cis-1,1,1,4,4,4-hexafluoro-2-butene) considered as working fluid and have almost 99% low GWP and relatively same thermophysical properties used as a replacement of HFC-245fa (1,1,1,3,3-pentafluoro-propane). The environmental, safety, healthy and thermophysical properties of HFO-1336mzz(Z) are needed to use it in a practical system. In this paper, Molecular dynamics simulations were used to investigate the Homogeneous condensation, thermophysical and structural properties of HFO-1336mzz(Z) and HFC-245fa. The effect of various temperatures and pressures on thermophysical properties and condensation was extensively investigated. The liquid densities and isobaric heat capacities of this refrigerant was simulated at 273.15K to 353.15K temperatures and pressure0.5-4.0MPa. The simulation outcomes were compared with experimental data to validate our simulation method. The mean square displacement for different temperatures was investigated for dynamical analysis. The variations in potential energies and condensation rate were simulated to get insight into the condensation process. The radial distribution function was simulated at the micro level for structural analysis and revealed that the phase transition of HFO-1336mzz(Z) did not affect the intramolecular structure.

Keywords: homogenous condensation, refrigerants, molecular dynamics simulations, organic rankine cycle

Procedia PDF Downloads 138

Authors: Strahinja Kovačević, Sanja Podunavac-Kuzmanović, Lidija Jevrić, Milica Karadžić

Abstract:

The present study is based on the quantitative structure-activity relationship (QSAR) analysis of eighteen compounds with anti-prion activity. The structures and anti-prion activities (expressed in response units, RU%) of the analyzed compounds are taken from CHEMBL database. In the first step of analysis 85 molecular descriptors were calculated and based on them the hierarchical cluster analysis (HCA) and principal component analysis (PCA) were carried out in order to detect potential significant similarities or dissimilarities among the studied compounds. The calculated molecular descriptors were physicochemical, lipophilicity and ADMET (absorption, distribution, metabolism, excretion and toxicity) descriptors. The first stage of the QSAR analysis was simple linear regression modeling. It resulted in one acceptable model that correlates Henry's law constant with RU% units. The obtained 2D-QSAR model was validated by cross-validation as an internal validation method. The validation procedure confirmed the model’s quality and therefore it can be used for prediction of anti-prion activity. The next stage of the analysis of anti-prion activity will include 3D-QSAR and molecular docking approaches in order to select the most promising compounds in treatment of prion diseases. These results are the part of the project No. 114-451-268/2016-02 financially supported by the Provincial Secretariat for Science and Technological Development of AP Vojvodina.

Keywords: anti-prion activity, chemometrics, molecular modeling, QSAR

Procedia PDF Downloads 281

Authors: Sajjad Seifoori

Abstract:

Impact behavior of striker on graphene sheet and carbon nanotube is investigated based on molecular dynamics (MD) simulations. A MD simulation is conducted to obtain the maximum dynamic deflections of a square and rectangular single-layered graphene sheets (SLGSs) with various values of side-length and striker parameter. Effect of (i) chirality, (ii) graphene side-length and nanotube length, (iii) striker mass on the maximum dynamic deflections of graphene and nanotube are investigated. The effect of different types of boundary condition on the maximum dynamic deflections is studied for zigzag and armchair SWCNTs with various aspect ratios (Length/Diameter).

Keywords: impact, molecular dynamic, graphene, spring mass

Procedia PDF Downloads 311

Authors: Anik Barua, Rabiul Hossain, Labonno Barua, Rashadul Hossain, Nurul Absar

Abstract:

Background: Globally, the burden of cancer is increasing consistently. Modern cancer therapies include lots of toxicity in the non-targeted organs reducing the life expectancy of the patients. Hence, scientists are trying to seek noble compounds from natural sources to treat cancer. Objectives: The objectives of the present study are to evaluate the phytochemicals, in vitro antioxidants, and in vivo and in silico anticancer study of various solvent fractions of Tinospora cordifolia (Willd.). Methodology: In this experiment, standard quantitative and qualitative assay methods were used to analyze the phytochemicals. The antioxidant activity was measured using the DPPH and ABTS scavenging methods. The in vivo antitumor activity is evaluated against Ehrlich ascites carcinoma (EAC) cell bearing in Swiss albino mice. In-silico ADME/T and molecular docking study were performed to assess the potential of stated phytochemicals against Transcription Factor STAT3b/DNA Complex of adenocarcinoma. Findings: Phytochemical screening confirmed the presence of flavonoids, alkaloids, glycosides, tannins, and carbohydrates. A significant amount of phenolic (20.19±0.3 mg/g GAE) and flavonoids (9.46±0.18 mg/g GAE) were found in methanolic extract in quantitative screening. Tinospora cordifolia methanolic extract showed promising DPPH and ABTS scavenging activity with the IC50 value of 1222.99 µg/mL and 1534.34 µg/mL, respectively, which was concentration dependent. In vivo anticancer activity in EAC cell-bearing mice showed significant (P < 0.05) percent inhibition of cell growth (60.12±1.22) was found at the highest dose compared with standard drug 5-Fluorouracil (81.18±1.28). Forty-two phytochemicals exhibit notable pharmacokinetics properties and passed drug-likeness screening tests in silico. In molecular docking study, (25S)-3Beta-acetoxy-5-alpha-22-beta-spirost-9(11)-en-12-beta-ol showed docking score (-8.5 kJ/mol) with significant non-bonding interactions with target enzyme. Conclusions: The results were found to be significant and confirmed that the methanolic extract of Tinospora cordifolia has remarkable antitumor activity with antioxidant potential. The Tinospora cordifolia methanolic extract may be considered a potent anticancer agent for advanced research.

Keywords: anticancer, antioxidant, Tinospora cordifolia, EAC cell

Procedia PDF Downloads 96

Authors: Elisa Y. M. Ang, Teng Yong Ng, Jingjie Yeo, Rongming Lin, Zishun Liu, K. R. Geethalakshmi

Abstract:

By stacking carbon nanotubes (CNT) one on top of another, single layer CNT arrays can perform water-salt separation with ultra-high permeability and selectivity. Such outer-wall CNT slit membrane is named as the transverse flow CNT membrane. By adjusting the slit size between neighboring CNTs, the membrane can be configured to sieve out different solutes, right down to the separation of monovalent salt ions from water. Molecular dynamics (MD) simulation results show that the permeability of transverse flow CNT membrane is more than two times that of conventional axial-flow CNT membranes, and orders of magnitude higher than current reverse osmosis membrane. In addition, by carrying out MD simulations with different CNT size, it was observed that the variance in desalination performance with CNT size is small. This insensitivity of the transverse flow CNT membrane’s performance to CNT size is a distinct advantage over axial flow CNT membrane designs. Not only does the membrane operate well under constant pressure desalination operation, but MD simulations further indicate that oscillatory operation can further enhance the membrane’s desalination performance, making it suitable for operation such as electrodialysis reversal. While there are still challenges that need to be overcome, particularly on the physical fabrication of such membrane, it is hope that this versatile membrane design can bring the idea of using low dimensional structures for desalination closer to reality.

Keywords: carbon nanotubes, membrane desalination, transverse flow carbon nanotube membrane, molecular dynamics

Procedia PDF Downloads 177

Authors: Mohammad Faheem, M. Tabish Rehman, Mohd Danishuddin, Asad U. Khan

Abstract:

The worldwide dissemination of CTX-M type β-lactamases is a threat to human health. Previously, we have reported the spread of blaCTX-M-15 gene in different clinical strains of Enterobacteriaceae from the hospital settings of Aligarh in north India. In view of the varying resistance pattern against cephalosporins and other β-lactam antibiotics, we intended to understand the correlation between MICs and catalytic activity of CTX-M-15. In this study, steady-state kinetic parameters and MICs were determined on E. coli DH5α transformed with blaCTX-M-15 gene that was cloned from Enterobacter cloacae (EC-15) strain of clinical background. The effect of conventional β-lactamase inhibitors (clavulanic acid, sulbactam and tazobactam) on CTX-M-15 was also studied. We have found that tazobactam is the best among these inhibitors against CTX-M-15. The inhibition characteristic of tazobactam is defined by its very low IC50 value (6 nM), high affinity (Ki = 0.017 µM) and better acylation efficiency (k+2/K9 = 0.44 µM-1s-1). It forms an acyl-enzyme covalent complex, which is quite stable (k+3 = 0.0057 s-1). Since increasing resistance has been reported against conventional b-lactam antibiotic-inhibitor combinations, we aspire to design a non-b-lactam core containing b-lactamase inhibitor. For this, we screened ZINC database and performed molecular docking to identify a potential non-β-lactam based inhibitor (ZINC03787097). The MICs of cephalosporin antibiotics in combination with this inhibitor gave promising results. Steady-state kinetics and molecular docking studies showed that ZINC03787097 is a reversible inhibitor which binds non-covalently to the active site of the enzyme through hydrogen bonds and hydrophobic interactions. Though, it’s IC50 (180 nM) is much higher than tazobactam, it has good affinity for CTX-M-15 (Ki = 0.388 µM). This study concludes that ZINC03787097 compound can be used as seed molecule to design more efficient non-b-lactam containing b-lactamase inhibitor that could evade pre-existing bacterial resistance mechanisms.

Keywords: ESBL, non-b-lactam-b-lactamase inhibitor, bioinformatics, biomedicine

Procedia PDF Downloads 222

Authors: R. O. Ocaya, J. J. Terblans

Abstract:

The laws of Newtonian mechanics allow ab-initio molecular dynamics to model and simulate particle trajectories in material science by defining a differentiable potential function. This paper discusses some considerations for the coding of ab-initio programs for simulation on a standalone computer and illustrates the approach by C language codes in the context of embedded metallic atoms in the face-centred cubic structure. The algorithms use velocity-time integration to determine particle parameter evolution for up to several thousands of particles in a thermodynamical ensemble. Such functions are reusable and can be placed in a redistributable header library file. While there are both commercial and free packages available, their heuristic nature prevents dissection. In addition, developing own codes has the obvious advantage of teaching techniques applicable to new problems.

Keywords: C language, molecular dynamics, simulation, embedded atom method

Procedia PDF Downloads 288

Authors: A. Rajabpour, A. Hadizadeh Kheirkhah

Abstract:

Heat transfer is of great importance to biological systems in order to function properly. In the present study, specific heat capacity as one of the most important heat transfer properties is calculated for a soluble in water Lysozyme protein. Using equilibrium molecular dynamics (MD) simulation, specific heat capacities of pure water, dry lysozyme, and lysozyme-water solution are calculated at 300K for different weight fractions. It is found that MD results are in good agreement with ideal binary mixing rule at small weight fractions. Results of all simulations have been validated with experimental data.

Keywords: specific heat capacity, molecular dynamics simulation, lysozyme protein, equilibrium

Procedia PDF Downloads 287

Authors: Abhishek Kumar Gupta

Abstract:

Atomistic Molecular Dynamics (MD) Simulations have been performed to study the effect of monovalent salt i.e. NaCl concentration (Cs) and chain degree of neutralization (f) on the structure and dynamics of anionic poly(methacrylic acid) (PMA) in dilute aqueous solutions. In the present study, the attention is to unveil the conformational structure, hydrogen-bonding, local polyion-counterion structure, h-bond dynamics, chain dynamics and thermodynamic enthalpy of solvation of a-PMA in dilute aqueous solutions as a function of salt concentration, Cs and f. The results have revealed that at low salt concentration, the conformational radius of gyration (Rg) increases and then decreases reaching a maximum in agreement with the reported light scattering experimental results. The Rg at f = 1 shows a continual decrease and acquire a plateau value at higher salt concentration in agreement with results obtained by light scattering experiments. The radial distribution functions between PMA, salt and water atoms has been computed with respect to atom and centre-of-mass to understand the intermolecular structure in detail. The results pertaining to PMA chain conformations and hydrogen bond autocorrelation function showcasing the h-bond dynamics will be presented. The results pertaining to chain dynamics will be presented. The results pertaining to counterion condensation on the PMA chain shows greater condensation of Na+ ions on to the carboxylate ions with increase in salt concentration. Moreover, the solvation enthalpy of the system as a function of salt concentration will be presented.

Keywords: conformations, molecular dynamics simulations, NaCl concentration, radial distribution functions

Procedia PDF Downloads 98

Authors: Juan Antonio Mondragon Sanchez, Ruben Santamaria

Abstract:

The DNA G-quadruplex is a four-stranded DNA structure formed by stacked planes of four base paired guanines (G-quartet). Guanine rich DNA sequences appear in many sites of genomic DNA and can potential form G-quadruplexes, such as those occurring at 3'-terminus of the human telomeric DNA. The formation and stabilization of a G-quadruplex by small ligands at the telomeric region can inhibit the telomerase activity. In turn, the ligands can be used to down regulate oncogene expression making G-quadruplex an attractive target for anticancer therapy. Many G-quadruplex ligands have been proposed with a planar core to facilitate the pi–pi stacking and electrostatic interactions with the G-quartets. However, many drug candidates are impossibilitated to discriminate a G-quadruplex from a double helix DNA structure. In this context, it is important to investigate the site topology for the interaction of a G-quadruplex with a ligand. In this work, we determine the free energy surface of a G-quadruplex-ligand to study the binding modes of the G-quadruplex (TG4T) with the daunomycin (DM) drug. The complex TG4T-DM is studied using classical molecular dynamics in combination with metadynamics simulations. The metadynamics simulations permit an enhanced sampling of the conformational space with a modest computational cost and obtain free energy surfaces in terms of the collective variables (CV). The free energy surfaces of TG4T-DM exhibit other local minima, indicating the presence of additional binding modes of daunomycin that are not observed in short MD simulations without the metadynamics approach. The results are compared with similar calculations on a different structure (the mutated mu-G4T-DM where the 5' thymines on TG4T-DM have been deleted). The results should be of help to design new G-quadruplex drugs, and understand the differences in the recognition topology sites of the duplex and quadruplex DNA structures in their interaction with ligands.

Keywords: g-quadruplex, cancer, molecular dynamics, metadynamics

Procedia PDF Downloads 442

Authors: Tripti Mishra, Shipra Shukla, Sanjeev Meena, , Ruchi Singh, Mahesh Pal, D. K. Upreti, Dipak Datta

Abstract:

Roccella montagnei belongs to lichen family Roccelleceae growing luxuriantly along the coastal regions of India. As Roccella has been shown to be bioactive, we prepared methanolic extract and assessed its anticancer potential. The methanolic extract showed significant in vitro cytotoxic activity against four human cancer cell lines such as Colon (DLD-1, SW-620), Breast (MCF-7), Head and Neck (FaDu). This prompted us to isolate bioactive compounds through column chromatography. Two compounds Roccellic acid and Everninic acid have been isolated, out of which Everninic acid is reported for the first time. Both the compounds have been tested for in vitro cytotoxic activity in which Roccellic acid showed strong anticancer activity as compared to the Everninic acid. CDK-10 (Cyclin-dependent kinase) contributes to proliferation of cancer cells, and aberrant activity of these kinases has been reported in a wide variety of human cancers. These kinases, therefore, constitute biomarkers of proliferation and attractive pharmacological targets for the development of anticancer therapeutics. Therefore both the isolated compounds were tested for in silico molecular docking study against CDK-10 isomer enzyme to support the cytotoxic activity.

Keywords: cytotoxic activity, everninic acid, roccellic acid, R. montagnei

Procedia PDF Downloads 306

Authors: Fathi Soliman

Abstract:

With the ultimate goal of developing the separation of n-paraffin as phase change material (PCM) by means of molecular dynamic simulations, we attempt to predict the solubility of aromatic n-paraffin in two organic solvents: Butyl Acetate (BA) and Methyl Iso Butyl Ketone (MIBK). A simple model of aromatic paraffin: 2-hexadecylantharacene with amorphous molecular structure and periodic boundary conditions was constructed. The results showed that MIBK is the best solvent to separate ultra-pure phase change materials and this data was compatible with experimental data done to separate ultra-pure n-paraffin from waste petroleum aromatic paraffin wax, the separated n-paraffin was characterized by XRD, TGA, GC and DSC, moreover; data revealed that the n-paraffin separated by using MIBK is better as PCM than that separated using BA.

Keywords: molecular dynamics simulation, n-paraffin, organic solvents, phase change materials, solvent extraction

Procedia PDF Downloads 174

Authors: Jineetkumar Gawad, Chandrakant Bonde

Abstract:

Tuberculosis (TB) is a major worldwide concern whose control has been exacerbated by HIV, the rise of multidrug-resistance (MDR-TB) and extensively drug resistance (XDR-TB) strains of Mycobacterium tuberculosis. The interest for newer and faster acting antitubercular drugs are more remarkable than any time. To search potent compounds is need and challenge for researchers. Here, we tried to design lead for inhibition of Decaprenyl phosphoryl-β-D-ribose 2′-epimerase (DprE1) enzyme. Arabinose is an essential constituent of mycobacterial cell wall. DprE1 is a flavoenzyme that converts decaprenylphosphoryl-D-ribose into decaprenylphosphoryl-2-keto-ribose, which is intermediate in biosynthetic pathway of arabinose. Latter, DprE2 converts keto-ribose into decaprenylphosphoryl-D-arabinose. We had a selection of 23 compounds from azaindole series for computational study, and they were drawn using marvisketch. Ligands were prepared using Maestro molecular modeling interface, Schrodinger, v10.5. Common pharmacophore hypotheses were developed by applying dataset thresholds to yield active and inactive set of compounds. There were 326 hypotheses were developed. On the basis of survival score, ADRRR (Survival Score: 5.453) was selected. Selected pharmacophore hypotheses were subjected to virtual screening results into 1000 hits. Hits were prepared and docked with protein 4KW5 (oxydoreductase inhibitor) was downloaded in .pdb format from RCSB Protein Data Bank. Protein was prepared using protein preparation wizard. Protein was preprocessed, the workspace was analyzed using force field OPLS 2005. Glide grid was generated by picking single atom in molecule. Prepared ligands were docked with prepared protein 4KW5 using Glide docking. After docking, on the basis of glide score top-five compounds were selected, (5223, 5812, 0661, 0662, and 2945) and the glide docking score (-8.928, -8.534, -8.412, -8.411, -8.351) respectively. There were interactions of ligand and protein, specifically HIS 132, LYS 418, TRY 230, ASN 385. Pi-pi stacking was observed in few compounds with basic Imidazo [4,5-b] pyridine ring. We had basic azaindole ring in parent compounds, but after glide docking, we received compounds with Imidazo [4,5-b] pyridine as a basic ring. That might be the new lead in the process of drug discovery.

Keywords: DprE1 inhibitors, in silico drug designing, imidazo [4, 5-b] pyridine, lead, tuberculosis

Procedia PDF Downloads 132

Authors: S. Reshma Reghu, Shiburaj Sugathan, T. G. Nandu, K. B. Ramesh Kumar, Mathew Dan

Abstract:

Bacterial diseases are considered to be one of the most prevalent health hazards in the developing world and many bacteria are becoming resistant to existing antibiotics making the treatment ineffective. Thus, it is necessary to find novel targets and develop new antibacterial drugs with a novel mechanism of action. The process of bacterial cell division is a novel and attractive target for new antibacterial drug discovery. FtsZ, a homolog of eukaryotic tubulin, is the major protein of the bacterial cell division machinery and is considered as an important antibacterial drug target. Zingiberaceae, the Ginger family consists of aromatic herbs with creeping rhizomes. Many of these plants have antimicrobial properties.This study aimed to determine the anti-bacterial activity of selected Zingiberaceous plants by targeting bacterial cell division protein, FtsZ. Essential oils and methanol extracts of Amomum ghaticum, Alpinia galanga, Kaempferia galanga, K. rotunda, and Zingiber officinale were tested to find its antibacterial efficiency using disc diffusion method against authentic bacterial strains obtained from MTCC (India). Essential oil isolated from A.galanga and Z.officinale were further assayed for FtsZ inhibition assay following non-radioactive malachite green-phosphomolybdate assay using E. coli FtsZ protein obtained from Cytoskelton Inc., USA. Z.officinale essential oil possess FtsZ inhibitory property. A molecular docking study was conducted with the known bioactive compounds of Z. officinale as ligands with the E. coli FtsZ protein homology model. Some of the major constituents of this plant like catechin, epicatechin, and gingerol possess agreeable docking scores. The results of this study revealed that several chemical constituents in Ginger plants can be utilised as potential source of antibacterial activity and it can warrant further investigation through drug discovery studies.

Keywords: antibacterial, FtsZ, zingiberaceae, docking

Procedia PDF Downloads 453

Authors: Wen-Shyong Tzou, Chih-Ching Huang, Chin-Hwa Hu, Ying-Tsang Lo, Tun-Wen Pai, Chia-Yin Chiang, Chung-Hao Li, Hong-Jyuan Jian

Abstract:

Plasmin plays an important role in the human circulatory system owing to its catalytic ability of fibrinolysis. The immediate injection of plasmin in patients of strokes has intrigued many scientists to design vectors that can transport plasmin to the desired location in human body. Here we predict the structure of human plasmin and investigate the interaction of plasmin with the gold-nanoparticle. Because the crystal structure of plasminogen has been solved, we deleted N-terminal domain (Pan-apple domain) of plasminogen and generate a mimic of the active form of this enzyme (plasmin). We conducted a simulated annealing process on plasmin and discovered a very large conformation occurs. Kringle domains 1, 4 and 5 had been observed to leave its original location relative to the main body of the enzyme and the original doughnut shape of this enzyme has been transformed to a V-shaped by opening its two arms. This observation of conformational change is consistent with the experimental results of neutron scattering and centrifugation. We subsequently docked the plasmin on the simulated gold surface to predict their interaction. The V-shaped plasmin could utilize its Kringle domain and catalytic domain to contact the gold surface. Our findings not only reveal the flexibility of plasmin structure but also provide a guide for the design of a plasmin-gold nanoparticle.

Keywords: docking, gold nanoparticle, molecular simulation, plasmin

Procedia PDF Downloads 462

Authors: Amit Chaudhary, B. S. Yadav, P. K. Maurya, A. M., S. Srivastava, S. Singh, A. Mani

Abstract:

Cold shock protein A (CspA) is major cold inducible protein present in Escherichia coli. The protein is involved in stabilizing secondary structure of RNA by working as chaperone during cold temperature. Two RNA binding motifs play key role in the stabilizing activity. This study aimed to investigate implications of low temperature on structure and RNA binding activity of E. coli CspA. Molecular dynamics simulations were performed to compare the stability of the protein at 37°C and 10 °C. The protein was mutated at RNA binding motifs and docked with RNA to assess the stability of both complexes. Results suggest that CspA as well as CspA-RNA complex is more stable at low temperature. It was also confirmed that RNP1 and RNP2 play key role in RNA binding.

Keywords: CspA, homology modelling, mutation, molecular dynamics simulation

Procedia PDF Downloads 358

Authors: R. Chanajaree, D. Luanwiset, K. Pongpratea

Abstract:

Dye removal is an environmental concern because the textile industries have been increasing by world population and industrialization. Adsorption is the technique to find adsorbents to remove dyes from wastewater. This method is low-cost and effective for dye removal. This work tries to develop effective adsorbents using the computational approach because it will be able to predict the possibility of the adsorbents for specific dyes in terms of binding free energies. The computational approach is faster and cheaper than the experimental approach in case of finding the best adsorbents. All starting structures of dyes and adsorbents are optimized by quantum calculation. The complexes between dyes and adsorbents are generated by the docking method. The obtained binding free energies from docking are compared to binding free energies from the experimental data. The calculated energies can be ranked as same as the experimental results. In addition, this work also shows the possible orientation of the complexes. This work used two experimental groups of the complexes of the dyes and adsorbents. In the first group, there are chitosan (adsorbent) and two dyes (reactive red (RR) and direct sun yellow (DY)). In the second group, there are poly(1,2-epoxy-3-phenoxy) propane (PEPP), which is the adsorbent, and 2 dyes of bromocresol green (BCG) and alizarin yellow (AY).

Keywords: dyes removal, binding free energies, quantum calculation, docking

Procedia PDF Downloads 133

Authors: Jiahui Song

Abstract:

The use of electric fields with high intensity (~ 100kV/cm or higher) and ultra short pulse durations (nanosecond range) has been a recent development. Most of the studies of electroporation have ignored possible thermal effects because of the small duration of the applied voltage pulses. However, it has been predicted membrane temperature gradients ranging from 0.2×109 to 109 K/m. This research focuses on thermal effects that drive for electroporative enhancements, even though the actual temperature values might not have changed appreciably from their equilibrium levels. The dynamics of pore formation with the application of an externally applied electric field is studied on the basis of molecular dynamics (MD) simulations using the GROMACS package. MD simulations of a lipid layer with constant electric field strength of 0.5 V/nm at 25 °C and 47 °C are implemented to simulate the appropriate thermal effects. The GROMACS provides the force fields for the lipid membranes, which is taken to comprise of dipalmitoyl-phosphatidyl-choline (DPPC) molecules. The water model mimicks the aqueous environment surrounding the membrane. Velocities of water and membrane molecules are generated randomly at each simulation run according to a Maxwellian distribution. The high background electric field is typically used in MD simulations to probe electroporation. It serves as an accelerated test of the pore formation process since low electric fields would take inordinately long simulation time. MD simulation shows no pore is formed in a 1-ns snapshot for a DPPC membrane set at a temperature of 25°C after a 0.5 V/nm electric field is applied. A nano-sized pore is clearly seen in a 0.75-ns snapshot on the same geometry, but with the membrane surfaces kept at temperatures of 47°C. And the pore increases at 1 ns. The MD simulation results suggest the possibility that the increase in temperature can result in different degrees of electrically stimulated bio-effects. The results points to the role of thermal effects in facilitating and accelerating the electroporation process.

Keywords: high-intensity, ultra-short, electroporation, thermal effects, molecular dynamics

Procedia PDF Downloads 35

Authors: Narendra Nath Patra

Abstract:

Several recent studies reinforce the existence of a thick molecular disc in galaxies along with the dynamically cold thin disc. Assuming a two-component molecular disc, we model the disc of NGC 6946 as a four-component system consists of stars, HI, thin disc molecular gas, and thick disc molecular gas in vertical hydrostatic equilibrium. Following, we set up the joint Poisson-Boltzmann equation of hydrostatic equilibrium and solve it numerically to obtain a three-dimensional density distribution of different baryonic components. Using the density solutions and the observed rotation curve, we further build a three-dimensional dynamical model of the molecular disc and consecutively produce simulated CO spectral cube and spectral width profile. We find that the simulated spectral width profiles distinguishably differs for different assumed thick disc molecular gas fraction. Several CO spectral width profiles are then produced for different assumed thick disc molecular gas fractions and compared with the observed one to obtain the best fit thick disc molecular gas fraction profile. We find that the thick disc molecular gas fraction in NGC 6946 largely remains constant across its molecular disc with a mean value of 0.70 +/- 0.09. We also estimate the amount of extra-planar molecular gas in NGC 6946. We find 60% of the total molecular gas is extra-planar at the central region, whereas this fraction reduces to ~ 35% at the edge of the molecular disc. With our method, for the first time, we estimate the thick disc molecular gas fraction as a function of radius in an external galaxy with sub-kpc resolution.

Keywords: galaxies: kinematics and dynamic, galaxies: spiral, galaxies: structure , ISM: molecules, molecular data

Procedia PDF Downloads 121

Authors: Javier Varillas, Jorge Alcalá

Abstract:

We investigate dislocation avalanche phenomena in face-centered cubic (FCC) and body-centered cubic (BCC) crystals using massive, large-scale molecular dynamics (MD) simulations. The analysis is focused on the intermittent development of dense dislocation arrangements subjected to uniaxial tensile straining under displacement control. We employ a novel computational scheme that allows us to inject an entangled dislocation structure in periodic MD domains. We assess the emission of plastic bursts (or dislocation avalanches) in terms of the sharp stress drops detected in the stress-strain curve. The plastic activity corresponds to the sporadic operation of specific dislocation glide processes exhibiting quiescent periods between successive avalanche events. We find that the plastic intermittences in our simulations do not overlap in time under sufficiently low strain rates as dissipation operates faster than driving, where the dense dislocation networks evolve through the emission of dislocation avalanche events whose carried slip adheres to self-organized power-law distributions. These findings enable the extension of the slip distributions obtained from strict displacement-controlled micropillar compression experiments towards smaller values of slip size. Our results furnish further understanding upon the development of entangled dislocation networks in metal plasticity, including specific mechanisms of dislocation propagation and annihilation, along with the evolution of specific dislocation populations through dislocation density analyses.

Keywords: dislocations, intermittent plasticity, molecular dynamics, slip distributions

Procedia PDF Downloads 123

Authors: Rong Li, Brian D. Wirth, Bing Liu

Abstract:

Graphite is the matrix, and structural material in the high temperature gas-cooled reactor exhibits an irradiation response. It is of significant importance to analyze the defect production and evaluate the role of graphite under irradiation. A vast experimental literature exists for graphite on the dimensional change, mechanical properties, and thermal behavior. However, simulations have not been applied to the atomistic perspective. Remarkably few molecular dynamics simulations have been performed to study the irradiation response in graphite. In this paper, irradiation-induced damage cascades in graphite were investigated with molecular dynamics simulation. Statistical results of the graphite defects were obtained by sampling a wide energy range (1–30 KeV) and 10 different runs for every cascade simulation with different random number generator seeds to the velocity scaling thermostat function. The chemical bonding in carbon was described using the adaptive intermolecular reactive empirical bond-order potential (AIREBO) potential coupled with the standard Ziegler–Biersack–Littmack (ZBL) potential to describe close-range pair interactions. This study focused on analyzing the number of defects, the final cascade morphology and the distribution of defect clusters in space, the length-scale cascade properties such as the cascade length and the range of primary knock-on atom (PKA), and graphite mechanical properties’ variation. It can be concluded that the number of surviving Frenkel pairs increased remarkably with the increasing initial PKA energy but did not exhibit a thermal spike at slightly lower energies in this paper. The PKA range and cascade length approximately linearly with energy which indicated that increasing the PKA initial energy will come at expensive computation cost such as 30KeV in this study. The cascade morphology and the distribution of defect clusters in space mainly related to the PKA energy meanwhile the temperature effect was relatively negligible. The simulations are in agreement with known experimental results and the Kinchin-Pease model, which can help to understand the graphite damage cascades and lifetime span under irradiation and provide a direction to the designs of these kinds of structural materials in the future reactors.

Keywords: graphite damage cascade, molecular dynamics, cascade morphology, cascade distribution

Procedia PDF Downloads 141

Authors: Nidhi Chadha, A. K. Tiwari, Marilyn D. Milton, Anil K. Mishra

Abstract:

Translocator protein (18 kDa) TSPO is highly expressed during microglia activation in neuroinflammation. Although a number of PET ligands have been developed for the visualization of activated microglia, one of the advantageous approaches is to develop potential optical imaging (OI) probe. Our study involves computational screening, synthesis and evaluation of TSPO ligand through various imaging modalities namely PET/SPECT/Optical. The initial computational screening involves pharmacophore modeling from the library designing having oxo-benzooxazol-3-yl-N-phenyl-acetamide groups and synthesis for visualization of efficacy of these compounds as multimodal imaging probes. Structure modeling of monomer, Ala147Thr mutated, parallel and anti-parallel TSPO dimers was performed and docking analysis was performed for distinct binding sites. Computational analysis showed pattern of variable binding profile of known diagnostic ligands and NBMP via interactions with conserved residues along with TSPO’s natural polymorphism of Ala147→Thr, which showed alteration in the binding affinity due to considerable changes in tertiary structure. Preliminary in vitro binding studies shows binding affinity in the range of 1-5 nm and selectivity was also certified by blocking studies. In summary, this skeleton was found to be potential probe for TSPO imaging due to ease in synthesis, appropriate lipophilicity and reach to specific region of brain.

Keywords: TSPO, molecular modeling, imaging, docking

Procedia PDF Downloads 441

Authors: Apisaraporn Baicharoen, Prapasiri Pongprayoon

Abstract:

Betaine aldehyde dehydrogenase 2 (BADH2) is an enzyme that inhibits the accumulation of 2-acetyl-1-pyrroline (2AP), a potent flavor compound in rice fragrance. BADH2 contains three domains (NAD-binding, substrate-binding, and oligomerization domains). It catalyzes the oxidation of amino aldehydes. The lack of BADH2 results in the formation of 2AP and consequently an increase in rice fragrance. To date, inadequate data on BADH2 structure and function are available. An insight into the nature of BADH2 can serve as one of key starting points for the production of high quality fragrant rice. In this study, we therefore constructed the homology model of BADH2 and employed 500-ns Molecular Dynamics simulations (MD) to primarily understand the structural and dynamic properties of BADH2. Initially, Ramachandran plot confirms the good quality of modeled protein structure. Principle Component Analysis (PCA) was also calculated to capture the protein dynamics. Among 3 domains, the results show that NAD binding site is found to be more flexible. Moreover, interactions from key amino acids (N162, E260, C294, and Y419) that are crucial for function are investigated.

Keywords: betaine aldehyde dehydrogenase 2, fragrant rice, homology modeling, molecular dynamics simulations

Procedia PDF Downloads 193

Authors: R. Pilafkan, M. Kaffash Irzarahimi, S. F. Asbaghian Namin

Abstract:

The biaxial buckling behavior of single-layered graphene sheets (SLGSs) is studied in the present work. To consider the size-effects in the analysis, Eringen’s nonlocal elasticity equations are incorporated into classical plate theory (CLPT). A Generalized Differential Quadrature Method (GDQM) approach is utilized and numerical solutions for the critical buckling loads are obtained. Then, molecular dynamics (MD) simulations are performed for a series of zigzag SLGSs with different side-lengths and with various boundary conditions, the results of which are matched with those obtained by the nonlocal plate model to numerical the appropriate values of nonlocal parameter relevant to each type of boundary conditions.

Keywords: biaxial buckling, single-layered graphene sheets, nonlocal elasticity, molecular dynamics simulation, classical plate theory

Procedia PDF Downloads 263

Authors: Jeoungsu Na, Jaehawn Lee, Changil Hong, Suhee Kim

Abstract:

A molecular dynamics simulation in a converging-diverging nozzle was performed to study molecular collisions and their influence to average flow velocity according to a variety of vacuum levels. The static pressures and the dynamic pressure exerted by the molecule collision on the selected walls were compared to figure out the intensity variances of the directional flows. With pressure differences constant between the entrance and the exit of the nozzle, the numerical experiment was performed for molecular velocities and directional flows. The result shows that the velocities increased at the nozzle exit as the vacuum level gets higher in that area because less molecular collisions.

Keywords: cavitation, molecular collision, nozzle, vacuum, velocity increase

Procedia PDF Downloads 418

Authors: Mahmoud A. Abdulhamid

Abstract:

Poly(ether-ether-ketone) (PEEK) has received increased attention due to its outstanding performance in different membrane applications including gas and liquid separation. However, it suffers from a semi-crystalline morphology, bad solubility and low porosity. To fabricate membranes from PEEK, the usage of harsh acid such as sulfuric acid is essential, regardless its hazardous properties. In this work, we report the molecular design of poly(ether-ether-ketones) (iPEEKs) with intrinsic porosity character, by incorporating kinked units into PEEK backbone such as spirobisindane, Tröger's base, and triptycene. The porous polymers were used to fabricate stable membranes for organic solvent nanofiltration application. To better understand the mechanism, we conducted molecular dynamics simulations to evaluate the possible interactions between the polymers and the solvents. Notable enhancement in separation performance was observed confirming the importance of molecular engineering of high-performance polymers. The iPEEKs demonstrated good solubility in polar aprotic solvents, a high surface area of 205–250 m² g⁻¹, and excellent thermal stability. Mechanically flexible nanofiltration membranes were prepared from N-methyl-2-pyrrolidone dope solution at iPEEK concentrations of 19–35 wt%. The molecular weight cutoff of the membranes was fine-tuned in the range of 450–845 g mol⁻¹ displaying 2–6 fold higher permeance (3.57–11.09 L m⁻² h⁻¹ bar⁻¹) than previous reports. The long-term stabilities were demonstrated by a 7 day continuous cross-flow filtration.

Keywords: molecular engineering, polymer synthesis, membrane fabrication, liquid separation

Procedia PDF Downloads 82

Authors: Sundaram Arulmozhiraja, Kanade Shimizu, Yuta Yamamoto, Satoshi Ichikawa, Maenaka Katsumi, Hiroaki Tokiwa

Abstract:

Drug discovery is shifting from small molecule based drugs targeting local active site to middle molecules (MM) targeting large, flat, and groove-shaped binding sites, for example, protein-protein interface because at least half of all targets assumed to be involved in human disease have been classified as “difficult to drug” with traditional small molecules. Hence, MMs such as peptides, natural products, glycans, nucleic acids with various high potent bioactivities become important targets for drug discovery programs in the recent years as they could be used for ‘undruggable” intracellular targets. Cell membrane permeability is one of the key properties of pharmacodynamically active MM drug compounds and so evaluating this property for the potential MMs is crucial. Computational prediction for cell membrane permeability of molecules is very challenging; however, recent advancement in the molecular dynamics simulations help to solve this issue partially. It is expected that MMs with high membrane permeability will enable drug discovery research to expand its borders towards intracellular targets. Further to understand the chemistry behind the permeability of MMs, it is necessary to investigate their conformational changes during the permeation through membrane and for that their interactions with the membrane field should be studied reliably because these interactions involve various non-bonding interactions such as hydrogen bonding, -stacking, charge-transfer, polarization dispersion, and non-classical weak hydrogen bonding. Therefore, parameters-based classical mechanics calculations are hardly sufficient to investigate these interactions rather, quantum mechanical (QM) calculations are essential. Fragment molecular orbital (FMO) method could be used for such purpose as it performs ab initio QM calculations by dividing the system into fragments. The present work is aimed to study the cell permeability of middle molecules using molecular dynamics simulations and FMO-QM calculations. For this purpose, a natural compound syringolin and its analogues were considered in this study. Molecular simulations were performed using NAMD and Gromacs programs with CHARMM force field. FMO calculations were performed using the PAICS program at the correlated Resolution-of-Identity second-order Moller Plesset (RI-MP2) level with the cc-pVDZ basis set. The simulations clearly show that while syringolin could not permeate the membrane, its selected analogues go through the medium in nano second scale. These correlates well with the existing experimental evidences that these syringolin analogues are membrane-permeable compounds. Further analyses indicate that intramolecular -stacking interactions in the syringolin analogues influenced their permeability positively. These intramolecular interactions reduce the polarity of these analogues so that they could permeate the lipophilic cell membrane. Conclusively, the cell membrane permeability of various middle molecules with potent bioactivities is efficiently studied using molecular dynamics simulations. Insight of this behavior is thoroughly investigated using FMO-QM calculations. Results obtained in the present study indicate that non-bonding intramolecular interactions such as hydrogen-bonding and -stacking along with the conformational flexibility of MMs are essential for amicable membrane permeation. These results are interesting and are nice example for this theoretical calculation approach that could be used to study the permeability of other middle molecules. This work was supported by Japan Agency for Medical Research and Development (AMED) under Grant Number 18ae0101047.

Keywords: fragment molecular orbital theory, membrane permeability, middle molecules, molecular dynamics simulation

Procedia PDF Downloads 169

Authors: Hamza Javar Magnier, Leslie V. Woodcock

Abstract:

It has been reported that at temperatures above the critical there is no “continuity of liquid and gas”, as originally hypothesized by van der Waals. Rather, both gas and liquid phases, with characteristic properties as such, extend to supercritical temperatures. Each phase is bounded by the locus of a percolation transition, i.e. a higher-order thermodynamic phase change associated with percolation of gas clusters in a large void, or liquid interstitial vacancies in a large cluster. Between these two-phase bounds, it is reported there exists a mesophase that resembles an otherwise homogeneous dispersion of gas micro-bubbles in liquid (foam) and a dispersion of liquid micro-droplets in gas (mist). Such a colloidal-like state of a pure one-component fluid represents a hitherto unchartered equilibrium state of matter besides pure solid, liquid or gas. Here we provide compelling evidence, from molecular dynamics (MD) simulations, for the existence of this supercritical mesophase and its colloidal nature. We report preliminary results of computer simulations for a model fluid using a simplistic representation of atoms or molecules, i.e. a hard-core repulsion with an attraction so short that the atoms are referred to as “adhesive spheres”. Molecular clusters, and hence percolation transitions, are unambiguously defined. Graphics of color-coded clusters show colloidal characteristics of the supercritical mesophase.

Keywords: critical phenomena, mesophase, supercritical, square-well, critical parameters

Procedia PDF Downloads 413