key: cord-0942897-67c395b0 authors: Du, Jianbin; Yang, Chunmei; Ma, Xiangyun; Li, Qifeng title: Insights into the conformation changes of SARS-CoV-2 spike receptor-binding domain on graphene date: 2021-11-27 journal: Appl Surf Sci DOI: 10.1016/j.apsusc.2021.151934 sha: 32b6410281f3bc06744a9a8f4011393a65f14948 doc_id: 942897 cord_uid: 67c395b0 Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been widely spread in the world, causing more than two million deaths and seriously threatening human life. Effective protection measures are important to prevent the infection and spreading of the virus. To explore the effects of graphene on the virus adsorption and its biological properties, the adsorption process of the receptor binding domain (RBD) of SARS-CoV-2 on graphene has been investigated by molecular dynamics simulations in this paper. The results show that RBD can be quickly adsorbed onto the surface of graphene due to [Formula: see text] stacking and hydrophobic interactions. Residue PHE486 with benzene ring has stronger adsorption force and the maximum contact area with graphene. Graphene significantly affects the secondary structure of RBD area, especially on the three key sites of binding with human ACE2, GLY476, PHE486 and ASN487. The binding free energy of RBD and graphene shows that the adsorption is irreversible. Undoubtedly, these changes will inevitably affect the pathogenicity of the virus. Therefore, this study provides a theoretical basis for the application of graphene in the protection of SARS-CoV-2, and also provides a reference for the potential application of graphene in the biomedical field. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been widely spread all over the world [1] [2] [3] [4] [5] . The diameter of SARS-CoV-2 is about 50-200 nm. The inner single stranded RNA consists of 29811 nucleotides, which is encapsulated by N (nucleocapsid) protein. E (envelope) protein and M (membrane) protein form the viral envelope, and S (spike) protein is embedded in the outermost layer of the virus [1] [2] [3] . The high affinity between the receptor binding domain (RBD) of S protein and human angiotensin converting enzyme 2 (ACE2) is the key to the strong infection of the virus [1, 6, 7] . SARS-CoV-2 has high homology with SARS-CoV appeared in [2002] [2003] and MERS-Cov appeared in 2012 [8] , possessing strong infection in human. It has rapidly caused the worldwide pandemic of coronavirus disease 2019 (COVID-19) [1] . According to the report of the World Health Organization (WHO), more than 100 million have been infected by COVID-19 until January 28, 2021, with 2.15 million deaths. In addition, the impact on economy and trade, industrial production and cultural exchange is enormous. Considering that SARS-CoV-2 is in nanometer size, nanotechnology may provide an effective way to solve this crisis [3] . As two-dimensional nano material with sp 2 hybrid, graphene has been widely studied due to its good properties and potential application prospects [9] [10] [11] . It plays a key role in the field of medical and electronic engineering since its inception [4, [12] [13] [14] [15] [16] . In the past decades, graphene based nanomaterials have been the most attractive in biosensor design due to their high affinity, low cost and ease of manufacture [17] [18] [19] . Being a single-layer sp 2 hybrid carbon atom, graphene has high surface area to volume ratio [20, 21] , it can easily detect single biomolecule. Therefore, it becomes a perfect material for sensor manufacture and implantable devices [22] . The research of cytochrome on graphene revealed a wide application prospect of graphene in 553 C biosensors and green bio-photovoltaic [23, 24] . The interaction of graphene with biomolecules has become one of the research hotspots in biomedical area. Using molecular dynamics simulation, Zhao et al. found that HIV Vpr protein could be adsorbed onto graphene and the protein conformation changed significantly during the adsorption process [25] . Zou and his team used sum frequency generation(SFG)and molecular dynamics simulation to explore the interaction of cecropin P1 and MSI78 (C1) with graphene, demonstrating that molecular dynamics simulation is effective in investigating the adsorption of biomolecules on graphene [11] . Immunoglobulin G could be rapidly expanded with conformation changes in the adsorption process of graphene [26] . Peptides with larger volume were easier to be adsorbed onto graphene [27] . The properties (size, oxidation rate, curvature, etc.) of graphene and the tension of pulmonary surfactant have significant effects on the extraction kinetics [28] . In addition, researchers found that graphene and its derivatives have good surface integrity in capturing virus [4, 29, 30] . If the structural arrangement and spatial configuration of the protein are changed, its function will be lost and may lead to virus inactivation [3] . Therefore, graphene may be used to fight SARS-CoV-2. The droplets discharged from respiratory tract of SARS-CoV-2 infected people contain virus particles. The virus in this aerosol remains active for up to 3 hours in the air [31] . People inhaling this aerosol will be infected, which is also the main path of human transmission. Medical experts advocate the use of surgical masks to prevent the virus from being inhaled [32] . Some people proposed to add graphene nanofilms into the masks to more effectively block the spread of the virus. Byeongtaek Oh et al. studied the application of graphene in promoting neuronal stem cell differentiation [33] . Graphene can be made into magnetic nanocapsules, which is able to deliver oral drugs to the stomach site-selective [34] . Graphene has large surface area ratio and delocalized electronic structure. It can be used as a loading platform for chemotherapeutic drugs  in cancer treatment. This strategy shows excellent effects in vitro and in vivo [35, 36] . Graphene nanocapsules (AuNR@G) were injected into mice caudal vein at certain concentration, which didn't cause any significant damage to cell morphology and main organs [36] . Mice orally administered with appropriate concentration of magnetic graphene capsules did not cause significant damage to cell morphology and gastric mucosa, and most of the particles were excreted within 12 h [37] . These researches indicate that graphene is feasible for in vitro protection. Whether graphene is used as a sensor to detect virus or added into protective equipment, the process of virus adsorption on graphene surface should be understood. The binding of RBD on S protein of SARS-CoV-2 with ACE2 of human [1, 6, 7] is the key for people to be infected. RBD is also the area the virus contacts outside. The adsorption process of RBD on graphene is the process of virus being captured by graphene. Thus, we investigated the whole process of RBD adsorption on graphene by molecular dynamics simulation, to excavate a significant way for graphene in virus detection and defence. To build the simulation model, the crystal structure of RBD and ACE2 complex was downloaded from Protein Data Bank (PDB: 6M0J). RBD was separated and the missed residues were completed by homologous modeling. Then, a two-dimensional graphene with the size of 118 × 89 Å 2 is constructed and placed in the X-Y plane at the bottom of the box [25, 38] . The length of the simulation system along the Z axis is 88 Å. RBD was placed on top of graphene with the outermost residue closer to graphene, which bounds to human ACE2. Residues embedded inside the virus were farther away from graphene. The initial distance between them was about 9 Å. To eliminate the influence of periodicity on the system, the distance between RBD and the edge of the simulation box was at least 10 Å. Finally, 27988 water molecules were filled into the box and solvated with SPC/E model [25, 39] . Three Clwere added into the box and randomly distributed to keep the system electrically neutral. All simulations were exceuted with Gromacs package 5.1.4 [40, 41] . The interaction between RBD and graphene was simulated by CHARMM36, which is widely used to describe the molecular dynamics behavior of proteins [42] [43] [44] [45] [46] . The C atoms in graphene were regarded as uncharged Lennard-Jones particles. The nonbonding used to describe them, which had been successfully used in the study of graphene surface adsorption [27, 42, 47] . The parameters of the nonbonding interactions between graphene, RBD and water were calculated by Lorentz-Berthelot combination rule. Periodic boundary conditions were applied in x, y and z direction during the simulation. The cutoff values of Lennard-Jones and Coulomb interactions were 10 Å. The particle mesh Ewald (PME) method was used to calculate the electrostatic interaction [48, 49] . The bond length of hydrogen atom in the system was constrained by LINCS algorithm [50] . Firstly, the energy of the system was minimized with the steepest descent method in 50000 steps. Then, pre equilibrium was carried out in the canonical ensemble for 1ns to relax the position restriction of solvent molecules on protein atoms, with V-rescale temperature coupler controlling the system temperature [51] . Finally, molecular dynamics simulation was carried out in NVT ensemble of 298K for 180 ns. There was no restriction on the protein and solvent molecules in this process, with Nose-Hoover temperature coupler [52] controlling the temperature. (On the basis of 180 ns molecular dynamics simulation, many 20 ns simulations were carried out. The results showed that RBD was stably adsorbed on the graphene surface.) The trajectory of the system was analyzed by Gromacs, and visualized by visual molecular dynamics (VMD) [53] . MM-PBSA is a method to calculate the binding free energy with the MD trajectory. It has been successfully applied into many researches in biomedical field [54] [55] [56] [57] [58] [59] . The free energy of a molecule can be expressed as and represent the free energy of the protein, graphene and pro G gra G com G their complex separately. Molecular dynamics simulation provides a possible way to clarify the effect of graphene on RBD structure. Therefore, the adsorption process of RBD on graphene surface and resulted changes in the structure were thoroughly investigated. As the adsorption of RBD on graphene is time-related, the adsorption trajectory is discussed in detail. Locations of the main nodes of RBD during the simulation are shown in Fig. 2 . As shown, RBD quickly approaches graphene in the initial stage of the simulation, to a distance about 3 Å at 3.5 ns. In other words, the initial adsorption process appears in a very short time, indicating the strong adsorption effect of graphene on RBD. At the moment, residue PHE486 has been adsorbed onto the surface of graphene, and thereafter never desorbs in the whole simulation process. It illustrates that the    stacking interaction between benzene ring of PHE486 and graphene plays an important role in the protein adsorption [60] , which is consistent with previous research results [61] [62] [63] [64] . With the increase of simulation time, residue SER477,THR488 and PRO479 are closer on to graphene. At 7.9 ns, the three amino acid residues are also adsorbed onto graphene and never desorb thereafter. From 3.5 to 7.9 ns, the distance between protein and graphene keeps about 3 Å. As shown by Fig. 2 and Fig. 4 , ASN481 begins to oscillate near graphene at 7.9 ns. The protein atoms contact the surface of graphene to the maximum extent after a long time of structural adjustment and interface rearrangement. In addition, water molecules are repulsed on the graphene surface and form a dehumidification layer to promote adsorption [65, 66] . ASN481 isn't stably adsorbed onto graphene until 95 ns. The five adsorbed residues never desorbe in 95-180 ns. Because the charge of graphene atom is 0, there is no electrostatic interaction between the protein and graphene. Therefore, van der Waals (vdW) force is the main driving force for RBD adsorption and stabilization. Hydrophobic and stacking interactions also assist in the adsorption, which are    included in vdW interaction of typical nonpolarized molecular dynamics force field [62, 67] . After the RBD was stably absorbed, the last 30 ns of the simulated trajectory was used to calculate the binding free energy by the shell script gmx_mmpbsa [57, 68] , with one frame at 1 ns interval. The results are shown in Fig. 3a . The binding free energy of RBD and graphene is -11 kJ/mol, the molecular mechanics energy is -163 kJ/mol, the polar solvation energy is 34 kJ/mol, the solvent accessible surface (SASA) energy is -19 kJ/mol and energy contributed by entropy is 137 kJ/mol. Negative values of energy is conducive to the binding of RBD and graphene, while the positive is not. The positive polar solvation energy and energy contributed by entropy are counterbalanced by the negative molecular mechanics energy and SASA energy. The binding free energy of RBD and graphene is negative, which means RBD will be adsorbed to graphene surface without desorbing. The binding energy (BE) of each residue with graphene was also calculated using gmx_mmpbsa [57, 68] , where the energy contributed by entropy was not included. The results are shown in Fig. 3b . It can be seen that the binding energies of SER477, THR488, PRO479, ASN481 and PHE486 are negative, while the binding energies of the remaining residues are basically zero. The results indicate that only these five residues are adsorbed with graphene. Among the five residues adsorbed, PHE486 has the smallest adsorption energy of -18 kJ/mol, which indicates that stacking    interaction between PHE486 and graphene is strong. and are LJ parameters between atom i and j, given by the corresponding ij  ij  force field of GROMACS. is the distance between atom i and j. As Graphene has ij r no charge, the energy between RBD and graphene is mainly LJ potential, whose variation with time is shown in Fig. 3c . It can be seen that the LJ potential decreases rapidly in the first 7.9 ns, corresponding to the process of RBD approaching the surface of graphene. After some residues of RBD were adsorbed onto graphene (such as SER477,THR488, PRO479 and PHE486), the energy curve shows a stable period of about 87 ns. LJ potential begins to decrease again at 95 ns, due to ASN481 contacting graphene with the change of RBD structure and interface. For another 10 ns, the adsorption of ASN481 is stable. The LJ potential between RBD and graphene keeps at -231 kJ/mol thereafter. In order to reveal the structural changes of RBD during simulation, the root mean square deviation (RMSD) was calculated according to Eq. (5) where is the mass of the ith atom, is the total mass of all atoms in RBD, is That is, RMSD converges. In Fig. 4a , RMSD increases rapidly in the process of RBD approaching graphene at the beginning of the simulation. When the protein is process. One part is the residues at the end of protein peptide chain (residue 331 and 530). As only one end of the these residues is connected with others, the other end can fluctuate freely. The other part is residues 477-486, which is the adsorption area of RBD on the surface of graphene. It indicates that the hydrophobic interface of graphene has a strong effect on these residues and changes their structure obviously. This domain is also the main binding region between RBD and ACE2. The contact area of RBD and some of its residues with graphene is shown in Fig. 5 . The average contact area of RBD with graphene is 129 Å 2 in 3.5-7.9 ns, and 354 Å 2 in 7.9 ns-95 ns. After 95 ns, the average contact area rests around 488 Å 2 . As seen from To further study the effect of graphene on RBD, the secondary structure of RBD adsorption region was calculated using Define Secondary Structure of Proteins (DSSP) [44, 69] , as shown in Fig. 6 It is worth noting that GLY476, PHE486 and ASN487 are the key sites for SARS-CoV-2 virus binding to human ACE2 [6, 8] . These changes in the secondary structure are bound to influence the biological properties of the virus, and even make it lose the pathogenicity. In this paper, the adsorption process of the receptor binding domain (RBD) of SARS-CoV-2 S protein with graphene was studied by molecular dynamics simulations. The results show that graphene has strong adsorption on RBD. Within 3.5 ns, RBD can be approach to and contact graphene. Five amino acid residues are adsorbed onto the surface of graphene. Compared to the other four residues, PHE486 gets stronger adsorption and more contact area, because of the stacking interaction between    benzene ring and graphene. Due to the adsorption of graphene, the secondary structure of RBD adsorption region changes significantly, with Turn formed by residues 486-487 becoming Bend, and the Bend formed by residues 476-478 becoming a Turn and a Coil. Among them, GLY476, PHE486 and ASN487 are the key sites for SARS-CoV-2 binding to ACE2. These changes of the secondary structure will inevitably affect the pathogenicity of the virus. The simulations after stable adsorption show that RBD will not desorb. The binding free energy of -11 kJ/mol also shows that the adsorption is irreversible. This is very important to prevent the spread of SARS-CoV-2. Therefore, it provides a theoretical basis for the application of graphene in defending SARS-CoV-2 and preventing its spread, as well as in the biomedical field. Computational Design of ACE2-Based Peptide Inhibitors of SARS-CoV-2 Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study Nanotechnology Responses to COVID-19 Role of graphene in biosensor and protective textile against viruses Nanomaterials and Nanotechnology-Associated Innovations against Viral Infections with a Focus on Coronaviruses Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor COVID-19: A Call for Physical Scientists and Engineers Structural and Functional Basis of SARS-CoV-2 Entry by Using Human ACE2 The electronic properties of graphene Graphene: Status and Prospects Molecular Interactions between Graphene and Biological Molecules Noncovalently functionalized monolayer graphene for sensitivity enhancement of surface plasmon resonance immunosensors Label-Free Biosensors Based on Aptamer-Modified Graphene Field-Effect Transistors Biomedical Applications of Graphene and Graphene Oxide Antiviral Activity of Graphene Oxide: How Sharp Edged Structure and Charge Matter Nanoscale graphene oxide (nGO) as artificial receptors: implications for biomolecular interactions and sensing New opportunities for emerging 2D materials in bioelectronics and biosensors Graphene biosensors for bacterial and viral pathogens Human virus detection with graphene-based materials Carbon nanotube-and graphene-reinforced multiphase polymeric composites: review on their properties and applications A review of the mechanical and thermal properties of graphene and its hybrid polymer nanocomposites for structural applications Graphene and its sensor-based applications: A review Assessing the Charge Transfer at the Cytochrome c553/Graphene Interface: A Multiscale Investigation Orientation of photosystem I on graphene through cytochrome c553 leads to improvement in photocurrent generation Molecular dynamics simulations of conformation changes of HIV-1 regulatory protein on graphene Assessing the Accuracy of Different Solvation Models To Describe Protein Adsorption Modulation of solid-waterpeptide interfacial properties towards surface adsorption/bioresistance Extracting pulmonary surfactants to form inverse micelles on suspended graphene nanosheets Carbon-based antiviral nanomaterials: graphene, C-dots, and fullerenes. A perspective, Chemical science Are graphene and graphene-derived products capable of preventing COVID-19 infection? Aerosol and Surface Stability of SARS-CoV-2 as Compared with SARS-CoV-1 It Is Time to Address Airborne Transmission of Coronavirus Disease 2019 (COVID-19) Modulating the Electrical and Mechanical Microenvironment to Guide Neuronal Stem Cell Differentiation, Advanced science In situ pepsin-assisted needle assembly of magnetic-graphitic-nanocapsules for enhanced gastric retention and mucus penetration Stable gold graphitic nanocapsule doped hydrogels for efficient photothermal antibacterial applications Simultaneous Application of Photothermal Therapy and an Anti-inflammatory Prodrug using Pyrene-Aspirin-Loaded Gold Nanorod Graphitic Nanocapsules In situ targeted MRI detection of Helicobacter pylori with stable magnetic graphitic nanocapsules Lipase adsorption on different nanomaterials: a multi-scale simulation study The Missing Term in Effective Pair Potentials Immobilization of laccase on a graphene interface: Direct electron transfer and molecular dynamics study GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers Entropy-Enthalpy Compensation in Peptide Adsorption on Solid Surfaces: Dependence on Surface Hydration Quantitatively Identifying the Roles of Interfacial Water and Solid Surface in Governing Peptide Adsorption Protein WW domain denaturation on defective graphene reveals the significance of nanomaterial defects in nanotoxicity Accelerated Molecular Dynamics Study of the Effects of Probing nano-patterned peptide self-organisation at the aqueous graphene interface Molecular dynamics simulations of the adsorption of bone morphogenetic protein-2 on surfaces with medical relevance Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems A smooth particle mesh Ewald method LINCS: A Linear Constraint Solver forMolecular Simulations Canonical sampling through velocity rescaling A unified formulation of the constant temperature molecular dynamics methods VMD: Visual Molecular Dynamics Fast Identification of Possible Drug Treatment of Coronavirus Disease-19 (COVID-19) through Computational Drug Repurposing Study Open Source Drug Discovery, A. Lynn, g_mmpbsa--a GROMACS tool for high-throughput MM-PBSA calculations py: An Efficient Program for End-State Free Energy Calculations Accurate Evaluation on the Interactions of SARS-CoV-2 with Its Receptor ACE2 and Antibodies CR3022/CB6*, Chinese Physics Letters Recent Developments and Applications of the MMPBSA Method An in-silico evaluation of different bioactive molecules of tea for their inhibition potency against non structural protein-15 of SARS-CoV-2 Binding of blood proteins to carbon nanotubes reduces cytotoxicity Probing carbon nanotube-amino acid interactions in aqueous solution with molecular dynamics simulations Assessment of standard force field models against high-quality ab initio potential curves for prototypes of pi-pi, CH/pi, and SH/pi interactions Modeling the Binding Affinity of Peptides for Graphitic Surfaces. Influences of Aromatic Content and Interfacial Shape Noncovalent interaction of carbon nanostructures Surface induced collapse of Abeta1-42 with the F19A replacement following adsorption on a single walled carbon nanotube Dewetting transition assisted clearance of (NFGAILS) amyloid fibrils from cell membranes by graphene DABCO-Metalloporphyrin Binding: Ternary Complexes, Host-Guest Chemistry, and the Measurement of π-π Interactions Dictionary of Protein Secondary Structure: Pattern Recognition of Hydrogen-Bonded and Geometrical Features The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.