Structural and Evolutionary Dynamics of Coronavirus Spike Protein: Integrating Cryo-EM, Molecular Dynamics, and Phylogenetic Surveillance

Introduction

Coronavirus spike (S) glycoproteins mediate host cell attachment and membrane fusion, constituting the primary determinant of host range, tissue tropism, and immune recognition across diverse coronaviruses infecting mammalian and avian species [1, 2]. In veterinary medicine, spike-driven entry mechanisms are central to the pathogenesis of porcine epidemic diarrhea virus (PEDV), feline infectious peritonitis virus (FIPV), canine respiratory coronavirus, bovine coronavirus, and avian infectious bronchitis virus (IBV), among others [3]. The structural biology of the spike protein, particularly its receptor-binding domain (RBD) and S1/S2 cleavage site, has been extensively characterized in SARS-CoV-2 as a model system, providing a framework for understanding analogous processes in animal coronaviruses [4, 5]. This article reviews the integration of cryo-electron microscopy (cryo-EM), molecular dynamics (MD) simulations, and phylogenetic surveillance in elucidating the structural and evolutionary dynamics of coronavirus spike proteins, with emphasis on computational methodologies applicable to veterinary virology.

For a detailed discussion of cryo-EM reconstruction workflows, see the companion article Cryo-EM Density Map Interpretation and Computational Structure Fitting. The evolutionary dynamics of RNA viruses are further explored in Evolutionary Dynamics of RNA Viruses.

Cryo-Electron Microscopy of Spike Protein Architecture

Cryo-EM has resolved the prefusion trimeric structure of coronavirus spikes at near-atomic resolution, revealing distinct domains: the N-terminal domain (NTD), the RBD in the S1 subunit, and the S2 fusion machinery including the fusion peptide and heptad repeats [6, 7]. The RBD undergoes hinge-like conformational changes between a "closed" (receptor-inaccessible) and "open" (receptor-accessible) state, a feature conserved across coronaviruses [1, 8]. Comparative structural analysis of the NTD across wild-type and emerging variants highlights the role of this domain in immune evasion through glycan shielding and loop insertions [6]. Glycosylation at conserved sites such as N343 in the RBD modulates co-receptor binding and antibody recognition across variants of concern [7].

The S1/S2 cleavage site, typically containing a polybasic motif in highly pathogenic coronaviruses, is a critical hotspot for furin-like protease processing, which primes the spike for membrane fusion after receptor engagement [5]. Structural constraints at this interface limit the space for viral adaptation, as demonstrated by evolutionary coupling analyses [1]. These static snapshots, while informative, do not capture the conformational ensembles that govern spike dynamics under physiological conditions.

Molecular Dynamics Simulations of Conformational Flexibility

All-atom and coarse-grained MD simulations complement cryo-EM by providing time-resolved trajectories of spike conformational changes at atomic resolution [9, 10, 11]. Simulations of the RBD-ACE2 complex have quantified the binding free energy contributions of individual mutations, revealing that affinity enhancement is often coupled with immune escape [9, 2]. For example, integrative MD analysis of mutation-driven adaptation in the RBD shows that substitutions such as N501Y increase hydrophobic packing with ACE2, while E484K alters electrostatic complementarity [12, 9, 8].

Extended simulations of the full-length spike ectodomain have elucidated the allosteric coupling between RBD opening and S2 domain pre-stabilization [13]. Markov state models constructed from microsecond-scale trajectories identify metastable states along the RBD opening pathway and hidden allosteric pockets that can be targeted by small molecules or antibodies [13]. The interplay of dynamics and convergent evolution modulates allostery, as observed in Omicron sublineages where compensatory mutations restore fitness losses [14, 10].

For protocols on setting up and analyzing protein-water systems, refer to GROMACS Molecular Dynamics: Setting Up, Simulating, and Analyzing Protein-Water Systems. A broader overview of force fields is provided in Molecular Dynamics Simulations of Proteins and Force Fields.

Integrative Modeling: Cryo-EM Restraints and Ensemble Refinement

Integrative approaches combine cryo-EM density maps as spatial restraints in MD simulations, producing ensembles that satisfy both experimental data and physical force fields [15, 16]. This hybrid strategy resolves heterogeneous conformations within a single cryo-EM dataset, revealing distinct dynamic signatures for antibody-bound versus unbound spikes [15]. For broadly neutralizing antibodies targeting conserved epitopes, the energy landscape of binding is characterized by frustrated interfaces that drive adaptive evolution [17, 18].

Multiscale modeling further bridges atomic detail with coarse-grained representations of the lipid membrane environment, enabling simulations of membrane fusion intermediates [19, 20]. These calculations compute viral fitness as a function of binding affinity and escape potential, providing a quantitative framework for predicting variant emergence [19, 21]. The integration of mutational profiling with ensemble-based network analysis identifies epistatic couplings that control ACE2 affinity across XBB lineages [11].

Phylogenetic Surveillance of Spike Evolution

Phylogenetic analysis of spike gene sequences from global surveillance databases (e.g., GISAID) tracks the emergence and spread of variants at regional and global scales [22, 23, 24]. Large-scale sequencing of PEDV S genes in China revealed extensive recombination and antigenic diversity, underscoring the need for continuous monitoring in swine populations [3]. For SARS-CoV-2, genomic epidemiology workflows have reconstructed transmission clusters and mutation dynamics during successive waves [24, 25]. In Morocco, whole-genome sequencing from 2020 to 2024 documented the replacement of early lineages by successive Omicron sublineages, driven by spike mutations conferring immune evasion [23].

The micro-evolution of spike within individual hosts, particularly in immunocompromised individuals, can generate novel mutations that seed new variants [26]. Defective viral genomes under selection pressure further contribute to spike diversity [27]. Bayesian walker algorithms coupled with computational workflows predict likely future mutations by extrapolating current evolutionary trajectories [28].

For an overview of global surveillance platforms, see The World Health Organization (WHO) and Global Genomic Surveillance. The specific article on Computational Modeling of Viral Quasispecies Diversity and Evolutionary Fitness Landscapes addresses intra-host diversity.

Binding Free Energy Calculations and Cross-Species Risk

Free energy perturbation and molecular mechanics generalized Born surface area (MM/GBSA) methods quantify the impact of individual mutations on ACE2 binding affinity and antibody escape [17, 21, 29]. These calculations have been applied to rank the fitness of emerging variants and predict zoonotic potential [30, 12]. The conservation of bacterial lipopolysaccharide binding by the spike across major variants suggests an additional layer of host interaction that may influence pathogenesis [30].

Computational alanine scanning and deep mutational scanning datasets are used to train models that predict which residues are essential for binding versus immune evasion [4, 29]. For example, the epitope landscape of the RBD can be classified into distinct escape categories, enabling proactive vaccine design [4]. Short functional peptides designed to inhibit ACE2-spike interaction, such as boomerang-shaped peptides, have been explored as therapeutic leads [31].

The Structural Bioinformatics of Viral Envelope Proteins and Entry Mechanisms article provides context for receptor engagement. Cross-species risk assessment is further discussed in Deep Learning for Predicting Viral Host-Range Transitions and Zoonotic Potential.

Allostery and Epistasis in Spike Evolution

Beyond additive mutation effects, epistatic interactions between spatially distant residues shape spike evolvability [32, 14, 33]. Allosteric communication networks within the spike trimer, identified through normal mode analysis and residue interaction networks, reveal that mutations in the NTD can exert long-range effects on RBD conformation and antibody binding [16, 33]. The integration of candidate adaptive polymorphisms with protein dynamics shows that evolutionary adaptations often exploit pre-existing dynamic modes [33].

Class I and Class 4/1 neutralizing antibodies overcome steric limitations through allosteric mechanisms, an insight derived from multimodal computational approaches combining docking, MD, and network analysis [15, 16, 20]. The balance of evolutionary adaptability and dynamic constraints determines the molecular determinants of immune escape [18]. These findings are directly relevant to veterinary vaccine design, where conserved spike epitopes are targeted for broad protection against diverse coronaviruses such as those affecting swine and poultry.

Machine Learning for Variant Effect Prediction

Machine learning models, including deep neural networks trained on large-scale mutagenesis data, predict the impact of unseen spike mutations on binding affinity, stability, and antibody escape [4, 28]. AlphaFold2-based predictions of conformational ensembles for Omicron sublineages have been combined with MD simulations to capture epistatic couplings that control ACE2 affinity [14, 10]. These hybrid approaches are increasingly used to screen emerging variants in near-real time during surveillance [22, 34].

The biological foundation models discussed in Biological Foundation Models in Veterinary Virology: From ESMFold to Genomic Surveillance represent the frontier of sequence-to-structure prediction. The Machine Learning for Variant Effect Prediction on Protein Stability article details methods directly applicable to spike mutagenesis.

Integrative Workflow Diagram

The following Mermaid diagram summarizes the integrated computational pipeline from structural determination to evolutionary prediction.

flowchart TD
    A[SARS-CoV-2 Spike Sequences\n(GISAID/Nextstrain)], > B[Phylogenetic Reconstruction\nand Variant Clustering]
    A, > C[AlphaFold2 Homology Modeling\nof Spike Trimer]
    C, > D[Cryo-EM 3D Reconstruction\n(Closed/Open States)]
    D, > E[All-Atom MD Simulations\n(GROMACS/AMBER)]
    E, > F[Markov State Models\nof Conformational Ensembles]
    E, > G[Binding Free Energy Calculations\n(MM/GBSA, FEP)]
    G, > H[Mutational Profiling\nand Epistatic Network Analysis]
    H, > I[Prediction of Binding Affinity\nand Immune Escape Hotspots]
    B, > J[Bayesian Evolutionary Analysis\nand Mutation Rate Estimation]
    J, > K[Forecasting of Emerging Variants]
    K, > I
    I, > L[Vaccine Antigen Design\nand Cross-Species Risk Assessment]

Figure 1. Integrated computational workflow combining cryo-EM, molecular dynamics, free energy calculations, and phylogenetic surveillance to study coronavirus spike protein structural and evolutionary dynamics.

Table of Computational Methods and Applications

Conclusions

The integration of cryo-EM, molecular dynamics, and phylogenetic surveillance provides a powerful framework for understanding coronavirus spike protein structure, dynamics, and evolution. While most detailed studies have focused on SARS-CoV-2 as a model, the computational methodologies apply directly to veterinary coronaviruses such as PEDV, IBV, and feline coronavirus. Structural constraints and allosteric communication networks limit the evolutionary space of spike, yet immune pressure drives continued diversification. Free energy calculations and machine learning models now enable near-real-time risk assessment of emerging variants. For the veterinary field, adopting these integrated approaches will enhance surveillance, inform vaccine design, and improve preparedness for cross-species transmission events involving coronaviruses of livestock and companion animals.

References

[1] Herzig JC, Magwira ML, Lovell SC. Structural Constraints Acting on the SARS-CoV-2 Spike Protein Reveal Limited Space for Viral Adaptation. Genome Biol Evol. 2026. URL: https://pubmed.ncbi.nlm.nih.gov/41876430/

[2] Ataya F, Alamro A, Alghamdi A, et al. SARS-CoV-2 spike mutations alter structure and energetics to modulate ACE2 binding immune evasion and viral adaptation. Sci Rep. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/41152301/

[3] Fu Y, Wang Y, Dai L, et al. Evolutionary dynamics and antigenic diversity of porcine epidemic diarrhea virus (PEDV) in China: phylogenetic and recombination analyses based on large-scale S gene sequences. BMC Vet Res. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/40604964/

[4] Teruel NFB, Crown M, Rajsbaum R, et al. Comprehensive analysis of SARS-CoV-2 Spike evolution: epitope classification and immune escape prediction. Virus Evol. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/40510757/

[5] Sinha A, Sangeet S, Roy S. Evolution of Sequence and Structure of SARS-CoV-2 Spike Protein: A Dynamic Perspective. ACS Omega. 2023. URL: https://pubmed.ncbi.nlm.nih.gov/37426203/ *** Disclaimer: This article is for educational and informational purposes only. It is not intended to substitute for professional veterinary advice, diagnosis, treatment, or regulatory guidance. Always consult a licensed veterinarian or qualified specialist regarding animal health, disease diagnosis, and therapeutic decisions.

[6] Quaranta M, Via A, Pascarella S. Structural Analysis of the SARS-CoV-2 Spike N-Terminal Domain Across Wild-Type and Recent Variants: A Comparative Study. Proteins. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/40485545/

[7] Ives CM, Nguyen L, Fogarty CA, et al. Role of N343 glycosylation on the SARS-CoV-2 S RBD structure and co-receptor binding across variants of concern. Elife. 2024. URL: https://pubmed.ncbi.nlm.nih.gov/38864493/

[8] Hsiao YW, Bray DJ, Taddese T, et al. Structure adaptation in Omicron SARS-CoV-2/hACE2: Biophysical origins of evolutionary driving forces. Biophys J. 2023. URL: https://pubmed.ncbi.nlm.nih.gov/37717145/

[9] Truong Hoai L, Ilham B, Sompornpisut T, et al. Mutation-driven adaptation of ACE2-RBD binding revealed by integrative molecular dynamics analysis. J Mol Graph Model. 2026. URL: https://pubmed.ncbi.nlm.nih.gov/41839468/

[10] Raisinghani N, Alshahrani M, Gupta G, et al. AlphaFold2 Predictions of Conformational Ensembles and Atomistic Simulations of the SARS-CoV-2 Spike XBB Lineages Reveal Epistatic Couplings between Convergent Mutational Hotspots that Control ACE2 Affinity. J Phys Chem B. 2024. URL: https://pubmed.ncbi.nlm.nih.gov/38696745/

[11] Raisinghani N, Alshahrani M, Gupta G, et al. Ensemble-Based Mutational Profiling and Network Analysis of the SARS-CoV-2 Spike Omicron XBB Lineages for Interactions with the ACE2 Receptor and Antibodies: Cooperation of Binding Hotspots in Mediating Epistatic Couplings Underlies Binding Mechanism and Immune Escape. Int J Mol Sci. 2024. URL: https://pubmed.ncbi.nlm.nih.gov/38673865/

[12] Soliman OA, Shahine Y, Baecker D, et al. Beyond the Mutation Abyss: Revisiting SARS-CoV-2 Receptor-Binding Domain Evolution from ACE2 Binding Optimization to Immune Epitope Remodeling. Pathogens. 2026. URL: https://pubmed.ncbi.nlm.nih.gov/41901725/

[13] Xiao S, Alshahrani M, Gupta G, et al. Markov State Models and Perturbation-Based Approaches Reveal Distinct Dynamic Signatures and Hidden Allosteric Pockets in the Emerging SARS-Cov-2 Spike Omicron Variant Complexes with the Host Receptor: The Interplay of Dynamics and Convergent Evolution Modulates Allostery and Functional Mechanisms. J Chem Inf Model. 2023. URL: https://pubmed.ncbi.nlm.nih.gov/37549201/

[14] Raisinghani N, Alshahrani M, Gupta G, et al. AlphaFold2 Modeling and Molecular Dynamics Simulations of the Conformational Ensembles for the SARS-CoV-2 Spike Omicron JN.1, KP.2 and KP.3 Variants: Mutational Profiling of Binding Energetics Reveals Epistatic Drivers of the ACE2 Affinity and Escape Hotspots of Antibody Resistance. Viruses. 2024. URL: https://pubmed.ncbi.nlm.nih.gov/39339934/

[15] Alshahrani M, Parikh V, Foley B, et al. Dissecting binding and immune evasion mechanisms for ultrapotent Class I and Class 4/1 neutralizing antibodies of SARS-CoV-2 spike protein using a multi-pronged computational approach: neutral frustration architecture of binding interfaces and immune escape hotspots drives adaptive evolution. Phys Chem Chem Phys. 2026. URL: https://pubmed.ncbi.nlm.nih.gov/41623222/

[16] Alshahrani M, Parikh V, Foley B, et al. Allosteric Control Overcomes Steric Limitations for Neutralizing Antibodies Targeting Conserved Binding Epitopes of the SARS-CoV-2 Spike Protein: Exploring the Intersection of Binding, Allostery, and Immune Escape with a Multimodal Computational Approach. Biomolecules. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/41008647/

[17] Alshahrani M, Gatlin W, Ludwick M, et al. Mechanisms of Binding and Immune Escape Resistance for Broadly Neutralizing Antibodies Targeting Distinct Conserved SARS-CoV-2 Spike Epitopes: A Hierarchical Approach Integrating Mutational Profiling and Energy Landscape Analysis. Int J Mol Sci. 2026. URL: https://pubmed.ncbi.nlm.nih.gov/42123600/

[18] Alshahrani M, Parikh V, Foley B, et al. Integrative Computational Modeling of Distinct Binding Mechanisms for Broadly Neutralizing Antibodies Targeting SARS-CoV-2 Spike Omicron Variants: Balance of Evolutionary and Dynamic Adaptability in Shaping Molecular Determinants of Immune Escape. Viruses. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/40573332/

[19] Mehra R, Thakur S. The structure-based approaches to computing viral fitness. Adv Protein Chem Struct Biol. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/40973410/

[20] Alshahrani M, Parikh V, Foley B, et al. Multiscale Modeling and Dynamic Mutational Profiling of Binding Energetics and Immune Escape for Class I Antibodies with SARS-CoV-2 Spike Protein: Dissecting Mechanisms of High Resistance to Viral Escape Against Emerging Variants. Viruses. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/40872744/

[21] Alshahrani M, Parikh V, Foley B, et al. Mutational Scanning and Binding Free Energy Computations of the SARS-CoV-2 Spike Complexes with Distinct Groups of Neutralizing Antibodies: Energetic Drivers of Convergent Evolution of Binding Affinity and Immune Escape Hotspots. Int J Mol Sci. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/40003970/

[22] Yan J, Liu F, Hu S, et al. Regional dynamics and mechanisms behind SARS-CoV-2 XDV.1 prevalence in Chongqing via genomic surveillance and molecular insights. Virus Res. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/40139569/

[23] Ghammaz H, Melloul M, Mbarki A, et al. Genomic evolution of SARS-CoV-2 in Morocco: Insights from whole genome sequences collected from 2020 to 2024. Virus Res. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/39864629/

[24] Razzaq A, Disoma C, Iqbal S, et al. Genomic epidemiology and evolutionary dynamics of the Omicron variant of SARS-CoV-2 during the fifth wave of COVID-19 in Pakistan. Front Cell Infect Microbiol. 2024. URL: https://pubmed.ncbi.nlm.nih.gov/39502171/

[25] Bi D, Luo X, Chen Z, et al. Genomic epidemiology reveals early transmission of SARS-CoV-2 and mutational dynamics in Nanning, China. Heliyon. 2023. URL: https://pubmed.ncbi.nlm.nih.gov/38125422/

[26] Huang W, Yin C, Briley KP, et al. Dynamic Evolution of SARS-CoV-2 in a Patient on Chemotherapy. Viruses. 2023. URL: https://pubmed.ncbi.nlm.nih.gov/37632101/

[27] Lin CH, Tam HM, Yang CY, et al. Evolution of the coronavirus spike protein in the full-length genome and defective viral genome under diverse selection pressures. J Gen Virol. 2023. URL: https://pubmed.ncbi.nlm.nih.gov/37997889/

[28] Ben Geoffrey AS, Gracia J. A Bayesian walker coupled with a computational workflow that generates the micro-evolution of SARS-CoV-2 and makes predictions of new mutations that can emerge. J Biomol Struct Dyn. 2024. URL: https://pubmed.ncbi.nlm.nih.gov/37771150/

[29] Alshahrani M, Parikh V, Foley B, et al. Quantitative Characterization and Prediction of the Binding Determinants and Immune Escape Hotspots for Groups of Broadly Neutralizing Antibodies Against Omicron Variants: Atomistic Modeling of the SARS-CoV-2 Spike Complexes with Antibodies. Biomolecules. 2025. URL: https://pubmed.ncbi.nlm.nih.gov/40001552/

[30] Samsudin F, Petruk G, Rui L, et al. Conservation of Bacterial Lipopolysaccharide Binding by SARS-CoV-2 Spike across Major Viral Variants. Comput Struct Biotechnol J. 2026. URL: https://pubmed.ncbi.nlm.nih.gov/41993878/

[31] Wei Y, Liu Z, Zhang M, et al. Inhibition of ACE2-S Protein Interaction by a Short Functional Peptide with a Boomerang Structure. Molecules. 2024. URL: https://pubmed.ncbi.nlm.nih.gov/38998974/

[32] Uzer F, Erendor F, Sanlioglu S. SARS-CoV-2 Variants and Immune Evasion: Mapping the Future of Vaccine Design. Rev Med Virol. 2026. URL: https://pubmed.ncbi.nlm.nih.gov/42047168/

[33] Ose NJ, Campitelli P, Modi T, et al. Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics. Elife. 2024. URL: https://pubmed.ncbi.nlm.nih.gov/38713502/

[34] Tsujino S, Deguchi S, Nomai T, et al. Virological characteristics of the SARS-CoV-2 Omicron EG.5.1 variant. Microbiol Immunol. 2024. URL: https://pubmed.ncbi.nlm.nih.gov/38961765/

[35] Chaudhuri A, Das S, Chakrabarti S. Mutational and evolutionary dynamics of non-structural and spike proteins from variants of concern (VOC) of SARS-CoV-2 in India. Int J Biol Macromol. 2024. URL: https://pubmed.ncbi.nlm.nih.gov/39488303/