Dynamical nonequilibrium molecular dynamics simulations identify allosteric sites and positions associated with drug resistance in the SARS-CoV-2 main protease

The SARS-CoV-2 main protease (Mpro) plays an essential role in the coronavirus lifecycle by catalysing hydrolysis of the viral polyproteins at specific sites. Mpro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how Mpro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of Mpro to the presence and absence of a substrate. The results highlight communication between the Mpro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behaviour of Mpro. More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance.

Molecular dynamics of spike variants in the locked conformation: RBD interfaces, fatty acid binding and furin cleavage sites.

Since December 2019 the SARS-CoV-2 virus has infected billions of people around the world and caused millions of deaths. The ability for this RNA virus to mutate has produced variants that have been responsible for waves of infections across the globe. The spike protein on the surface of the SARS-CoV-2 virion is responsible for cell entry in the infection process. Here we have studied the spike proteins from the Original, Alpha (B.1.1.7), Delta (B1.617.2), Delta-plus (B1.617.2-AY1), Omicron BA.1 and Omicron BA.2 variants. Using models built from cryo-EM structures with linoleate bound (6BZ5.pdb) and the N-terminal domain from 7JJI.pdb, each is built from the first residue, with missing loops modelled and 45 disulphides per trimer. Each spike variant was modified from the same Original model framework to maximise comparability. Three replicate, 200 ns atomistic molecular dynamics simulations were performed for each case. (These data also provide the basis for further, non-equilibrium molecular dynamics simulations, published elsewhere.) The analysis of our equilibrium molecular dynamics reveals that sequence variation at the closed receptor binding domain interface particularly for Omicron BA.2 has implications for the avidity of the locked conformation, with potential effects on Omicron BA.1 and Delta-plus. Linoleate binding has a mildly stabilizing effect on furin cleavage site motions in the Original and Alpha variants, but has no effect in Delta, Delta-plus and slightly increases motions at this site for Omicron BA.1, but not BA.2, under these simulation conditions.

SARS-CoV-2 spike variants differ in their allosteric response to linoleic acid

The SARS-CoV-2 spike protein contains a fatty acid binding site, also found in some other coronaviruses (e.g. SARS-CoV), which binds linoleic acid and is functionally important. When occupied by linoleic acid, it reduces infectivity, by locking the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the response of spike variants to linoleic acid removal. These simulations show that the fatty acid site is coupled to functional regions of the protein, some of them far from the site (e.g. in the receptor-binding motif, N-terminal domain, the furin cleavage site located in position 679-685 and the fusion peptide-surrounding regions) and identify the allosteric networks involved in these connections. Comparison of the response of the original ( Wuhan) spike with four variants: Alpha, Delta, Delta plus and Omicron BA.1 show that the variants differ significantly in their response to linoleic acid removal. The allosteric connections to the fatty acid site on Alpha are generally similar to the original protein, except for the receptor-binding motif and S71-R78 region which show a weaker link to the FA site. In contrast, Omicron is the most affected variant exhibiting significant differences in the receptor-binding motif, N-terminal domain, V622-L629 and the furin cleavage site. These differences in allosteric modulation may be of functional relevance, e.g. in differences in transmissibility and virulence. Experimental comparison of the effects of linoleic acid on different variants is warranted.

Discovery of SARS-CoV-2 Mpro Peptide Inhibitors from Modelling Substrate and Ligand Binding

The main protease (Mpro) of SARS-CoV-2 is central to its viral lifecycle and is a promising drug target, but little is known concerning structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of classical molecular mechanics and quantum mechanical techniques, including automated docking, molecular dynamics (MD) simulations, linear-scaling DFT, QM/MM, and interactive MD in virtual reality, to investigate the molecular features underlying recognition of the natural Mpro substrates. Analyses of the subsite interactions of modelled 11-residue cleavage site peptides, ligands from high-throughput crystallography, and designed covalently binding inhibitors were performed. Modelling studies reveal remarkable conservation of hydrogen bonding patterns of the natural Mpro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular at the P2/S2 sites. The binding modes of the natural substrates, together with extensive interaction analyses of inhibitor and fragment binding to Mpro, reveal new opportunities for inhibition. Building on our initial Mpro-substrate models, computational mutagenesis scanning was employed to design peptides with improved affinity and which inhibit Mpro competitively. The combined results provide new insight useful for the development of Mpro inhibitors.

The fatty acid site is coupled to functional motifs in the SARS-CoV-2 spike protein and modulates spike allosteric behaviour

The SARS-CoV-2 spike protein is the first contact point between the SARS-CoV-2 virus and host cells and mediates membrane fusion. Recently, a fatty acid binding site was identified in the spike (Toelzer et al. Science 2020). The presence of linoleic acid at this site modulates binding of the spike to the human ACE2 receptor, stabilizing a locked conformation of the protein. Here, dynamical-nonequilibrium molecular dynamics simulations reveal that this fatty acid site is coupled to functionally relevant regions of the spike, some of them far from the fatty acid binding pocket. Removal of a ligand from the fatty acid binding site significantly affects the dynamics of distant, functionally important regions of the spike, including the receptor-binding motif, furin cleavage site and fusion-peptide-adjacent regions. The results also show significant differences in behaviour between clinical variants of the spike: e.g. the D614G mutation shows a significantly different conformational response for some structural motifs relevant for binding and fusion. The simulations identify structural networks through which changes at the fatty acid binding site are transmitted within the protein. These communication networks significantly involve positions that are prone to mutation, indicating that observed genetic variation in the spike may alter its response to linoleate binding and associated allosteric communication.

Structure and mechanism of SARS-CoV-2 Spike N679-V687 deletion variant elucidate cell-type specific evolution of viral fitness

As the global burden of SARS-CoV-2 infections escalates, so does the evolution of viral variants which is of particular concern due to their potential for increased transmissibility and pathology. In addition to this entrenched variant diversity in circulation, RNA viruses can also display genetic diversity within single infected hosts with co-existing viral variants evolving differently in distinct cell types. The BriS{Delta} variant, originally identified as a viral subpopulation by passaging SARS-CoV-2 isolate hCoV-19/England/02/2020, comprises in the spike glycoprotein an eight amino-acid deletion encompassing the furin recognition motif and S1/S2 cleavage site. Here, we elucidate the structure, function and molecular dynamics of this variant spike providing mechanistic insight into how the deletion correlates to viral cell tropism, ACE2 receptor binding and infectivity of this SARS-CoV-2 variant. Moreover, our study reveals long-range allosteric communication between functional regions within the spike that differ in wild-type and deletion variant. Our results support a view of SARS-CoV-2 probing multiple evolutionary trajectories in distinct cell types within the same infected host.

Simulations support the interaction of the SARS-CoV-2 spike protein with nicotinic acetylcholine receptors and suggest subtype specificity

Changeux et al. recently suggested that the SARS-CoV-2 spike (S) protein may interact with nicotinic acetylcholine receptors (nAChRs). Such interactions may be involved in pathology and infectivity. Here, we use molecular simulations of validated atomically detailed structures of nAChRs, and of the S protein, to investigate this nicotinic hypothesis. We examine the binding of the Y674-R685 loop of the S protein to three nAChRs, namely the human 4{beta}2 and 7 subtypes and the muscle-like {beta}{gamma}d receptor from Tetronarce californica. Our results indicate that Y674-R685 has affinity for nAChRs and the region responsible for binding contains the PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. In particular, R682 has a key role in the stabilisation of the complexes as it forms interactions with loops A, B and C in the receptors binding pocket. The conformational behaviour of the bound Y674-R685 region is highly dependent on the receptor subtype, adopting extended conformations in the 4{beta}2 and 7 complexes and more compact ones when bound to the muscle-like receptor. In the 4{beta}2 and {beta}{gamma}d complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation similar to other known nAChR antagonists. In contrast, in the 7 complex, Y674-R685 penetrates deeply into the binding pocket where it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1 and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of the simulations of the full-length S protein show that the Y674-R685 region is accessible for binding, and suggest a potential binding orientation of the S protein with nAChRs.