Molecular Dynamics Simulations of Viral Glycoproteins: Predicting Host Receptor Binding and Immune Escape

Molecular dynamics (MD) simulations have become an indispensable tool in computational virology for dissecting the biophysical properties of viral glycoproteins at atomic resolution. These simulations enable the prediction of host receptor binding affinity, the identification of conformational epitopes, and the characterization of escape mutations that arise under immune pressure [1, 2]. By integrating structural biology data with all-atom force fields, researchers can model the dynamic behavior of glycoproteins from diverse viral families, including influenza A hemagglutinin (HA), paramyxovirus fusion (F) and attachment (G) proteins, coronavirus spike (S) proteins, and retrovirus envelope (Env) trimers [3, 4]. This article provides a technical review of the MD simulation workflow, its application to predicting receptor tropism and immune evasion, and its relevance to veterinary medicine and zoonotic risk assessment.

Computational Workflow for MD Simulations of Viral Glycoproteins

The typical MD simulation pipeline for viral glycoproteins begins with preparation of a three-dimensional structure, sourced either from experimental methods such as X-ray crystallography or cryo-electron microscopy, or from de novo prediction using tools such as AlphaFold2 [5, 6]. The structure is then placed in a solvated simulation box with explicit water molecules and ions, parameterized with a molecular mechanics force field such as CHARMM or AMBER [7, 8]. All-atom MD simulations are run for trajectories spanning hundreds of nanoseconds to several microseconds, using integration time steps of 1 to 2 femtoseconds under periodic boundary conditions [9]. Equilibration phases involve energy minimization, heating, and pressure coupling to mimic physiological conditions [4, 10].

Post-simulation analysis focuses on metrics that quantify conformational stability and binding energetics. Root mean square deviation (RMSD) and root mean square fluctuation (RMSF) are used to assess overall structural drift and per-residue flexibility, respectively [11, 12]. Principal component analysis (PCA) is applied to capture essential dynamics and identify dominant motions linked to receptor recognition or antibody evasion [13]. Binding free energies are estimated using the molecular mechanics generalized Born surface area (MM/GBSA) or the molecular mechanics Poisson Boltzmann surface area (MM/PBSA) approaches [14, 15]. Solvated interaction energy (SIE) calculations further prioritize ligand binding at protein interfaces [16].

Predicting Host Receptor Binding and Tropism

One of the most powerful applications of MD simulations is the prediction of host receptor binding specificity and tropism. For coronaviruses, the interaction between the receptor-binding domain (RBD) of the spike protein and the host angiotensin-converting enzyme 2 (ACE2) receptor is a critical determinant of cross-species transmission [17, 18, 19]. MD simulations have revealed how mutations in the RBD, such as N481K in SARS-CoV-2, alter hydrogen bonding networks and electrostatic complementarity at the interface [11]. Similarly, the D614G substitution in the spike protein reshapes allosteric networks and modulates the opening mechanism of the RBD, increasing receptor accessibility [20].

For influenza A viruses, MD simulations of HA in complex with sialic acid receptors distinguish avian-type alpha-2,3 linkages from mammalian-type alpha-2,6 linkages, providing a computational correlate of host range [7, 21]. The convergent evolution of the N156K mutation in pandemic H1N1 HA has been shown through MD to alter antigenic drift while maintaining receptor binding affinity [7]. In paramyxoviruses such as Nipah virus, simulations of the attachment glycoprotein (G) in complex with ephrin receptors reveal key residue contacts that mediate bat-to-human spillover [10, 22]. For the tick-borne encephalitis virus, MD models of the E protein suggest molecular mimicry with host proteins, potentially linking viral entry to autoimmune sequelae [6].

Non-enveloped viruses also benefit from MD insights. For white spot syndrome virus, simulations of envelope proteins VP28, VP26, and VP24 have identified multi-target inhibitors that block viral attachment in shrimp [23]. The Oropouche virus glycoprotein and RNA-dependent RNA polymerase have been targeted in immunoinformatic vaccine design workflows that rely on MD validation of epitope stability [3].

Identification of Immune Escape Mutations

MD simulations are uniquely suited to characterize the conformational basis of antibody escape. The HIV-1 envelope trimer exhibits high conformational variability, which facilitates immune evasion by cloaking conserved epitopes [1]. Long-timescale MD simulations of the HIV-1 Env trimer on the virion surface have shown that N-glycan moieties modulate the tilting of the glycoprotein spikes, thereby altering accessibility to broadly neutralizing antibodies [24]. For SARS-CoV-2, hierarchical mutational profiling combined with energy landscape analysis has identified distinct mechanisms of resistance against antibodies that target the RBD and the N-terminal domain [21]. The glycosylation state of the spike protein plays a central role in immunogenicity; MD simulations demonstrate that glycan shielding can sterically hinder antibody binding, while specific glycosylation sites can be engineered to enhance vaccine responses [19, 25].

MD simulations also inform the design of escape mutation prediction pipelines. For influenza A, deep mutational scanning data coupled with MD free energy calculations can rank HA mutations by their likelihood of escaping polyclonal sera [7]. In the context of dengue virus, consensus physics-based and artificial intelligence driven screening of the pre-fusion envelope protein has identified potent inhibitors that stabilize the closed conformation, thereby reducing exposure of fusion loop epitopes [26]. For Ebola virus, fragment-based drug design approaches guided by quantitative structure-activity relationships have identified small molecules that bind the glycoprotein and block entry [27]. Similarly, natural products have been computationally screened against the Ebola glycoprotein to identify entry inhibitors [8].

Implications for Veterinary Vaccine Design and Antiviral Development

The veterinary applications of MD simulations of viral glycoproteins are extensive. For Newcastle disease virus (NDV), computational and experimental studies have demonstrated that polyene macrolides can disrupt the F protein mediated membrane fusion, offering a basis for antiviral therapy in poultry [5]. The measles virus, which also infects non-human primates, has been studied with long-timescale MD to identify cannabichromevarin as a stabilizer of the prefusion F protein, a strategy that could be applied to paramyxovirus vaccine antigen design for livestock [4]. For the Middle East respiratory syndrome coronavirus (MERS-CoV), in silico exploration of the heptad repeat 2 domain of the spike fusion machinery has led to the design of antiviral peptides that inhibit membrane fusion [28].

MD simulations are also integrated with immunoinformatics to design multi-epitope vaccines. For infectious hematopoietic necrosis virus (IHNV) in fish, a multiple-epitope vaccine was designed using MD to confirm the structural stability of the epitopes displayed on a carrier protein [22]. For the monkeypox virus, a multi-epitope peptide vaccine targeting the A35R glycoprotein and E8L membrane protein was refined through MD to ensure proper folding and antigen presentation [29]. The same approach has been applied to the Oropouche virus, where the glycoprotein and RdRp were used as targets for a multi-epitope precision vaccine [3].

The role of pH in glycoprotein dynamics is also relevant for veterinary receptors. MD simulations have shown that microenvironmental pH influences the structural architecture of junctional adhesion molecules, which serve as receptors for reoviruses and other pathogens, suggesting that tissue specific pH changes can affect viral entry efficiency [18]. For the SARS-like bat coronavirus spike, structural dynamics and allosteric communication have been mapped using MD to predict which bat species maintain spikes with high human ACE2 binding affinity [30]. These models are critical for assessing zoonotic spillover risk from bat reservoirs [17].

Integrated Workflow Diagram

A typical integrated workflow for using MD simulations to predict receptor binding and immune escape is depicted in the following diagram:

flowchart TD
    A[Structural data: X-ray, Cryo-EM, AlphaFold2], > B[Structure preparation: solvation, ionization, force field assignment]
    B, > C[All-atom MD simulation: equilibration, production runs]
    C, > D[Trajectory analysis: RMSD, RMSF, PCA]
    D, > E[Binding free energy calculation: MM/GBSA, SIE]
    D, > F[Mutation scanning: in silico mutagenesis, free energy perturbation]
    E, > G[Rank receptor binding affinity predictions]
    F, > H[Identify escape mutations by antibody pressure]
    G, > I[Predict host tropism and spillover risk]
    H, > J[Design stable vaccine antigens or escape resistant therapeutics]
    I, > K[Experimental validation: surface plasmon resonance, pseudovirus entry assays]
    J, > K

Key Studies Using MD Simulations of Viral Glycoproteins

Future Directions

The continued development of coarse-grained MD models will enable simulations of larger glycoprotein complexes, such as entire virion surfaces, over physiologically relevant timescales [9, 12]. Integration of MD with deep mutational scanning and machine learning will further automate the prediction of escape mutations under defined antibody pressure [21, 16]. For veterinary virology, these approaches can be directly applied to pathogens of livestock, poultry, and wildlife to inform surveillance and vaccine antigen design [5, 22]. The combination of MD with structural prediction methods like AlphaFold2 will allow rapid modeling of glycoproteins from newly emerged viruses without experimental structures [6, 30]. As computational power increases, routine MD screening of variant effects on receptor binding and immune escape will become a standard component of outbreak preparedness programs [13, 18].

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.

References

[1] Cao Y, Im W. Conformational variability of HIV-1 Env trimer and viral vulnerability. Elife. URL: https://pubmed.ncbi.nlm.nih.gov/42360802/

[2] Almaghrabi M, Alturki MS. AI-Driven Discovery of Prototype CLEC4M Inhibitors Targeting Marburg Virus Entry via Integrated Machine Learning and Molecular Modeling. Int J Mol Sci. URL: https://pubmed.ncbi.nlm.nih.gov/42353044/

[3] Alissa M, Suleman M. Immunoinformatic based Development of a Multi-Epitope Precision Vaccine Targeting Glycoprotein and RdRp of Oropouche Virus: An Innovative Approach to Counter Emerging Public Health Threats. Saudi Med J. URL: https://pubmed.ncbi.nlm.nih.gov/42293716/

[4] Abbou H, Gaouzi Z, Zegrari R et al. Identification of cannabichromevarin as a potent stabilizer of the measles virus prefusion F protein: structural insights from long-timescale molecular dynamics. Sci Rep. URL: https://pubmed.ncbi.nlm.nih.gov/42286054/

[5] Mukhametkaliyev A, Bogoyavlenskiy A, Alexyuk P et al. Antiviral Activity of Polyene Macrolides Against Newcastle Disease Virus: Computational and Experimental Insights. Molecules. URL: https://pubmed.ncbi.nlm.nih.gov/42280215/

[6] Timofeeva AM, Aulova KS, Ulyanova YS et al. Molecular Mimicry by the Tick-Borne Encephalitis Virus E Protein: A Hidden Link to Autoimmunity. Int J Mol Sci. URL: https://pubmed.ncbi.nlm.nih.gov/42278275/

[7] Tian Y, Jiang D, Ma W et al. Convergent evolution of the N156K mutation in A(H1N1)pdm09 hemagglutinin contributes to antigenic drift and cluster transition. Emerg Microbes Infect. URL: https://pubmed.ncbi.nlm.nih.gov/42267384/

[8] Ait Lahcen N, Yang L, Chen W et al. Natural product-based Ebola virus entry inhibitors targeting the viral glycoprotein: A combined computational and experimental study. Antiviral Res. URL: https://pubmed.ncbi.nlm.nih.gov/42263866/

[9] Urano R, Yoshimoto S, Hori K et al. Synergistic Protein-Protein and Protein-Lipid Interactions Drive SARS-CoV-2 Envelope Assembly. J Chem Inf Model. URL: https://pubmed.ncbi.nlm.nih.gov/42262910/

[10] de Oliveira VM, Marinho MM, Dos Santos HS et al. Molecular Docking and Dynamics Simulations of Natural Limonoids Interacting with the Nipah Virus Attachment Glycoprotein. Curr Drug Targets. URL: https://pubmed.ncbi.nlm.nih.gov/42261155/

[11] Smatti MK, Al-Khatib H, Suleman M et al. Functional and structural characterization of the SARS-CoV-2 spike N481K mutation. Arch Virol. URL: https://pubmed.ncbi.nlm.nih.gov/42243594/

[12] Majumder A, Dutta M, Cherek L et al. Structure and Dynamics of the HIV-1 Envelope Protein on the Virion Envelope. bioRxiv. URL: https://pubmed.ncbi.nlm.nih.gov/42239368/

[13] Altaf M, Kin Cheng MT, González-Vázquez LD et al. SARS-CoV-2 intra-host recombination promotes epistatic spike interactions and temperature-dependent adaptation. Cell Rep. URL: https://pubmed.ncbi.nlm.nih.gov/42213788/

[14] Navhaya LT, Monama MZ, Matsebatlela TM et al. In-Depth Molecular Dynamics Simulations Reveal Ligand-Induced Modulations of the HSPA8-SARS-CoV-2 Spike Protein Interaction. Int J Mol Sci. URL: https://pubmed.ncbi.nlm.nih.gov/42196270/

[15] Delgado-Maldonado T, Vargas-Salas G, Bandyopadhyay D et al. Ligand-Based Pharmacophore Modeling and Structure-Based Virtual Screening for Identifying Potential Therapeutic Agents Targeting the SARS-CoV-2 Spike Protein-ACE2 Receptor Interaction. ChemMedChem. URL: https://pubmed.ncbi.nlm.nih.gov/42174382/

[16] Kongtaworn N, Sinsulpsiri S, Hanpaibool C et al. Molecular Dynamics and Solvated Interaction Energy Prioritize Cannabidiol and Cannabinol as Variant-Spanning SARS-CoV-2 RBD-ACE2 Interface Blockers. Molecules. URL: https://pubmed.ncbi.nlm.nih.gov/42075932/

[17] Kori L, Chandra A, Mukhopadhyay L et al. E-protein variability in Zika virus strains: A possible new O-glycosylation site and its implications. Indian J Med Res. URL: https://pubmed.ncbi.nlm.nih.gov/42165713/

[18] Karagöl T, Karagöl A. Microenvironmental pH and host-pathogen co-evolution potentially influence the structural architecture of junctional adhesion molecules. Mol Genet Genomics. URL: https://pubmed.ncbi.nlm.nih.gov/42154072/

[19] Wang H, Li X, He P et al. Investigating the effects and underlying mechanisms of glycosylation sites on the immunogenicity of COVID-19 vaccines. Signal Transduct Target Ther. URL: https://pubmed.ncbi.nlm.nih.gov/42129141/

[20] Kearns FL, Bogetti AT, Calvó-Tusell C et al. D614G reshapes allosteric networks and opening mechanisms of SARS-CoV-2 spikes. Proc Natl Acad Sci U S A. URL: https://pubmed.ncbi.nlm.nih.gov/42101997/

[21] 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. URL: https://pubmed.ncbi.nlm.nih.gov/42123600/

[22] Abdellaoui N, Lim JW, Han HJ et al. Immunoinformatics-Based Design of a Multiple-Epitope Vaccine Against IHNV. Mar Biotechnol (NY). URL: https://pubmed.ncbi.nlm.nih.gov/42113364/

[23] Hossen Z, Sajjad Hossain M, Naim Uddin Forhad M et al. Multi-target inhibitors of white spot syndrome virus envelope proteins (VP28, VP26, and VP24) to protect shrimp (Penaeus monodon): An integrated virtual screening, pharmacokinetic, and molecular dynamics study. Microb Pathog. URL: https://pubmed.ncbi.nlm.nih.gov/42055173/

[24] Shehata M, Casalino L, Duquette M et al. N-Glycans modulate tilting of HIV-1 envelope glycoprotein. Nat Commun. URL: https://pubmed.ncbi.nlm.nih.gov/41986331/

[25] Zhou X, Du H, Wang C et al. Identification of a Potential Dual-Target Candidate Against RSV F Protein and 15-LOX from TCMSP: Integrating Virtual Screening, Molecular Dynamics, and Experimental Evaluation. Int J Mol Sci. URL: https://pubmed.ncbi.nlm.nih.gov/42074092/

[26] Pravin MA, Khan MA, Singh SK. Identification of potent inhibitors targeting the pre-fusion DENV envelope protein: A consensus physics-and AI-based in silico multi-tier screening approach. Biochem Biophys Res Commun. URL: https://pubmed.ncbi.nlm.nih.gov/42008872/

[27] Ait Lahcen N, Liman W, Zekri S et al. QSAR-Guided and Fragment-Based Drug Design of Monoterpenoid Inhibitors Targeting Ebola Virus Glycoprotein. Int J Mol Sci. URL: https://pubmed.ncbi.nlm.nih.gov/41977174/

[28] Alotaiq N, Dermawan D, Chtita S. Targeting Middle East Respiratory Syndrome Coronavirus Spike Fusion Machinery With Antiviral Peptides: In Silico Exploration of the Heptad Repeat 2 Domain. Microbiologyopen. URL: https://pubmed.ncbi.nlm.nih.gov/42082899/

[29] Attique L, Jamal SB, Gulistan T et al. Multiepitope-Based Peptide Vaccine Against A35R Glycoprotein and E8L Membrane Protein of Monkeypox Virus Using an Immunoinformatics Approach. Biology (Basel). URL: https://pubmed.ncbi.nlm.nih.gov/41972525/

[30] Balogun TA, Kearns FL, Calvó-Tusell C et al. Structural dynamics and allosteric communication of a SARS-like bat coronavirus spike glycoprotein. Biophys J. URL: https://pubmed.ncbi.nlm.nih.gov/42026866/

[31] Boršová V, Zajaček D, Bucinsky L et al. Integrated computational protocol for sampling molecular databases towards allosteric inhibition of SARS-CoV-2 spike activation. J Biomol Struct Dyn. URL: https://pubmed.ncbi.nlm.nih.gov/42101152/

[32] Tariq SS, Zia K, Ali N et al. Immunoinformatics driven multi-epitope vaccine design targeting clade IIb Mpox variant. Sci Rep. URL: https://pubmed.ncbi.nlm.nih.gov/42092088/

[33] Bogoyavlenskiy A, Manakbayeva A, Kerimov T et al. Structural Basis of Anthocyanin-Mediated Modulation of IL-2, IL-17, and TNF-α: A Docking and Molecular Dynamics Study. Int J Mol Sci. URL: https://pubmed.ncbi.nlm.nih.gov/42074123/

[34] Dubey A, Kumar M, Tufail A. Suppressing viral assembly in human metapneumovirus by targeting fusion protein with natural compounds: a structural dynamics and energetics study. J Mol Model. URL: https://pubmed.ncbi.nlm.nih.gov/42043606/

[35] Ušjak L, Skalicka-Woźniak K, Kulinowski Ł et al. Integrated Phytochemical Profiling and Bioactivity Evaluation of Micromeria nervosa, with Emphasis on Antimicrobial and Antiviral Properties. Antibiotics (Basel). URL: https://pubmed.ncbi.nlm.nih.gov/42041337/