key: cord-0999764-hr49rxtr authors: Karimipour, Aliakbar; Amini, Ali; Nouri, Mohammad; D’Orazio, Annunziata; Sabetvand, Roozbeh; Hekmatifar, Maboud; Marjani, Azam; Bach, Quang-vu title: Molecular dynamics performance for coronavirus simulation by C, N, O, and S atoms implementation dreiding force field: drug delivery atomic interaction in contact with metallic Fe, Al, and steel date: 2020-11-17 journal: Comput Part Mech DOI: 10.1007/s40571-020-00367-w sha: b8b2471d739e954f7de104adc554dd52308ed993 doc_id: 999764 cord_uid: hr49rxtr Coronavirus causes some illnesses to include cold, COVID-19, MERS, and SARS. This virus can be transmitted through contact with different atomic matrix between humans. So, this atomic is essential in medical cases. In this work, we describe the atomic manner of this virus in contact with various metallic matrix such as Fe, Al, and steel with equilibrium molecular dynamic method. For this purpose, we reported physical properties such as temperature, total energy, distance and angle of structures, mutual energy, and volume variation of coronavirus. In this approach, coronavirus is precisely simulated by O, C, S, and N atoms and they are implemented dreiding force field. Our simulation shows that virus interaction with steel matrix causes the maximum removing of the virus from the surfaces. After 1 ns, the atomic distance between these two structures increases from 45 to 75 Å. Furthermore, the volume of coronavirus 14.62% increases after interaction with steel matrix. This atomic manner shows that coronavirus removes and destroyed with steel surface, and this metallic structure can be a promising material for use in medical applications. Coronaviruses are a group of connected RNA viruses that origin illnesses in birds and mammals. In humans, these viruses cause some problem in respiratory system infections that can range from mild to lethal. Mild diseases contain some cases of cold and fatal varieties can cause diseases such as MERS, COVID-19, and SARS. There are any vaccines or antiviral drugs to stop or remedy human coronaviruses infections yet. Coronaviruses constitute the subfamily Orthocoronavirinae, in the family Coronaviridae, order Nidovirales, and realm Riboviria [1, 2] . Historically, coronavirus was recognized in the 1930s for the first time [3] . This virus depicted in Fig. 1 . Researchers reported a new respiratory infection of chickens in North Dakota in 1931. The infection of chicks was distinguished by listlessness and gasping. The death rate of the chickens was up to 90% [5, 6] . Transmissible gastroen- 8 Sustainable Management of Natural Resources and Environment Research Group, Faculty of Environment and Labour Safety, Ton Duc Thang University, Ho Chi Minh City, Vietnam Fig. 1 a The common structure of coronavirus recognized by doing some experiment researches. b All-atom depiction of the coronavirus simulated using the MD approach [4] teritis and mouse hepatitis virus were detected in the 1940s [7] . Chemically, these structures are enveloped viruses with a nucleocapsid of helical symmetry and a single-stranded RNA genome which located in an icosahedral protein shell [8] . This atomic structure is one of the biggest thorugh RNA viruses. The genome size of coronaviruses ranges from about 26 to 32 kilobases [4, 9] . Today, numerous researches have been done on the characteristics of this virus. Masters addressed the present state of comprehension of coronavirus genome packaging, which has mainly been studied in 2 prototype species, transmissible gastroenteritis virus, and mouse hepatitis virus [10] . Lu er al. [11] phylogenetic analysis shows that bats, an animal sold at the seafood market in Wuhan, might be the original host of this virus in human. Furthermore, the structural analysis of this research group suggested that coronavirus might be able to bind to the angiotensin-converting enzyme 2 receptors in humans. Then researchers introduce evolutionary features of the coronaviruses and successfully predicted new coronavirus outbreak by pointing out that novel pathogenic variants will readily appear from very various severe acute respiratory syndrome-related coronaviruses of the bat origin through their close coexistence and high genetic recombination potential [12] . Contrary to many studies, it is not known how the coronavirus interacts with the human environment, such as metallic structures. In this work, we study the coronavirus atomic interaction with metallic surfaces with MD [13] [14] [15] [16] [17] [18] . Investigations are provided for the coronavirus and metallic Table 2 The r 0 and θ 0 rates of coronavirus structure in our MD simulations [52] Bond/angle r 0 (Å) θ 0 (°) CC In this computational study, we used MD simulation to calculate the atomic behavior of coronavirus in contact with metallic surfaces such as Fe, Al, and steel. Molecular dynamics method is a way for calculating the dynamical evolution of atoms and molecules. In MD method, atoms and molecules are allowed to atomic interaction and give a view of the position and velocity changes of the total system. In the research, all MD simulations were done via LAMMPS [19] [20] [21] [22] . This computational package designing began in the 1990s by Sandia and LLNL laboratories. To use this package to calculate the atomic behavior of coronavirus, this atomic structure, H 2 O molecules, and metallic matrix were simulated. For this, coronavirus structure was situated in the middle region of MD simulation package, while the other regions of MD simulation box was filled by H 2 O structure. This initial atomic positions prepared via Packmol package [23] . [24] . Particle base methods such as MD, LBM [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] , or nanoparticles dispersed in the base fluid can be estimated as the various approaches of investigation through such that problem [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] . In this molecular dynamics simulations, periodic boundary conditions were used in x, y, and z directions. For preparing the initial temperature of simulated structures, the thermostat was used in the MD simulation box to equilibrate this physical parameter at 300 K, 325 K and 350 K with 1 femtosecond time step. After 1,000,000 MD simulation time steps, at the selected temperature, the simulated atomic structures were equilibrated; then, micro-canonical ensemble (NVE) was implemented to describe the virus and metallic surfaces atomic interaction. To simulate the coronavirus atomic structure, we use the dreiding force field [52] . This interatomic force field is the appropriate choice to biological mixtures MD simulation. In dreiding force field, the potential of various atoms was described as a non-bonded and bonded forces. Non-bond force between particles in dreiding force field described by the Lennard-Jones (LJ) formalism. This formalism is a mathematically simple relation about the interatomic force between a pair of particles. This simple In Eq. (1), σ is the distance at which the function is zero; ε is the depth of the potential well, and r ij is the distance between the two atoms. In MD simulations, both σ and ε constants related to the kind of the particles in the MD simulation package. The length scale parameter, cutoff radius, and energy for various atoms in coronavirus simulation are classified in Table 1 [52] . The bonded forces consist of bond angle bend, bond strength, and dihedral angle torsion terms. The bond and Computational Particle Mechanics (2) and (3) formulas: In these formulas, K θ and K r are harmonic oscillator constants. θ 0 is the equilibrium value of angles, and r 0 is the atomic bond length. K θ and K r constants in coronavirus MD study were selected at 100 (kcal/mol)/degree 2 and 300 (kcal/mol)/A 2 , respectively. Furthermore, the r 0 and θ 0 constants in our MD simulations are classified in Table 2 [52] . Dihedral term in atomic interaction described with Eq. (4) and its coefficients are chosen from a dreiding force field [52] : In Eq. (4), K is an oscillator constant with + 1 or − 1 rates and the integer number is n [52] . Today, various atomic models such as SPC, TIP4P, and TIP3P are used to molecular dynamics simulation of water atomic structure. In SPC model, three sites were used for the electrostatic interactions and the positive charges on the hydrogen atoms were neutralized by a negative charge. The interaction between H 2 O molecules were simulated using a LJ potential. These interaction parameters for atomic water structure are reported in [54, 55] : Furthermore, metallic structures in our MD simulations described by inserted atom model (EAM) force field [56] . EAM force field is represented by Eq. (5): In this equation, r ij is the distance between i and j particles, φ αβ is a pair-wise potential function, ρ αβ is the contribution to the electron density from particle j of type β at the location of particle i, and F α is an embedding function that stats the energy required to place particle i of type α into the electron distribution. After atomic modeling and force field implementing, the atomic manner of coronavirus in the vicinity of metallic matrix estimated. For computations of these atomic evolutions, Newton's second law [13, 14] , Which leads to [13, 14] : Temperature from the Gaussian distribution [13, 14] : Association of Eq. (6) is provided by velocity Verlet algorithm. (9) In these Eqs. (9) and (10), r (t + δt), v(t + δt) is the final coordination and velocity of particles and v(t) and r(t) are the rate of these physical parameters at t 0. Finally, we can say that MD simulations of this work in 2 steps [57] [58] [59] : Step A Coronavirus and metallic matrix were simulated at 300, 325, and 350 K with 1 fs time step. Step B In the second step of our molecular dynamics simulation, the atomic interaction between the virus and metallic matrix was carried out for 1,000,000 MD simulation time steps. In the first step of this molecular dynamics simulations, the atomic structure of the coronavirus and metal surfaces was studied. Figure 3 shows the atomic arrangement of structures in our simulation box [52] . Numerically, structures atomic stability are described by reporting of temperature and potential energy of them at 300 K, 325 K, and 350 K. Temperature variation of various structures is depicted in Fig. 3 . From this figure, we can say that all structures thermally equilibrated after 1,000,000 MD simulation time steps passing. Figure 4 displays the potential energy. We can say that the simulated structures potential energy converged after 1,000,000-time steps. This result shows that dreiding force field has good ability in bio structures simulations. Numerically, the potential energy of the coronavirus structure with Fe matrix decreases from − 28,545 to − 25,665 eV by temperature increasing from 300 to 350 K. This atomic manner shows the kinetic energy increasing in simulated structures which this parameter increasing cause the mean distance of atomic structure rises. It can be said that simulated structures with Fe matrix have better stability rather than Al and Steel matrix and by temperature increasing this stability decreases (see Fig. 4 ). Furthermore, coronavirus atomic stability reaches to minimum rate in the presence of steel matrix, while this atomic behavior of the coronavirus makes the steel matrix work best in medical applications. The COM of the atomic structures is one of the geometrical. The interatomic force between the coronavirus and metallic surface is repulsive (Fig. 5) . Numerically, the distance of coronavirus and Fe atoms varies from 45 to 57 Å at 300 K. This physical parameter increases by temperature increasing from 300 to 350 K (see Fig. 6 and Table 3 ). As reported in Table 3 , steel and aluminum matrix interaction with coronavirus shows a similar manner and steel matrix shows maximum atomic distance from the coronavirus after 1,000,000-time steps. Furthermore, as seen in Table 4 , the angle of the coronavirus with all matrix surface fluctuates with temperature increasing. This geometrical parameter has 92°, 95°, and 88°at 300 K, 325 K and 350 K for Fe matrix Fig. 8 Variation of coronavirus volume with time steps: a n 0, b n 250,000, n 500,000, n 750,000, and n 1,000,000 surface, respectively. So by increasing temperature, the angle of the coronavirus does not follow a logical relation. Mutual energy of two groups of simulated structures describes their atomic interaction in MD simulations. In this step of our calculations, we report the mutual energy of coronavirus and metallic matrix. Figure 7 shows that the mutual energy between the coronavirus and Fe matrix varies from 32,261/30,569/29,777 to 0 eV by temperature increasing. Zero rates of mutual energy rate show that the coronavirus structure and Fe matrix distance are bigger than the cutoff radius of them at 350 K and these structures departed from each other after 1,000,000-time steps. Coronavirus with Al and steel matrix interacts with a similar manner, and mutual energy of these structures decreases by MD simulation time steps passing (Table 5) . Between metallic matrix, steel one has minimum mutual interaction with the coronavirus. That implies the antivirus property of this structure which nominated this metal to antibacterial applications. Larger simulation time steps, volume of virus structure rises in the presence of steel atoms and so, the atomic distance increases, too (see Fig. 8 ). More the distance of atoms in the coronavirus structure, the stability of the virus decreases. From Fig. 9 by the time evolution from 0 to 1,000,000-time steps, the volume of coronaviruses increases from 179,091 to 205,283 Å 3 . Coronavirus volume does not change drastically by simulation time evolution in the presence of Fe and Al (Table 7) . Figure 10 shows the volume variation of the coronavirus structure in the presence of various metallic matrix. From these results, we can say that coronavirus atomic stability do not reduce in the presence of Fe and Al matrix. The present study investigates the atomic interaction between the coronavirus and metallic matrix such as Fe, Al, and steel with MD simulations. In our calculations, the coronavirus is shown by S,O, N, and C atoms. Result of simulations shows that steel matrix has good properties to prevent coronavirus transmission, which can be used for medical purposes. Generally, the results are as following: • Dreiding force field is the appropriate interatomic potential for MD simulation of coronavirus. Numerically, the total energy of our simulated structures converged to − 28,545 eV, − 22,835 eV, and − 17,124 eV after 1 ns at 300 K in the presence of Fe, Al, and steel matrix. • Center of the mass distance of coronavirus and metallic structures increases by the passage of MD simulations time from 45 to 75 Å where the largest rate of this physical parameter occurs for steel matrix. • Generally, by increasing the simulation temperature, the atomic repulsion between the metallic matrix and coronavirus rises. • By increasing the temperature of simulated structures, the angle of the coronavirus does not follow a logical relation • The mutual energy of coronavirus and metallic structures varies from 19,965 to 32,261 eV after 1 ns. The highest rate of this physical parameter calculated for the Fe matrix, which describes maximum interaction between the coronavirus and these metallic structures in our MD simulations. The volume of the coronavirus increases by MD simulation time passing in the presence of steel matrix. Numerically, coronavirus volume increases from 179,091 to 253,244 after atomic interaction with steel matrix. International Committee on Taxonomy of Viruses, International Union of Microbiological Societies Feigin and Cherry's textbook of pediatric infectious diseases Coronaviruses, a new group of animal RNA viruses Virology: coronaviruses The early history of infectious bronchitis Gammacoronavirus ‡: coronaviridae. The Springer Index of Viruses Coronaviruses: a comparative review Bat coronaviruses in China Coronavirus genomics and bioinformatics analysis Coronavirus genomic RNA packaging Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding Origin and evolution of pathogenic coronaviruses Studies in molecular dynamics. I. General method Correlations in the motion of atoms in liquid argon Studying temperature effects on electronic and optical properties of cubic CH 3 NH 3 SnI 3 perovskite Investigation of the oxidation mechanism of dopamine functionalization in an AZ31 magnesium alloy for biomedical applications Investigation of thermal properties of DNA structure with precise atomic arrangement via equilibrium and non-equilibrium molecular dynamics approaches Calculation of the thermal conductivity of Human Serum Albumin (HSA) with equilibrium/nonequilibrium molecular dynamics approaches Implementing molecular dynamics on hybrid high performance computers-short range forces Implementing molecular dynamics on hybrid high performance computers-particle-particle particle-mesh Fast parallel algorithms for short-range molecular dynamics Substructured molecular dynamics using multibody dynamics algorithms PACK-MOL: a package for building initial configurations for molecular dynamics simulations Visualization and analysis of atomistic simulation data with OVITO-the open visualization tool A modified two-phase mixture model of nanofluid flow and heat transfer in a 3-D curved microtube Investigation of rib's height effect on heat transfer and flow parameters of laminar water-Al 2 O 3 nanofluid in a rib-microchannel Efficacy of hybrid nanofluid in a new microchannel heat sink equipped with both secondary channels and ribs Comparison of the finite volume and lattice Boltzmann methods for solving natural convection heat transfer problems inside cavities and enclosures Heat transfer improvement of water/single-wall carbon nanotubes (SWCNT) nanofluid in a novel design of a truncated double-layered microchannel heat sink Analysis of heat transfer and nanofluid fluid flow in microchannels with trapezoidal, rectangular and triangular shaped ribs The effect of attack angle of triangular ribs on heat transfer of nanofluids in a microchannel Application of nanofluid to improve the thermal performance of horizontal spiral coil utilized in solar ponds: geometric study Flow and heat transfer in non-Newtonian nanofluids over porous surfaces Heat transfer and fluid flow of pseudo-plastic nanofluid over a moving permeable plate with viscous dissipation and heat absorption/generation Assessment of thermal conductivity enhancement of nano-antifreeze containing single-walled carbon nanotubes: optimal artificial neural network and curve-fitting A smoothed particle hydrodynamics approach for numerical simulation of nano-fluid flows: application to forced convection heat transfer over a horizontal cylinder Evaluating the effect of temperature and concentration on the thermal conductivity of ZnO-TiO 2 /EG hybrid nanofluid using artificial neural network and curve fitting on experimental data Effect of employing a new biological nanofluid containing functionalized graphene nanoplatelets on thermal and hydraulic characteristics of a spiral heat exchanger Natural convection heat transfer enhancement in new designs of plate-fin based heat sinks MHD mixed convection in a vertical annulus filled with Al 2 O 3 -water nanofluid considering nanoparticle migration Numerical study on mixed convection of a non-Newtonian nanofluid with porous media in a two lid-driven square cavity An experimental study on stability and thermal conductivity of water/silica nanofluid: eco-friendly production of nanoparticles An experimental study on heat transfer and pressure drop of water/graphene oxide nanofluid in a copper tube under air crossflow: applicable as a heat exchanger Empirical analysis of heat transfer and friction factor of water/graphene oxide nanofluid flow in turbulent regime through an isothermal pipe Entropy generation of boehmite alumina nanofluid flow through a minichannel heat exchanger considering nanoparticle shape effect Effects of graphene oxide-silicon oxide hybrid nanomaterials on rheological behavior of water at various time durations and temperatures: synthesis, preparation and stability Adaptive particle method based on moments for simulating the mass transport in natural flows Explicit incompressible smoothed particle hydrodynamics in a multi-GPU environment for large-scale simulations Unique solvability and stability analysis for incompressible smoothed particle hydrodynamics method Convergence study and optimal weight functions of an explicit particle method for the incompressible Navier-Stokes equations Investigation of the mechanical responses of copper nanowires based on molecular dynamics and coarse-grained molecular dynamics DREIDING: a generic force field for molecular simulations On the determination of molecular fields Molecular-dynamics study of atomic motions in water The missing term in effective pair potentials Embedded-atom method: derivation and application to impurities, surfaces, and other defects in metals The art of molecular dynamics simulation A unified formulation of the constant temperature molecular-dynamics methods Canonical dynamics: equilibrium phase-space distributions Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Conflict of interest There is no conflict of interest.