Search results for: reactive molecular dynamics

Authors: Ali Rahnamoun, Adri van Duin

Abstract:

The concept of high velocity particle impact on the spacecraft surface materials has been one of the important issues in the design of such materials. Among these particles, water clusters might be the most abundant and the most important particles to be studied. The importance of water clusters is that upon impact on the surface of the materials, they can cause damage to the material and also if they are sub-cooled water clusters, they can attach to the surface of the materials and cause ice accumulation on the surface which is very problematic in spacecraft and also aircraft operations. The dynamics of the collisions between amorphous silica structures and water clusters with impact velocities of 1 km/s to 10 km/s are studied using the ReaxFF reactive molecular dynamics simulation method. The initial water clusters include 150 water molecules and the water clusters are collided on the surface of amorphous fully oxidized and suboxide silica structures. These simulations show that the most abundant molecules observed on the silica surfaces, other than reflecting water molecules, are H3O+ and OH- for the water cluster impacts on suboxide and fully oxidized silica structures, respectively. The effect of impact velocity on the change of silica mass is studied. At high impact velocities the water molecules attach to the silica surface through a chemisorption process meaning that water molecule dissociates through the interaction with silica surface. However, at low impact velocities, physisorbed water molecules are also observed, which means water molecule attaches and accumulates on the silica surface. The amount of physisorbed waters molecules at low velocities is higher on the suboxide silica surfaces. The evolution of the temperatures of the water clusters during the collisions indicates that the possibility of electron excitement at impact velocities less than 10 km/s is minimal and ReaxFF reactive molecular dynamics simulation can predict the chemistry of these hypervelocity impacts. However, at impact velocities close to 10 km/s the average temperature of the impacting water clusters increase to about 2000K, with individual molecules oocasionally reaching temperatures of over 8000K and thus will be prudent to consider the concept of electron excitation at these higher impact velocities which goes beyond the current ReaxFF ability.

Keywords: spacecraft materials, hypervelocity impact, reactive molecular dynamics simulation, amorphous silica

Procedia PDF Downloads 381

Authors: R. Mohan, V. Jadhav, A. Ahmed, J. Rivas, A. Kelkar

Abstract:

Cementitious materials are an excellent example of a composite material with complex hierarchical features and random features that range from nanometer (nm) to millimeter (mm) scale. Multi-scale modeling of complex material systems requires starting from fundamental building blocks to capture the scale relevant features through associated computational models. In this paper, molecular dynamics (MD) modeling is employed to predict the effect of plasticizer additive on the mechanical properties of key hydrated cement constituent calcium-silicate-hydrate (CSH) at the molecular, nanometer scale level. Due to complexity, still unknown molecular configuration of CSH, a representative configuration widely accepted in the field of mineral Jennite is employed. The effectiveness of the Molecular Dynamics modeling to understand the predictive influence of material chemistry changes based on molecular/nanoscale models is demonstrated.

Keywords: cement composite, mechanical properties, molecular dynamics, plasticizer additives

Procedia PDF Downloads 415

Authors: K. Veluraja

Abstract:

Glycan of glycolipids and glycoproteins is playing a significant role in living systems particularly in molecular recognition processes. Molecular recognition processes are attributed to their occurrence on the surface of the cell, sequential arrangement and type of sugar molecules present in the oligosaccharide structure and glyosidic linkage diversity (glycoinformatics) and conformational diversity (glycoconformatics). Molecular Dynamics Simulation study is a theoretical-cum-computational tool successfully utilized to establish glycoconformatics of glycan. The study on various oligosaccharides of glycan clearly indicates that oligosaccharides do exist in multiple conformational states and these conformational states arise due to the flexibility associated with a glycosidic torsional angle (φ,ψ) . As an example: a single disaccharide structure NeuNacα(2-3) Gal exists in three different conformational states due to the differences in the preferential value of glycosidic torsional angles (φ,ψ). Hence establishing three dimensional structural and conformational models for glycan (cartesian coordinates of every individual atoms of an oligosaccharide structure in a preferred conformation) is quite crucial to understand various molecular recognition processes such as glycan-toxin interaction and glycan-virus interaction. The gycoconformatics models obtained for various glycan through Molecular Dynamics Simulation stored in our 3DSDSCAR (3DSDSCAR.ORG) a public domain database and its utility value in understanding the molecular recognition processes and in drug design venture will be discussed.

Keywords: glycan, glycoconformatics, molecular dynamics simulation, oligosaccharide

Procedia PDF Downloads 97

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

Abstract:

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

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

Procedia PDF Downloads 404

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 144

Authors: T. Rajamäki

Abstract:

In clean room environment molecular contamination in very small concentrations can cause significant harm for the components and processes. This is commonly referred as airborne molecular contamination (AMC). There is a shortage of high sensitivity continuous measurement data for existence and behavior of several of these contaminants. Accordingly, in most cases correlation between concentration of harmful molecules and their effect on processes is not known. In addition, the formation and distribution of contaminating molecules are unclear. In this work sensitive optical techniques are applied in clean room facilities for investigation of concentrations, forming mechanisms and effects of contaminating molecules. Special emphasis is on reactive acid and base gases ammonia (NH3) and hydrogen fluoride (HF). They are the key chemicals in several operations taking place in clean room processes.

Keywords: AMC, clean room, concentration, reactive gas

Procedia PDF Downloads 251

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 36

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 365

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 410

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 269

Authors: A. Rajabpour, A. Hadizadeh Kheirkhah

Abstract:

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

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

Procedia PDF Downloads 273

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 325

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 81

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 291

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

Abstract:

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

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

Procedia PDF Downloads 343

Authors: A. S. Walkey, N. P. Patidar

Abstract:

It is not always economical to provide reactive power using synchronous alternators. The cost of reactive power can be minimized by optimal placing of FACTS devices in power systems. In this paper a Particle Swarm Optimization- Sequential Quadratic Programming (PSO-SQP) algorithm is applied to minimize the cost of reactive power generation along with real power generation to alleviate the bus voltage violations. The effectiveness of proposed approach tested on IEEE-14 bus systems. In this paper in addition to synchronous generators, an opportunity of FACTS devices are also proposed to procure the reactive power demands in the power system.

Keywords: reactive power, reactive power cost, voltage security margins, capability curve, FACTS devices

Procedia PDF Downloads 472

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 80

Authors: Mahdi Golriz

Abstract:

The aim of this work is to estimate the change in molecular structure of linear low-density polyethylene (LLDPE) during peroxide modification can be detected by a simple rheological method. For this purpose a commercial grade LLDPE (Exxon MobileTM LL4004EL) was reacted with different doses of dicumyl peroxide (DCP). The samples were analyzed by size-exclusion chromatography coupled with a light scattering detector. The dynamic shear oscillatory measurements showed a deviation of the δ-׀G ׀٭curve from that of the linear LLDPE, which can be attributed to the presence of long-chain branching (LCB). By the use of a simple rheological method that utilizes melt rheology, transformations in molecular architecture induced on an originally linear low density polyethylene during the early stages of reactive modification were indicated. Reasonable and consistent estimates are obtained, concerning the degree of LCB, the volume fraction of the various molecular species produced in peroxide modification of LLDPE.

Keywords: linear low-density polyethylene, peroxide modification, long-chain branching, rheological method

Procedia PDF Downloads 124

Authors: Kaining Zhang, Lang Chen, Danyang Liu, Jianying Lu, Kun Yang, Junying Wu

Abstract:

In order to solve the problem of a large amount of simulation and limited simulation scale in the first-principle molecular dynamics simulation of energetic material detonation reaction, we established an artificial intelligence model for analyzing and predicting the detonation reaction mechanism of CL-20 based on the first-principle molecular dynamics simulation of the multiscale shock technique (MSST). We employed principal component analysis to identify the dominant charge features governing molecular reactions. We adopted the K-means clustering algorithm to cluster the reaction paths and screen out the key reactions. We introduced the neural network algorithm to construct the mapping relationship between the charge characteristics of the molecular structure and the key reaction characteristics so as to establish a calculation method for predicting detonation reactions based on the charge characteristics of CL-20 and realize the rapid analysis of the reaction mechanism of energetic materials.

Keywords: energetic material detonation reaction, first-principle molecular dynamics simulation of multiscale shock technique, neural network, CL-20

Procedia PDF Downloads 68

Authors: Maedeh Rahimnejad, Bahman Vahidi, Bahman Ebrahimi Hoseinzadeh, Fatemeh Yazdian, Puria Motamed Fath, Roghieh Jamjah

Abstract:

Introduction: The intensive labor and high cost of developing new vehicles for controlled drug delivery highlights the need for a change in their discovery process. Computational models can be used to accelerate experimental steps and control the high cost of experiments. Methods: In this work, to better understand the interaction of anti-cancer drug and the nanocarrier with the cell membrane, we have done molecular dynamics simulation using NAMD. We have chosen paclitaxel for the drug molecule and dipalmitoylphosphatidylcholine (DPPC) as a natural phospholipid nanocarrier. Results: Next, center of mass (COM) between molecules and the van der Waals interaction energy close to the cell membrane has been analyzed. Furthermore, the simulation results of the paclitaxel interaction with the cell membrane and the interaction of DPPC as a nanocarrier loaded by the drug with the cell membrane have been compared. Discussion: Analysis by molecular dynamics (MD) showed that not only the energy between the nanocarrier and the cell membrane is low, but also the center of mass amount decreases in the nanocarrier and the cell membrane system during the interaction; therefore they show significantly better interaction in comparison to the individual drug with the cell membrane.

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

Procedia PDF Downloads 262

Authors: J. S. Savier

Abstract:

Real power is the component power which is converted into useful energy whereas reactive power is the component of power which cannot be converted to useful energy but it is required for the magnetization of various electrical machineries. If the reactive power is compensated at the consumer end, the need for reactive power flow from generators to the load can be avoided and hence the overall power loss can be reduced. In this scenario, this paper presents a succinct method called JSS method for allocation of reactive power losses to consumers connected to radial distribution networks in a deregulated environment. The proposed method has the advantage that no assumptions are made while deriving the reactive power loss allocation method.

Keywords: deregulation, reactive power loss allocation, radial distribution systems, succinct method

Procedia PDF Downloads 340

Authors: Christy Rosaline, Rathankar Roa, Waheeta Hopper

Abstract:

Persistence of tuberculosis, emergence of multidrug-resistance and extensively drug-resistant forms of the disease, has increased the interest in developing new antitubercular drugs. Developing inhibitors for 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase from Mycobacterium tuberculosis (MtbDAH7Ps), an enzyme involved in shikimate pathway, gives a selective target for antitubercular agents. MtbDAH7Ps was screened against ZINC database, and shortlisted compounds were subjected to induce fit docking. Prime/Molecular Mechanics Generalized Born Surface Area calculation was used to validate the binding energy of ligand-protein complex. Molecular Dynamics analysis for of the lead compounds–MtbDAH7Ps complexes showed that the backbone of MtbDAH7Ps in their complexes were stable. These results suggest that the shortlisted lead compounds ZINC04097114, ZINC15163225, ZINC16857013, ZINC06275603, and ZINC05331260 could be developed into novel drug leads to inhibit DAH7Ps in Mycobacterium tuberculosis.

Keywords: MtbDAH7Ps, Mycobacterium tuberculosis, HTVS, molecular dynamics

Procedia PDF Downloads 144

Authors: M. Dehestani, M. Ghasemi-Kooch

Abstract:

In this work, adsorption of chlorophylls a and b pigments in aqueous solution on the inner and outer surfaces of single-walled carbon nanotube (SWCNT) has been studied using molecular dynamics simulation. The linear interaction energy algorithm has been used to calculate the binding free energy. The results show that the adsorption of two pigments is fine on the both positions. Although there is the close similarity between these two pigments, their interaction with the nanotube is different. This result is useful to separate these pigments from one another. According to interaction energy between the pigments and carbon nanotube, interaction between these pigments-SWCNT on the inner surface is stronger than the outer surface. The interaction of SWCNT with chlorophylls phytol tail is stronger than the interaction of SWCNT with porphyrin ring of chlorophylls.

Keywords: adsorption, chlorophyll, interaction, molecular dynamics simulation, nanotube

Procedia PDF Downloads 201

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 271

Authors: M. Dehestani, F. Kamali, M. Klantari Pour, L. Zeidabadi-Nejad

Abstract:

Chemical bonding between polyethylene glycol (PEG) with pharmaceutical proteins called pegylation is one of the most effective methods of improving the pharmacological properties. The covalent attachment of polyethylene glycol (PEG) to proteins will increase their pharmacologic properties. For the formation of a combination of pegylated protein should first be activated PEG and connected to the protein. Interferons(IFNs) are a family of cytokines which show antiviral effects in front of the biological and are responsible for setting safety system. In this study, the nature and properties of the interaction between active positions of IFNs and polyethylene glycol have been investigated using molecular dynamics simulation. The main aspect of this theoretical work focuses on the achievement of valuable data on the reaction pathways of PEG-IFNs and the transition state energy. Our results provide a new perspective on the interactions, chemical properties and reaction pathways between IFNs and PEG.

Keywords: interaction, interferons, molecular dynamics simulation, polyethylene glycol

Procedia PDF Downloads 199

Authors: Sumit Sharma, Rakesh Chandra, Pramod Kumar, Navin Kumar

Abstract:

Molecular dynamics (MD) simulation has been used to study the effect of carbon nanofiber (CNF) volume fraction (Vf) and aspect ratio (l/d) on mechanical properties of CNF reinforced polypropylene (PP) composites. Materials Studio 5.5 has been used as a tool for finding the modulus and damping in composites. CNF composition in PP was varied by volume from 0 to 16%. Aspect ratio of CNF was varied from l/d=5 to l/d=100. To the best of the knowledge of the authors, till date there is no study, either experimental or analytical, which predict damping for CNF-PP composites at the nanoscale. Hence, this will be a valuable addition in the area of nanocomposites. Results show that with only 2% addition by volume of CNF in PP, E11 increases 748%. Increase in E22 is very less in comparison to the increase in E11. With increase in CNF aspect ratio (l/d) till l/d=60, the longitudinal loss factor (η11) decreases rapidly. Results of this study have been compared with those available in literature.

Keywords: carbon nanofiber, elasticity, mechanical properties, molecular dynamics

Procedia PDF Downloads 455

Authors: G. V. Novikov, V. S. Sivozhelezov, S. S. Kolesnikov, K. V. Shaitan

Abstract:

The study reports about the influence of binding of orthosteric ligands as well as point mutations on the conformational dynamics of β-2-adrenoreceptor. Using molecular dynamics simulation we found that there was a little fraction of active states of the receptor in its apo (ligand free) ensemble corresponded to its constitutive activity. Analysis of MD trajectories indicated that such spontaneous activation of the receptor is accompanied by the motion in intracellular part of its alpha-helices. Thus receptor’s constitutive activity directly results from its conformational dynamics. On the other hand the binding of a full agonist resulted in a significant shift of the initial equilibrium towards its active state. Finally, the binding of the inverse agonist stabilized the receptor in its inactive state. It is likely that the binding of inverse agonists might be a universal way of constitutive activity inhibition in vivo. Our results indicate that ligand binding redistribute pre-existing conformational degrees of freedom (in accordance to the Monod-Wyman-Changeux-Model) of the receptor rather than cause induced fit in it. Therefore, the ensemble of biologically relevant receptor conformations is encoded in its spatial structure, and individual conformations from that ensemble might be used by the cell in conformity with the physiological behaviour.

Keywords: seven-transmembrane receptors, constitutive activity, activation, x-ray crystallography, principal component analysis, molecular dynamics simulation

Procedia PDF Downloads 225

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

Abstract:

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

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

Procedia PDF Downloads 126

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

Abstract:

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

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

Procedia PDF Downloads 248

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 474