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.

Nanomechanical collective vibration of SARS‐CoV‐2 spike proteins

Physical-Chemical Regulation of Membrane Receptors Dynamics in Viral Invasion and Immune Defense

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

Mass Spectrometry Analysis of SARS-CoV-2 Nucleocapsid Protein Reveals Camouflaging Glycans and Unique Post-Translational Modifications

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

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

Mechanical activation of spike fosters SARS-CoV-2 viral infection

Computational analysis of affinity dynamics between the variants of SARS-CoV-2 spike protein (RBD) and human ACE-2 receptor

Probing the Allosteric Inhibition Mechanism of a Spike Protein Using Molecular Dynamics Simulations and Active Compound Identifications

Probing Nanomechanical Interactions of SARS-CoV-2 Variants Omicron and XBB with Common Surfaces.

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

Single-molecule imaging reveals allosteric stimulation of SARS-CoV-2 spike receptor binding domain by host sialic acid

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

Nanomechanical sonification of the 2019-nCoV coronavirus spike protein through a materiomusical approach

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

Molecular Dynamics Studies on the Structural Characteristics for the Stability Prediction of SARS-CoV-2

Understanding the Excitation Wavelength Dependence and Thermal Stability of the SARS-CoV-2 Receptor-Binding Domain Using Surface-Enhanced Raman Scattering and Machine Learning

Activation Pathways and Free Energy Landscapes of the SARS-CoV-2 Spike Protein

Unraveling the binding mechanisms of SARS-CoV-2 variants through molecular simulations

Effect of receptors on the resonant and transient harmonic vibrations of Coronavirus

Nanomechanical footprint of SARS-CoV-2 variants in complex with a potent nanobody by molecular simulations