Abstract:Molecular dynamics (MD) simulations produce a substantial volume of high-dimensional data, and traditional methods for analyzing these data pose significant computational demands. Advances in MD simulation analysis combined with deep learning-based approaches have led to the understanding of specific structural changes observed in MD trajectories, including those induced by mutations. In this study, we model the trajectories resulting from MD simulations of the SARS-CoV-2 spike protein-ACE2, specifically the receptor-binding domain (RBD), as interresidue distance maps, and use deep convolutional neural networks to predict the functional impact of point mutations, related to the virus's infectivity and immunogenicity. Our model was successful in predicting mutant types that increase the affinity of the S protein for human receptors and reduce its immunogenicity, both based on MD trajectories (precision = 0.718; recall = 0.800; = 0.757; MCC = 0.488; AUC = 0.800) and their centroids. In an additional analysis, we also obtained a strong positive Pearson's correlation coefficient equal to 0.776, indicating a significant relationship between the average sigmoid probability for the MD trajectories and binding free energy (BFE) changes. Furthermore, we obtained a coefficient of determination of 0.602. Our 2D-RMSD analysis also corroborated predictions for more infectious and immune-evading mutants and revealed fluctuating regions within the receptor-binding motif (RBM), especially in the loop. This region presented a significant standard deviation for mutations that enable SARS-CoV-2 to evade the immune response, with RMSD values of 5Å in the simulation. This methodology offers an efficient alternative to identify potential strains of SARS-CoV-2, which may be potentially linked to more infectious and immune-evading mutations. Using clustering and deep learning techniques, our approach leverages information from the ensemble of MD trajectories to recognize a broad spectrum of multiple conformational patterns characteristic of mutant types. This represents a strategic advantage in identifying emerging variants, bypassing the need for long MD simulations. Furthermore, the present work tends to contribute substantially to the field of computational biology and virology, particularly to accelerate the design and optimization of new therapeutic agents and vaccines, offering a proactive stance against the constantly evolving threat of COVID-19 and potential future pandemics.

Enhanced sampling protocol to elucidate fusion peptide opening of SARS-CoV-2 spike protein

Cryo-electron Microscopy Structures of the SARS-CoV Spike Glycoprotein Reveal a Prerequisite Conformational State for Receptor Binding

Distinct conformational states of SARS-CoV-2 spike protein

Conformational dynamics and allosteric modulation of the SARS-CoV-2 spike

Conformational dynamics of SARS-CoV-2 Omicron spike trimers during fusion activation at single molecule resolution

Structure, Dynamics, Receptor Binding, and Antibody Binding of the Fully Glycosylated Full-Length SARS-CoV-2 Spike Protein in a Viral Membrane

Study of fusion peptide release for the spike protein of SARS-CoV-2

Simulation-Driven Design of Stabilized SARS-CoV-2 Spike S2 Immunogens

Deciphering the protein motion of S1 subunit in SARS-CoV-2 spike glycoprotein through integrated computational methods

Highly Conserved Homotrimer Cavity Formed by the SARS-CoV-2 Spike Glycoprotein: A Novel Binding Site

Nanometer resolution in situ structure of SARS-CoV-2 post-fusion spike

Modulation of SARS-CoV-2 spike binding to ACE2 through conformational selection

Host Cell Membrane Capture by the SARS-CoV-2 Spike Protein Fusion Intermediate

Deep learning for discriminating non-trivial conformational changes in molecular dynamics simulations of SARS-CoV-2 spike-ACE2

Exploring the Regulatory Function of the N-terminal Domain of SARS-CoV-2 Spike Protein Through Molecular Dynamics Simulation

In Silico Discovery of Small Molecule Modulators Targeting the Achilles’ Heel of SARS-CoV-2 Spike Protein

Nanometer-resolution in situ structure of the SARS-CoV-2 postfusion spike protein

Computational epitope map of SARS-CoV-2 spike protein

Structures and distributions of SARS-CoV-2 spike proteins on intact virions

Capture of fusion-intermediate conformations of SARS-CoV-2 spike requires receptor binding and cleavage at either the S1/S2 or S2' site

SARS-CoV-2 spike engagement of ACE2 primes S2′ site cleavage and fusion initiation