> top > projects > LitCovid-sentences > docs > PMC:7605337 > annotations
LitCovid-sentences  

PMC:7605337 JSONTXT 17 Projects

Annnotations TAB TSV DIC JSON TextAE-old TextAE

Id Subject Object Predicate Lexical cue
T1 0-129 Sentence denotes Critical Sequence Hotspots for Binding of Novel Coronavirus to Angiotensin Converter Enzyme as Evaluated by Molecular Simulations
T2 131-139 Sentence denotes Abstract
T3 140-369 Sentence denotes The novel coronavirus (nCOV-2019) outbreak has put the world on edge, causing millions of cases and hundreds of thousands of deaths all around the world, as of June 2020, let alone the societal and economic impacts of the crisis.
T4 370-601 Sentence denotes The spike protein of nCOV-2019 resides on the virion’s surface mediating coronavirus entry into host cells by binding its receptor binding domain (RBD) to the host cell surface receptor protein, angiotensin converter enzyme (ACE2).
T5 602-870 Sentence denotes Our goal is to provide a detailed structural mechanism of how nCOV-2019 recognizes and establishes contacts with ACE2 and its difference with an earlier severe acute respiratory syndrome coronavirus (SARS-COV) in 2002 via extensive molecular dynamics (MD) simulations.
T6 871-996 Sentence denotes Numerous mutations have been identified in the RBD of nCOV-2019 strains isolated from humans in different parts of the world.
T7 997-1139 Sentence denotes In this study, we investigated the effect of these mutations as well as other Ala-scanning mutations on the stability of the RBD/ACE2 complex.
T8 1140-1306 Sentence denotes It is found that most of the naturally occurring mutations to the RBD either slightly strengthen or have the same binding affinity to ACE2 as the wild-type nCOV-2019.
T9 1307-1412 Sentence denotes This means that the virus had sufficient binding affinity to its receptor at the beginning of the crisis.
T10 1413-1532 Sentence denotes This also has implications for any vaccine design endeavors since these mutations could act as antibody escape mutants.
T11 1533-1759 Sentence denotes Furthermore, in silico Ala-scanning and long-timescale MD simulations highlight the crucial role of the residues at the interface of RBD and ACE2 that may be used as potential pharmacophores for any drug development endeavors.
T12 1760-1939 Sentence denotes From an evolutional perspective, this study also identifies how the virus has evolved from its predecessor SARS-COV and how it could further evolve to become even more infectious.
T13 1941-1953 Sentence denotes Introduction
T14 1954-2089 Sentence denotes The novel coronavirus (nCOV-2019) outbreak emerging from China has become a global pandemic and a major threat for human public health.
T15 2090-2442 Sentence denotes According to the World Health Organization as of August 28th 2020, there has been about 25 million confirmed cases and approaching 900,000 deaths because of coronavirus in the world.1,2 Much of the human population including the United States of America were under lockdown or official stay-at-home orders to minimize the continued spread of the virus.
T16 2443-2511 Sentence denotes Coronaviruses are a family of single-stranded enveloped RNA viruses.
T17 2512-3063 Sentence denotes Phylogenetic analysis of coronavirus genome has shown that nCOV-2019 belongs to the beta-coronavirus family, which also includes Middle East respiratory syndrome coronavirus (MERS-COV), severe acute respiratory syndrome coronavirus (SARS-COV), and bat-SARS-related coronaviruses.3,4 It is worth mentioning that SARS-COV, which was widespread in 2002 caused more than 8000 cases and about 800 deaths and MERS-COV in 2012 also spread in more than 25 countries, causing about 2500 cases and more than 850 deaths. (www.who.int/health-topics/coronavirus).5
T18 3064-3348 Sentence denotes In all coronaviruses, a homotrimeric spike glycoprotein on the virion’s envelope mediates coronavirus entry into host cells through a mechanism of receptor binding followed by fusion of viral and host membranes.3,6 Coronavirus spike protein contains two functional subunits S1 and S2.
T19 3349-3810 Sentence denotes The S1 subunit is responsible for binding to host cell receptor, and the S2 subunit is responsible for fusion of viral and host cell membranes.3,7 The spike protein in nCOV-2019 exists in a meta-stable pre-fusion conformation that undergoes a substantial conformational rearrangement to fuse the viral membrane with the host cell membrane.7,8 nCOV-2019 is closely related to bat coronavirus RaTG13 with about 93.1% sequence similarity in the spike protein gene.
T20 3811-4032 Sentence denotes The sequence similarity of nCOV-2019 and SARS-COV is less than 80% in the spike sequence.2 The S1 subunit in the spike protein includes a receptor binding domain (RBD) that recognizes and binds to the host cells receptor.
T21 4033-4736 Sentence denotes The RBD of nCOV-2019 shares 72.8% sequence identity to SARS-COV RBD and the root mean squared deviation (RMSD) for the structure between the two proteins is 1.2Å, which shows the high structural similarity.4,8,9 Experimental binding affinity measurements using surface plasmon resonance (SPR) have shown that nCOV-2019 spike protein binds its receptor human angiotensin converter enzyme (ACE2) with 10 to 20 fold higher affinity than SARS-COV binding to ACE2.7 Based on the sequence similarity between RBD of nCOV-2019 and SARS-COV and also the tight binding between the RBD of nCOV-2019 and ACE2, it is most probable that nCOV-2019 uses this receptor on human cells to gain entry into the body.3,6,7,10
T22 4737-4849 Sentence denotes The spike protein and specifically the RBD in coronaviruses have been a major target for therapeutic antibodies.
T23 4850-5083 Sentence denotes However, no monoclonal antibodies targeted to RBD have been able to bind efficiently and neutralize nCOV-2019.7,11 The core of nCOV-2019 RBD is a five-stranded antiparallel β-sheet with connected short α-helices and loops (Figure 1).
T24 5084-5183 Sentence denotes The binding interface of nCOV-2019 and SARS-COV with ACE2 is very similar with less than 1.3Å RMSD.
T25 5184-5342 Sentence denotes An extended insertion inside the core containing short strands, α-helices, and loops called the receptor binding motif (RBM) makes all the contacts with ACE2.
T26 5343-5472 Sentence denotes In nCOV-2019 RBD, the RBM forms a concave surface with a ridge loop on one side and it binds to a convex exposed surface of ACE2.
T27 5473-5542 Sentence denotes The overlay of SARS and nCOV-2019 RBD proteins is shown in Figure 1A.
T28 5543-5653 Sentence denotes The binding interface in nCOV-2019 contains loops L1 to L4 and short β-strands β5 and β6 and a short helix α5.
T29 5654-5762 Sentence denotes The location of RBM in nCOV-2019 RBD as well as different helices, strands, and loops is shown in Figure 1B.
T30 5763-5940 Sentence denotes Figure 1 (A) Superposition of the RBD of SARS-COV (yellow) and nCOV-2019 (red). (B) Different regions in the binding domain of nCOV-2019 defining the extended loop (nonyellow).
T31 5941-6073 Sentence denotes The sequence alignment between SARS-COV in human, SARS civet, Bat RaTG13 coronavirus, and nCOV-2019 in the RBM is shown in Figure 2.
T32 6074-6151 Sentence denotes There is a 50% sequence similarity between the RBM of nCOV-2019 and SARS-COV.
T33 6152-6323 Sentence denotes RBM mutations played an important role in the SARS epidemic in 2002.3,12 Two mutations in the RBM of SARS-2002 from SARS-Civet were observed from strains of these viruses.
T34 6324-6365 Sentence denotes These two mutations were K479N and S487T.
T35 6366-6470 Sentence denotes These two residues are close to the virus binding hotspots in ACE2 including hotspot-31 and hotspot-353.
T36 6471-6599 Sentence denotes Hotspot-31 centers on the salt-bridge between K31-E35 and hotspot-353 are centered on the salt-bridge between K353-E358 on ACE2.
T37 6600-7287 Sentence denotes Residues K479 and S487 in SARS-Civet are in close proximity with these hotspots and mutations at these residues caused SARS to bind ACE2 with significantly higher affinity than SARS-civet and played a major role in civet-to-human and human-to-human transmission of SARS coronavirus in 2002.3,13−15 Numerous mutations in the interface of SARS-COV RBD and ACE2 from different strains of SARS isolated from humans in 2002 have been identified and the effect of these mutations on binding ACE2 has been investigated by SPR.14,16 Two identified RBD mutations (Y442F and L472F) increased the binding affinity of SARS-COV to ACE2 and two mutations (N479K, T487S) decreased the binding affinity.
T38 7288-8239 Sentence denotes It was demonstrated that these mutations were viral adaptations to either human or civet ACE2.14,16 A pseudotyped viral infection assay of the interaction between different spike proteins and ACE2 confirmed the correlation between high affinity mutants and their high infection.16 Further investigation of RBD residues in binding of SARS-COV and ACE2 was performed through ala-scanning mutagenesis, which resulted in identification of residues that reduce binding affinity to ACE2 upon mutation to alanine.17 RBD mutations have also been identified in MERS-COV, which affected their affinity to receptor (DPP4) on human cells.14 Multiple monoclonal antibodies have been developed for SARS since 2002 that neutralized the spike glycoprotein on the SARS-COV surface.18−22 However, multiple escape mutations exist in the RBD of SARS-COV that affect neutralization with antibodies, which led to the use of a cocktail of antibodies as a robust treatment.23
T39 8240-8333 Sentence denotes Figure 2 Sequence comparison of the RBM in SARS-2002, SARS-civet, Bat RaTG13, and nCOV-2019.
T40 8334-8393 Sentence denotes Mutations from SARS-2002 to nCOV-2019 are marked with blue.
T41 8394-8444 Sentence denotes Important mutations in RBM are marked with yellow.
T42 8445-8550 Sentence denotes Red color shows the three-residue motif in SARS and civet and four-residue motif in RaTG13 and nCOV-2019.
T43 8551-8709 Sentence denotes Full genome analysis of nCOV-2019 in different countries and the receptor binding surveillance have shown multiple mutations in the RBD of glycosylated spike.
T44 8710-8834 Sentence denotes The GISAID database24 (www.gisaid.org/) contains genomes on nCOV-2019 from researchers across the world since December 2019.
T45 8835-9104 Sentence denotes The latest report by the GISAID database on June 2020 has shown 25 different variants of RBD from strains of nCOV-2019 collected from different countries along with the number of occurrences in these regions which is listed below for the seven most occurring mutations:
T46 9105-9413 Sentence denotes 213x N439K (211 Scotland, England, and Romania), 65x T478I in England, 30x V483A (26 USA/WA, 2 USA/UN, USA/CT, and England), 10x G476S (8 USA/WA, USA/OR, and Belgium), 7x S494P (3 USA/MI, England, Spain, India, and Sweden), 5x V483F (4x Spain and England), and 4x A475V (2 USA/AZ, USA/NY, and Australia/NSW).
T47 9414-9513 Sentence denotes It is not known whether these mutations are linked to the severity of coronavirus in these regions.
T48 9514-9737 Sentence denotes Starr and co-workers25 performed a deep mutational scanning of nCOV-2019 RBD and used flow cytometry to measure the effect of single mutations on the expression of the folded protein as well as its binding affinity to ACE2.
T49 9738-9878 Sentence denotes They showed that RBD is very tolerant to these mutations to maintain its expression level as well as binding affinity to ACE2 in most cases.
T50 9879-9994 Sentence denotes According to their results, most natural mutations exert similar binding affinities to ACE2 as wild-type nCOV-2019.
T51 9995-10228 Sentence denotes Furthermore, they showed that mutations at critical positions at the RBD-ACE2 interface at nCOV-2019 such as residues Q493 and Q498 do not reduce the binding affinity to ACE2 which shows the substantial plasticity of the interface.25
T52 10229-10418 Sentence denotes Different groups have computationally studied the binding of nCOV-2019 RBD with ACE2.25−29 All these studies point to the higher binding affinity of nCOV-2019 RBD than SARS-COV RBD to ACE2.
T53 10419-10765 Sentence denotes Interestingly, the role of water-mediated interactions has been pointed out to be a driving force which is shown to be similar for both SARS-COV and nCOV-2019 RBD.27 Spinello and co-workers30 studied the binding of nCOV-2019 and SARS-COV RBD to ACE2 and found that the former binds its receptor with 30 kcal/mol higher affinity than SARS-COV RBD.
T54 10766-10932 Sentence denotes Gao et al.31 used free energy perturbation (FEP) and showed that most amino acid mutations at the RBM from SARS-COV to nCOV-2019 increase the affinity of RBD to ACE2.
T55 10933-11215 Sentence denotes The focus of this article is to elucidate the differences between the interface of SARS-COV and nCOV-2019 with ACE2 to understand with atomic resolution the interaction mechanism and hotspot residues at the RBD/ACE2 interface using long-timescale molecular dynamics (MD) simulation.
T56 11216-11404 Sentence denotes An alanine-scanning mutagenesis in the RBM of nCOV-2019 helped to identify the key residues in the interaction, which could be used as potential pharmacophores for future drug development.
T57 11405-11585 Sentence denotes Furthermore, we performed molecular simulations on the seven most common mutations found from the surveillance of RBD mutations N439K, T478I, V483A, G476S, S494P, V483F, and A475V.
T58 11586-11731 Sentence denotes From an evolutionary perspective this study shows the residues in which the virus might further evolve to be even more dangerous to human health.
T59 11733-11740 Sentence denotes Methods
T60 11742-11784 Sentence denotes Sequence Comparison and Mutant Preparation
T61 11785-12453 Sentence denotes nCOV-2019 shares 76% sequence similarity with SARS-2002 spike protein, 73% sequence identity for RBD and 50% for the RBM.1 Bat coronavirus RaTG13 seems to be the closest relative of nCOV-2019 sharing about 93% sequence identity in the spike protein.6 The sequence alignment of SARS-2002 (accession number: AFR58742), SARS-civet (accession number: AY304486.1), Bat RaTG13 (accession number: MN996532.1), and nCOV-2019 (accession number: MN908947.1) is shown in Figure 2.6 To investigate the roles of critical mutations on the complex stability of nCOV-2019 with ACE2, mCSM-PPI2 webserver32 was used to find the residues in nCOV-2019 that are at the interface with ACE2.
T62 12454-12724 Sentence denotes This method uses a graph-based signature framework and predicts the effect of alanine substitution at interface residues on the binding energy of the complex, and 21 different residues were identified to be in contact with ACE2 and were chosen for further MD simulation.
T63 12725-13001 Sentence denotes On the other hand, mutations are also observed in the RBD domain from full genome analysis of different nCOV-2019 variants collected from different countries compiled in the GISAID database.24 The mutations selected are listed in Table S1 along with their location in the RBD.
T64 13003-13017 Sentence denotes MD Simulations
T65 13018-13151 Sentence denotes The crystal structure of nCOV-2019 in the complex with hACE2 (pdb id:6M0J)17 as well as the SARS-COV complex with human ACE2 (pdb id:
T66 13152-13308 Sentence denotes 6ACJ)33 were obtained from RCSB (www.rcsb.org).34 The RBD domain of nCOV-2019 comprises 194 residues (333–526) and SARS-COV includes 190 residues (323–512).
T67 13309-13395 Sentence denotes ACE2 protein contains 597 residues (19–615) in the complex structure for both systems.
T68 13396-13651 Sentence denotes All the structures including nCOV-2019, SARS-COV, and all 21 alanine substitutions of nCOV-2019 were prepared and solvated in GROMACS.35 A TIP3P water model was used for the solvent and Param99SB-ILDN AMBER force field36,37 was used for all the complexes.
T69 13652-13845 Sentence denotes A few counter ions were added to each system to neutralize the charges on the RBD and ACE2 as the PBSA method for binding energy calculation is known to be problematic with high ionic strength.
T70 13846-13888 Sentence denotes Each system contained about 260,000 atoms.
T71 13889-13981 Sentence denotes It is important to note that, none of the RBD/ACE2 complexes studied here were glycosylated.
T72 13982-14101 Sentence denotes The glycosylation sites on RBD are far from the binding interface and do not interfere with the binding of RBD to ACE2.
T73 14102-14282 Sentence denotes Moreover, there are nine Cys residues at the RBD of nCOV-2019 and eight of them form four pairs of disulfide bonds (Cys336-Cys361, Cys379-Cys432, Cys391-Cys525, and Cys480-Cys488).
T74 14283-14375 Sentence denotes In total, 5000 steps of energy minimizations were done using the steepest descent algorithm.
T75 14376-14530 Sentence denotes In all steps, the LINCS algorithm was used to constrain all bonds containing hydrogen atoms and a time step of 2 fs was used as the integration time step.
T76 14531-14589 Sentence denotes Equilibration of all systems was performed in three steps.
T77 14590-14966 Sentence denotes In the first step, 100,000 steps of simulation were performed using a velocity-rescaling thermostat to maintain the temperature at 310 K with a 0.1 ps coupling constant in an NVT ensemble under periodic boundary conditions and harmonic restraints on the backbone and side-chain atoms of the complex.38 A velocity rescaling thermostat was used in all other steps of simulation.
T78 14967-15287 Sentence denotes In the next step, we performed 300,000 steps in the isothermal-isobaric NPT ensemble at a temperature of 310 K and pressure of 1 bar using a Berendsen barostat.39 This was done through decreasing the force constant of the restraint on the backbone and side-chain atoms of the complex from 1000 to 100 and finally to 10 .
T79 15288-15408 Sentence denotes The Berendsen barostat was only used for the equilibration step because of its usefulness in rapidly correcting density.
T80 15409-15546 Sentence denotes In the next step, the restraints were removed, and the systems were subjected to 1,000,000 steps of MD simulation under the NPT ensemble.
T81 15547-15817 Sentence denotes In the production run, harmonic restraints were removed and all the systems were simulated using a NPT ensemble where the pressure was maintained at 1 bar using the Parrinello-Rahman barostat40 with a compressibility of 4.5 × 10–5bar–1 and a coupling constant of 0.5 ps.
T82 15818-16320 Sentence denotes It is important to note that the Berendsen barostat was only used for the equilibration step as it was shown that this barostat can cause unrealistic temperature gradients.41 The production run lasted for 500 ns for SARS-COV and nCOV-2019 complexes and 300 ns for all the mutants with a 2 fs timestep and the particle-mesh Ewald42 for long range electrostatic interactions using the GROMACS 2018.3 package.43 All mutant systems were constructed as described before and ran for 300 ns of production run.
T83 16321-16428 Sentence denotes In addition, the simulation time for a few mutants (Y449A, T478I, Y489A, and S494P) was extended to 500 ns.
T84 16430-16470 Sentence denotes Gibbs Free Energy and Correlated Motions
T85 16471-16877 Sentence denotes The last 400 ns of simulation was used to explore the dominant motions in SARS-COV, nCOV-2019 and the mutations with extended simulation, and last 200 ns for all other mutants using principal component analysis (PCA) as part of the quasiharmonic analysis method.44 For this method the rotational and translational motions of RBD of all systems were eliminated by fitting to a reference (crystal) structure.
T86 16878-17121 Sentence denotes Next, 4000 snapshots from the last 400 ns of SARS-COV, nCOV-2019 and mutations with extended simulation time, and 2000 snapshots from last 200 ns of all other mutant systems were taken to generate the covariance matrix between Cα atoms of RBD.
T87 17122-17208 Sentence denotes In the mutant systems with production run, the last 400 ns was used for this analysis.
T88 17209-17542 Sentence denotes Diagonalization of this matrix resulted in a diagonal matrix of eigenvalues and their corresponding eigenvectors.43,45 The first eigenvector which indicate the first principal component was used to visualize the dominant global motions of all complexes through porcupine plots using the (PorcupinePlot.tcl) script in visual MD (VMD).
T89 17543-17645 Sentence denotes The principal components were used to calculate and plot the approximate free energy landscape (aFEL).
T90 17646-18006 Sentence denotes We refer to the FEL produced by this approach to be approximate in that the ensemble with respect to the first few PC’s (lowest frequency quasiharmonic modes) is not close to convergence, but the analysis can still provide valuable information and insight. g_sham, g_covar, and g_anaeig functions in GROMACS35 were used to obtain principal components and aFEL.
T91 18007-18258 Sentence denotes The dynamic cross-correlation maps (DCCM) were obtained using the MD_TASK package to identify the correlated motions of RBD residues.46 In DCCM, the cross-correlation matrix Cij is obtained from displacement of backbone Cα atoms at a time interval Δt.
T92 18259-18436 Sentence denotes The DCCM was constructed using the last 400 ns of SARS-COV and nCOV-2019 and the extended mutant systems and last 200 ns of all other mutant systems with a 100 ps time interval.
T93 18437-18568 Sentence denotes Hydrogen bonds were analyzed in VMD where the distance cutoff was 3.2Å and the angle cutoff between the donor and acceptor was 30°.
T94 18570-18652 Sentence denotes Binding Free Energy from Molecular Mechanics Poisson-Boltzmann Surface Area Method
T95 18653-19028 Sentence denotes The molecular mechanics Poisson-Boltzmann surface area (MMPBSA) method was applied to calculate the binding energy between RBD and ACE2 in all complexes.47,48 For SARS-COV and nCOV-2019, 200 snapshots of the last 400 ns and for the mutant systems and 100 snapshots of the last 200 ns simulation were used for the calculation of binding free energies with an interval of 2 ns.
T96 19029-19226 Sentence denotes Simulation for a few mutant systems (Y449A, T478I, Y489A, and S494P) was extended to 400 ns, and the binding energies were calculated for the last 400 ns to assess the convergence of free energies.
T97 19227-19310 Sentence denotes The binding free energy of a ligand–receptor complex can be calculated as:1 2 3 4 5
T98 19311-19518 Sentence denotes In these equations, ΔEMM, ΔGbind, solv, and −TΔS are calculated in the gas phase. ΔEMM is the gas phase molecular mechanical energy changes which includes covalent, electrostatic, and van der Waals energies.
T99 19519-19723 Sentence denotes Based on previous studies, the entropy change during binding is small and neglected in these calculations.48−50 ΔGbind, solv is the solvation free energy which comprises the polar and nonpolar components.
T100 19724-19857 Sentence denotes The polar solvation is calculated using the MMPBSA method by setting a value of 80 and 2 for solvent and solute dielectric constants.
T101 19858-19979 Sentence denotes The nonpolar free energy is simply estimated from the solvent accessible surface area (SASA) of the solute from the eq 5.
T102 19981-19988 Sentence denotes Results
T103 19990-20009 Sentence denotes Structural Dynamics
T104 20010-20243 Sentence denotes To compute the RMSD of systems, the rotational and translational movements were removed by first fitting the Cα atoms of the RBD to the crystal structure and then computing the RMSD with respect to the Cα atoms of RBD in each system.
T105 20244-20333 Sentence denotes Figure 3 shows the RMSD plot in the RBD of SARS-COV, nCOV-2019, and some of its variants.
T106 20334-20454 Sentence denotes Comparison of the RMSD of SARS and nCOV-2019 RBD shows that SARS-COV has a larger RMSD throughout the 500 ns simulation.
T107 20455-20624 Sentence denotes In nCOV-2019, the RMSD is very stable with a value of about 1.5 Å, whereas in SARS-COV, the RMSD increases up to ∼4 Å after 100 ns and then fluctuates between 3 and 4 Å.
T108 20625-20730 Sentence denotes The change in RMSD of SARS is partially related to the motion in the C-terminal which is a flexible loop.
T109 20731-20812 Sentence denotes Figure 3 Cα RMSD plots for nCOV-2019 and SARS-COV and a few nCOV-2019 mutations.
T110 20813-20893 Sentence denotes The RMSD plots for the nCOV-2019 mutants show similar behaviors to nCOV-2019-wt.
T111 20894-21035 Sentence denotes In most of the variants, the RMSD is very stable during the 300 ns simulation which shows the great tolerance of the interface for mutations.
T112 21036-21087 Sentence denotes However, a few mutations showed some RMSD variance.
T113 21088-21175 Sentence denotes In mutation Y489A, the RMSD increases from 1.37 ± 0.21 Å to1.88 ± 0.16 Å after 2000 ns.
T114 21176-21288 Sentence denotes Mutation Y505A resulted in an increase in RMSD up to 100 ns to a value to 1.98 ± 0.20 Å and decreased afterward.
T115 21289-21376 Sentence denotes The RMSD for mutation N487A shows an increasing behavior with a value of 2.10 ± 0.23 Å.
T116 21377-21442 Sentence denotes Mutations N439K, V483A, and V483F showed a stable RMSD of ∼1.5 Å.
T117 21443-21519 Sentence denotes For mutations T478I, G476S, S494P, and A475V, the RMSD increases up to ∼2 Å.
T118 21520-21622 Sentence denotes These variations in the backbone RMSD show the involvement of these residues in the complex stability.
T119 21623-21677 Sentence denotes RMSD plots for other mutations are shown in Figure S1.
T120 21678-21879 Sentence denotes Since the extended loop (residues 449 to 510 shown in Figure 1B from α4 to α5) of the RBD makes all contacts with ACE2, the RMSD was computed by also fitting to the Cα atoms of this region (Figure S2).
T121 21880-22050 Sentence denotes The extended loop in nCOV-2019 is very stable with less deviation (RMSD = 0.86 ± 0.017 Å) from the crystal structure compared to SARS-COV having an RMSD of 2.79 ± 0.05 Å.
T122 22051-22126 Sentence denotes Few of the mutants show an increase in the loop RMSD during the simulation.
T123 22127-22274 Sentence denotes In mutants N487A and Y449A the loop RMSD jumps to a value of about 1.95 ± 0.12 Å and 1.94 ± 0.24 Å, respectively, after about 200 ns of simulation.
T124 22275-22395 Sentence denotes Mutants G447A and E484A show a loop RMSD values of 2.22 ± 0.03 Å and 1.96 ± 0.02 Å during the last 100 ns of simulation.
T125 22396-22482 Sentence denotes Other mutants showed a stable extended loop (observed in loop RMSD) during simulation.
T126 22483-22583 Sentence denotes The stability of extended loop for mutant systems confirms the high tolerance of this region of RBD.
T127 22584-22720 Sentence denotes To characterize the dynamic behavior for each amino acid in the RBD, we analyzed the root mean square fluctuation (RMSF) of all systems.
T128 22721-22864 Sentence denotes The RMSF plots for nCOV-2019, SARS-COV, and four other mutations are shown in Figure 4. nCOV-2019 shows less fluctuations compared to SARS-COV.
T129 22865-23007 Sentence denotes L3 in nCOV-2019 corresponding to residues 476 to 487 (shown in red in Figure 4) has smaller RMSF (1.5 Å) than SARS-COV L3 residues 463 to 474.
T130 23008-23090 Sentence denotes L1 in both nCOV-2019 and SARS-COV (green) has small fluctuation (less than 1.5 Å).
T131 23091-23186 Sentence denotes Moreover, the C-terminal residues of SARS-COV show high fluctuation (Figure 4 shown in orange).
T132 23187-23233 Sentence denotes Few mutants show higher fluctuation in the L1.
T133 23234-23331 Sentence denotes Mutants Y505A and S494A had a RMSF of 2.5 Å and mutation N487A had a RMSF of about 4 Å in the L1.
T134 23332-23388 Sentence denotes Mutation Y449A has a higher RMSF of about 3 Å in the L3.
T135 23389-23470 Sentence denotes Mutants G496A, E484A, and G447A show a high fluctuation of about 4.5 Å in the L3.
T136 23471-23520 Sentence denotes The RMSF of other variants is shown in Figure S3.
T137 23521-23601 Sentence denotes Figure 4 RMSF plots for nCOV-2019-wt, SARS-COV, Y505A, N487A, G496A, and E484A.
T138 23602-23708 Sentence denotes The red shaded region shows the fluctuation in L1 and the green shaded region shows the fluctuation in L3.
T139 23709-23786 Sentence denotes The orange shaded region in SARS-COV shows the fluctuation in the C-terminal.
T140 23787-23860 Sentence denotes For comparison, the RMSF of nCOV-2019-wt is shown as cyan in other plots.
T141 23862-23874 Sentence denotes PCA and aFEL
T142 23875-23975 Sentence denotes To identify the dominant motions in the nCOV-2019, SARS-COV, and all the mutants, PCA was performed.
T143 23976-24181 Sentence denotes Most of the combined motions were captured by the first ten eigenvectors generated from the last 400 ns for SARS-COV, nCOV-2019, and extended mutant systems and the last 200 ns for other nCOV-2019 mutants.
T144 24182-24296 Sentence denotes The percentage of the motions captured by the first three eigenvectors was 51% for nCOV-2019 and 68% for SARS-COV.
T145 24297-24390 Sentence denotes In all mutations, more than 50% of the motions were captured by the first three eigenvectors.
T146 24391-24551 Sentence denotes The first few PC’s describe the highest motions in a protein which are related to a functional motion such as binding or unbinding of protein from the receptor.
T147 24552-24751 Sentence denotes The first three eigenvectors were used to calculate the aFEL using the last 400 ns of simulation for nCOV-2019 and SARS-COV as shown in Figure 5, which displays the variance in conformational motion.
T148 24752-24850 Sentence denotes SARS-COV showed two distinct low free energy states shown as blue separated by a metastable state.
T149 24851-24952 Sentence denotes There is a clear separation between the two regions by a free energy barrier of about 6–7.5 kcal/mol.
T150 24953-25068 Sentence denotes These two states correspond to the loop motions in the L3 as well as the motion in C-terminal residues of SARS-COV.
T151 25069-25217 Sentence denotes The L3 motion in nCOV-2019 is stabilized by the H-bond between N487 on RBD and Y83 on ACE2 as well as a π-stacking interaction between F486 and Y83.
T152 25218-25414 Sentence denotes It is evident that the nCOV-2019 RBD is more stable than SASR-COV RBD and exists in one conformation whereas the SARS-COV interface fluctuates and the aFEL is separated into two different regions.
T153 25415-25636 Sentence denotes The first two eigenvectors were used to calculate and plot the aFEL as a function of first two principal components using the last 200 ns of the simulation for mutant systems. aFEL for other systems is shown in Figure S4.
T154 25637-25796 Sentence denotes Figure 5 Mapping of the principal components of the RBD for the aFEL from the last 400 ns of simulations for SARS-COV (top row) and nCOV-2019-wt (bottom row).
T155 25797-25855 Sentence denotes The color bar is relative to the lowest free energy state.
T156 25856-26134 Sentence denotes In each system, the first eigenvector was used to construct the porcupine plots to visualize the most dominant movements (Figure S5). nCOV-2019 showed a small motion in L3 and the core region and the extended loop region are very rigid showing small cones in the porcupine plot.
T157 26135-26240 Sentence denotes The core structure of the RBD remains dormant as the cones are blue in most of the regions (Figure S5-A).
T158 26241-26310 Sentence denotes In SARS-COV, the C-terminal region shows large motions (Figure S5-B).
T159 26311-26368 Sentence denotes Mutation N487A showed a large motion in L1 (Figure S5-C).
T160 26369-26447 Sentence denotes Mutations Y449A, G447A, and E484A demonstrate large motions in L3 (Figure S5).
T161 26448-26557 Sentence denotes Overall, these plots show the involvement of these residues in the dynamic stability of the RBD/ACE2 complex.
T162 26559-26563 Sentence denotes DCCM
T163 26564-26810 Sentence denotes The correlated motions of RBD atoms were also analyzed with the DCCM based on the Cα atoms of RBD from the last 400 ns of simulation for nCOV-2019, SARS-COV, and extended mutant systems and the last 200 ns for the other mutant systems (Figure 6).
T164 26811-27010 Sentence denotes The DCCM for nCOV-2019 showed a correlation between residues 490–505 (containing α5, L4 and β5 regions) and residues 440–455 (containing α4, L1 and β5 regions) shown in the red rectangle in Figure 6.
T165 27011-27098 Sentence denotes This correlation showed the coordination of these regions for binding ACE2 effectively.
T166 27099-27217 Sentence denotes Another important correlation that appears in the DCCM of nCOV-2019 is between residues 473–481 with residues 482–491.
T167 27218-27319 Sentence denotes These residues are in L3 and β6 regions and their correlation in nCOV-2019 is stronger than SARS-COV.
T168 27320-27431 Sentence denotes This is due to the presence of the β6 strand in nCOV-2019, whereas in SARS-COV these residues all belong to L3.
T169 27432-27621 Sentence denotes This indicates that L3 in nCOV-2019 has evolved from SARS-COV to adopt a new secondary structure, which causes strong correlation and makes the loop act as a recognition region for binding.
T170 27622-27685 Sentence denotes The correlation in L3 is shown as a blue rectangle in Figure 6.
T171 27686-27782 Sentence denotes Some of the mutations disrupted the patterns of correlation and anticorrelation in nCOV-2019-wt.
T172 27783-27871 Sentence denotes Mutation N487A showed a stronger correlation in L3 and β6 strand than the wild-type RBD.
T173 27872-27943 Sentence denotes In mutation E484A, correlation in L3 is stronger than the nCOV-2019-wt.
T174 27944-27990 Sentence denotes DCCM for other mutants are shown in Figure S6.
T175 27991-28162 Sentence denotes It is worth mentioning that mutation F486A disrupts the DCCM of nCOV-2019 by introducing strong correlations in the core region of RBD as well as the extended loop region.
T176 28163-28312 Sentence denotes Residue F486 resides in L3 and plays a crucial role in stabilizing the recognition loop by making a π-stacking interaction with residue Y83A on ACE2.
T177 28313-28404 Sentence denotes Figure 6 DCCM for nCOV-2019, SARS-COV, and mutants with residue numbers of the RBD domain.
T178 28405-28494 Sentence denotes Red boxes show the correlation between α5, L4, and β5 regions and α4, L1, and β5 regions.
T179 28495-28548 Sentence denotes Blue boxes show the correlation in L3 and β6 regions.
T180 28550-28571 Sentence denotes Binding Free Energies
T181 28572-29076 Sentence denotes The binding energetics between ACE2 and the RBD of SARS-COV, nCOV-2019, and all its mutant complexes were investigated by the MMPBSA method.48 The binding energy was partitioned into its individual components including: electrostatic, van der Waals, polar solvation, and SASA to identify important factors affecting the interface of RBD and ACE2 in all complexes. nCOV-2019 has a total binding energy of −50.22 ± 1.93 kcal/mol, whereas SARS-COV has a much higher binding energy of −18.79 ± 1.53 kcal/mol.
T182 29077-29319 Sentence denotes Decomposition of binding energy to its components show that the most striking difference between nCOV-2019 and SARS-COV is the electrostatic contribution which is −746.69 ± 2.66 kcal/mol for nCOV-2019 and −600.14 ± 7.65 kcal/mol for SARS-COV.
T183 29320-29624 Sentence denotes This high electrostatic contribution is compensated by a large polar solvation free energy which is 797.30 ± 3.12kcal/mol for nCOV-2019 and 659.61 ± 8.98 kcal/mol for SARS-COV. nCOV-2019 also possess a higher van der Waals (vdw) contribution (−89.93 ± 0.46 kcal/mol than SARS-COV (−70.07 ± 1.22 kcal/mol.
T184 29625-29744 Sentence denotes Furthermore, the SASA contribution to binding for SARS-COV was −8.30 ± 0.15 kcal/mol and −10.58 kcal/mol for nCOV-2019.
T185 29745-29876 Sentence denotes Both hydrophobic and electrostatic interactions play major roles in the higher affinity of nCOV-2019 RBD than SARS-COV RBD to ACE2.
T186 29877-30133 Sentence denotes The binding free energies for nCOV-2019 and SARS-COV were decomposed into a per-residue based binding energy to find the residues that contribute strongly to the binding and are responsible for higher binding affinity of nCOV-2019 than SARS-COV (Figure 7).
T187 30134-30253 Sentence denotes Most of the residues in the RBM of nCOV-2019 had more favorable contribution to the total binding energy than SARS-COV.
T188 30254-30404 Sentence denotes Residues Q498, Y505, N501, Q493, and K417 in nCOV-2019 RBM contributed more than 5 kcal/mol to binding affinity and are crucial for complex formation.
T189 30405-30494 Sentence denotes A few residues such as E484 and S494 contributed unfavorably to the total binding energy.
T190 30495-30615 Sentence denotes Among all the interface residues K417 had the highest contribution to the total binding energy (−12.34 ± 0.23 kcal/mol).
T191 30616-30786 Sentence denotes The corresponding residue in SARS-COV is V404 only had a −0.02 ± 0.01 kcal/mol contribution, which points to the importance of this residue for nCOV-2019 binding to ACE2.
T192 30787-30945 Sentence denotes Residue Q498 contributed −6.72 ± 0.18 kcal/mol and its corresponding residue in SARS-COV is a Y484 that contributed to total binding by −1.83 ± 0.06 kcal/mol.
T193 30946-31125 Sentence denotes Other important residues Y505 and N501 have more negative contribution to total binding energy than their counterparts in SARS-COV residues Y491 and T487, respectively (Figure 7).
T194 31126-31346 Sentence denotes Residue D480 in SARS-COV contributed positively to binding energy by 6.2 ± 0.15 kcal/mol and the corresponding residue in nCOV-2019 which is a S494 residue lowered this positive contribution to only 1.17 ± 0.06 kcal/mol.
T195 31347-31755 Sentence denotes Mutation D480A/G appeared to be a dominant mutation in SARS-COV in 2002–2003.51 This mutation was reported to escape neutralization by antibody 80R.52 To investigate the effect of this point mutation on binding of SARS-COV RBD to ACE2 we performed an additional simulation and calculated the binding affinity for this mutant in SARS-COV RBD with the same approach for other mutation in this study (Figure 8).
T196 31756-31887 Sentence denotes D480A mutation showed a binding affinity of 23.46 ± 3.07 kcal/mol which is about 5 kcal/mol higher than the wild-type SARS-COV RBD.
T197 31888-32095 Sentence denotes In SARS-COV, residue R426 had the highest contribution to the total binding energy (−6.27 ± 0.22 kcal/mol although the corresponding residue in nCOV-2019 is N439 with a contribution of −0.32 ± 0.02 kcal/mol.
T198 32096-32241 Sentence denotes These important mutations on RBM of nCOV-2019 from SARS-COV caused RBD of nCOV-2019 to bind ACE2 with much stronger (about 30 kcal/mol) affinity.
T199 32242-32331 Sentence denotes Figure 7 Binding energy decomposition per residue for the RBM of nCOV-2019 and SARS-COV.
T200 32332-32404 Sentence denotes Figure 8 Total free binding energy of SARS-COV, nCOV-2019, and mutants.
T201 32405-32454 Sentence denotes Natural mutants are shown with X at the bar base.
T202 32455-32561 Sentence denotes Binding free energy decomposition to its individual components for all mutants is represented in Table S2.
T203 32562-32748 Sentence denotes In all complexes, a large positive polar solvation free energy disfavors the binding and complex formation, which is compensated by a large negative electrostatic free energy of binding.
T204 32749-32788 Sentence denotes All variants had similar SASA energies.
T205 32789-32902 Sentence denotes The vdw free energy of binding ranged from −84.68 ± 0.68 kcal/mol for Q493A to −103.85 ± 0.66 kcal/mol for Y489A.
T206 32903-33140 Sentence denotes Mutant K417A had the lowest electrostatic contribution to binding −415.67 ± 5.07 kcal/mol and mutants N439K and E484A had the highest electrostatic binding contribution of −989.80 ± 5.6 kcal/mol and −941.20 ± 3.95 kcal/mol, respectively.
T207 33141-33303 Sentence denotes Most alanine substitutions exhibited similar or lower total binding affinities to nCOV-2019, however a few mutants had higher binding affinity than the wild type.
T208 33304-33434 Sentence denotes Mutant Y489A had a total binding energy of −61.78 ± 2.59 kcal/mol which was about 11 kcal/mol lower than wild type binding energy.
T209 33435-33532 Sentence denotes Mutants G446A, G447A, and T478I also demonstrated higher total binding affinities than nCOV-2019.
T210 33533-33618 Sentence denotes Other alanine substitutions had similar or lower total binding energy than nCOV-2019.
T211 33619-33734 Sentence denotes Mutant G502A has the lowest binding affinity among all the mutants with a binding energy of −24.31 ± 2.98 kcal/mol.
T212 33735-33880 Sentence denotes Mutant systems K417A, L455A, T500A, and N501A are the other mutants with total binding affinities significantly lower than the wild type complex.
T213 33881-33989 Sentence denotes The electrostatic component of binding contributes the most to the low binding affinities for these mutants.
T214 33990-34116 Sentence denotes The contribution of RBM residues to binding with ACE2 for nCOV-2019 was mapped to the RBD structure and is shown in Figure 9B.
T215 34117-34262 Sentence denotes Figure 9 (A) H bonds between RBD of nCOV-2019 and SARS-COV. (B) Mapping contribution of interface residues to structure in the RBD of nCOV-2019.
T216 34263-34304 Sentence denotes The RBD is purple and the ACE2 is yellow.
T217 34305-34482 Sentence denotes The RBD in contact with AC2 is rendered in a surface format with more red being a favorable contribution to binding (more negative) and blue unfavorable contribution (positive).
T218 34483-34595 Sentence denotes Most natural mutants exhibited similar binding affinities compared to wild-type nCOV-2019 with a few exceptions.
T219 34596-34767 Sentence denotes Mutation T478I which is one of the most frequent mutations based on the GISAID database has a binding affinity which is about 6 kcal/mol higher than that of the wild-type.
T220 34768-34852 Sentence denotes S494P and A475V showed a slightly lower binding affinity than the wild-type complex.
T221 34853-34926 Sentence denotes Other natural mutants showed binding affinities similar to wild-type RBD.
T222 34927-35126 Sentence denotes N439K demonstrated a high electrostatic energy which is compensated by large polar solvation energy and this mutant has a total binding energy of −48.27 ± 3.07 kcal/mol which is similar to nCOV-2019.
T223 35128-35188 Sentence denotes Hydrogen Bond, Salt-Bridge, and Hydrophobic Contact Analysis
T224 35189-35493 Sentence denotes Important hydrogen bonds (H-bonds) and salt bridges between nCOV-2019 RBD or SARS-COV RBD and ACE2 for the last 400 ns of trajectory are shown in Table 1. nCOV-2019 RBD makes 10 H-bonds/1 salt bridge with ACE2, whereas SARS-COV makes only 5 H-bonds/1 salt bridge with ACE2 with more than 30% persistence.
T225 35494-35599 Sentence denotes Table 1 H-Bonds and Salt-Bridges between nCOV-2019 and ACE2 and SARS-COV and ACE2 that Persist for >30%a
T226 35600-35660 Sentence denotes # nCOV-2019 ACE2 % occupancy SARS-COV ACE2 % occupancy
T227 35661-35693 Sentence denotes 1 G502 K353 89 Y436 D38 96
T228 35694-35729 Sentence denotes 2 Q493 E35 83 R426 E329 87
T229 35730-35762 Sentence denotes 3 N487 Y83 80 T486 D355 83
T230 35763-35795 Sentence denotes 4 Q498 D38 73 G488 K353 80
T231 35796-35830 Sentence denotes 5 K417 D30 55 N479 K31 52
T232 35831-35863 Sentence denotes 6 T500 D355 53 Y440 H34 47
T233 35864-35889 Sentence denotes 7 Y505 E37 52      
T234 35890-35916 Sentence denotes 8 Q498 K353 49      
T235 35917-35942 Sentence denotes 9 Y449 D38 45      
T236 35943-35970 Sentence denotes 10 G496 K353 37      
T237 35971-35997 Sentence denotes 11 Q493 K31 32      
T238 35998-36030 Sentence denotes a Salt bridge is shown as bold.
T239 36031-36145 Sentence denotes The evolution of the coronavirus from SARS-COV to nCOV-2019 has reshaped the interfacial hydrogen bonds with ACE2.
T240 36146-36214 Sentence denotes G502 in nCOV-2019 has a persistent H-bond with residue K353 on ACE2.
T241 36215-36296 Sentence denotes This residue was G488 in SARS-COV, which also makes the H-bond with K353 on ACE2.
T242 36297-36373 Sentence denotes Q493 in nCOV-2019 makes H-bond with E35 and another H-bond with K31 on ACE2.
T243 36374-36457 Sentence denotes This residue was an N479 in SARS-COV, which only makes one H-bond with K31 on ACE2.
T244 36458-36551 Sentence denotes An important mutation from SARS-COV to nCOV-2019 is residue Q498, which was Y484 in SARS-COV.
T245 36552-36662 Sentence denotes Q498 makes two H-bonds with residues D38 and K353 on ACE2, whereas Y484 in SARS-COV does not make any H-bonds.
T246 36663-36807 Sentence denotes Importantly, a salt bridge between K417 and D30 in the nCOV-2019/ACE2 complex contributes to the total binding energy by −12.34 ± 0.23 kcal/mol.
T247 36808-36918 Sentence denotes This residue is V404 in SARS-COV which is not able to make any salt-bridge and does not make H-bond with ACE2.
T248 36919-37064 Sentence denotes Gao et al.27 used a FEP approach and showed that mutation V404 to K417 lowers the binding energy of nCOV-2019 RBD to ACE2 by −2.2 ± 0.9 kcal/mol.
T249 37065-37156 Sentence denotes A salt bridge between R426 on RBD and E329 on ACE2 stabilizes the complex in SARS-COV/ACE2.
T250 37157-37250 Sentence denotes This residue is N439 in nCOV-2019 which is unable to make salt-bridge with ACE2 residue E329.
T251 37251-37422 Sentence denotes One of the most observed mutations in nCOV-2019 according to the GISAID database is N439K which recovers some of the electrostatic interactions with ACE2 at this position.
T252 37423-37497 Sentence denotes Y436 in SARS-COV and Y449 in nCOV-2019 both make H-bonds with D38 on ACE2.
T253 37498-37624 Sentence denotes The unchanged T486 in SARS-COV corresponds to T500 in nCOV-2019, both of which make consistent H-bonds with ACE2 residue D355.
T254 37625-37727 Sentence denotes Hydrophobic interactions also play an important role in stabilizing the RBD/ACE2 complex in nCOV-2019.
T255 37728-37848 Sentence denotes An important interaction between nCOV-2019 RBD and ACE2 is the π-stacking interaction between F486 (RBD) and Y83 (ACE2).
T256 37849-37951 Sentence denotes This interaction helps in stabilizing L3 in nCOV-2019 compared to SARS-COV where this residue is L472.
T257 37952-38098 Sentence denotes It was observed by Gao et al.26 that mutation L472 to F486 in nCOV-2019 results in a net change in the binding free energy of −1.2 ± 0.2 kcal/mol.
T258 38099-38238 Sentence denotes Other interfacial residues in nCOV-2019 RBD that participate in the hydrophobic interaction with ACE2 are L455, F456, Y473, A475, and Y489.
T259 38239-38328 Sentence denotes It is interesting to note that all these residues except Y489 have mutated from SARS-COV.
T260 38329-38581 Sentence denotes Spinello and co-workers30 performed long-timescale (1μs) simulation of nCOV-2019/ACE2 and SARS-COV ACE2 and found that L3 in nCOV-2019 is more stable due to presence of the β6 strand and existence of two H-bonds in L3 (H-bonds G485-C488 and Q474-G476).
T261 38582-38757 Sentence denotes Importantly, an amino acid insertion in L3 makes this loop longer than L3 in SARS-COV and enables it to act like a recognition loop and make more persistent H-bonds with ACE2.
T262 38758-38918 Sentence denotes L455 in nCOV-2019 RBD is important for hydrophobic interaction with ACE2 and mutation L455A lowers the vdw contribution of binding affinity by about 5 kcal/mol.
T263 38919-38992 Sentence denotes The H-bonds between RBD of nCOV-2019 and SARS-COV are shown in Figure 9A.
T264 38993-39116 Sentence denotes The structural details discussed here are in agreement with other structural studies of the nCOV-2019 RBD/ACE2 complex.4,53
T265 39117-39256 Sentence denotes H-bond analysis was also performed for the mutant systems and the results for H-bonds with more than 40% consistency are shown in Table S3.
T266 39257-39364 Sentence denotes Few of the alanine substitutions increase the number of interfacial H-bonds between nCOV-2019 RBD and ACE2.
T267 39365-39467 Sentence denotes Interestingly, the ala-substitution at Y489A increased the number of H-bonds in the wild-type complex.
T268 39468-39644 Sentence denotes Mutation in some of the residues having consistent H-bonds in the wild type complex such as Q498A and Q493A, stunningly maintain the number of H-bonds in the wild-type complex.
T269 39645-39812 Sentence denotes This indicates that the plasticity in the network of H-bonds in RBM of nCOV-2019 can reshape the network and strengthen other H-bonds upon mutation in these locations.
T270 39813-39894 Sentence denotes However, few mutations decrease the number of H-bonds from the wild-type complex.
T271 39895-40014 Sentence denotes Alanine substitution at residue G502 has a significant effect on the network of H-bonds between nCOV-2019 and SARS-COV.
T272 40015-40106 Sentence denotes This residue locates at the end of L4 loop near two other important residues Q498 and T500.
T273 40107-40158 Sentence denotes This mutation breaks the H-bonds at these residues.
T274 40159-40259 Sentence denotes Mutation K417A decreases the number of H-bonds to only 5 where the H-bond at residue Q498 is broken.
T275 40260-40394 Sentence denotes This indicates the delicate nature of the H-bond from residue Q498 which can easily be broken upon ala-substitution at other residues.
T276 40395-40490 Sentence denotes Furthermore, mutation N487 also decreases the number of H-bonds by breaking the H-bond at Q498.
T277 40492-40502 Sentence denotes Discussion
T278 40503-40656 Sentence denotes In this work, we preformed MD simulations to unveil the detailed molecular mechanism for the receptor binding of nCOV-2019 and compared it with SARS-COV.
T279 40657-40769 Sentence denotes The role of key residues at the interface of nCOV-2019 with ACE2 was investigated by computational ala-scanning.
T280 40770-40936 Sentence denotes A rigorous 500 ns MD simulation was performed for nCOV-2019, SARS-COV, and few mutants (Y449, T478I, Y489A, and S494P) as well as 300 ns MD simulation on each mutant.
T281 40937-41154 Sentence denotes These simulations aid in understanding the dynamic role of RBD/ACE2 interface residues and estimating the binding free energy of these variants, which shed light on crucial residues for the RBD/ACE2 complex stability.
T282 41155-41467 Sentence denotes Moreover, numerous mutations have been identified in the RBD of different nCOV-2019 strains from all over the world not known to be critical for infection.54 The effect of these mutations on the stability of the RBD/ACE2 complex was investigated to shed light on their role in the viral infection of coronavirus.
T283 41468-41727 Sentence denotes Changes in the RBD structure of nCOV-2019, SARS-COV, and mutants from their crystal structure were analyzed by RMSD and RMSF. nCOV-2019 showed a stable structure with a RMSD =1.5 Å, whereas SARS-COV had a larger RMSD value between 3–4 Å during the simulation.
T284 41728-41802 Sentence denotes Most mutations of nCOV-2019 maintained similar stability to the wild-type.
T285 41803-41934 Sentence denotes However, a few nCOV-2019 mutations resulted in larger deviations (>2 Å), i.e., Y489A, F456A, Y505A, N487A, K417A, Y473A, and Y449A.
T286 41935-42179 Sentence denotes We further investigated the structure of the extended loop domain (Figure 1B) and discovered that nCOV-2019 is stable with an RMSD of less than 1 Å, whereas the extended loop in SARS-COV shows an RMSD of about 3 Å during simulation (Figure S2).
T287 42180-42225 Sentence denotes Some mutants showed high RMSD in this region.
T288 42226-42305 Sentence denotes Alanine-substitution at residue N487 increased the extended loop RMSD to 2.5 Å.
T289 42306-42399 Sentence denotes Other mutations that increased the extended loop RMSD (>2 Å) include Y449A, G477A, and E484A.
T290 42400-42490 Sentence denotes The dynamic behavior of RBD was further investigated by analyzing the RMSF of all systems.
T291 42491-42566 Sentence denotes As shown in Figure 4, nCOV-2019 shows less fluctuation in L3 than SARS-COV.
T292 42567-42792 Sentence denotes This is due to the presence of a four-residue motif (GQTQ) in nCOV-2019 L3, which forces the loop to adopt a compact structure by making two H-bonds (G485-C488 and Q474-G476) and thereby reducing the fluctuations in the loop.
T293 42793-42939 Sentence denotes Residues F486 and N487 play major roles in stabilizing the recognition loop by making π-stacking and H-bond interactions with residue Y83 on ACE2.
T294 42940-42999 Sentence denotes Alanine substitution at N487 introduced a large RMSF to L1.
T295 43000-43354 Sentence denotes Mutation L472 to F486 in SARS-COV was shown to favor binding by −1.2 ± 0.2 kcal/mol using FEP.26 In addition, this mutation was shown to be among the five mutations that produce a super affinity ACE2 binder based on SARS-COV RBD.6 Alanine mutations at residues Y449, G447, and E484 increased the motion in L3 characterized by a large RMSF in this region.
T296 43355-43593 Sentence denotes Using PCA, the aFEL for nCOV-2019 and SARS-COV demonstrated that the former occupies only one low energy state whereas the latter forms two distinct low energy basins separated by a metastable state with a barrier of about 6–7.5 kcal/mol.
T297 43594-43709 Sentence denotes This confirms that the level of binding for the RBD domain is weaker in SARS-COV due to the presence of two basins.
T298 43710-43844 Sentence denotes Similarly, alanine-substitution for a few residues caused the FEL to degenerate into separate multiple low energy regions (Figure S4).
T299 43845-44023 Sentence denotes Dominant motions in the RBD are visualized in Figure S5 using the first eigenvector of the PCA. nCOV-2019 and SARS-COV did not show any strong motion in the extended loop region.
T300 44024-44193 Sentence denotes Porcupine plots of alanine-mutants demonstrated that mutant N487A shows large motion in the L1 region and Y449A, G447A, and E484A showed large motions in L3 (Figure S5).
T301 44194-44404 Sentence denotes To better characterize the functional motions of RBD, DCCM for all systems are constructed and shown in Figure 6 and Figure S6. nCOV-2019 showed a large correlation between the α4-L1- β5 and α5- L4- β5 regions.
T302 44405-44495 Sentence denotes This correlation was stronger in SARS-COV and few mutants such as Y449A, G447A, and E484A.
T303 44496-44559 Sentence denotes Another important correlation in nCOV-2019 is inside L3 and β6.
T304 44560-44691 Sentence denotes This correlation is stronger in nCOV-2019 than SARS-COV due to the presence of β6 which makes the loop to adopt correlated motions.
T305 44692-44756 Sentence denotes Few mutants impact the correlation in this region such as N487A.
T306 44757-45012 Sentence denotes Interestingly, mutant F486A which is in L3 and participates in binding by π-stacking interaction with Y83 on ACE2, disrupts the DCCM of wild-type nCOV-2019 and introduces strong correlation in the extended loop region as well as the core structure of RBD.
T307 45013-45285 Sentence denotes The details of hydrogen bond and salt-bridge pattern in nCOV-2019 and SARS-COV to ACE2 (Table 1) are key to the virus attachment to the host. nCOV-2019 residues participate in 10 H-bonds/1 salt bridge with ACE2, whereas SARS-COV only has 5 H-bonds/1 salt bridge with ACE2.
T308 45286-45403 Sentence denotes This significantly contributes to ∼30 kcal/mol difference in the total binding free energy of nCOV-2019 and SARS-COV.
T309 45404-45780 Sentence denotes The binding energies calculated here for nCOV-2019 and SARS-COV (−50.22 ± 1.93 and −18.79 ± 1.53 kcal/mol, respectively) are in good agreement with the binding energies calculated using the generalized Born method by Spinello et al.30 Moreover, the patterns of H-bonds between nCOV-2019 and ACE2 were also already characterized by other groups26,30 which agrees with our work.
T310 45781-45868 Sentence denotes An important H-bond between nCOV-2019 and ACE2 is between G502 on RBD and K353 of ACE2.
T311 45869-45948 Sentence denotes G502 is in the L4 region, which is populated by 5 H-bonds between RBD and ACE2.
T312 45949-46239 Sentence denotes The contribution of this residue to the total binding energy is −2.03 ± 0.04 kcal/mol and the Ala-substitution at G502 has the highest effect on the binding energy among all the residues by lowering the total binding affinity to 24.31 ± 2.98 kcal/mol, which is the lowest among all mutants.
T313 46240-46329 Sentence denotes This mutation breaks the other H-bonds in L4 such as H-bonds from residues Q498 and T500.
T314 46330-46453 Sentence denotes This residue is preserved and corresponds to residue G488 in SARS-COV, which also makes a H-bond with residue K353 on ACE2.
T315 46454-46565 Sentence denotes Residue Q493 in nCOV-2019 participates in binding ACE2 by making two H-bonds with residues E35 and K31 on ACE2.
T316 46566-46667 Sentence denotes Q493 corresponds to residue N479 in SARS-COV, which only makes one H-bond with residue Lys31 on ACE2.
T317 46668-46754 Sentence denotes This caused Q493 to have more contribution to total binding than its counterpart N479.
T318 46755-46969 Sentence denotes However, alanine substitution at Q493 did not affect the total binding energy and this mutant had a total binding energy similar to the wild-type complex as it maintains the number H-bonds in the wild-type complex.
T319 46970-47087 Sentence denotes Residues Q498 and T500 in nCOV-2019 are crucial for binding by making H-bonds with ACE2 residues D38, D355, and K353.
T320 47088-47201 Sentence denotes Residue Q498 corresponds to residue Y484 in SARS-COV which does not make any H-bond in the SARS-COV/ACE2 complex.
T321 47202-47335 Sentence denotes Q498 contributes to binding by −6.72 ± 0.18 kcal/mol which is more than the contribution of Y484 in SARS-COV (−1.83 ± 0.06 kcal/mol).
T322 47336-47415 Sentence denotes Ala-substitution at Q498 did not show large impact on the total binding energy.
T323 47416-47520 Sentence denotes Residue T500 is conserved and corresponds to residue T486 which also makes a H-bond with Asp355 on ACE2.
T324 47521-47598 Sentence denotes Mutation of T500 to Alanine lowers the binding affinity by about 10 kcal/mol.
T325 47599-47736 Sentence denotes Residue N487 in nCOV-2019 locates in L3 and plays a crucial role in stabilizing the recognition loop by making a H-bond with Y83 on ACE2.
T326 47737-47952 Sentence denotes This residue contributes to the total binding energy of nCOV-2019 by −1.52 ± 0.06 kcal/mol, whereas its corresponding residue in SARS-COV does not show any contribution to the binding energy (−0.02 ± 0.05 kcal/mol).
T327 47953-48093 Sentence denotes This demonstrates that L3 in SARS-COV has evolved to be an important recognition loop in nCOV-2019, which participates in binding with ACE2.
T328 48094-48247 Sentence denotes Residue K417 in nCOV-2019 has the most contribution to the total binding energy (−12.34 ± 0.23 kcal/mol by making a salt-bridge with residue D30 on ACE2.
T329 48248-48387 Sentence denotes This residue is crucial for the binding of RBD and ACE2 and alanine substitution lowers the total binding energy to −29.56 ± 2.95 kcal/mol.
T330 48388-48582 Sentence denotes This salt-bridge is found to be important for the stability of the crystal structure of the RBD/ACE2 complex in nCOV-2019.4 K417 is Val404 in SARS-COV which does not participate in binding ACE2.
T331 48583-48684 Sentence denotes Another important residue in nCOV-2019 is L455 which contributes to binding by −1.86 ± 0.03 kcal/mol.
T332 48685-48831 Sentence denotes This residue is important for hydrophobic interaction with ACE2 and mutating it to alanine lowers the total binding affinity by about 17 kcal/mol.
T333 48832-48982 Sentence denotes The hydrophobic residue F456 in nCOV-2019 also has a favorable contribution to the binding energy and F456A lowers the binding affinity by 5 kcal/mol.
T334 48983-49346 Sentence denotes These results are in fair agreement with experimental binding measurements with deep mutational scanning of RBD in nCOV-2019 where they used flow cytometry for different ACE2 concentrations to measure the dissociation constant KD.25 It was shown that mutations at K417, N487, T500, and G502 are detrimental for binding to ACE2, which agrees with the results here.
T335 49347-49618 Sentence denotes These experiments showed that mutations at Q493 and Q498 do not impact the binding affinity of RBD to ACE2 which demonstrates the high plasticity of the network of H-bonds at the interface where upon mutation at these residues the network can reshape to form new H-bonds.
T336 49619-49729 Sentence denotes Mutations at hydrophobic residues L455 and F456 are shown to reduce the binding affinity in these experiments.
T337 49730-49936 Sentence denotes The total binding energy calculation of all the variants showed that mutation Y489A has the highest binding affinity among all systems which is about 11 kcal/mol stronger than that of the nCOV-2019 complex.
T338 49937-50035 Sentence denotes This residue is located in β6, which is part of the recognition region of RBD for binding to ACE2.
T339 50036-50210 Sentence denotes Removal of this bulky hydrophobic residue at the interface with ACE2 caused the extended loop to move closer to the ACE2 interface and make more H-bonds with ACE2 (Table S3).
T340 50211-50338 Sentence denotes A high electrostatic interaction energy is the reason for the higher binding energy of mutant Y489A than the wild-type complex.
T341 50339-50530 Sentence denotes It is interesting to note that among the five residues L455, F456, Y473, A475, and Y489 that make hydrophobic interactions with ACE2, Y489 is the only residue that is conserved from SARS-COV.
T342 50531-50692 Sentence denotes However, the experimental binding affinity measurements using deep mutational scanning showed that mutations at this position lower the binding affinity to ACE2.
T343 50693-50774 Sentence denotes Other alanine substitutions that increase the binding energy are G446A and G447A.
T344 50775-50909 Sentence denotes Residues G446 and G447 reside in L1 and mutation to alanine can make L1 take a more rigid form as shown in the RMSF plot in Figure S3.
T345 50910-51243 Sentence denotes However, experiment showed that these mutations have similar or lower binding affinities to ACE2 than the wild-type RBD and care must be taken when interpreting these results.25 This discrepancy could be due to force-field inaccuracy and the deficiencies in the PBSA method for the treatment of solvent in binding energy calculation.
T346 51244-51353 Sentence denotes Further studies are needed to investigate whether these mutations will increase the binding affinity to ACE2.
T347 51354-51701 Sentence denotes Deep mutational scanning using flow cytometry is a qualitative method to measure the impact of a large number of mutations of protein–protein interactions and further experiments such as SPR or isothermal titration calorimetry which are conventional methods for measuring binding affinities needed to study the effect of these mutations in detail.
T348 51702-51816 Sentence denotes Important mutations found in naturally occurring nCOV-2019 appear to influence to some extent the binding to ACE2.
T349 51817-51978 Sentence denotes Mutation T478I which is one of the most frequent mutations according to GISAID database, increases the binding affinity of nCOV-2019 to ACE2 by about 6 kcal/mol.
T350 51979-52160 Sentence denotes Mutation N439K has the highest occurrence among all strains of coronavirus in the GISAID database which demonstrated the highest electrostatic interaction among all studied systems.
T351 52161-52263 Sentence denotes This residue corresponds to R426 in SARS-COV which exerts a salt-bridge interaction with E329 on ACE2.
T352 52264-52390 Sentence denotes Mutation N439K recovers some of this ACE2 interaction; however, it exerts a binding affinity similar to that of wild-type RBD.
T353 52391-52533 Sentence denotes Contribution of important interface residues to binding affinity was compared for mutations T478I, N439K, and wild-type nCOV-2019 (Figure S7).
T354 52534-52715 Sentence denotes The most striking differences between wild-type RBD and mutation T478I are residues Y449 and Q498 which have significantly higher contribution to binding than the wild type residue.
T355 52716-52802 Sentence denotes Most other residues at the interface have similar binding affinities to the nCOV-2019.
T356 52803-52989 Sentence denotes A higher H-bond persistence is also seen for these two residues Y449 and Q498 compared to the wild type RBD which is the reason for their higher contribution to the total binding energy.
T357 52990-53074 Sentence denotes Mutation N439K has a slightly lower binding affinity to ACE2 than the wild type RBD.
T358 53075-53342 Sentence denotes Per residue binding energy decomposition showed that K439 in this system has a favorable contribution of −1.80 ± 0.15 kcal/mol to the total binding energy which is balanced by a lower contribution of K417 which resulted in a binding affinity similar to wild-type RBD.
T359 53343-53483 Sentence denotes Mutant E484A, which is also one of the observed mutations based on GISAID database, demonstrates a high electrostatic interaction with ACE2.
T360 53484-53656 Sentence denotes E484 contributes to binding by 3.56 ± 0.15 kcal/mol whereas the corresponding residue in SARS-COV, P469 contributes to binding of SARS-COV to ACE2 by −0.27 ± 0.01 kcal/mol.
T361 53657-53744 Sentence denotes This residue is close to D30 on ACE2 and has electrostatic repulsion with this residue.
T362 53745-54088 Sentence denotes Most natural mutants including N439K, A475V, G476S, V483A, V483F, E484A, and S494P showed similar or slightly lower binding affinities to ACE2 compared to wild-type complex which agrees with experimental binding measurements.25 However, the experimental binding affinity for T478I also showed similar binding affinity to the wild-type complex.
T363 54089-54287 Sentence denotes This difference could be due to the use of MMPBSA approach for calculation of polar solvation and further studies are needed to study the effect of this mutation on viral infectivity of coronavirus.
T364 54288-54378 Sentence denotes Additional sequence differences between nCOV-2019 and SARS-COV influence RBD/ACE2 binding.
T365 54379-54579 Sentence denotes Residue D480 in SARS-COV contributes negatively to total binding energy (6.25 ± 0.14 kcal/mol) and mutating this residue to S494 in nCOV-2019 lowers this negative contribution to 1.17 ± 0.06 kcal/mol.
T366 54580-54684 Sentence denotes D480 in SARS-COV is located in a region of high negative charge from residues E35, E37, and D38 on ACE2.
T367 54685-54860 Sentence denotes Electrostatic repulsion between D480 on SARS-COV and the acidic residues on ACE2 is the reason for highly negative contribution of this residue to binding of SARS-COV to ACE2.
T368 54861-54937 Sentence denotes Mutation to S494 in this location removes this highly negative contribution.
T369 54938-55100 Sentence denotes Gao and co-workers26 computed the relative free energies of binding because of mutations from the RBD-ACE2 of SARS-COV to the corresponding residues in nCOV-2019.
T370 55101-55259 Sentence denotes They used a FEP approach and showed that mutation D480S in SARS-COV changed the binding free energy by −1.9 ± 0.8 kcal/mol which is consistent with our study.
T371 55260-55481 Sentence denotes Furthermore, we performed an additional simulation on D480A mutant in SARS-COV and found that this mutation has a binding affinity of 23.46 ± 3.07 kcal/mol which is about 5 kcal/mol higher than the wild-type SARS-COV RBD.
T372 55482-55667 Sentence denotes In addition, experimental binding affinity measurements showed that mutations of S494 to an acidic residue highly reduce the binding affinity to ACE2 which confirms the hypothesis here.
T373 55668-55794 Sentence denotes To our knowledge this is first detailed molecular simulation study on the effect of mutations on binding of nCOV-2019 to ACE2.
T374 55795-56368 Sentence denotes Previous computational studies have found that nCOV-2019 binds to ACE2 with a total binding affinity which was about 30 kcal/mol stronger than SARS-COV and is in fair agreement with the results here.56 The critical role of interface residues is computationally investigated here and in other articles and the results of all the studies indicate the importance of these residues for the stability of the complex and finding hotspot residues for the interaction with receptor ACE2.26,30,55,56 It is interesting to note the role of L3 in the stability of the RBD/ACE2 complex.
T375 56369-56499 Sentence denotes The amino acid insertions in L3 for nCOV-2019 have converted an unessential part of RBD in SARS to a functional domain of the RBD.
T376 56500-56655 Sentence denotes This loop participates in binding ACE2 by making H-bond as well as π-stacking interactions with ACE2, which makes this region to act as a recognition loop.
T377 56656-56955 Sentence denotes Previous studies on SARS-COV have shown that there is a correlation between the higher binding affinity to the receptor and higher infection rate by coronavirus.6,13,57 The higher binding affinity of nCOV-2019 for ACE2 than SARS-COV to ACE2 is suggested to be the reason for its high infection rate.
T378 56956-57132 Sentence denotes Most natural mutations showed similar binding affinities to wild-type ACE2 which indicates that the virus was already effective at the beginning of the crisis for binding ACE2.
T379 57133-57224 Sentence denotes A few mutations such as N489A and T478I are shown to increase the binding affinity to ACE2.
T380 57225-57313 Sentence denotes However, more studies are needed to investigate the effect of these mutations in detail.
T381 57314-57504 Sentence denotes Mutations of nCOV-2019 RBD that do not change the binding affinity and complex stability could have implications for antibody design purposes since they could act as antibody escape mutants.
T382 57505-57691 Sentence denotes Escape from monoclonal antibodies is observed for mutations of SARS-COV in 2002 and these mutations should be considered for any antibody design endeavors against these escape mutations.
T383 57693-57704 Sentence denotes Conclusions
T384 57705-57964 Sentence denotes In conclusion, this study unraveled key molecular traits underlying the higher affinity of nCOV-2019 for ACE2 compared to SARS-COV and unveiled critical residues for the interaction by in silico alanine scanning mutations and binding free energy calculations.
T385 57965-58108 Sentence denotes The higher affinity of nCOV-2019 to binding with ACE2 correlates with higher human-to-human transmissibility of nCOV-2019 compared to SARS-COV.
T386 58109-58352 Sentence denotes Ala-scanning mutagenesis of the interface residues of nCOV-2019 RBM has shed light on the crucial interface residues and helped obtain an atomic-level understanding of the interaction between coronavirus and the receptor ACE2 on the host cell.
T387 58353-58555 Sentence denotes MD simulations on RBD mutations found in strains of nCOV-2019 from different countries aid in the understanding of how these mutations can play an important role in viral infection with ACE2 attachment.
T388 58556-58746 Sentence denotes In addition to previously reported residues, it was found that residue F486 locating in L3 plays a crucial role in the dynamic stability of the complex by a π-stacking interaction with ACE2.
T389 58747-58868 Sentence denotes Per-residue free energy decomposition pinpoints the critical role of residues K417, Y505, Q498, and Q493 in binding ACE2.
T390 58869-59067 Sentence denotes Alanine scanning of interface residues in nCOV-2019 RBD showed that alanine substitution at some residues such as G502, K417, and L455 can significantly decrease the binding affinity of the complex.
T391 59068-59233 Sentence denotes Moreover, mutation T478I, which is one of the most probable mutations in RBD of nCOV-2019 is found to bind ACE2 with about 7 kcal/mol higher affinity than wild-type.
T392 59234-59458 Sentence denotes It is also alerting that some of the alanine substitutions at residues G446, G447, and Y489 substantially increased the binding affinity that may lead to a strong RBD attachment to ACE2 and influence the infection virulence.
T393 59459-59580 Sentence denotes However, details of interaction between these mutants and ACE2 should be carefully studied using experimental techniques.
T394 59581-59805 Sentence denotes On the other hand, most mutations are found not to impact the binding affinity of RBD with ACE2 in nCOV-2019 which could have implications for vaccine design endeavors as these mutations could act as antibody escape mutants.
T395 59806-60047 Sentence denotes Receptor recognition is the first line of attack for coronavirus and this study gives novel insights into key structural features of interface residues for the advancement of effective therapeutic strategies to stop the coronavirus pandemic.
T396 60049-60081 Sentence denotes Supporting Information Available
T397 60082-60553 Sentence denotes The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpcb.0c05994.Location of residues picked for mutation; binding free energies of mutants; details of H-bonds for all mutations; RMSD of mutants; loop RMSD of all systems; RMSF plots for all mutants; FELs using first two principal components; porcupine plots for structural fluctuations; DCCM plots for mutants; and binding energy decomposition for T478I and N439K mutants (PDF)
T398 60555-60577 Sentence denotes Supplementary Material
T399 60578-60598 Sentence denotes jp0c05994_si_001.pdf
T400 60600-60652 Sentence denotes The authors declare no competing financial interest.
T401 60654-60669 Sentence denotes Acknowledgments
T402 60670-60829 Sentence denotes This work was supported by the Intramural Research Program of the National Heart, Lung, and Blood Institute (NHLBI) at the National Institutes of Health (NIH).
T403 60830-60963 Sentence denotes The authors acknowledge the Biowulf high performance Linux cluster at the NIH for providing computational resources for this project.
T404 60964-61162 Sentence denotes The authors would like to dedicate this article to the doctors and nurses who sacrificed their time, health, and even their lives to fight COVID-19, particularly those in Iran and the United States.
T405 61163-61310 Sentence denotes J.B.K. would also like to dedicate this work to family friend Joe Kaplan (Silver Spring, MD) who passed away because of COVID-19 on April 22, 2020.