Evolution of Stronger SARS-CoV-2 Variants as Revealed Through the Lens of Molecular Dynamics Simulations

Writing a literature review as a class project in an upper-level undergraduate biochemistry course

A literature review is an important part of conducting academic research. Knowing how to conduct a literature search and write a high-quality literature review is a valuable skill. Herein, the authors describe the method of introducing a literature review writing exercise in an upper-level biochemistry course. Since 2020, authors have collaborated with numerous undergraduates writing literature reviews on topics in biochemistry that resulted in peer-reviewed publications. Authors believe that this unique idea of providing a course-based undergraduate research experience (CURE) to many undergraduates, especially those who otherwise do not receive collaborative research experience through traditional research paths, must be shared with other instructors.

Insights from in silico study of receptor energetics of SARS-CoV-2 variants

The emergence of new variants of the novel coronavirus SARS-CoV-2 with increased infectivity, superior virulence, high transmissibility, and unmatched immune escape has demonstrated the adaptability and evolutionary fitness of the virus. The subject of relative order of the binding affinity of SARS-CoV-2 variants with the human ACE2 (hACE2) receptor is hotly debated and its resolution has implications for drug design and development. In this work, we have investigated the energetics of the binding of receptor binding domain (RBD) of SARS-CoV-2 variants of concern (VOCs): Beta (B.1.351), Delta (B.1.617.2), Omicron (B.1.1.529), variant of interest (VOI): Kappa (B.1.617.1), and Delta Plus (B.1.617.2.1) variant with the human ACE2 receptor by using the umbrella sampling (US) method. Our work indicates that Delta and Delta Plus variants have greater values of the US binding free energy than Wild-type (WT), whereas Beta, Kappa, and Omicron variants have lower values. Further analysis of hydrogen bonding, salt bridges, non-bonded interaction energy, and contact surface area at the RBD-hACE2 interface establish Delta as the variant with the highest binding affinity among these variants.

Identification of key mutations responsible for the enhancement of receptor-binding affinity and immune escape of SARS-CoV-2 Omicron variant

The Omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has raised concerns worldwide due to its enhanced transmissibility and immune escapability. The first dominant Omicron BA.1 subvariant harbors more than 30 mutations in the spike protein from the prototype virus, of which 15 mutations are located at the receptor binding domain (RBD). These mutations in the RBD region attracted significant attention, which potentially enhance the binding of the receptor human angiotensin-converting enzyme 2 (hACE2) and decrease the potency of neutralizing antibodies/nanobodies. This study applied the molecular dynamics simulations combined with the molecular mechanics-generalized Born surface area (MMGBSA) method, to investigate the molecular mechanism behind the impact of the mutations acquired by Omicron on the binding affinity between RBD and hACE2. Our results indicate that five key mutations, i.e., N440K, T478K, E484A, Q493R, and G496S, contributed significantly to the enhancement of the binding affinity by increasing the electrostatic interactions of the RBD-hACE2 complex. Moreover, fourteen neutralizing antibodies/nanobodies complexed with RBD were used to explore the effects of the mutations in Omicron RBD on their binding affinities. The calculation results indicate that the key mutations E484A and Y505H reduce the binding affinities to RBD for most of the studied neutralizing antibodies/nanobodies, mainly attributed to the elimination of the original favorable gas-phase electrostatic and hydrophobic interactions between them, respectively. Our results provide valuable information for developing effective vaccines and antibody/nanobody drugs.