Search results for: ab initio molecular dynamics

Authors: Nikhil Agrawal, Adam A. Skelton

Abstract:

Alzheimer’s disease (AD) is one of the most common forms of dementia, which is caused by misfolding and aggregation of amyloid beta (Aβ) peptides into amyloid-β fibrils (Aβ fibrils). To disrupt the remodeling of Aβ fibrils, a number of candidate molecules have been proposed. To study the molecular mechanisms of Aβ fibrils remodeling we performed a series of all-atom molecular dynamics simulations, a total time of 3µs, in explicit solvent. Several previously undiscovered candidate molecule-Aβ fibrils binding modes are unraveled; one of which shows the direct conformational change of the Aβ fibril by understanding the physicochemical factors responsible for binding and subsequent remodeling of Aβ fibrils by the candidate molecule, open avenues into structure-based drug design for AD can be opened.

Keywords: alzheimer’s disease, amyloid, MD simulations, misfolded protein

Procedia PDF Downloads 347

Authors: Weirong Wang, Shenghong Huang, Xisheng Luo, Zhenyu Li

Abstract:

Material mixing process and related dynamic issues at extreme compressing conditions have gained more and more concerns in last ten years because of the engineering appealings in inertial confinement fusion (ICF) and hypervelocity aircraft developments. However, there lacks models and methods that can handle fully coupled turbulent material mixing and complex fluid evolution under conditions of high energy density regime up to now. In aspects of macro hydrodynamics, three numerical methods such as direct numerical simulation (DNS), large eddy simulation (LES) and Reynolds-averaged Navier–Stokes equations (RANS) has obtained relative acceptable consensus under the conditions of low energy density regime. However, under the conditions of high energy density regime, they can not be applied directly due to occurrence of dissociation, ionization, dramatic change of equation of state, thermodynamic properties etc., which may make the governing equations invalid in some coupled situations. However, in view of micro/meso scale regime, the methods based on Molecular Dynamics (MD) as well as Monte Carlo (MC) model are proved to be promising and effective ways to investigate such issues. In this study, both classical MD and first-principle based electron force field MD (eFF-MD) methods are applied to investigate Richtmyer-Meshkov Instability of metal Lithium and gas Hydrogen (Li-H2) interface mixing at different shock loading speed ranging from 3 km/s to 30 km/s. It is found that: 1) Classical MD method based on predefined potential functions has some limits in application to extreme conditions, since it cannot simulate the ionization process and its potential functions are not suitable to all conditions, while the eFF-MD method can correctly simulate the ionization process due to its ‘ab initio’ feature; 2) Due to computational cost, the eFF-MD results are also influenced by simulation domain dimensions, boundary conditions and relaxation time choices, etc., in computations. Series of tests have been conducted to determine the optimized parameters. 3) Ionization induced by strong shock compression has important effects on Li-H2 interface evolutions of RMI, indicating a new micromechanism of RMI under conditions of high energy density regime.

Keywords: first-principle, ionization, molecular dynamics, material mixture, Richtmyer-Meshkov instability

Procedia PDF Downloads 225

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 329

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 115

Authors: Maria Pintea, Nigel Mason

Abstract:

Induced chemistry is one of the newest pathways in the nanotechnology field with applications in the focused electron beam induced processes for deposition of nm scale structures. Si(OPr)₄ and Ti(OEt)₄ are two of the precursors that have not been so extensively researched, though highly sought for semiconductor and medical applications fields, the two compounds make good candidates for FEBIP and are the subject of velocity slice map imaging analysis for deposition purposes, offering information on kinetic energies, fragmentation channels, and angular distributions. The velocity slice map imaging technique is a method used for the characterization of molecular dynamics of the molecule and the fragmentation channels as a result of induced chemistry. To support the gas-phase analysis, Meso-Bio-Nano simulations of irradiation dynamics studies are employed with final results on Fe(CO)₅ deposited on various substrates. The software is capable of running large scale simulations for complex biomolecular, nano- and mesoscopic systems with applications to thermos-mechanical DNA damage, complex materials, gases, nanoparticles for cancer research and deposition applications for nanotechnology, using a large library of classical potentials, many-body force fields, molecular force fields involved in the classical molecular dynamics.

Keywords: focused electron beam induced deposition, FEBID, induced chemistry, molecular dynamics, velocity map slice imaging

Procedia PDF Downloads 109

Authors: Ruqaiya Khalil, Saman Usmani, Zaheer Ul-Haq

Abstract:

The accurate characterization of ligand binding affinity is indispensable for designing molecules with optimized binding affinity. Computational tools help in many directions to predict quantitative correlations between protein-ligand structure and their binding affinities. Molecular dynamics (MD) simulation is a modern state-of-the-art technique to evaluate the underlying basis of ligand-protein interactions by characterizing dynamic and energetic properties during the event. Autoimmune diseases arise from an abnormal immune response of the body against own tissues. The current regimen for the described condition is limited to immune-modulators having compromised pharmacodynamics and pharmacokinetics profiles. One of the key player mediating immunity and tolerance, thus invoking autoimmunity is Interleukin-2; a cytokine influencing the growth of T cells. Molecular dynamics simulation techniques are applied to seek insight into the inhibitory mechanisms of newly synthesized compounds that manifested immunosuppressant potentials during in silico pipeline. In addition to estimation of free energies associated with ligand binding, MD simulation yielded us a great deal of information about ligand-macromolecule interactions to evaluate the pattern of interactions and the molecular basis of inhibition. The present study is a continuum of our efforts to identify interleukin-2 inhibitors of both natural and synthetic origin. Herein, we report molecular dynamics simulation studies of Interluekin-2 complexed with different antagonists previously reported by our group. The study of protein-ligand dynamics enabled us to gain a better understanding of the contribution of different active site residues in ligand binding. The results of the study will be used as the guide to rationalize the fragment based synthesis of drug-like interleukin-2 inhibitors as immune-modulators.

Keywords: immuno-modulators, MD simulation, protein-ligand interaction, structure-based drug design

Procedia PDF Downloads 262

Authors: Said Belaaouad

Abstract:

This study aimed to explore the quantitative structure-activity relationship (QSAR) of 1,2,4-Triazine-3(2H)-one derivative as a potential anticancer agent against breast cancer. The electronic descriptors were obtained using the Density Functional Theory (DFT) method, and a multiple linear regression techniques was employed to construct the QSAR model. The model exhibited favorable statistical parameters, including R2=0.849, R2adj=0.656, MSE=0.056, R2test=0.710, and Q2cv=0.542, indicating its reliability. Among the descriptors analyzed, absolute electronegativity (χ), total energy (TE), number of hydrogen bond donors (NHD), water solubility (LogS), and shape coefficient (I) were identified as influential factors. Furthermore, leveraging the validated QSAR model, new derivatives of 1,2,4-Triazine-3(2H)-one were designed, and their activity and pharmacokinetic properties were estimated. Subsequently, molecular docking (MD) and molecular dynamics (MD) simulations were employed to assess the binding affinity of the designed molecules. The Tubulin colchicine binding site, which plays a crucial role in cancer treatment, was chosen as the target protein. Through the simulation trajectory spanning 100 ns, the binding affinity was calculated using the MMPBSA script. As a result, fourteen novel Tubulin-colchicine inhibitors with promising pharmacokinetic characteristics were identified. Overall, this study provides valuable insights into the QSAR of 1,2,4-Triazine-3(2H)-one derivative as potential anticancer agent, along with the design of new compounds and their assessment through molecular docking and dynamics simulations targeting the Tubulin-colchicine binding site.

Keywords: QSAR, molecular docking, ADMET, 1, 2, 4-triazin-3(2H)-ones, breast cancer, anticancer, molecular dynamic simulations, MMPBSA calculation

Procedia PDF Downloads 97

Authors: Giorgio Visentin, Alexei A. Buchachenko

Abstract:

Recent proposals to use Yb dimer as an optical clock and as a sensor for non-Newtonian gravity imply the knowledge of its interaction potential. Here, the ground-state Born-Oppenheimer Yb₂ potential energy curve is represented by a semi-analytical function, consisting of short- and long-range contributions. For the former, the systematic ab initio all-electron exact 2-component scalar-relativistic CCSD(T) calculations are carried out. Special care is taken to saturate diffuse basis set component with the atom- and bond-centered primitives and reach the complete basis set limit through n = D, T, Q sequence of the correlation-consistent polarized n-zeta basis sets. Similar approaches are used to the long-range dipole and quadrupole dispersion terms by implementing the CCSD(3) polarization propagator method for dynamic polarizabilities. Dispersion coefficients are then computed through Casimir-Polder integration. The semiclassical constraint on the number of the bound vibrational levels known for the ¹⁷⁴Yb isotope is used to scale the potential function. The scaling, based on the most accurate ab initio results, bounds the interaction energy of two Yb atoms within the narrow 734 ± 4 cm⁻¹ range, in reasonable agreement with the previous ab initio-based estimations. The resulting potentials can be used as the reference for more sophisticated models that go beyond the Born-Oppenheimer approximation and provide the means of their uncertainty estimations. The work is supported by Russian Science Foundation grant # 17-13-01466.

Keywords: ab initio coupled cluster methods, interaction potential, semi-analytical function, ytterbium dimer

Procedia PDF Downloads 152

Authors: Shagufta Tabassum

Abstract:

The study of dielectric properties in a binary mixture of liquids is very useful to understand the liquid structure, molecular interaction, dynamics, and kinematics of the mixture. Time-domain reflectometry (TDR) is a powerful tool for studying the cooperation and molecular dynamics of the H-bonded system. In this paper, we discuss the basic calibration and normalization procedure for time-domain reflectometry measurements. Our approach is to explain the different types of error occur during TDR measurements and how these errors can be eliminated or minimized.

Keywords: time domain reflectometry measurement techinque, cable and connector loss, oscilloscope loss, and normalization technique

Procedia PDF Downloads 206

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 188

Authors: Ncoza Dlova, Tivani Mashamba-Thompson

Abstract:

Garcinia livingstonei (Clusiaceae) extracts, morelloflavone- 7″- sulphate and sargaol were shown to be effective against hyper-pigmentation through inhibition of tyrosinase enzyme, in vitro . The aim of this study is to elucidate the structural mechanism through which morelloflavone- 7″- sulphate and sargaol binds human tyrosinase. Implementing a homology model to construct a tyrosinase model using the crystal structure of a functional unit from Octopus hemocyanin (PDB: 1JS8) as a reference template enabled us to create a human tyrosinase model. Molecular dynamics and binding free energy calculations were optimized to enable molecular dynamics simulation of the copper dependent inhibitors. Results show the importance of the hydrogen bond formation morelloflavone- 7″- sulphate and sargaol between compound and active site residues. Both complexes demonstrated the metallic coordination between compound and arginine residue as well as copper ions within the active site. The comprehensive molecular insight gained from this study should be vital in understanding the binding mechanism morelloflavone- 7″- sulphate and sargaol. Moreover, these results will assist in the design of novel of metal ion dependent enzyme inhibitors as potential anti-hyper-pigmentation disorder therapies.

Keywords: hyper-pigmentation disorders, dyschromia African skin, morelloflavone- 7″- sulphate, sagoal

Procedia PDF Downloads 406

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 195

Authors: Zeeshan Ahmed, Ajinkya Sarode, Pratik Basarkar, Atul Bhargav, Debjyoti Banerjee

Abstract:

The use of CO₂ in oil recovery and in CO₂ capture and storage is gaining traction in recent years. These applications involve heat transfer between CO₂ and the base fluid, and hence, there arises a need to improve the thermal conductivity of CO₂ to increase the process efficiency and reduce cost. One way to improve the thermal conductivity is through nanoparticle addition in the base fluid. The nanofluid model in this study consisted of copper (Cu) nanoparticles in varying concentrations with CO₂ as a base fluid. No experimental data are available on thermal conductivity of CO₂ based nanofluid. Molecular dynamics (MD) simulations are an increasingly adopted tool to perform preliminary assessments of nanoparticle (NP) fluid interactions. In this study, the effect of the formation of a nanolayer (or molecular layering) at the gas-solid interface on thermal conductivity is investigated using equilibrium MD simulations by varying NP diameter and keeping the volume fraction (1.413%) of nanofluid constant to check the diameter effect of NP on the nanolayer and thermal conductivity. A dense semi-solid fluid layer was seen to be formed at the NP-gas interface, and the thickness increases with increase in particle diameter, which also moves with the NP Brownian motion. Density distribution has been done to see the effect of nanolayer, and its thickness around the NP. These findings are extremely beneficial, especially to industries employed in oil recovery as increased thermal conductivity of CO₂ will lead to enhanced oil recovery and thermal energy storage.

Keywords: copper-CO2 nanofluid, molecular dynamics simulation, molecular interfacial layer, thermal conductivity

Procedia PDF Downloads 337

Authors: Jingjing Huang, Nengwei Li, Guanghua Wei, Jiabin You, Chao Wang, Junliang Zhang

Abstract:

In this work, molecular dynamics (MD) simulation is applied to analyze the mass transport process in the cathode of proton exchange membrane fuel cell (PEMFC), of which all types of molecules situated in the cathode is considered. a reasonable and effective MD simulation process is provided, and models were built and compared using both Materials Studio and LAMMPS. The mass transport is one of the key issues in the study of proton exchange membrane fuel cells (PEMFCs). In this report, molecular dynamics (MD) simulation is applied to analyze the influence of Nafion ionomer distribution and Pt nano-particle size on mass transport process in the cathode. It is indicated by the diffusion coefficients calculation that a larger quantity of Nafion, as well as a higher equivalent weight (EW) value, will hinder the transport of oxygen. In addition, medium-sized Pt nano-particles (1.5~2nm) are more advantageous in terms of proton transport compared with other particle sizes (0.94~2.55nm) when the center-to-center distance between two Pt nano-particles is around 5 nm. Then mass transport channels are found to be formed between the hydrophobic backbone and the hydrophilic side chains of Nafion ionomer according to the radial distribution function (RDF) curves. And the morphology of these channels affected by the Pt size is believed to influence the transport of hydronium ions and, consequently the performance of PEMFC.

Keywords: cathode catalytic layer, mass transport, molecular dynamics, proton exchange membrane fuel cell

Procedia PDF Downloads 243

Authors: M. Ibrir, S. Berri, S. Lakel, D. Maouche And Y. Medkour

Abstract:

We have performed first-principle calculations of the structural, electronic and magnetic properties of KFeF3, KCoF3, KMnF3, KVF3, using full-potential linearized augmented plane-wave (FP-LAPW) scheme within GGA. Features such as the lattice constant, bulk modulus and its pressure derivative are reported. Also, we have presented our results of the band structure and the density of states. The magnetic moments of KFeF3, KCoF3, KMnF3, KVF3 compounds are in most came from the exchange-splitting of X-3d orbital.

Keywords: Ab initio calculations, electronic structure, magnetic materials

Procedia PDF Downloads 420

Authors: Jeong-Won Kang

Abstract:

Graphene, a single-atom sheet, has been considered as the most promising material for making future nanoelectromechanical systems as well as purely electrical switching with graphene transistors. Graphene-based devices have advantages in scaled-up device fabrication due to the recent progress in large area graphene growth and lithographic patterning of graphene nanostructures. Here we investigated its mechanical responses of circular graphene nanoflake under the nanoindentation using classical molecular dynamics simulations. A correlation between the load and the indentation depth was constructed. The nanoindented force in this work was applied to the center point of the circular graphene nanoflake and then, the resonance frequency could be tuned by a nanoindented depth. We found the hardening or the softening of the graphene nanoflake during its nanoindented-deflections, and such properties were recognized by the shift of the resonance frequency. The calculated mechanical parameters in the force vs deflection plot were in good agreement with previous experimental and theoretical works. This proposed schematics can detect the pressure via the deflection change or/and the resonance frequency shift, and also have great potential for versatile applications in nanoelectromechanical systems.

Keywords: graphene, pressure sensor, circular graphene nanoflake, molecular dynamics

Procedia PDF Downloads 387

Authors: M. Zemouli, K. Amara, M. Elkeurti, Y. Benallou

Abstract:

We present the molecular dynamic simulations results of the structural and dynamical properties of the zinc-blende GaBi over a wide range of temperature (300-1000) K. Our simulation where performed in the framework of the three-body Tersoff potential, which accurately reproduces the lattice constants and elastic constants of the GaBi. A good agreement was found between our calculated results and the available theoretical data of the lattice constant, the bulk modulus and the cohesive energy. Our study allows us to predict the thermodynamic properties such as the specific heat and the lattice thermal expansion. In addition, this method allows us to check its ability to predict the phase transition of this compound. In particular, the transition pressure to the rock-salt phase is calculated and the results are compared with other available works.

Keywords: Gallium compounds, molecular dynamics simulations, interatomic potential thermodynamic properties, structural phase transition

Procedia PDF Downloads 445

Authors: Mouna Baassi, Mohamed Moussaoui, Hatim Soufi, Sanchaita RajkhowaI, Ashwani Sharma, Subrata Sinha, Said Belaaouad

Abstract:

Human Immunodeficiency Virus type 1 protease (HIV-1 PR) is one of the most challenging targets of antiretroviral therapy used in the treatment of AIDS-infected people. The performance of protease inhibitors (PIs) is limited by the development of protease mutations that can promote resistance to the treatment. The current study was carried out using statistics and bioinformatics tools. A series of thirty-three compounds with known enzymatic inhibitory activities against HIV-1 protease was used in this paper to build a mathematical model relating the structure to the biological activity. These compounds were designed by software; their descriptors were computed using various tools, such as Gaussian, Chem3D, ChemSketch and MarvinSketch. Computational methods generated the best model based on its statistical parameters. The model’s applicability domain (AD) was elaborated. Furthermore, one compound has been proposed as efficient against HIV-1 protease with comparable biological activity to the existing ones; this drug candidate was evaluated using ADMET properties and Lipinski’s rule. Molecular Docking performed on Wild Type and Mutant Type HIV-1 proteases allowed the investigation of the interaction types displayed between the proteases and the ligands, Darunavir (DRV) and the new drug (ND). Molecular dynamics simulation was also used in order to investigate the complexes’ stability, allowing a comparative study of the performance of both ligands (DRV & ND). Our study suggested that the new molecule showed comparable results to that of Darunavir and may be used for further experimental studies. Our study may also be used as a pipeline to search and design new potential inhibitors of HIV-1 proteases.

Keywords: QSAR, ADMET properties, molecular docking, molecular dynamics simulation.

Procedia PDF Downloads 39

Authors: Samy Madbouly

Abstract:

The relaxation behavior of biorenewable castor oil-based polyurethane-lignin thin films synthesized in homogenous waterborne dispersions was investigated as a function of concentration at different temperatures and frequencies using broadband dielectric relaxation spectroscopy (BDRS). The molecular dynamics of the glass relaxation process and the local relaxation process of the PU-LS thin films were studied over a wide range of temperatures (-70 to 30 ℃) and frequencies (5 × 10−2 to 0.5 × 107 Hz) for different lignin concentration. Four relaxation processes have been observed namely; ?-, β-, γ-relaxations and ionic conductivity for pure castor oil-based PU and castor oil-lignin-based PU thin films at different temperatures and frequencies ranges. The Vogel-Fulcher-Tammann equation was found to be well described the temperature dependence of the characteristic relaxation times of the ?-relaxation process. However, on the other hand, the molecular dynamics of both β- and γ-relaxation processes were given by the Arrhenius equation. The incorporation of lignin into the castor oil-based PU significantly increased the glass transition temperature and primitivity of the thin films. In addition, the broadness, intensity, and molecular dynamics of the only observed ?-relaxation process were found to be strongly dependent on lignin concentration.

Keywords: castor oil, lignin, polyurethane, dielectric, dispersions

Procedia PDF Downloads 204

Authors: M. Rahimnejad, B. Vahidi, B. Ebrahimi Hoseinzadeh, F. Yazdian, P. Motamed Fath, R. Jamjah

Abstract:

In this study, we developed and simulated nano-drug delivery systems efficacy in compare to free drug prescription. Computational models can be utilized to accelerate experimental steps and control the experiments high cost. Molecular dynamics simulation (MDS), in particular NAMD was utilized to better understand the anti-cancer drug interaction with cell membrane model. Paclitaxel (PTX) and dipalmitoylphosphatidylcholine (DPPC) were selected for the drug molecule and as a natural phospholipid nanocarrier, respectively. This work focused on two important interaction parameters between molecules in terms of center of mass (COM) and van der Waals interaction energy. Furthermore, we compared the simulation results of the PTX interaction with the cell membrane and the interaction of DPPC as a nanocarrier loaded by the drug with the cell membrane. The molecular dynamic analysis resulted in low energy between the nanocarrier and the cell membrane as well as significant decrease of COM amount in the nanocarrier and the cell membrane system during the interaction. Thus, the drug vehicle showed notably better interaction with the cell membrane in compared to free drug interaction with the cell membrane.

Keywords: anti-cancer drug, center of mass, interaction energy, molecular dynamics simulation, nanocarrier

Procedia PDF Downloads 341

Authors: Salah Belaidi

Abstract:

We have carried out extensive and in-depth research to search for bioactive compounds based on Algerian plants. A selection of 50 ligands from Algerian medicinal plants. Several compounds used in herbal medicine have been drawn using Marvin Sketch software. We determined the three-dimensional structures of the ligands with the MMFF94 force field in order to prepare these ligands for molecular docking. The 3D protein structure of the SARS-CoV-2 main protease was taken from the Protein Data Bank. We used AutoDockVina software to apply molecular docking. The hydrogen atoms were added during the molecular docking process, and all the twist bonds of the ligands were added using the (ligand) module in the AutoDock software. The COVID-19 main protease (Mpro) is a key enzyme that plays a vital role in viral transcription and mediating replication, so it is a very attractive drug target for SARS-CoV-2. In this work, an evaluation was carried out on the biologically active compounds present in these selected medicinal plants as effective inhibitors of the protease enzyme of COVID-19, with an in-depth computational calculation of the molecular docking using the Autodock Vina software. The top 7 ligands: Phloroglucinol, Afzelin, Myricetin-3-O- rutinosidTricin 7-neohesperidoside, Silybin, Silychristinthat and Kaempferol are selected among the 50 molecules studied which are Algerian medicinal plants, whose selection is based on the best binding energy which is relatively low compared to the reference molecule with binding affinities of -9.3, -9.3, -9, -8.9, -8 .5, 8.3 and -8.3 kcal mol-1 respectively. Then, we analyzed the ADME properties of the best7 ligands using the web server SwissADME. Two ligands (Silybin, Silychristin) were found to be potential candidates for the discovery and design of novel drug inhibitors of the protease enzyme of SARS-CoV-2. The stability of the two ligands in complexing with the Mpro protease was validated by molecular dynamics simulation; they revealed a stable trajectory in both techniques, RMSD and RMSF, by showing molecular properties with coherent interactions in molecular dynamics simulations. Finally, we conclude that the Silybin ligand forms a more stable complex with the Mpro protease compared to the Silychristin ligand.

Keywords: COVID-19, medicinal plants, molecular docking, ADME properties, molecular dynamics

Procedia PDF Downloads 34

Authors: Tsu-Hsu Yen

Abstract:

The interfacial gas adsorption presents a frequent challenge and opportunity for micro-/nano-fluidic operation. In this study, we investigate the wettability, gas accumulation, and nanobubble formation on various homogeneous surface conditions by using MD simulation, including a series of 3D and quasi-2D argon-water-solid systems simulation. To precisely determine the wettability on various substrates, several indicators were calculated. Among these wettability indicators, the water PMF (potential of mean force) has the most correlation tendency with interfacial water molecular orientation than depletion layer width and droplet contact angle. The results reveal that the aggregation of argon molecules on substrates not only depending on the level of hydrophobicity but also determined by the competition between gas-solid and water-solid interaction as well as water molecular structure near the surface. In addition, the surface nanobubble is always observed coexisted with the gas enrichment layer. The water structure adjacent to water-gas and water-solid interfaces also plays an important factor in gas out-flux and gas aggregation, respectively. The quasi-2D simulation shows that only a slight difference in the curved argon-water interface from the plane interface which suggests no noticeable obstructing effect on gas outflux from the gas-water interfacial water networks.

Keywords: gas aggregation, interfacial nanobubble, molecular dynamics simulation, wettability

Procedia PDF Downloads 115

Authors: Hao-Wen Chang, Santhanamoorthi Nachimuthu, Jyh-Chiang Jiang

Abstract:

All-solid-state lithium batteries (ASSLBs) are increasingly recognized as a safer and more reliable alternative to conventional lithium-ion batteries due to their non-flammable nature and enhanced safety performance. ASSLBs utilize a range of solid-state electrolytes, including solid polymer electrolytes (SPEs), inorganic solid electrolytes (ISEs), and composite polymer electrolytes (CPEs). SPEs are particularly valued for their flexibility, ease of processing, and excellent interfacial compatibility with electrodes, though their ionic conductivity remains a significant limitation. ISEs, on the other hand, provide high ionic conductivity, broad electrochemical windows, and strong mechanical properties but often face poor interfacial contact with electrodes, impeding performance. CPEs, which merge the strengths of SPEs and ISEs, represent a compelling solution for next-generation ASSLBs by addressing both electrochemical and mechanical challenges. Despite their potential, the mechanisms governing lithium-ion transport within these systems remain insufficiently understood. In this study, we designed CPEs based on argyrodite-type Li₆PS₅Cl (LPSC) combined with two distinct polymer matrices: poly(ethylene oxide) (PEO) with 24.5 wt% lithium bis(trifluoromethane)sulfonimide (LiTFSI) and polycaprolactone (PCL) with 25.7 wt% LiTFSI. Through density functional theory (DFT) calculations, we investigated the interfacial chemistry of these materials, revealing critical insights into their stability and interactions. Additionally, ab initio molecular dynamics (AIMD) simulations of lithium electrodes interfaced with LPSC layers containing polymers and LiTFSI demonstrated that the polymer matrix significantly mitigates LPSC decomposition, compared to systems with only a lithium electrode and LPSC layers. These findings underscore the pivotal role of CPEs in improving the performance and longevity of ASSLBs, offering a promising path forward for next-generation energy storage technologies.

Keywords: all-solid-state lithium-ion batteries, composite solid electrolytes, DFT calculations, Li-ion transport

Procedia PDF Downloads 20

Authors: Leila Momenzadeh, Alexander V. Evteev, Elena V. Levchenko, Tanvir Ahmed, Irina Belova, Graeme Murch

Abstract:

In this work, the phonon thermal conductivity of f.c.c. Al is investigated in detail in the temperature range 100 – 900 K within the framework of equilibrium molecular dynamics simulations making use of the Green-Kubo formalism and one of the most reliable embedded-atom method potentials. It is found that the heat current auto-correlation function of the f.c.c. Al model demonstrates a two-stage temporal decay similar to the previously observed for f.c.c Cu model. After the first stage of decay, the heat current auto-correlation function of the f.c.c. Al model demonstrates a peak in the temperature range 100-800 K. The intensity of the peak decreases as the temperature increases. At 900 K, it transforms to a shoulder. To describe the observed two-stage decay of the heat current auto-correlation function of the f.c.c. Al model, we employ decomposition model recently developed for phonon-mediated thermal transport in a monoatomic lattice. We found that the electronic contribution to the total thermal conductivity of f.c.c. Al dominates over the whole studied temperature range. However, the phonon contribution to the total thermal conductivity of f.c.c. Al increases as temperature decreases. It is about 1.05% at 900 K and about 12.5% at 100 K.

Keywords: aluminum, gGreen-Kubo formalism, molecular dynamics, phonon thermal conductivity

Procedia PDF Downloads 413

Authors: Seyedmajid Mehrnia, Maximilan Kuhr, Peter F. Pelz

Abstract:

Molecular Dynamics (MD) simulation is a proven method to inspect behaviours of lubricant oils in nano-scale gaps. However, most MD simulations on tribology have been performed with atomically smooth walls to determine wall slip and friction properties. This study will investigate the effect of porosity, specifically nano-porous walls, on wall slip properties of hydrocarbon oils confined between two walls in a Couette flow. Different pore geometries will be modelled to investigate the effect on wall slip and bulk shear. In this paper, the Polyalphaolefin (PAO) molecules are confined to a stationary and a moving wall. A hybrid force field consisting of different potential energy functions was employed in this MD simulation. Newton’s law defines how those forces will influence the atoms' movements. The interactions among surface atoms were simulated with an Embedded Atom Method (EAM) potential function which can represent the characteristics of metallic arrangements very strongly. We implemented NERD forcefield for intramolecular potential energy function. Also, Lennard-Jones potential was employed for nonbonded intermolecular interaction.

Keywords: slip length, molecular dynamics, critical shear rate, Couette flow

Procedia PDF Downloads 131

Authors: Moulay Driss Mellaoui, Khadija Zaki, Khalid Abbiche, Abdallah Imjjad, Rachid Boutiddar, Abdelouahid Sbai, Aaziz Jmiai, Souad El Issami, Al Mokhtar Lamsabhi, Hanane Zejli

Abstract:

This study presents a comprehensive analysis of six isoxazolidine and isoxazoline derivatives, leveraging a multifaceted approach that combines Density Functional Theory (DFT), AdmetSAR analysis, and molecular docking simulations to explore their electronic, pharmacokinetic, and anticancer properties. Through DFT analysis, using the B3LYP-D3BJ functional and the 6-311++G(d,p) basis set, we optimized molecular geometries, analyzed vibrational frequencies, and mapped Molecular Electrostatic Potentials (MEP), identifying key sites for electrophilic attacks and hydrogen bonding. Frontier Molecular Orbital (FMO) analysis and Density of States (DOS) plots revealed varying stability levels among the compounds, with 1b, 2b, and 3b showing slightly higher stability. Chemical potential assessments indicated differences in binding affinities, suggesting stronger potential interactions for compounds 1b and 2b. AdmetSAR analysis predicted favorable human intestinal absorption (HIA) rates for all compounds, highlighting compound 3b superior oral effectiveness. Molecular docking and molecular dynamics simulations were conducted on isoxazolidine and 4-isoxazoline derivatives targeting the EGFR receptor (PDB: 1JU6). Molecular docking simulations confirmed the high affinity of these compounds towards the target protein 1JU6, particularly compound 3b, among the isoxazolidine derivatives, compound 3b exhibited the most favorable binding energy, with a g score of -8.50 kcal/mol. Molecular dynamics simulations over 100 nanoseconds demonstrated the stability and potential of compound 3b as a superior candidate for anticancer applications, further supported by structural analyses including RMSD, RMSF, Rg, and SASA values. This study underscores the promising role of compound 3b in anticancer treatments, providing a solid foundation for future drug development and optimization efforts.

Keywords: isoxazolines, DFT, molecular docking, molecular dynamic, ADMET, drugs.

Procedia PDF Downloads 47

Authors: Mohammad Reza Ghaani, Niall English

Abstract:

Hydrogen is looked upon as the next-generation clean-energy carrier; the search for an efficient material and method for storing hydrogen has been, and is, pursued relentlessly. Clathrate hydrates are inclusion compounds wherein guest gas molecules like hydrogen are trapped in a host water-lattice framework. These types of materials can be categorised as potentially attractive hosting environments for physical hydrogen storage (i.e., no chemical reaction upon storage). Non-equilibrium molecular dynamics (NEMD) simulations have been performed to investigate thermal-driven break-up of propane-hydrate interfaces with liquid water at 270-300 K, with the propane hydrate containing either one or no hydrogen molecule in each of its small cavities. In addition, two types of hydrate-surface water-lattice molecular termination were adopted, at the hydrate edge with water: a 001-direct surface cleavage and one with completed cages. The geometric hydrate-ice-liquid distinction criteria of Báez and Clancy were employed to distinguish between the hydrate, ice lattices, and liquid-phase. Consequently, the melting temperatures of interface were estimated, and dissociation rates were observed to be strongly dependent on temperature, with higher dissociation rates at larger over-temperatures vis-à-vis melting. The different hydrate-edge terminations for the hydrate-water interface led to statistically-significant differences in the observed melting point and dissociation profile: it was found that the clathrate with the planar interface melts at around 280 K, whilst the melting temperature of the cage-completed interface was determined to be circa 270 K.

Keywords: hydrogen storage, clathrate hydrate, molecular dynamics, thermal dissociation

Procedia PDF Downloads 275

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 152

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 196

Authors: Danish Idrees, Vijay Kumar, Samudrala Gourinath

Abstract:

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease caused by misfolding and aggregation of Cu, Zn superoxide dismutase (SOD1). The use of small molecules has been shown to stabilize the SOD1 dimer and preventing its dissociation and aggregation. In this study, we employed molecular docking, molecular dynamics simulation and surface plasmon resonance (SPR) to study the interactions between SOD1 and natural polyphenolic compounds. In order to explore the noncovalent interaction between SOD1 and natural polyphenolic compounds, molecular docking and molecular dynamic (MD) simulations were employed to gain insights into the binding modes and free energies of SOD1-polyphenolic compounds. MM/PBSA methods were used to calculate free energies from obtained MD trajectories. The compounds, Hesperidin, Ergosterol, and Rutin showed the excellent binding affinity in micromolar range with SOD1. Ergosterol and Hesperidin have the strongest binding affinity to SOD1 and was subjected to further characterization. Biophysical experiments using Circular Dichroism and Thioflavin T fluorescence spectroscopy results show that the binding of these two compounds can stabilize SOD1 dimer and inhibit the aggregation of SOD1. Molecular simulation results also suggest that these compounds reduce the dissociation of SOD1 dimers through direct interaction with the dimer interface. This study will be helpful to develop other drug-like molecules which may have the effect to reduce the aggregation of SOD1.

Keywords: amyotrophic lateral sclerosis, molecular dynamics simulation, surface plasmon resonance, superoxide dismutase

Procedia PDF Downloads 139