Search results for: molecular simulations

Authors: Jannes Quer, Amir Niknejad, Marcus Weber

Abstract:

Computational Drug Design is often based on Molecular Dynamics simulations of molecular systems. Molecular Dynamics can be used to simulate, e.g., the binding and unbinding event of a small drug-like molecule with regard to the active site of an enzyme or a receptor. However, the time-scale of the overall binding event is many orders of magnitude longer than the time-scale of simulation. Thus, there is a need to speed-up molecular simulations. In order to speed up simulations, the molecular dynamics trajectories have to be ”steared” out of local minimizers of the potential energy surface – the so-called metastabilities – of the molecular system. Increasing the kinetic energy (temperature) is one possibility to accelerate simulated processes. However, with temperature the entropy of the molecular system increases, too. But this kind ”stearing” is not directed enough to stear the molecule out of the minimum toward the saddle point. In this article, we give a new mathematical idea, how a potential energy surface can be changed in such a way, that entropy is kept under control while the trajectories are still steared out of the metastabilities. In order to compute the unsteared transition behaviour based on a steared simulation, we propose to use extrapolation methods. In the end we mathematically show, that our method accelerates the simulations along the direction, in which the curvature of the potential energy surface changes the most, i.e., from local minimizers towards saddle points.

Keywords: extrapolation, Eyring-Kramers, metastability, multilevel sampling

Procedia PDF Downloads 289

Authors: Hai Hoang, Thanh Xuan Nguyen Thi, Guillaume Galliero

Abstract:

Coarse-Grained (CG) molecular simulations have shown to be an efficient way to estimate thermophysical (static and dynamic) properties of fluids. Several strategies have been developed and reported in the literature for defining CG molecular models. Among them, those based on a top-down strategy (i.e. CG molecular models related to macroscopic observables), despite being heuristic, have increasingly gained attention. This is probably due to its simplicity in implementation and its ability to provide reasonable results for not only simple but also complex systems. Regarding simple Force-Fields associated with these CG molecular models, it has been found that the four parameters Mie chain model is one of the best compromises to describe thermophysical static properties (e.g. phase diagram, saturation pressure). However, parameterization procedures of these Mie-chain GC molecular models given in literature are generally insufficient to simultaneously provide static and dynamic (e.g. viscosity) properties. To deal with such situations, we have extended the corresponding states by using a quantity associated with the liquid viscosity. Results obtained from molecular simulations have shown that our approach is able to yield good estimates for both static and dynamic thermophysical properties for various real non-associating fluids. In addition, we will show that on simple (e.g. phase diagram, saturation pressure) and complex (e.g. thermodynamic response functions, thermodynamic energy potentials) static properties, results of our scheme generally provides improved results compared to existing approaches.

Keywords: coarse-grained model, mie potential, molecular simulations, thermophysical properties, phase equilibria

Procedia PDF Downloads 303

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 308

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 406

Authors: Jiahui Song

Abstract:

The use of high-intensity, nanosecond electric pulses has been a recent development in biomedical. High-intensity (∼100 kV/cm), nanosecond duration-pulsed electric fields have been shown to induce cellular electroporation. This will lead to an increase in transmembrane conductivity and diffusive permeability. These effects will also alter the electrical potential across the membrane. The applications include electrically triggered intracellular calcium release, shrinkage of tumors, and temporary blockage of the action potential in nerves. In this research, the dynamics of pore formation with the presence of an externally applied electric field is studied on the basis of molecular dynamics (MD) simulations using the GROMACS package. MD simulations show pore formation occurs for a pulse with the amplitude of 0.5V/nm at 1ns at temperature 316°K. Also increasing temperatures facilitate pore formation. When the temperature is increased to 323°K, pore forms at 0.75ns with the pulse amplitude of 0.5V/nm. For statistical significance, a total of eight MD simulations are carried out with different starting molecular velocities for each simulation. Also, actual experimental observations are compared against MD simulation results.

Keywords: molecular dynamics, high-intensity, nanosecond, electroporation

Procedia PDF Downloads 77

Authors: Seyyed Feisal Asbaghian Namin, Reza Pilafkan, Mahmood Kaffash Irzarahimi

Abstract:

TThis paper considers vibration of single-layered graphene sheets using molecular dynamics (MD) and nonlocal elasticity theory. Based on the MD simulations, Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), an open source software, is used to obtain fundamental frequencies. On the other hand, governing equations are derived using nonlocal elasticity and first order shear deformation theory (FSDT) and solved using generalized differential quadrature method (GDQ). The small-scale effect is applied in governing equations of motion by nonlocal parameter. The effect of different side lengths, boundary conditions and nonlocal parameter are inspected for aforementioned methods. Results are obtained from MD simulations is compared with those of the nonlocal elasticity theory to calculate appropriate values for the nonlocal parameter. The nonlocal parameter value is suggested for graphene sheets with various boundary conditions. Furthermore, it is shown that the nonlocal elasticity approach using classical plate theory (CLPT) assumptions overestimates the natural frequencies.

Keywords: graphene sheets, molecular dynamics simulations, fundamental frequencies, nonlocal elasticity theory, nonlocal parameter

Procedia PDF Downloads 472

Authors: Akarsh Verma, Avinash Parashar

Abstract:

The aim of this article is to study the effects of hydroxyl functional group on the mechanical strength and fracture toughness of graphene. This functional group forms the backbone of intrinsic atomic structure of graphene oxide (GO). Molecular dynamics-based simulations were performed in conjunction with reactive force field (ReaxFF) parameters to capture the mode-I fracture toughness of hydroxyl functionalised graphene. Moreover, these simulations helped in concluding that spatial distribution and concentration of hydroxyl functional group significantly affects the fracture morphology of graphene nanosheet. In contrast to literature investigations, atomistic simulations predicted a transition in the failure morphology of hydroxyl functionalised graphene from brittle to ductile as a function of its spatial distribution on graphene sheet.

Keywords: graphene, graphene oxide, ReaxFF, molecular dynamics

Procedia PDF Downloads 142

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 55

Authors: Adam Abate, Elham Rasti, Philip Romero

Abstract:

Beta-glucosidase is the key enzyme component present in cellulase and completes the final step during cellulose hydrolysis by converting the cellobiose to glucose. The regulatory properties of beta-glucosidases are most commonly found for the retaining and inverting enzymes. Hydrolysis of a glycoside typically occurs with general acid and general base assistance from two amino acid side chains, normally glutamic or aspartic acids. In order to obtain more detailed information on the dynamic events origination from the interaction with enzyme active site, we carried out molecular dynamics simulations of beta-glycosidase in protonated state (Glu-H178) and deprotonated state (Glu178). The theoretical models generated from our molecular dynamics simulations complement and advance the structural information currently available, leading to a more detailed understanding of Beta-glycosidase structure and function. This article presents the important role of Asn307 in enzyme activity of beta-glucosidase

Keywords: Beta-glucosidase, GROMACS, molecular dynamics simulation, structural parameters

Procedia PDF Downloads 364

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

Abstract:

We present a theoretical study of the structural, elastic and thermodynamic properties of the zinc-blende AlBi for a wide temperature range. The simulation calculation is performed in the framework of the molecular dynamics method using the three-body Tersoff potential which reproduces provide, with reasonable accuracy, the lattice constants and elastic constants. Our results for the lattice constant, the bulk modulus and cohesive energy are in good agreement with other theoretical available works. Other thermodynamic properties such as the specific heat and the lattice thermal expansion can also be predicted. 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: aluminium compounds, molecular dynamics simulations, interatomic potential, thermodynamic properties, structural phase transition

Procedia PDF Downloads 266

Authors: Betim Bahtiri, B. Arash, R. Rolfes

Abstract:

Epoxy-based materials are frequently exposed to high-humidity environments in many engineering applications. As a result, their material properties would be degraded by water absorption. A full characterization of the material properties under hygrothermal conditions requires time- and cost-consuming experimental tests. To gain insights into the physics of diffusion mechanisms, atomistic simulations have been shown to be effective tools. Concerning the diffusion of water in polymers, spatial trajectories of water molecules are obtained from molecular dynamics (MD) simulations allowing the interpretation of diffusion pathways at the nanoscale in a polymer network. Conventional MD simulations of water diffusion in amorphous polymers lead to discrepancies at low temperatures due to the short timescales of the simulations. In the proposed model, this issue is solved by using a combined scheme of autonomous basin climbing (ABC) with kinetic Monte Carlo and reactive MD simulations to investigate the diffusivity of water molecules in epoxy resins across a wide range of temperatures. It is shown that the proposed simulation framework estimates kinetic properties of water diffusion in epoxy resins that are consistent with experimental observations and provide a predictive tool for investigating the diffusion of small molecules in other amorphous polymers.

Keywords: epoxy resins, water diffusion, autonomous basin climbing, kinetic Monte Carlo, reactive molecular dynamics

Procedia PDF Downloads 32

Authors: Aibek Kukpayev, Dhawal Shah

Abstract:

Combustion of sour fuels containing high amount of sulfur leads to the formation of sulfur oxides, which adversely harm the environment and has a negative impact on human health. Considering this, several legislations have been imposed to bring down the sulfur content in fuel to less than 10 ppm. In recent years, novel deep eutectic solvents (DESs) have been developed to achieve deep desulfurization, particularly to extract thiophenic compounds from liquid fuels. These novel DESs, considered as analogous to ionic liquids are green, eco-friendly, inexpensive, and sustainable. We herein, using molecular dynamic simulation, analyze the interactions of metal-based DESs with model oil consisting of thiophenic compounds. The DES used consists of polyethylene glycol (PEG-200) as a hydrogen bond donor, choline chloride (ChCl) or tetrabutyl ammonium chloride (TBAC) as a hydrogen bond acceptor, and cobalt chloride (CoCl₂) as metal salt. In particular, the combination of ChCl: PEG-200:CoCl₂ at a ratio 1:2:1 and the combination of TBAC:PEG-200:CoCl₂ at a ratio 1:2:0.25 were simulated, separately, with model oil consisting of octane and thiophenes at 25ᵒC and 1 bar. The results of molecular dynamics simulations were analyzed in terms of interaction energies between different components. The simulations revealed a stronger interaction between DESs/thiophenes as compared with octane/thiophenes, suggestive of an efficient desulfurization process. In addition, our analysis suggests that the choice of hydrogen bond acceptor strongly influences the efficiency of the desulfurization process. Taken together, the results also show the importance of the metal ion, although present in small amount, in the process, and the role of the polymer in desulfurization of the model fuel.

Keywords: deep eutectic solvents, desulfurization, molecular dynamics simulations, thiophenes

Procedia PDF Downloads 110

Authors: Mohamed Moussaoui, Mouna Baassi, Soukayna Baammi, Hatim Soufi, Mohammed Salah, Rachid Daoud, Achraf EL Allali, Mohammed Elalaoui Belghiti, Said Belaaouad

Abstract:

The present study aims to investigate the quantitative structure-activity relationship (QSAR) of a series of Thiazole derivatives reported as anticancer agents (hepatocellular carcinoma), using principally the electronic descriptors calculated by the density functional theory (DFT) method and by applying the multiple linear regression method. The developed model showed good statistical parameters (R²= 0.725, R²ₐ𝒹ⱼ= 0.653, MSE = 0.060, R²ₜₑₛₜ= 0.827, Q²𝒸ᵥ = 0.536). The energy of the highest occupied molecular orbital (EHOMO) orbital, electronic energy (TE), shape coefficient (I), number of rotatable bonds (NROT), and index of refraction (n) were revealed to be the main descriptors influencing the anti-cancer activity. Additional Thiazole derivatives were then designed and their activities and pharmacokinetic properties were predicted using the validated QSAR model. These designed molecules underwent evaluation through molecular docking (MD) and molecular dynamic (MD) simulations, with binding affinity calculated using the MMPBSA script according to a 100 ns simulation trajectory. This process aimed to study both their affinity and stability towards Cyclin-Dependent Kinase 2 (CDK2), a target protein for cancer disease treatment. The research concluded by identifying four CDK2 inhibitors - A1, A3, A5, and A6 - displaying satisfactory pharmacokinetic properties. MDs results indicated that the designed compound A5 remained stable in the active center of the CDK2 protein, suggesting its potential as an effective inhibitor for the treatment of hepatocellular carcinoma. The findings of this study could contribute significantly to the development of effective CDK2 inhibitors.

Keywords: QSAR, ADMET, Thiazole, anticancer, molecular docking, molecular dynamic simulations, MMPBSA calculation

Procedia PDF Downloads 57

Authors: Titus A. Beu

Abstract:

Many modern gene-delivery protocols rely on condensed complexes of DNA with polycations to introduce the genetic payload into cells by endocytosis. In particular, polyethyleneimine (PEI) stands out by a high buffering capacity (enabling the efficient condensation of DNA) and relatively simple fabrication. Realistic computational studies can offer essential insights into the formation process of DNA-PEI polyplexes, providing hints on efficient designs and engineering routes. We present comprehensive computational investigations of solvated PEI and DNA-PEI polyplexes involving calculations at three levels: ab initio, all-atom (AA), and coarse-grained (CG) molecular mechanics. In the first stage, we developed a rigorous AA CHARMM (Chemistry at Harvard Macromolecular Mechanics) force field (FF) for PEI on the basis of accurate ab initio calculations on protonated model pentamers. We validated this atomistic FF by matching the results of extensive molecular dynamics (MD) simulations of structural and dynamical properties of PEI with experimental data. In a second stage, we developed a CG MARTINI FF for PEI by Boltzmann inversion techniques from bead-based probability distributions obtained from AA simulations and ensuring an optimal match between the AA and CG structural and dynamical properties. In a third stage, we combined the developed CG FF for PEI with the standard MARTINI FF for DNA and performed comprehensive CG simulations of DNA-PEI complex formation and condensation. Various technical aspects which are crucial for the realistic modeling of DNA-PEI polyplexes, such as options of treating electrostatics and the relevance of polarizable water models, are discussed in detail. Massive CG simulations (with up to 500 000 beads) shed light on the mechanism and provide time scales for DNA polyplex formation independence of PEI chain size and protonation pattern. The DNA-PEI condensation mechanism is shown to primarily rely on the formation of DNA bundles, rather than by changes of the DNA-strand curvature. The gained insights are expected to be of significant help for designing effective gene-delivery applications.

Keywords: DNA condensation, gene-delivery, polyethylene-imine, molecular dynamics.

Procedia PDF Downloads 87

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 299

Authors: Neslihan Demirci, Serdar Durdağı

Abstract:

Mycobacterium tuberculosis is still a worldwide disease-causing agent that, according to WHO, led to the death of 1.5 million people from tuberculosis (TB) in 2020. The bacteria reside in macrophages located specifically in the lung. There is a known quadruple drug therapy regimen for TB consisting of isoniazid (INH), rifampin (RIF), pyrazinamide (PZA), and ethambutol (EMB). Over the past 60 years, there have been great contributions to treatment options, such as recently approved delamanid (OPC67683) and bedaquiline (TMC207/R207910), targeting mycolic acid and ATP synthesis, respectively. Also, there are natural compounds that can block the tryptophanyl-tRNA synthetase (TrpRS) enzyme, chuangxinmycin, and indolmycin. Yet, already the drug resistance is reported for those agents. In this study, the newly released TrpRS enzyme structure is investigated for potential inhibitor drugs from already synthesized molecules to help the treatment of resistant cases and to propose an alternative drug for the quadruple drug therapy of tuberculosis. Maestro, Schrodinger is used for docking and molecular dynamic simulations. In-house library containing ~8000 compounds among FDA-approved indole-containing compounds, a total of 57 obtained from the ChemBL were used for both ATP and tryptophan binding pocket docking. Best of indole-containing 57 compounds were subjected to hit expansion and compared later with virtual screening workflow (VSW) results. After docking, VSW was done. Glide-XP docking algorithm was chosen. When compared, VSW alone performed better than the hit expansion module. Best scored compounds were kept for ten ns molecular dynamic simulations by Desmond. Further, 100 ns molecular dynamic simulation was performed for elected molecules according to Z-score. The top three MMGBSA-scored compounds were subjected to steered molecular dynamic (SMD) simulations by Gromacs. While SMD simulations are still being conducted, ponesimod (for multiple sclerosis), vilanterol (β₂ adrenoreceptor agonist), and silodosin (for benign prostatic hyperplasia) were found to have a significant affinity for tuberculosis TrpRS, which is the propulsive force for the urge to expand the research with in vitro studies. Interestingly, top-scored ponesimod has been reported to have a side effect that makes the patient prone to upper respiratory tract infections.

Keywords: drug repurposing, molecular dynamics, tryptophanyl-tRNA synthetase, tuberculosis

Procedia PDF Downloads 82

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 73

Authors: Catherine Vasnetsov, Victor Vasnetsov

Abstract:

Context: Thermo-responsive polymers are materials that undergo significant changes in their physical properties in response to temperature changes. These polymers have gained significant attention in research due to their potential applications in various industries and medicine. However, the molecular mechanisms underlying their behavior are not well understood, particularly in relation to cosolvency, which is crucial for practical applications. Research Aim: This study aimed to theoretically investigate the phenomenon of cosolvency in long-chain polymers using the Flory-Huggins statistical-mechanical framework. The main objective was to understand the interactions between the polymer, solvent, and cosolvent under different conditions. Methodology: The research employed a combination of Monte Carlo computer simulations and advanced machine-learning methods. The Flory-Huggins mean field theory was used as the basis for the simulations. Spinodal graphs and ternary plots were utilized to develop an initial computer model for predicting polymer behavior. Molecular dynamic simulations were conducted to mimic real-life polymer systems. Machine learning techniques were incorporated to enhance the accuracy and reliability of the simulations. Findings: The simulations revealed that the addition of very low or very high volumes of cosolvent molecules resulted in smaller radii of gyration for the polymer, indicating poor miscibility. However, intermediate volume fractions of cosolvent led to higher radii of gyration, suggesting improved miscibility. These findings provide a possible microscopic explanation for the cosolvency phenomenon in polymer systems. Theoretical Importance: This research contributes to a better understanding of the behavior of thermo-responsive polymers and the role of cosolvency. The findings provide insights into the molecular mechanisms underlying cosolvency and offer specific predictions for future experimental investigations. The study also presents a more rigorous analysis of the Flory-Huggins free energy theory in the context of polymer systems. Data Collection and Analysis Procedures: The data for this study was collected through Monte Carlo computer simulations and molecular dynamic simulations. The interactions between the polymer, solvent, and cosolvent were analyzed using the Flory-Huggins mean field theory. Machine learning techniques were employed to enhance the accuracy of the simulations. The collected data was then analyzed to determine the impact of cosolvent volume fractions on the radii of gyration of the polymer. Question Addressed: The research addressed the question of how cosolvency affects the behavior of long-chain polymers. Specifically, the study aimed to investigate the interactions between the polymer, solvent, and cosolvent under different volume fractions and understand the resulting changes in the radii of gyration. Conclusion: In conclusion, this study utilized theoretical modeling and computer simulations to investigate the phenomenon of cosolvency in long-chain polymers. The findings suggest that moderate cosolvent volume fractions can lead to improved miscibility, as indicated by higher radii of gyration. These insights contribute to a better understanding of the molecular mechanisms underlying cosolvency in polymer systems and provide predictions for future experimental studies. The research also enhances the theoretical analysis of the Flory-Huggins free energy theory.

Keywords: molecular modelling, flory-huggins, cosolvency, stimuli-responsive polymers

Procedia PDF Downloads 38

Authors: Eunae Lee, Oh-Kuen Kwon, Ki-Sub Kim, Jeong Won Kang

Abstract:

We investigated the dynamic properties of graphene-nanoribbon (GNR) memory encapsulating graphene-nanoflake (GNF) shuttle in the potential to be applicable as a non-volatile random access memory via molecular dynamics simulations. This work explicitly demonstrates that the GNR encapsulating the GNF shuttle can be applied to nonvolatile memory. The potential well was originated by the increase of the attractive vdW energy between the GNRs when the GNF approached the edges of the GNRs. So the bistable positions were located near the edges of the GNRs. Such a nanoelectromechanical non-volatile memory based on graphene is also applicable to the development of switches, sensors, and quantum computing.

Keywords: graphene nanoribbon, graphene nanoflake, shuttle memory, molecular dynamics

Procedia PDF Downloads 408

Authors: Aliyu Adamu, Fahrul Huyop, Roswanira Abdul Wahab, Mohd Shahir Shamsir

Abstract:

Halogenated organic compounds occur in huge amount in biosphere. This is attributable to the diverse use of halogen-based compounds in the synthesis of various industrially important products. Halogenated compound is toxic and may persist in the environment, thereby causing serious health and environmental pollution problems. L-2-haloacid dehalogenases (EC 3.8.1.2) catalyse the specific cleavage of carbon-halogen bond in L-isomers of halogenated compounds, which consequently reverse the effects of environmental halogen-associated pollution. To enhance the efficiency and utility of these enzymes, this study investigates the catalytic amino acid residues and the molecular functional mechanism of DehL, by classical molecular dynamic simulations, MM-PBSA and ab initio fragments molecular orbital (FMO) calculations. The results of the study will serve as the basis for the molecular engineering of the enzyme.

Keywords: DehL, Functional mechanism, Catalytic residues, L-2-haloacid dehalogenase

Procedia PDF Downloads 325

Authors: Jayanthi Sivaraman

Abstract:

Claudins are the important components of the tight junctions that play a key role in paracellular permeability. Among various members of Claudin family, Claudin 4 (CLDN4) is found to be overexpressed in ovarian, pancreatic carcinomas and other epithelial malignancies. Therefore, in this study, an attempt has been made to identify potent inhibitors for CLDN4 from the ZINC database using virtual screening, molecular docking and molecular dynamics simulations. A well refined molecular model of CLDN4 was built using Prime of Schrodinger v10.2(Template- PDB ID: 4P79). Approximately, 6 million compounds from ZINC database are subjected to high-throughput virtual screening (HTVS) against the active site of CLDN4. Molecular docking using GLIDE predicted ARG31, ASN142, ASP146 and ARG158 as critically important residues. Furthermore, three compounds from ZINC database (ZINC96331839, ZINC36533519 and ZINC75819394) showed highly promising ADME properties and binding affinity with stable conformation. The therapeutic efficiency of these lead compounds is evaluated and confirmed by in-vitro and in-vivo studies which leads to the development of novel anti-cancer drugs.

Keywords: ADME property, inhibitors, molecular docking, virtual screening

Procedia PDF Downloads 301

Authors: Sakina Fatima, Joseph-Patrick W. E. Clarke, Patricia A. Thibault, Subha Kalyaanamoorthy, Michael Levin, Aravindhan Ganesan

Abstract:

Heteronuclear ribonucleoprotein (hnRNPA1 or A1) is associated with the pathology of different diseases, including neurological disorders and cancers. In particular, the aggregation and dysfunction of A1 have been identified as a critical driver for neurodegeneration (NDG) in Multiple Sclerosis (MS). Structurally, A1 includes a low-complexity domain (LCD) and two RNA-recognition motifs (RRMs), and their interdomain coordination may play a crucial role in A1 aggregation. Previous studies propose that RNA-inhibitors or nucleoside analogs that bind to RRMs can potentially prevent A1 self-association. Therefore, molecular-level understanding of the structures, dynamics, and nucleotide interactions with A1 RRMs can be useful for developing therapeutics for NDG in MS. In this work, a combination of computational modelling and biochemical experiments were employed to analyze a set of RNA-A1 RRM complexes. Initially, the atomistic models of RNA-RRM complexes were constructed by modifying known crystal structures (e.g., PDBs: 4YOE and 5MPG), and through molecular docking calculations. The complexes were optimized using molecular dynamics simulations (200-400 ns), and their binding free energies were computed. The binding affinities of the selected complexes were validated using a thermal shift assay. Further, the most important molecular interactions that contributed to the overall stability of the RNA-A1 RRM complexes were deduced. The results highlight that adenine and guanine are the most suitable nucleotides for high-affinity binding with A1. These insights will be useful in the rational design of nucleotide-analogs for targeting A1 RRMs.

Keywords: hnRNPA1, molecular docking, molecular dynamics, RNA-binding proteins

Procedia PDF Downloads 78

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 324

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 114

Authors: Mouna Baassi, Mohamed Moussaoui, Sanchaita Rajkhowa, Hatim Soufi, Said Belaaouad

Abstract:

The objective of this paper is to address the challenging task of targeting Human Immunodeficiency Virus type 1 Reverse Transcriptase (HIV-1 RT) in the treatment of AIDS. Reverse Transcriptase inhibitors (RTIs) have limitations due to the development of Reverse Transcriptase mutations that lead to treatment resistance. In this study, a combination of statistical analysis and bioinformatics tools was adopted to develop a mathematical model that relates the structure of compounds to their inhibitory activities against HIV-1 Reverse Transcriptase. Our approach was based on a series of compounds recognized for their HIV-1 RT enzymatic inhibitory activities. These compounds were designed via software, with their descriptors computed using multiple tools. The most statistically promising model was chosen, and its domain of application was ascertained. Furthermore, compounds exhibiting comparable biological activity to existing drugs were identified as potential inhibitors of HIV-1 RT. The compounds underwent evaluation based on their chemical absorption, distribution, metabolism, excretion, toxicity properties, and adherence to Lipinski's rule. Molecular docking techniques were employed to examine the interaction between the Reverse Transcriptase (Wild Type and Mutant Type) and the ligands, including a known drug available in the market. Molecular dynamics simulations were also conducted to assess the stability of the RT-ligand complexes. Our results reveal some of the new compounds as promising candidates for effectively inhibiting HIV-1 Reverse Transcriptase, matching the potency of the established drug. This necessitates further experimental validation. This study, beyond its immediate results, provides a methodological foundation for future endeavors aiming to discover and design new inhibitors targeting HIV-1 Reverse Transcriptase.

Keywords: QSAR, ADMET properties, molecular docking, molecular dynamics simulation, reverse transcriptase inhibitors, HIV type 1

Procedia PDF Downloads 48

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 159

Authors: Tivani Phosa Mashamba-Thompson, Mahmoud E. S. Soliman

Abstract:

To date, the cause of neurodegeneration is not well understood and diseases that stem from neurodegeneration currently have no known cures. Cathepsin B (CB) enzyme is known to be involved in the production of peptide neurotransmitters and toxic peptides in neurodegenerative diseases (NDs). CA-074Me is a membrane-permeable irreversible selective cathepsin B (CB) inhibitor as confirmed by in vivo studies. Due to the lack of the crystal structure, the binding mode of CA-074Me with the human CB at molecular level has not been previously reported. The main aim of this study is to gain an insight into the binding mode of CB CA-074Me to human CB using various computational tools. Herein, molecular dynamics simulations, binding free energy calculations and per-residue energy decomposition analysis were employed to accomplish the aim of the study. Another objective was to identify novel CB inhibitors based on the structure of CA-074Me using fragment based drug design using scaffold hoping drug design approach. Results showed that two of the designed ligands (hit 1 and hit 2) were found to have better binding affinities than the prototype inhibitor, CA-074Me, by ~2-3 kcal/mol. Per-residue energy decomposition showed that amino acid residues Cys29, Gly196, His197 and Val174 contributed the most towards the binding. The Van der Waals binding forces were found to be the major component of the binding interactions. The findings of this study should assist medicinal chemist towards the design of potential irreversible CB inhibitors.

Keywords: cathepsin B, scaffold hopping, docking, molecular dynamics, binding-free energy, neurodegerative diseases

Procedia PDF Downloads 343

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 295

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 265

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 71