Rapidly spreading SARS-CoV-2 variants present not only an increased threat to human health due to the confirmed greater transmissibility of several of these new strains but, due to conformational changes induced by the mutations, may render first-wave SARS-CoV-2 convalescent sera, vaccine-induced antibodies, or recombinant neutralizing antibodies (nAbs) ineffective. To be able to assess the risk of viral escape from neutralization by first-wave antibodies, we leveraged our capability for Molecular Dynamic (MD) simulation of the spike receptor binding domain (S RBD) and its binding to human angiotensin-converting enzyme 2 (hACE2) to predict alterations in molecular interactions resulting from the presence of the E484K, K417N, and N501Y variants found in the South African 501Y.V2 strain – alone and in combination. We report here the combination of E484K, K417N and N501Y results in the highest degree of conformational alterations of S RBD when bound to hACE2, compared to either E484K or N501Y alone. Both E484K and N501Y increase affinity of S RBD for hACE2 and E484K in particular switches the charge on the flexible loop region of RBD which leads to the formation of novel favorable contacts. Enhanced affinity of S RBD for hACE2 very likely underpins the greater transmissibility conferred by the presence of either E484K or N501Y; while the induction of conformational changes may provide an explanation for evidence that the 501Y.V2 variant, distinguished from the B.1.1.7 UK variant by the presence of E484K, is able to escape neutralization by existing first-wave anti-SARS-CoV-2 antibodies and re-infect COVID-19 convalescent individuals.

Northwestern blot analysis in the presence of competitor RNA was used to examine the interaction between the mouse hepatitis virus (MHV) nucleocapsid protein (N) and virus-specific RNAs. Our accompanying article demonstrates that anti-N monoclonal antibodies immunoprecipitated all seven MHV-specific RNAs as well as the small leader-containing RNAs from infected cells. In this article we report that a Northwestern blotting protocol using radiolabeled viral RNAs in the presence of host cell competitor RNA can be used to demonstrate a high-affinity interaction between the MHV N protein and the virus-specific RNAs. Further, RNA probes prepared by in vitro transcription were used to define the sequences that participate in such high-affinity binding. A specific interaction occurs between the N protein and sequences contained with the leader RNA which is conserved at the 5' end of all MHV RNAs. We have further defined the binding sites to the area of nucleotides 56 to 65 at the 3' end of the leader RNA and suggest that this interaction may play an important role in the discontinuous nonprocessive RNA transcriptional process unique to coronaviruses.

Alchemical free energy (AFE) calculations based on molecular dynamics (MD) simulations are key tools in both improving our understanding of a wide variety of biological processes and accelerating the design and optimization of therapeutics for numerous diseases. Computing power and theory have, however, long been insufficient to enable AFE calculations to be routinely applied in early stage drug discovery. One of the major difficulties in performing AFE calculations is the length of time required for calculations to converge to an ensemble average. CPU implementations of MD-based free energy algorithms can effectively only reach tens of nanoseconds per day for systems on the order of 50,000 atoms, even running on massively parallel supercomputers. Therefore, converged free energy calculations on large numbers of potential lead compounds are often untenable, preventing researchers from gaining crucial insight into molecular recognition, potential druggability and other crucial areas of interest. Graphics Processing Units (GPUs) can help address this. We present here a seamless GPU implementation, within the PMEMD module of the AMBER molecular dynamics package, of thermodynamic integration (TI) capable of reaching speeds of >140 ns/day for a 44,907-atom system, with accuracy equivalent to the existing CPU implementation in AMBER. The implementation described here is currently part of the AMBER 18 beta code and will be an integral part of the upcoming version 18 release of AMBER. © 2018 Wiley Periodicals, Inc.

The SARS-CoV-2 variants replacing the first wave strain pose an increased threat by their potential ability to escape pre-existing humoral protection. An angiotensin converting enzyme 2 (ACE2) decoy that competes with endogenous ACE2 for binding of the SARS-CoV-2 spike receptor binding domain (S RBD) and inhibits infection may offer a therapeutic option with sustained efficacy against variants. Here, we used Molecular Dynamics (MD) simulation to predict ACE2 sequence substitutions that might increase its affinity for S RBD and screened candidate ACE2 decoys in vitro. The lead ACE2(T27Y/H34A)-IgG1FC fusion protein with enhanced S RBD affinity shows greater live SARS-CoV-2 virus neutralization capability than wild type ACE2. MD simulation was used to predict the effects of S RBD variant mutations on decoy affinity that was then confirmed by testing of an ACE2 Triple Decoy that included an additional enzyme activity-deactivating H374N substitution against mutated S RBD. The ACE2 Triple Decoy maintains high affinity for mutated S RBD, displays enhanced affinity for S RBD N501Y or L452R, and has the highest affinity for S RBD with both E484K and N501Y mutations, making it a viable therapeutic option for the prevention or treatment of SARS-CoV-2 infection with a high likelihood of efficacy against variants.

The Receptor Binding Domain (RBD) of the SARS-CoV-2 surface spike (S) protein interacts with host angiotensin converting enzyme 2 (ACE2) to gain entry to host cells and initiate infection1–3. Detailed, accurate understanding of key interactions between S RBD and ACE2 provides critical information that may be leveraged in the development of strategies for the prevention and treatment of COVID-19. Utilizing the published sequences and cryo-EM structures of both the viral S RBD and ACE24,5, we performed in silico molecular dynamics (MD) simulations of free S RBD and of its interaction with ACE2 over the exceptionally long durations of 2.9 and 2 milliseconds, respectively, to elucidate the nature and relative affinity of S RBD surface residues for the ACE2 binding region. Our findings reveal that free S RBD has assumed an optimized ACE2 binding-ready conformation, incurring little entropic penalty for binding, an evolutionary adaptation that contributes to its high affinity for the receptor6. We further identified high probability molecular binding interactions that inform both vaccine design and therapeutic development, which may include recombinant ACE2-based spike decoys7 and/or allosteric S RBD-ACE2 binding inhibitors8,9 to prevent or arrest infection and thus disease.