Atomistic Simulations and In Silico Mutational Profiling of Protein Stability and Binding in the SARS-CoV-2 Spike Protein Complexes with Nanobodies: Molecular Determinants of Mutational Escape Mechanisms

Interaction of SARS-CoV-2 with host cells and antibodies: experiment and simulation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the devastating global COVID-19 pandemic announced by WHO in March 2020. Through unprecedented scientific effort, several vaccines, drugs and antibodies have been developed, saving millions of lives, but the fight against COVID-19 continues as immune escape variants of concern such as Delta and Omicron emerge. To develop more effective treatments and to elucidate the side effects caused by vaccines and therapeutic agents, a deeper understanding of the molecular interactions of SARS-CoV-2 with them and human cells is required. With special interest in computational approaches, we will focus on the structure of SARS-CoV-2 and the interaction of its spike protein with human angiotensin-converting enzyme-2 (ACE2) as a prime entry point of the virus into host cells. In addition, other possible viral receptors will be considered. The fusion of viral and human membranes and the interaction of the spike protein with antibodies and nanobodies will be discussed, as well as the effect of SARS-CoV-2 on protein synthesis in host cells.

Binding of synthetic nanobodies to the SARS-CoV-2 receptor-binding domain: the importance of salt bridges

In this study, five different SARS-CoV-2 receptor-binding domain (RBD) models were created based on the crystal structures of RBD complexes with two synthetic nanobodies (Sb16 and Sb45). Microsecond all-atom MD simulations revealed that Sb16 and Sb45 substantially stabilized the flexible RBD loop (residues GLU471-SER494) due to the salt bridges and hydrogen bonding interactions between RBD and the synthetic nanobodies. However, the calculation of binding free energy displayed that Sb45 had a higher binding affinity to RBD than Sb16, in agreement with the experimental result. This is because Sb45 has stronger electrostatic attraction to RBD as compared to Sb16. In particular, the salt bridge GLU484-ARG33 in Sb45-RBD is stronger than the GLU484-LYS32 in Sb16-RBD. Furthermore, by comparing the binding affinity of Sb16 for two RBD mutants (E484K and K417N), we found that E484K mutation substantially reduced the binding affinity to Sb16, and K417N mutation had no significant effect, qualitatively in agreement with experimental studies. According to the binding free energy calculation, the strong electrostatic repulsion between LYS32 and LYS484 caused by E484K mutation destroys the salt bridge between LYS32 and GLU484 in the RBD wild type (WT). In contrast, the binding of the K417N mutant to Sb16 effectively maintains the salt bridge between LYS32 and GLU484. Therefore, our research suggests that the salt bridges between RBD and synthetic nanobodies are crucial for binding synthetic nanobodies to RBD, and a SARS-CoV-2 variant can escape neutralization from nanobodies by creating electrostatic repulsion between them.

Outsmarting Pathogens with Antibody Engineering

There is growing interest in identifying antibodies that protect against infectious diseases, especially for high-risk individuals and pathogens for which no vaccine is yet available. However, pathogens that manifest as opportunistic or latent infections express complex arrays of virulence-associated proteins and are adept at avoiding immune responses. Some pathogens have developed strategies to selectively destroy antibodies, whereas others create decoy epitopes that trick the host immune system into generating antibodies that are at best nonprotective and at worst enhance pathogenesis. Antibody engineering strategies can thwart these efforts by accessing conserved neutralizing epitopes, generating Fc domains that resist capture or degradation and even accessing pathogens hidden inside cells. Design of pathogen-resistant antibodies can enhance protection and guide development of vaccine immunogens against these complex pathogens. Here, we discuss general strategies for design of antibodies resistant to specific pathogen defense mechanisms.

Machine learning and protein allostery

The fundamental biological importance and complexity of allosterically regulated proteins stem from their central role in signal transduction and cellular processes. Recently, machine-learning approaches have been developed and actively deployed to facilitate theoretical and experimental studies of protein dynamics and allosteric mechanisms. In this review, we survey recent developments in applications of machine-learning methods for studies of allosteric mechanisms, prediction of allosteric effects and allostery-related physicochemical properties, and allosteric protein engineering. We also review the applications of machine-learning strategies for characterization of allosteric mechanisms and drug design targeting SARS-CoV-2. Continuous development and task-specific adaptation of machine-learning methods for protein allosteric mechanisms will have an increasingly important role in bridging a wide spectrum of data-intensive experimental and theoretical technologies.

General Trends of the Camelidae Antibody VHHs Domain Dynamics

Conformational flexibility plays an essential role in antibodies' functional and structural stability. They facilitate and determine the strength of antigen-antibody interactions. Camelidae express an interesting subtype of single-chain antibody, named Heavy Chain only Antibody. They have only one N-terminal Variable domain (V_HH) per chain, composed of Frameworks (FRs) and Complementarity Determining regions (CDRs) like their VH and VL counterparts in IgG. Even when expressed independently, V_HH domains display excellent solubility and (thermo)stability, which helps them to retain their impressive interaction capabilities. Sequence and structural features of V_HH domains contributing to these abilities have already been studied compared to classical antibodies. To have the broadest view and understand the changes in dynamics of these macromolecules, large-scale molecular dynamics simulations for a large number of non-redundant V_HH structures have been performed for the first time. This analysis reveals the most prevalent movements in these domains. It reveals the four main classes of V_HHs dynamics. Diverse local changes were observed in CDRs with various intensities. Similarly, different types of constraints were observed in CDRs, while FRs close to CDRs were sometimes primarily impacted. This study sheds light on the changes in flexibility in different regions of V_HH that may impact their in silico design.

Current Update on Rotavirus in-Silico Multiepitope Vaccine Design

Rotavirus gastroenteritis is one of the leading causes of pediatric morbidity and mortality worldwide in infants and under-five populations. The World Health Organization (WHO) recommended global incorporation of the rotavirus vaccine in national immunization programs to alleviate the burden of the disease. Implementation of the rotavirus vaccination in certain regions of the world brought about a significant and consistent reduction of rotavirus-associated hospitalizations. However, the efficacy of licensed vaccines remains suboptimal in low-income countries where the incidences of rotavirus gastroenteritis continue to happen unabated. The problem of low efficacy of currently licensed oral rotavirus vaccines in low-income countries necessitates continuous exploration, design, and development of new rotavirus vaccines. Traditional vaccine development is a complex, expensive, labor-intensive, and time-consuming process. Reverse vaccinology essentially utilizes the genome and proteome information on pathogens and has opened new avenues for in-silico multiepitope vaccine design for a plethora of pathogens, promising time reduction in the complete vaccine development pipeline by complementing the traditional vaccinology approach. A substantial number of reviews on licensed rotavirus vaccines and those under evaluation are already available in the literature. However, a collective account of rotavirus in-silico vaccines is lacking in the literature, and such an account may further fuel the interest of researchers to use reverse vaccinology to expedite the vaccine development process. Therefore, the main focus of this review is to summarize the research endeavors undertaken for the design and development of rotavirus vaccines by the reverse vaccinology approach utilizing the tools of immunoinformatics.

Structural and Computational Studies of the SARS-CoV-2 Spike Protein Binding Mechanisms with Nanobodies: From Structure and Dynamics to Avidity-Driven Nanobody Engineering

Nanobodies provide important advantages over traditional antibodies, including their smaller size and robust biochemical properties such as high thermal stability, high solubility, and the ability to be bioengineered into novel multivalent, multi-specific, and high-affinity molecules, making them a class of emerging powerful therapies against SARS-CoV-2. Recent research efforts on the design, protein engineering, and structure-functional characterization of nanobodies and their binding with SARS-CoV-2 S proteins reflected a growing realization that nanobody combinations can exploit distinct binding epitopes and leverage the intrinsic plasticity of the conformational landscape for the SARS-CoV-2 S protein to produce efficient neutralizing and mutation resistant characteristics. Structural and computational studies have also been instrumental in quantifying the structure, dynamics, and energetics of the SARS-CoV-2 spike protein binding with nanobodies. In this review, a comprehensive analysis of the current structural, biophysical, and computational biology investigations of SARS-CoV-2 S proteins and their complexes with distinct classes of nanobodies targeting different binding sites is presented. The analysis of computational studies is supplemented by an in-depth examination of mutational scanning simulations and identification of binding energy hotspots for distinct nanobody classes. The review is focused on the analysis of mechanisms underlying synergistic binding of multivalent nanobodies that can be superior to single nanobodies and conventional nanobody cocktails in combating escape mutations by effectively leveraging binding avidity and allosteric cooperativity. We discuss how structural insights and protein engineering approaches together with computational biology tools can aid in the rational design of synergistic combinations that exhibit superior binding and neutralization characteristics owing to avidity-mediated mechanisms.

Allosteric Determinants of the SARS-CoV-2 Spike Protein Binding with Nanobodies: Examining Mechanisms of Mutational Escape and Sensitivity of the Omicron Variant

Structural and biochemical studies have recently revealed a range of rationally engineered nanobodies with efficient neutralizing capacity against the SARS-CoV-2 virus and resilience against mutational escape. In this study, we performed a comprehensive computational analysis of the SARS-CoV-2 spike trimer complexes with single nanobodies Nb6, VHH E, and complex with VHH E/VHH V nanobody combination. We combined coarse-grained and all-atom molecular simulations and collective dynamics analysis with binding free energy scanning, perturbation-response scanning, and network centrality analysis to examine mechanisms of nanobody-induced allosteric modulation and cooperativity in the SARS-CoV-2 spike trimer complexes with these nanobodies. By quantifying energetic and allosteric determinants of the SARS-CoV-2 spike protein binding with nanobodies, we also examined nanobody-induced modulation of escaping mutations and the effect of the Omicron variant on nanobody binding. The mutational scanning analysis supported the notion that E484A mutation can have a significant detrimental effect on nanobody binding and result in Omicron-induced escape from nanobody neutralization. Our findings showed that SARS-CoV-2 spike protein might exploit the plasticity of specific allosteric hotspots to generate escape mutants that alter response to binding without compromising activity. The network analysis supported these findings showing that VHH E/VHH V nanobody binding can induce long-range couplings between the cryptic binding epitope and ACE2-binding site through a broader ensemble of communication paths that is less dependent on specific mediating centers and therefore may be less sensitive to mutational perturbations of functional residues. The results suggest that binding affinity and long-range communications of the SARS-CoV-2 complexes with nanobodies can be determined by structurally stable regulatory centers and conformationally adaptable hotspots that are allosterically coupled and collectively control resilience to mutational escape.