Next Article in Journal
An In Silico Methodology That Facilitates Decision Making in the Engineering of Nanoscale Protein Materials
Next Article in Special Issue
The 2-hydroxy-3-(4-aryl-1-piperazinyl)propyl Phthalimide Derivatives as Prodrugs—Spectroscopic and Theoretical Binding Studies with Plasma Proteins
Previous Article in Journal
Novel Antibacterial Copolymers Based on Quaternary Ammonium Urethane-Dimethacrylate Analogues and Triethylene Glycol Dimethacrylate
Previous Article in Special Issue
Design and Nonadiabatic Photoisomerization Dynamics Study of a Three-Stroke Light-Driven Molecular Rotary Motor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Dynamics Studies on the Structural Stability Prediction of SARS-CoV-2 Variants Including Multiple Mutants

1
Graduate School of New Drug Discovery and Development, Chungnam National University, 99 Daehak-ro, Yuseong-gu, Daejeon 34134, Korea
2
Division of Emerging Infectious Diseases, Bureau of Infectious Disease Diagnosis Control, Korea Disease Control and Prevention Agency, 187 Osongsaengmyeong 2-ro, Osong-eup, Heungdeok-gu, Cheongju-si 28159, Korea
3
Bureau of Infectious Disease Diagnosis Control, Korea Disease Control and Prevention Agency, 187 Osongsaengmyeong 2-ro, Osong-eup, Heungdeok-gu, Cheongju-si 28159, Korea
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(9), 4956; https://doi.org/10.3390/ijms23094956
Submission received: 24 March 2022 / Revised: 21 April 2022 / Accepted: 26 April 2022 / Published: 29 April 2022
(This article belongs to the Collection Computational Studies of Biomolecules)

Abstract

:
Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has caused the Coronavirus Disease (COVID-19) pandemic worldwide. The spike protein in SARS-CoV-2 fuses with and invades cells in the host respiratory system by binding to angiotensin-converting enzyme 2 (ACE2). The spike protein, however, undergoes continuous mutation from a D614G single mutant to an omicron variant, including multiple mutants. In this study, variants, including multiple mutants (double, triple mutants, B.1.620, delta, alpha, delta_E484Q, mu, and omicron) were investigated in patients. The 3D structure of the full-length spike protein was used in conformational analysis depending on the SARS-CoV-2 variants. The structural stability of the variant types was analyzed based on the distance between the receptor-binding domain (RBD) of each chain in the spike protein and the binding free energy between the spike protein and bound ACE2 in the one-, two-, and three-open-complex forms using molecular dynamics (MD) simulation. Omicron variants, the most prevalent in the recent history of the global pandemic, which consist of 32 mutations, showed higher stability in all open-complex forms compared with that of the wild type and other variants. We suggest that the conformational stability of the spike protein is the one of the important determinants for the differences in viral infectivity among variants, including multiple mutants.

1. Introduction

Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2), an RNA virus, was initially discovered during an outbreak in Wuhan, China and has since become a global pandemic [1]. SARS-CoV-2 invades respiratory organs with symptoms, including cough, followed by fever, myalgia, dyspnea, headache, and death in extreme cases [2]. The number of COVID-19-infected patients has increased continuously, with over 500,000 deaths in the USA (2021) [3]. The four major structural proteins of coronaviruses (CoVs) include spike (S) glycoproteins, envelope (E) proteins, membrane (M) proteins, and nucleocapsid (N) proteins, along with 16 non-structural proteins. The S protein, located in the membrane, plays an important role in viral interactions with the host as it is involved in the first contact with the host cells [4,5].
The S protein cleavage by host proteases leads to the production of the S1 subunit at the N-terminus and the S2 subunit at the C-terminus [6]. The biological mechanism for S protein attachment to host cells is caused by the receptor-binding domain (RBD) interaction of the S1 subunit with human angiotensin-converting enzyme 2 (ACE2), including proteases, such as transmembrane serine protease 2, protein cleavage enzyme furin, and extracellular matrix metalloproteinase inducer CD 147 [7].
The S1 subunit consists of various domains, namely the N-terminal, RBD, and C-terminal (CT). Depending on the RBD position state, the S protein structure changes into multiple conformational states, such as closed, one-, two-, or full three-open forms [8]. Three RBDs in full open forms can interact with and convert a maximum of three ACE2 molecules from the closed form to the open form. The biological mechanism between ACE2 and RBD has been a research topic for antigenicity under different conditions, such as temperature [9], environment [10], and species [11].
The S protein is continuously variable in protein residues. Multiple amino acid residues of the S protein have mutated evolutionally from a single major D614G mutant to various mutants, including the alpha (from England), beta, delta, mu (from Colombia) [12,13], and omicron variants [14]. In addition, the omicron and delta variants (from India) are the dominant strains globally [15]. However, structural considerations have not been actively researched. Previously, our group reported that single mutants showed higher stability compared with the wild type in all open-complex forms [16]. This study focused on the conformational stability of the various mutants.
Computational methods using released 3D structures have been conducted to investigate the antigenic evolution of viruses [17,18]. Further, these approaches can be used to study the mechanisms of the activity of proteins produced by the virus. Nelson et al. simulated the S protein RBD at the long millisecond scale to study the interaction intensity of ACE2 with the flexible loop region of RBD prior to post-fusion [19]. Deganutti et al. calculated the binding affinity for small molecules to the RBD–ACE2 complex structure [20]. Peng et al. conducted a binding free energy calculation in the system to identify differences in binding affinity of ACE2 for SARS-CoV-2 and SARS-CoV-1 resulting from the 2003 pandemic [21].
The S proteins are frequently variable to avoid the antibodies produced by the host. From April 2020, the single mutant D614G of the S protein exhibited major variations, and D614G became the dominant mutant form of the S protein during the pandemic [22]. Over time, multiple variants, such as alpha, delta, and mu, have been identified. In particular, alpha [23], and delta [24] variants include deletion mutants, whereas omicron variants include both deletion and insertion mutants unlike previous single variants for point mutations.
Although genomics was actively conducted in the pandemic-struck population, structural research has not been conducted on multiple variants. In this study, we focused on the conformational stability of variants, including multiple mutants, using computational approaches with full-length S protein-ACE2 complex structure in one-, two-, and three-open complexes. Until now, full-length structure-based analysis for S protein-ACE2 complex has not been actively studied unlike only RBD-based structural research.
In this study, the full-length SARS-CoV-2 S protein–human ACE2 complex structure, including the RBD domain, was selected to analyze the structural stability differences according to variant types, including multiple mutations, such as D614G/E484K double, D614G/E484K/N440K triple mutant, B.1.620, delta lacking E484Q mutation, alpha, delta_E484Q including E484Q mutation, mu, and omicron variants. We expect that this study of the structural stability of SARS-CoV-2 variants, including multiple mutants using the full-length S protein-ACE2 complex, may provide insights explaining the differences in viral infectivity among variant types.

2. Results

2.1. SARS-CoV-2 S Protein Genome/Protein Analysis in Patients with COVID-19

We identified multiple variants (D614G/E484K double mutant, D614G/E484K/N440K triple mutant, alpha variant, delta variant, mu variant, omicron variant, and B.1.620 lineage) related to the structure of the SARS-CoV-2 S protein using whole-genome sequence analysis from COVID-19 patients in South Korea. E484Q variants were also used, although they have not yet been identified in South Korea.

2.2. The Stability Analysis for RBDs in Full-Length Trimeric S Protein Chains

In our previous study, we investigated the structural stability of single mutants. In this study, we focused on multiple variants, and the Protein Data Bank (PDB) data for SARS-CoV-2 S proteins complexed with ACE2 were used. The MD simulation was conducted in the gas phase with a microcanonical (NVE) ensemble as an isolated system with particles (N) and energy (E) in a volume (V) because SARS-CoV-2 infects hosts from the outer air to inner respiratory systems (as does influenza virus [25]), thereby, characterizing the infection process in a variable environment for temperature and pressure.
We analyzed the stability among chains consisting of trimeric S proteins using the distance and standard deviation (SD) for the distances between two chains (AB and AC, AC and BC, and AC and BC) in each chain for V503 and N501 residues. In order to identify the conformational relationship in each chain of trimeric full-length S proteins, V503 and N501 residues were selected for analysis as in the previous study. V503 residues were involved in the closest distance below 1.0 nm with symmetry among RBDs of chains in closed form (PDB ID 6VXX) of whole trimeric S proteins.
N501 residues involved in the interaction with ACE2 were confirmed as 8.7 nm with the symmetry distance among RBDs of chains in full-open-complex form (PDB ID 7A98) [16]. These analytical methods were used to quantify the stability between the chains of full-length structures considering the V503 residues with three symmetries in the closed form and the N501 residues with three symmetries interacting with the ACE2 interphase in the open form. The mutation sites for each variant are shown in Figure 1. Additionally, the omicron variant with 32 mutation sites (A67V, ΔH69/ΔV70, T95I, G142D/ΔV143-ΔY145, ΔN211/L212I, ins214EPE, G339D, S371L, S373P, S375F, K417N, N440K, G446S, S477N, T478K, E484A, Q493K, G496S, Q498R, N501Y, Y505H, T547K, D614G, H655Y, N679K, P681H, N764K, D796Y, N856K, Q954H, N969K, and L981F) is shown in Figure S1.
Depending on the MD trajectories, the distance and SD of the V503 and N501 residues among the three chains bound to ACE2 were analyzed (Table 1) in all open-complex forms (one-, two-, and three-open). In the one-open-complex form, the SD value for distance calculations showed no differences caused by the different mutation types (<0.5 nm) (Figure 2). In the two-open-complex form, all variants (i.e., double, triple, B.1.620, alpha, delta, delta_E484Q, and omicron) showed commonly lower distance and SD values compared with the wild type. The SD values among the chains were approximately 5 and 4 nm in the wild type and variants, respectively (Figure 3). In the three-open-complex form, all variants had lower SD values compared with the wild type, and delta_E484Q had the lowest SD value among the chains. The mu variant showed also lower distance and SD values (Figure 4).
The raw data for the distance and SD for the distances between two chains (AB and AC, AC and BC, and AC and BC) in each chain for the V503 and N501 residues depending on all MD trajectories are presented in Supplementary Materials. The distance fluctuations between N501 residues showed a pattern similar to that of the V503 residues. In the one-open-complex form, there were no differences caused by the different variants.
The SD values for the distances among each chain were below 0.5 nm. In the two-open-complex form, the SD values among chains were approximately 4.5 and 4 nm in the wild type and variants (i.e., double, triple, B.1.620, alpha, delta, delta_E484Q, and omicron), respectively. The quantification of the RBDs of the full-length S protein indicated that the variants were more stable than the wild type. From the results for each open-complex form, the delta and mu variants appeared to be more stable compared with the other variants. The summaries of the quantification of V503 residue and N501 residue distances are presented in Table 1 and Table S1, respectively.

2.3. Binding Free Energy Analysis on Stability between the S Protein and ACE2

Almost all interaction energies of all RBD–ACE2 complex structures have been studied previously. Our research focused on the energy calculation for all atoms between the full-length S protein and ACE2 using MM/PBSA calculations. The results between the two groups are presented as the van der Waals energy, electrostatic energy, and total energy. Variants consisting of only point mutations and lacking deletion residues were additionally calculated at deleted residue conditions with alpha and delta variant scales. Commonly, in all open-complex forms, the variants were more stable than the wild type.
In the one-open-complex form, the delta and mu variants had lower MM/PBSA values (Figure 5 and Table 2). In the two-open-complex form, the MM/PBSA was calculated between the S protein and ACE2 for the A–D and B–E chains in the two sets. The omicron variant had the lowest stable energy (Figure S8 and Table 3). The dela_E484Q variant remained stable in both steps. In the three-open-complex form, the MM/PBSA was calculated between the S protein and ACE2 for the A–D, B–E, and C–F chains in three sets. The Delta_E484Q variant, except for omicron, showed energetic stability among various variants (Figure S9 and Table 4).
Other variants, including triple mutants and the B.1.620 variant, showed higher stability but decreased in the delta scale. The omicron variant commonly exhibited the lowest energy in all complex forms.

3. Discussion

SARS-CoV-2 is an RNA virus that is currently causing severe pandemic disease. The S protein of SARS-CoV-2 is frequently mutated. The Korea Disease Control and Prevention Agency (KDCA) investigated the S protein variants of SARS-CoV-2 in Korean patients using RT-PCR analysis and identified variants (double, triple, B.1.620, alpha, delta, mu, and omicron) in addition to the previous five major single mutants (D614G, D614A, L455F, F456L, and Q787H) depending on the domain region. In our previous study [16], we showed that the D614G mutant in the CT2 domain has the highest stability among all tested single mutants. In this study, variants, including multiple mutants, were used for the stability analysis.
Spreading variants, including double, triple, B.1.620, alpha, delta, mu, and omicron variants, were mutated in various domains comprising RBD, CT1, and CT2. The wild type commonly showed lower stability compared with the mutant types through the RBD distance and MM/PBSA values. In particular, the double, triple, B.1.620, delta_E484Q, and mu variants, including both D614G in CT2 and E484K in RBD, showed higher stability compared with the wild type, alpha, and Korean delta variants lacking E484K mutation. In the distance and SD analysis of the distances between two chains (AB and AC, AC and BC, and AC and BC) for the structural stability analysis using V503 and N501 among RBDs, the results for the wild type were similar compared with the variants in one-open-complex form but showed larger SD values in the two- and three-open-complex forms.
In the binding free energy analysis for stability between the S protein and ACE2 using MM/PBSA, the energy value for the wild type was slightly unstable in the one-open-complex form as expected according to the results from the before-mentioned analysis. The difference in the binding free energies between the wild type and variant types increased significantly in the two- and three-open-complex forms. These results show that the S protein variants, including multiple mutants, have a stronger interaction with ACE2 prior to membrane fusion compared with the wild type.
The large differences in the analysis for the distance and SD of distances between two chains among the RBDs of the S protein and in the analysis for the binding free energy between the S protein and ACE2 were conspicuous, as these progressed from the one-open form toward the three-open-complex form. These tendencies were clear in variants (double, triple mutants, B.1.620, delta_E484Q, and mu variants), including major mutations, such as D614G and E484K. Therefore, these mutations are expected to play a critical role in the conformational stability for the interaction with ACE2 in the pre-fusion state prior to the metastable state and post-fusion.
In the metastable state, which requires higher energy than the pre-fusion state, S proteins undergo a rearrangement process to enter membrane fusion between SARS-CoV-2 and the host cell [26]. The conformational stability of the S protein is required for strong interaction with ACE2 in the pre-fusion state prior to the dynamic conformational change to post-fusion [27]. These differences in structural stability in variant types of SARS-CoV-2 were related to the differences in biological function for the interaction with ACE2 in the pre-fusion state from our research results. The alpha and Korean delta variants, with a non-mutation in the E484 residue, showed higher energy than the other variants, including a mutation in the E484 residue, although this showed lower energy than the wild type.
From these structural perspective, we suggest that the mutation in E484—in addition to D614G as deduced in our previous study [16]—is one of the major determinants in conformational stability. Therefore, although the E484 residue is located in the RBD, and thus its interaction is weak in the interphase between the S protein and ACE2, our study showed that mutation in the E484 residue is related to the conformational stability considering the whole structure using full-length residues in one-, two-, and three-open-complex forms.
In addition to the importance of mutation in the E484 residue, we investigated a relationship between the deletion residues contained in variants, such as delta and alpha and the conformational stability. Deletion residues in the delta variant were reported as avoiding antibody recognition and attack [28]. The binding free energy values in double, triple mutants, B.1.620, and mu variants were decreased in the delta scale (deletion of F157 and R158 residues) but not in the alpha scale (deletion of H69, H70, and Y144). As well as the mutation in E484 residue, deletions in the F157 and R158 residues were slightly involved in conformational stability.
The recently identified omicron variant has not been studied in detail compared to other variants but has been found to be able to avoid the host immune response [29]. The studies by Koleya et al. for the omicron variant using the RBD-ACE2 complex structure resulted in a slightly strong interaction with ACE2 due to showing a lower binding free energy than the wild type [30]. Eileen et al. analyzed the interaction changes of residues in only RBD-ACE2 interphase for wild type, delta, and omicron variants to identify the altered interaction pattern, and two salt-bridge formations in RBD-ACE2 interphase for omicron variants, unlike wild type and delta variants, were found [31].
The analysis using only RBD-ACE2, however, considers no conformational stability because the S protein adopts a trimeric whole structure. In our study, we calculated the conformational stability for the omicron variant with 32 mutation sites, including both deletion and insertion in a full-length S protein. The omicron variant showed the largest conformational stability between the S protein and ACE2 in all-open-complex forms based on our results for the binding free energy value using MM/PBSA calculations.
Therefore, although the recently identified omicron variant has not been actively studied compared with other mutant types, the highest infectivity of the omicron variant among variant types is related to various causes, such as antibody evasion [29], the binding affinity in the RBD-ACE2 interphase [30,31], and the conformational stability considering the whole structure in this study. SARS-CoV-2 variants have evolved into variants, including multiple mutants, such as the omicron variant, in order to enhance the infectivity. As a result of the evolution, the omicron variant appears to have the most stable structure among variant types in connection with increased infectivity.
It is important to identify and target a large number of mutation sites with conformational stability, such as omicron variants in addition to other variants, including mutations in the D614 and E484 residues for applications in therapeutic development, such as antibody design. Omicron variants, the most prevalent in the recent global pandemic, which consist of 32 mutations, showed higher stability in all open-complex forms compared with the wild type and other variants. We suggest that the conformational stability of the spike protein is an important determinant for the differences in viral infectivity among variants, including multiple mutants.

4. Materials and Methods

4.1. Ethical Considerations

This study was approved by the Institutional Review Board at the Korea Disease Control and Prevention Agency (2020-03-01-P-A) and is considered to be a public health act to the outbreak. Thus, the board has waived the requirement for written consent as outlined in the Title Laboratory Respondence to COVID-19. All the methods presented in this study were conducted in accordance with the relevant guidelines and regulations.

4.2. Genome/Protein Sequence Investigation of Korean Patients with COVID-19

Nasopharyngeal and oropharyngeal swab specimens were collected from symptomatic patients to detect Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) by real-time reverse transcriptase-polymerase chain reaction (RT-PCR). RNA was extracted from the specimens using a Qiagen viral RNA mini kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. Next, real-time RT-PCR was performed on the cycle threshold value of the SARS-CoV-2 target gene (ORF 1b and E) [32]. For whole-genome sequencing, the cDNA was amplified using the QIAseq SARS-CoV-2 Primer Panel and QIAseq FX DNA Library UDI Kit (QIAGEN, Hilden, Germany).
DNA libraries were extracted using the Nextera DNA Flex Library Prep Kit (Illumina, San Diego, CA, USA), and sequencing was performed on the MiSeq instrument using a MiSeq reagent kit V2 (Illumina, San Diego, CA, USA) to obtain an average genome coverage greater than 1000× for all the isolates. The reads were trimmed and mapped to reference genome MN908947.3 using CLC Genomics Workbench version 20.0.3 (CLC Bio, Aarhus, Denmark). Using the reference genome, single-nucleotide variants (SNVs) were called using the BioNumerics version 7.6 SARS-CoV-2 plugin (Applied Maths, Sint-Martens-Latem, Belgium) [33].

4.3. Dataset Preparation for SARS-CoV-2 S Protein Structure

The 3D structures of the full-length SARS-CoV-2 S protein complexed with ACE2 were selected prior to the analysis. The PDB structures 7A94, 7A97, and 7A98 were selected as examples of one-, two-, and three-open complexes with ACE2. Variants investigated in Korean patients with COVID-19 were used as double, triple, B.1.620, delta, alpha, mu, and omicron variants; additionally, the delta_E484Q variant identified in other countries was investigated in this study. 3D structures for the variants, including point mutations, were built using the ‘Build and Edit Protein’ module Discovery Studio 2021 (BIOVIA, San Diego, CA, USA). Overall, 24 structures (except for the previous wild type) were built as the inputs to the MD simulation calculation.

4.4. Self-Homology Modeling for Deletion Mutant

For the modeling of deletion or insertion mutants, a full-length S protein sequence, including deletion mutants (ΔF157 and ΔR158 for delta variant, and ΔV69, ΔH70, and ΔY144 for alpha variant) and omicron variants (ΔV69, ΔH70, ΔV143, ΔY144, ΔY145, ΔN211, and Ins214EPE) was used as the input. Each open-complex PDB was used as the template structure in SWISS-MODEL [34]. The remaining point mutation residues were built using the ‘Build and Edit Protein’ module in Discovery Studio 2021 (BIOVIA, San Diego, CA, USA).

4.5. MD Simulation

All MD simulations for the above 24-input system were performed using the GROMACS 5.1.3 package [35] with the CHARMM27 all-atom force field [36]. The PDB entries 7A94 (one-open-complex form), 7A97 (two-open-complex form), and 7A98 (three-open-complex form) were used as the starting point for the simulations. MD simulations were performed at 298 K in a cubic box with molecules. The system was locally minimized using the steepest descent algorithm with 50,000 steps to ensure that the system had no steric clashes, inappropriate geometry, or structural distortions. MD production was conducted under NVE ensembles in a vacuum phase.
The Linear Constraint Solver (LINCS) algorithm was used to constrain the bond lengths. A cutoff of 10 Å was used for short-range interactions, long-range electrostatic interactions were treated using the particle-mesh Ewald (PME) summation method [37], and production runs were performed for 10 ns. The integration time-step with leap-frog was set to 1 fs, and the trajectory coordinates and energies were saved at 1 ps intervals. The built-in programs presented as gmx modules of the GROMACS software package were used for the analysis. The re-built 3D structures with self-homology modeling for delta and alpha variants due to deletion mutants conducted more minimization and additional 10 ns MD production in a total of 20 ns.

4.6. V503 and N501 Residue Distance Calculation

V503 and N501 residues in RBDs of trimeric S protein chains were selected using the “gmx index” module. The distance calculation between V503 residues and between N501 residues among A–B, A–C, and B–C chains in the S protein were conducted as the minimum distance of residues consisting of atoms in the “gmx mindist” module. Standard deviations (SDs) between V503 residues and between N501 residues in two chains and among three chains, such as A–B and A–C, A–B and B–C, A–C and B–C, and A–B–C, were calculated from the distance values depending on trajectories. The unit for the distance value was represented as nanometers (nm).

4.7. Binding Free Energy Analysis Using MM/PBSA Calculation

In order to calculate the binding free energy between the full-length S protein and ACE2 complexes, the molecular mechanics/Poisson–Boltzmann surface area (MM/PBSA) calculation was conducted using the “g_mmpbsa” module [38] in each MD trajectory file. All atom interaction energies between the whole domains of the full-length S protein (not only the RBD) and ACE2 were calculated.
The structure of 7A94 (one-open-complex form) was used for the MM/PBSA calculation between the A and D chain, and the structure of 7A97 (two-open-complex form) was used for the MM/PBSA calculation between the B and E chain in addition to calculation between the A and D chain. The structure of 7A98 (three-open-complex form) was used for the MM/PBSA calculation between the A–D, B–E, and C–F chains. Due to the long computation time, the MM/PBSA was calculated at 10 ps intervals.
The van der Waals, electrostatic, and total energies were represented as kJ/mol in the vacuum phase system. All variants in the dataset preparation were used in this binding free energy calculation. In order to correct an effect for deletion residues, such as alpha and delta variants, the binding free energy for the wild type and the variant types consisting of only point mutants, such as double, triple mutants, B.1.620, and mu variants, was also calculated with deletion condition (ΔF157 and ΔR158 for delta variant, and ΔV69, ΔH70, and ΔY144 for alpha variant) in the final trajectory.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23094956/s1.

Author Contributions

K.-E.C., J.-M.K., J.E.R., A.K.P., E.-J.K., and C.K.Y. designed the work, performed the simulations, analyzed the results, and wrote the manuscript. N.S.K. guided and approved the research work, and contributed to supervision, project administration, and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Korea Disease Control and Prevention Agency Grant (6300-6331-301) and the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (NRF2020R1A2C100691511).

Institutional Review Board Statement

2020-03-01-PA (Korea Disease Control and Prevention Agency).

Data Availability Statement

The SARS-CoV-2 whole-genome sequences described are available in GISAID (Accession IDs EPI_ISL_738142, 2361153, 2361182, 2887353, 3024839, and 8885887).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, Y.; Yang, C.; Xu, X.-f.; Xu, W.; Liu, S.-w. Structural and functional properties of SARS-CoV-2 spike protein: Potential antivirus drug development for COVID-19. Acta Pharmacol. Sin. 2020, 41, 1141–1149. [Google Scholar] [CrossRef] [PubMed]
  2. Sharma, A.; Tiwari, S.; Deb, M.K.; Marty, J.L. Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2): A global pandemic and treatment strategies. Int. J. Antimicrob. Agents 2020, 56, 106054. [Google Scholar] [CrossRef] [PubMed]
  3. Sivakumar, B.; Deepthi, B. Complexity of COVID-19 Dynamics. Entropy 2022, 24, 50. [Google Scholar] [CrossRef] [PubMed]
  4. Duan, L.; Zheng, Q.; Zhang, H.; Niu, Y.; Lou, Y.; Wang, H. The SARS-CoV-2 spike glycoprotein biosynthesis, structure, function, and antigenicity: Implications for the design of spike-based vaccine immunogens. Front. Immunol. 2020, 11, 2593. [Google Scholar] [CrossRef]
  5. Yadav, R.; Chaudhary, J.K.; Jain, N.; Chaudhary, P.K.; Khanra, S.; Dhamija, P.; Sharma, A.; Kumar, A.; Handu, S. Role of Structural and Non-Structural Proteins and Therapeutic Targets of SARS-CoV-2 for COVID-19. Cells 2021, 10, 821. [Google Scholar] [CrossRef]
  6. Wu, Y.; Zhao, S. Furin cleavage sites naturally occur in coronaviruses. Stem Cell Res. 2021, 50, 102115. [Google Scholar] [CrossRef]
  7. Radzikowska, U.; Ding, M.; Tan, G.; Zhakparov, D.; Peng, Y.; Wawrzyniak, P.; Wang, M.; Li, S.; Morita, H.; Altunbulakli, C. Distribution of ACE2, CD147, CD26, and other SARS-CoV-2 associated molecules in tissues and immune cells in health and in asthma, COPD, obesity, hypertension, and COVID-19 risk factors. Allergy 2020, 75, 2829–2845. [Google Scholar] [CrossRef]
  8. Ray, D.; Le, L.; Andricioaei, I. Distant residues modulate conformational opening in SARS-CoV-2 spike protein. BioRxiv 2021, 118, 43. [Google Scholar] [CrossRef]
  9. Rath, S.L.; Kumar, K. Investigation of the effect of temperature on the structure of SARS-CoV-2 Spike Protein by Molecular Dynamics Simulations. Front. Mol. Biosci. 2020, 7, 583523. [Google Scholar] [CrossRef]
  10. Hozhabri, H.; Piceci Sparascio, F.; Sohrabi, H.; Mousavifar, L.; Roy, R.; Scribano, D.; De Luca, A.; Ambrosi, C.; Sarshar, M. The global emergency of novel coronavirus (SARS-CoV-2): An update of the current status and forecasting. Int. J. Environ. Res. Public Health 2020, 17, 5648. [Google Scholar] [CrossRef]
  11. Piplani, S.; Singh, P.K.; Winkler, D.A.; Petrovsky, N. In silico comparison of spike protein-ACE2 binding affinities across species; significance for the possible origin of the SARS-CoV-2 virus. arXiv 2020, arXiv:2005.06199. [Google Scholar]
  12. López-Cortés, G.I.; Palacios-Pérez, M.; Zamudio, G.S.; Veledíaz, H.F.; Ortega, E.; José, M.V. Neutral evolution test of the spike protein of SARS-CoV-2 and its implications in the binding to ACE2. Sci. Rep. 2021, 11, 18847. [Google Scholar] [CrossRef]
  13. Xie, X.; Han, J.-B.; Ma, G.; Feng, X.-L.; Li, X.; Zou, Q.-C.; Deng, Z.-H.; Zeng, J. Emerging SARS-CoV-2 B. 1.621/Mu variant is prominently resistant to inactivated vaccine-elicited antibodies. Zool. Res. 2021, 42, 789. [Google Scholar] [CrossRef] [PubMed]
  14. He, X.; Hong, W.; Pan, X.; Lu, G.; Wei, X. SARS-CoV-2 Omicron variant: Characteristics and prevention. MedComm 2021, 2, 838–845. [Google Scholar] [CrossRef] [PubMed]
  15. Tian, D.; Sun, Y.; Zhou, J.; Ye, Q. The Global Epidemic of the SARS-CoV-2 Delta Variant, Key Spike Mutations and Immune Escape. Front. Immunol. 2021, 12, 751778. [Google Scholar] [CrossRef] [PubMed]
  16. Choi, K.-E.; Kim, J.-M.; Rhee, J.; Park, A.K.; Kim, E.-J.; Kang, N.S. Molecular Dynamics Studies on the Structural Characteristics for the Stability Prediction of SARS-CoV-2. Int. J. Mol. Sci. 2021, 22, 8714. [Google Scholar] [CrossRef]
  17. Martí, D.; Alsina, M.; Alemán, C.; Bertran, O.; Turon, P.; Torras, J. Unravelling the molecular interactions between the SARS-CoV-2 RBD spike protein and various specific monoclonal antibodies. Biochimie 2021, 193, 90–102. [Google Scholar] [CrossRef]
  18. Ali, A.; Vijayan, R. Dynamics of the ACE2–SARS-CoV-2/SARS-CoV spike protein interface reveal unique mechanisms. Sci. Rep. 2020, 10, 1–12. [Google Scholar] [CrossRef]
  19. Nelson, G.; Buzko, O.; Bassett, A.; Spilman, P.; Niazi, K.; Rabizadeh, S.; Soon-Shiong, P. Millisecond-scale molecular dynamics simulation of spike RBD structure reveals evolutionary adaption of SARS-CoV-2 to stably bind ACE2. bioRxiv 2020. [Google Scholar] [CrossRef]
  20. Deganutti, G.; Prischi, F.; Reynolds, C.A. Supervised molecular dynamics for exploring the druggability of the SARS-CoV-2 spike protein. J. Comput. Aided Mol. Des. 2021, 35, 195–207. [Google Scholar] [CrossRef]
  21. Peng, C.; Zhu, Z.; Shi, Y.; Wang, X.; Mu, K.; Yang, Y.; Zhang, X.; Xu, Z. Exploring the binding mechanism and accessible angle of SARS-CoV-2 spike and ACE2 by molecular dynamics simulation and free energy calculation. ChemRxiv 2021. [Google Scholar] [CrossRef]
  22. Korber, B.; Fischer, W.M.; Gnanakaran, S.; Yoon, H.; Theiler, J.; Abfalterer, W.; Hengartner, N.; Giorgi, E.E.; Bhattacharya, T.; Foley, B. Tracking changes in SARS-CoV-2 spike: Evidence that D614G increases infectivity of the COVID-19 virus. Cell 2020, 182, 812–827.e819. [Google Scholar] [CrossRef] [PubMed]
  23. Meng, B.; Kemp, S.A.; Papa, G.; Datir, R.; Ferreira, I.A.; Marelli, S.; Harvey, W.T.; Lytras, S.; Mohamed, A.; Gallo, G. Recurrent emergence of SARS-CoV-2 spike deletion H69/V70 and its role in the Alpha variant B. 1.1. 7. Cell Rep. 2021, 35, 109292. [Google Scholar] [CrossRef] [PubMed]
  24. Hamill, V.; Noll, L.; Lu, N.; Tsui, W.N.T.; Porter, E.P.; Gray, M.; Sebhatu, T.; Goerl, K.; Brown, S.; Palinski, R. Molecular detection of SARS-CoV-2 strains and differentiation of Delta variant strains. Transbound. Emerg. Dis. 2021. ahead of print. [Google Scholar] [CrossRef]
  25. Miner, J.C.; Lappala, A.; Fenimore, P.W.; Fischer, W.M.; McMahon, B.H.; Hengartner, N.W.; Sanbonmatsu, K.Y.; Tung, C.-S. Modeling the Influenza A NP-vRNA-Polymerase Complex in Atomic Detail. Biomolecules 2021, 11, 124. [Google Scholar] [CrossRef]
  26. Krumm, Z.A.; Lloyd, G.M.; Francis, C.P.; Nasif, L.H.; Mitchell, D.A.; Golde, T.E.; Giasson, B.I.; Xia, Y. Precision therapeutic targets for COVID-19. Virol. J. 2021, 18, 1–22. [Google Scholar] [CrossRef]
  27. Giron, C.C.; Laaksonen, A.; Barroso da Silva, F.L. Up state of the SARS-CoV-2 spike homotrimer favors an increased virulence for new variants. Front. Med. Technol. 2021, 3, 694347. [Google Scholar] [CrossRef]
  28. Harvey, W.T.; Carabelli, A.M.; Jackson, B.; Gupta, R.K.; Thomson, E.C.; Harrison, E.M.; Ludden, C.; Reeve, R.; Rambaut, A.; Peacock, S.J. SARS-CoV-2 variants, spike mutations and immune escape. Nat. Rev. Microbiol. 2021, 19, 409–424. [Google Scholar] [CrossRef]
  29. Zhang, X.; Wu, S.; Wu, B.; Yang, Q.; Chen, A.; Li, Y.; Zhang, Y.; Pan, T.; Zhang, H.; He, X. SARS-CoV-2 Omicron strain exhibits potent capabilities for immune evasion and viral entrance. Signal Transduct. Target. Ther. 2021, 6, 1–3. [Google Scholar] [CrossRef]
  30. Koleya, T.; Kumara, M.; Goswami, A.; Ethayathulla, A.S.; Hariprasada, G. Structural modeling of Omicron spike protein and its complex with human ACE-2 receptor: Molecular basis for high transmissibility of the virus. Biochem. Biophys. Res. Commun. 2022, 592, 51–53. [Google Scholar] [CrossRef]
  31. Socher, E.; Heger, L.; Paulsen, F.; Zunke, F.; Arnold, P. Molecular dynamics simulations of the delta and omicron SARS-CoV-2 spike–ACE2 complexes reveal distinct changes between both variants. Comput. Struct. Biotechnol. J. 2022. ahead of print. [Google Scholar] [CrossRef] [PubMed]
  32. Kim, J.-M.; Chung, Y.-S.; Jo, H.J.; Lee, N.-J.; Kim, M.S.; Woo, S.H.; Park, S.; Kim, J.W.; Kim, H.M.; Han, M.-G. Identification of coronavirus isolated from a patient in Korea with COVID-19. Osong Public Health Res. Perspect. 2020, 11, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Park, A.K.; Kim, I.-H.; Kim, J.; Kim, J.-M.; Kim, H.M.; Young Lee, C.; Han, M.-G.; Rhie, G.-E.; Kwon, D.; Nam, J.-G. Genomic surveillance of SARS-CoV-2: Distribution of clades in the Republic of Korea in 2020. Osong Public Health Res. Perspect. 2021, 12, 37. [Google Scholar] [CrossRef] [PubMed]
  34. Socher, E.; Conrad, M.; Heger, L.; Paulsen, F.; Sticht, H.; Zunke, F.; Arnold, P. Mutations in the B. 1.1. 7 SARS-CoV-2 spike protein reduce receptor-binding affinity and induce a flexible link to the fusion peptide. Biomedicines 2021, 9, 525. [Google Scholar] [CrossRef]
  35. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1, 19–25. [Google Scholar] [CrossRef] [Green Version]
  36. Bjelkmar, P.; Larsson, P.; Cuendet, M.A.; Hess, B.; Lindahl, E. Implementation of the CHARMM force field in GROMACS: Analysis of protein stability effects from correction maps, virtual interaction sites, and water models. J. Chem. Theory Comput. 2010, 6, 459–466. [Google Scholar] [CrossRef]
  37. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N log (N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  38. Kumari, R.; Kumar, R.; Consortium, O.S.D.D.; Lynn, A. g_mmpbsa—A GROMACS tool for high-throughput MM-PBSA calculations. J. Chem. Inf. Modeling 2014, 54, 1951–1962. [Google Scholar] [CrossRef]
Figure 1. Mutation sites for S protein multiple variants. Trimeric RBDs in each chain are indicated in pink and ACE2 in red. Yellow, blue, and green colors indicate other domains except for RBD in each chain of the S protein. The black circular mark is the mutation site. (a) A double/triple mutant, (b) B.1.620, (c) alpha, (d) delta (E484Q mutation lacking in the Korean delta variant), and (e) mu variant.
Figure 1. Mutation sites for S protein multiple variants. Trimeric RBDs in each chain are indicated in pink and ACE2 in red. Yellow, blue, and green colors indicate other domains except for RBD in each chain of the S protein. The black circular mark is the mutation site. (a) A double/triple mutant, (b) B.1.620, (c) alpha, (d) delta (E484Q mutation lacking in the Korean delta variant), and (e) mu variant.
Ijms 23 04956 g001
Figure 2. The standard deviation for the distance between V503 residues in each chain for the final MD trajectory in the one-open-complex form (7A94). The X axis denotes variant types, and the Y axis denotes the standard deviation for distance (nm). Blue indicates the standard deviation (SD) between the AB and BC chains, orange indicates the SD between the AB and AC chains, gray indicates the SD between the AC and BC chains, and yellow indicates the SD among the A, B, and C chains.
Figure 2. The standard deviation for the distance between V503 residues in each chain for the final MD trajectory in the one-open-complex form (7A94). The X axis denotes variant types, and the Y axis denotes the standard deviation for distance (nm). Blue indicates the standard deviation (SD) between the AB and BC chains, orange indicates the SD between the AB and AC chains, gray indicates the SD between the AC and BC chains, and yellow indicates the SD among the A, B, and C chains.
Ijms 23 04956 g002
Figure 3. The standard deviation for the distance between V503 residues in each chain for the final MD trajectory in the two-open-complex form (7A97). The X axis denotes mutant types, and the Y axis denotes the standard deviation for distance (nm). Blue indicates the standard deviation between the AB and BC chains, orange indicates the SD between the AB and AC chains, gray indicates the SD between the AC and BC chains, and yellow indicates the SD among the A, B, and C chains.
Figure 3. The standard deviation for the distance between V503 residues in each chain for the final MD trajectory in the two-open-complex form (7A97). The X axis denotes mutant types, and the Y axis denotes the standard deviation for distance (nm). Blue indicates the standard deviation between the AB and BC chains, orange indicates the SD between the AB and AC chains, gray indicates the SD between the AC and BC chains, and yellow indicates the SD among the A, B, and C chains.
Ijms 23 04956 g003
Figure 4. The standard deviation for the distance between V503 residues in each chain for the final MD trajectory in the three-open-complex form (7A98). The X axis denotes mutant types, and the Y axis denotes the standard deviation for distance (nm). Blue indicates the standard deviation between the AB and BC chains, orange indicates the SD between the AB and AC chains, gray indicates the SD between the AC and BC chains, and yellow indicates the SD among the A, B, and C chains.
Figure 4. The standard deviation for the distance between V503 residues in each chain for the final MD trajectory in the three-open-complex form (7A98). The X axis denotes mutant types, and the Y axis denotes the standard deviation for distance (nm). Blue indicates the standard deviation between the AB and BC chains, orange indicates the SD between the AB and AC chains, gray indicates the SD between the AC and BC chains, and yellow indicates the SD among the A, B, and C chains.
Ijms 23 04956 g004
Figure 5. MM/PBSA results depending on the MD trajectories in the one-open-complex form (7A94). (a) Wild type and (bi) variants. (b) double, (c) triple, (d) B.1.620, (e) delta, (f) alpha, (g) delta_E484Q, (h) mu, and (i) omicron. The X axis denotes the trajectory (ps), and the Y axis denotes the MM/PBSA energy value (kJ/mol). Blue indicates the van der Waals energy between the S protein and ACE2, orange indicates electrostatic energy, and gray indicates the total energy.
Figure 5. MM/PBSA results depending on the MD trajectories in the one-open-complex form (7A94). (a) Wild type and (bi) variants. (b) double, (c) triple, (d) B.1.620, (e) delta, (f) alpha, (g) delta_E484Q, (h) mu, and (i) omicron. The X axis denotes the trajectory (ps), and the Y axis denotes the MM/PBSA energy value (kJ/mol). Blue indicates the van der Waals energy between the S protein and ACE2, orange indicates electrostatic energy, and gray indicates the total energy.
Ijms 23 04956 g005
Table 1. Summary of distance and SD between V503 residues (nm). The values were measured in the final trajectory (10 ns).
Table 1. Summary of distance and SD between V503 residues (nm). The values were measured in the final trajectory (10 ns).
PDBVariant TypeA–BA–CB–CS.D.
(A–B) & (A–C)
S.D.
(A–B) & (B–C)
S.D.
(A–C) & (B–C)
S.D.
(A–B–C)
7A94
(one-open-complex form)
Wild3.003.010.250.011.951.951.59
D614G/E484K2.782.510.240.191.791.611.39
D614G/E484K/N440K2.782.750.460.021.641.621.33
B.1.6202.552.720.780.121.251.251.07
Delta2.392.830.400.311.411.721.29
Alpha2.962.580.220.271.931.671.48
Delta (+E484Q)3.082.660.580.301.771.471.34
Mu2.952.390.400.401.811.411.34
Omicron3.333.280.910.041.721.681.38
7A97
(two-open-complex form)
Wild13.963.5710.637.352.365.005.31
D614G/E484K11.243.577.235.432.832.593.84
D614G/E484K/N440K9.523.166.164.492.372.123.18
B.1.6209.794.235.353.933.143.142.94
Delta11.133.597.455.332.602.733.77
Alpha10.654.236.094.553.231.323.31
Delta (+E484Q)10.593.926.664.722.781.943.35
Mu8.204.003.862.973.070.102.47
Omicron9.813.346.734.572.182.403.24
7A98
(three-open-complex form)
Wild8.2310.329.751.481.080.401.08
D614G/E484K7.729.819.061.480.950.531.06
D614G/E484K/N440K7.638.728.950.770.930.160.70
B.1.6208.718.029.590.490.620.620.79
Delta9.1810.009.880.580.500.080.44
Alpha8.809.1010.620.211.291.070.98
Delta (+E484Q)9.268.998.890.190.260.070.19
Mu9.609.609.140.000.330.330.27
Omicron10.328.819.891.060.310.760.78
Table 2. Summary of the MM/PBSA between the S protein and ACE2 in the one-open-complex form. VdW indicates van der Waals energy, E denotes electrostatic energy, and T is the total energy (kJ/mol). The values were measured in the final trajectory (10 ns).
Table 2. Summary of the MM/PBSA between the S protein and ACE2 in the one-open-complex form. VdW indicates van der Waals energy, E denotes electrostatic energy, and T is the total energy (kJ/mol). The values were measured in the final trajectory (10 ns).
Variant TypeVET
Wild−1461.5−10,311.8−11,773.3
Wild_deltascale (ΔF157/ΔR158)−1456.7−9227.7−10,684.4
Wild_ alphascale (ΔH69/ΔV70/ΔY144)−1457.6−10152.2−11609.8
D614G/E484K−1283.4−12,360.2−13,643.7
D614G/E484K_deltascale (ΔF157/ΔR158)−1280.2−11,380.0−12,660.2
D614G/E484K_alphascale (ΔH69/ΔV70/ΔY144)−1283.1−12,359.2−13,642.3
D614G/E484K/N440K−1116.2−14,581.0−15,697.2
D614G/E484K/N440K deltascale (ΔF157/ΔR158)−1114.5−13,526.4−14,640.9
D614G/E484K/N440K_alphascale (ΔH69/ΔV70/ΔY144)−1115.0−14,570.2−15,685.2
B.1.620 (S477N/E484K/D614G/P681H)−1318.0−13,832.5−15,150.5
B.1.620_deltascale (ΔF157/ΔR158)−1316.3−12,772.9−14,089.2
B.1.620_alphascale (ΔH69/ΔV70/ΔY144)−1317.7−13,821.9−15,139.6
Delta (ΔF157/ΔR158/L452R/T478K/D614G/P681R)−968.9−11,150.0−12,118.9
Alpha (ΔH69/ΔV70/ΔY144/N501Y/A570D/D614G/P681H)−1157.5−11,690.1−12,847.7
Delta_E484Q (VF157/ΔR158/L452R/T478K/D614G/P681R/E484Q)−1046.3−13,738.3−14,784.6
Mu (T95I/Y144S/Y145N/R346K/E484K/N501Y/D614G/P681H/D950N)−1290.4−14,776.8−16,067.2
Mu_deltascale (ΔF157/ΔR158)−1280.3−13,360.1−14,640.4
Mu_alphascale (ΔH69/ΔV70/ΔY144)−1290.0−14,782.2−16,072.2
Omicron−921.2−16,059.5−16,980.7
Table 3. Summary of the MM/PBSA between the S protein and ACE2 in the two-open-complex form. V indicates van der Waals energy, E denotes electrostatic energy, and T is the total energy (kJ/mol). The values were measured in the final trajectory (10 ns).
Table 3. Summary of the MM/PBSA between the S protein and ACE2 in the two-open-complex form. V indicates van der Waals energy, E denotes electrostatic energy, and T is the total energy (kJ/mol). The values were measured in the final trajectory (10 ns).
Variant Type7A97 AD7A97 BE7A97 Total
VETVETVET
Wild−1326.4−13,200.4−14,526.7−677.3−5553.0−6230.3−2003.7−18,753.4−20,757.0
Wild_deltascale (ΔF157/ΔR158)−1305.0−11,646.9−12,951.9−677.2−5195.1−5872.3−1982.3−16,842.0−18,824.3
Wild_alphascale (ΔH69/ΔV70/ΔY144)−1311.7−12,967.2−14,278.9−677.2−5550.9−6228.1−1989.0−18,518.1−20,507.1
D614G/E484K−1265.9−13,766.8−15,032.8−725.2−8535.5−9260.7−2003.2−20,445.4−22,448.6
D614G/E484K_deltascale (ΔF157/ΔR158)−1278.0−12,305.2−13,583.2−725.2−8140.3−8865.5−2003.2−20,445.6−22,448.8
D614G/E484K_alphascale (ΔH69/ΔV70/ΔY144)−1250.0−13,611.7−14,861.7−725.2−8538.6−9263.8−1975.3−22,150.4−24,125.6
D614G/E484K/N440K−874.3−13,307.9−14,182.2−626.7−10,773.4−11,400.1−1504.7−22,321.5−23,826.2
D614G/E484K/N440K_deltascale (ΔF157/ΔR158)−878.0−11,944.6−12,822.6−626.7−10,376.9−11,003.6−1504.7−22,321.5−23,826.2
D614G/E484K/N440K_alphascale (ΔH69/ΔV70/ΔY144)−873.4−13,308.3−14,181.8−626.7−10,774.4−11,401.1−1500.1−24,082.7−25,582.9
B.1.620 (S477N/E484K/D614G/P681H)−1696.4−14,874.1−16,570.6−676.0−10,506.1−11,182.1−2368.8−23,361.0−25,729.8
B.1.620_deltascale (ΔF157/ΔR158)−1692.8−13,415.8−15,108.6−676.0−9945.1−10,621.1−2368.8−23,360.9−25,729.7
B.1.620_alphascale (ΔH69/ΔV70/ΔY144)−1647.6−14,698.5−16,346.1−676.0−10,500.0−11,176.0−2323.6−25,198.5−27,522.1
Delta (ΔF157/ΔR158/L452R/T478K/D614G/P681R)−1344.8−13,867.5−15,212.3−534.5−8634.5−9169.0−1879.3−22,502.0−24,831.3
Alpha (ΔH69/ΔV70/ΔY144/N501Y/A570D/D614G/P681H)−1273.8−11,310.1−12,583.9−917.2−6604.0−7521.2−2191.0−17,914.1−20,105.1
Delta_E484Q (ΔF157/ΔR158/L452R/T478K/D614G/P681R/E484Q)−1385.1−13,658.2−15,043.3−676.6−10,319.8−10,996.4−2061.7−23,978.0−26,039.7
Mu (T95I/Y144S/Y145N/R346K/E484K/N501Y/D614G/P681H/D950N)−1084.8−13,097.6−14,182.3−1356.9−14,208.0−15,564.9−2441.7−27,305.6−29,747.2
Mu_deltascale (ΔF157/ΔR158)−1093.6−11,809.3−12,902.9−1356.4−13,230.0−14,586.4−2450.0−25,039.3−27,489.3
Mu_alphascale (ΔH69/ΔV70/ΔY144)−1074.1−13,035.9−14,110.0−1355.7−14,169.2−15,524.9−2429.8−27,205.1−29,634.9
Omicron−1343.9−16,994.6−18,338.5−1215.5−17,733.8−18,949.3−2559.5−34,728.4−37,287.8
Table 4. Summary of the MM/PBSA between the S protein and ACE2 in the three-open-complex form. V indicates van der Waals energy, E denotes electrostatic energy, and T is the total energy (kJ/mol). The values were measured in the final trajectory (10 ns).
Table 4. Summary of the MM/PBSA between the S protein and ACE2 in the three-open-complex form. V indicates van der Waals energy, E denotes electrostatic energy, and T is the total energy (kJ/mol). The values were measured in the final trajectory (10 ns).
Variant Type7A98 AD7A98 BE7A98 CF7A98 Total
VETVETVETVET
Wild−535.4−5278.5−5813.8−527.7−4872.3−5400.0−583.0−5587.4−6170.4−1646.1−15,738.2−17,384.2
Wild_deltascale(ΔF157/ΔR158)−535.4−4883.9−5419.3−527.7−4404.4−4932.1−583.0−5176.7−5759.7−1646.0−14,465.0−16,111.0
Wild_ alphascale (ΔH69/ΔV70/ΔY144)−535.4−5276.2−5811.6−527.7−4873.3−5401.0−583.0−5586.5−6169.5−1646.0−15,736.1−17,382.1
D614G/E484K−1236.3−12,616.4−13,852.7−931.4−13,292.4−14,223.9−710.6−10,229.1−10,939.7−2878.3−36,137.9−39,016.3
D614G/E484K_deltascale (ΔF157/ΔR158)−1236.0−11,740.8−12,976.9−931.1−12,398.8−13,329.9−710.3−9402.2−10,112.6−2877.4−33,541.8−36,419.4
D614G/E484K_ alphascale (ΔH69/ΔV70/ΔY144)−1209.9−12,601.0−13,810.9−909.4−13,278.9−14,188.3−717.8−10,124.4−10,842.2−2837.1−36,004.3−38,841.4
D614G/E484K/N440K−1203.1−12,512.8−13,716.0−1351.0−15,489.1−16,840.2−1010.8−12,006.6−13,017.4−3565.0−40,008.5−43,573.6
D614G/E484K/N440K_deltascale (ΔF157/ΔR158)−1202.8−11,658.0−12,860.9−1350.7−14,612.7−15,963.4−1009.5−11,050.2−12,059.8−3563.0−37,320.9−40,884.1
D614G/E484K/N440K_alphascale (ΔH69/ΔV70/ΔY144)−1171.4−12,430.8−13,602.2−1276.9−15,413.7−16,690.5−1000.6−11,966.5−12,967.1−3448.8−39,811.0−43,259.8
B.1.620 (S477N/E484K/D614G/P681H)−866.3−11,861.2−12,727.5−733.9−10,837.8−11,571.8−1275.6−13,000.2−14,275.8−2875.8−35,699.2−38,575.1
B.1.620_deltascale (ΔF157/ΔR158)−866.0−10,973.5−11,839.5−733.7−10,017.6−10,751.3−1275.4−12,173.6−13,449.0−2875.1−33,164.7−36,039.8
B.1.620_alphascale (ΔH69/ΔV70/ΔY144)−864.6−11,860.4−12,725.0−731.7−10,830.2−11,561.9−1237.9−12,914.8−14,152.7−2834.2−35,605.4−38,439.6
Delta (ΔF157/ΔR158/L452R/T478K/D614G/P681R)−724.0−8939.1−9663.1−753.4−9047.1−9800.4−687.1−8021.5−8708.7−2164.5−26,007.7−28,172.2
Alpha(ΔH69/ΔV70/ΔY144/N501Y/A570D/D614G/P681H)−939.2−9293.1−10,232.4−950.8−8220.7−9171.4−684.4−7538.3−8222.7−2574.4−25,052.1−27,626.5
Delta_E484Q(ΔF157/ΔR158/L452R/T478K/D614G/P681R/E484Q)−1017.0−12,096.5−13,113.4−1048.5−10,457.3−11,505.8−1022.2−12,888.7−13,910.9−3087.6−35,442.5−38,530.1
Mu(T95I/Y144S/Y145N/R346K/E484K/N501Y/D614G/P681H/D950N)−1168.3−12,901.2−14,069.5−632.7−11,080.4−11,713.1−1258.5−13,583.7−14,842.2−3059.5−37,565.3−40,624.8
Mu_deltascale (ΔF157/ΔR158)−1168.1−12,064.9−13,233.1−632.3−10,248.1−10,880.4−1258.0−12,631.4−13,889.4−3058.4−34,944.4−38,002.9
Mu_alphascale (ΔH69/ΔV70/ΔY144)−1147.2−12,842.0−13,989.3−630.7−11,062.9−11,693.6−1209.4−13,442.6−14,652.0−2987.3−37,347.5−40,334.9
Omicron−872.7−14,757.3−15,630.0−570.9−13,743.9−14,314.8−1081.3−14,247.0−15,328.3−2524.9−42,748.2−45,273.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Choi, K.-E.; Kim, J.-M.; Rhee, J.E.; Park, A.K.; Kim, E.-J.; Yoo, C.K.; Kang, N.S. Molecular Dynamics Studies on the Structural Stability Prediction of SARS-CoV-2 Variants Including Multiple Mutants. Int. J. Mol. Sci. 2022, 23, 4956. https://doi.org/10.3390/ijms23094956

AMA Style

Choi K-E, Kim J-M, Rhee JE, Park AK, Kim E-J, Yoo CK, Kang NS. Molecular Dynamics Studies on the Structural Stability Prediction of SARS-CoV-2 Variants Including Multiple Mutants. International Journal of Molecular Sciences. 2022; 23(9):4956. https://doi.org/10.3390/ijms23094956

Chicago/Turabian Style

Choi, Kwang-Eun, Jeong-Min Kim, Jee Eun Rhee, Ae Kyung Park, Eun-Jin Kim, Cheon Kwon Yoo, and Nam Sook Kang. 2022. "Molecular Dynamics Studies on the Structural Stability Prediction of SARS-CoV-2 Variants Including Multiple Mutants" International Journal of Molecular Sciences 23, no. 9: 4956. https://doi.org/10.3390/ijms23094956

APA Style

Choi, K. -E., Kim, J. -M., Rhee, J. E., Park, A. K., Kim, E. -J., Yoo, C. K., & Kang, N. S. (2022). Molecular Dynamics Studies on the Structural Stability Prediction of SARS-CoV-2 Variants Including Multiple Mutants. International Journal of Molecular Sciences, 23(9), 4956. https://doi.org/10.3390/ijms23094956

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop