Reconstruction of the unbinding pathways of new inhibitors of the SARS-CoV-2 papain-like protease using molecular dynamics simulation

Developing novel antiviral drugs against the SARS-CoV-2 virus and COVID-19 disease is imperative as the vaccines may not offer absolute protection. PLpro plays a crucial role in the viral life cycle, making it an attractive target for drug development. Several PLpro inhibitors have been developed, and their 3D structures in complex with PLpro are available. In this work, we employed Supervised Molecular Dynamics (SuMD), a specific Unbiased Molecular Dynamics (UMD) method, to investigate unbinding pathways of the novel inhibitors of PLpro (PDB IDs: 7LBR, 7RZC, 7SDR and 7E35) and GRL0617 (PDB ID: 7JRN) as a reference. We conducted three simulations for each ligand and achieved unbinding events in the nanosecond timescale in all simulations. We found that unbinding events are commonly affected by altering the conformation of the BL2 loop, which is caused by the natural fluctuations of the loop that are required to trap the substrate and throw out the product. BL2 loop is crucial for keeping the ligand and unbinding and acts as a double-edged sword. Any inhibitor designed to be effective must prevent the loop's natural fluctuations. We perceived that increasing ligands interactions with the binding pocket interior and the BL2 loop will help prevent natural fluctuation of the BL2 loop, Although the interactions with the binding pocket's inner side are more critical than the BL2 loop. These findings may be helpful in developing more potent inhibitors against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

HT-SuMD: making molecular dynamics simulations suitable for fragment-based screening. A comparative study with NMR

Fragment-based lead discovery (FBLD) is one of the most efficient methods to develop new drugs. We present here a new computational protocol called High-Throughput Supervised Molecular Dynamics (HT-SuMD), which makes it possible to automatically screen up to thousands of fragments, representing therefore a new valuable resource to prioritise fragments in FBLD campaigns. The protocol was applied to Bcl-XL, an oncological protein target involved in the regulation of apoptosis through protein-protein interactions. Initially, HT-SuMD performances were validated against a robust NMR-based screening, using the same set of 100 fragments. These independent results showed a remarkable agreement between the two methods. Then, a virtual screening on a larger library of additional 300 fragments was carried out and the best hits were validated by NMR. Remarkably, all the in silico selected fragments were confirmed as Bcl-XL binders. This represents, to date, the largest computational fragments screening entirely based on MD.

Supervised Molecular Dynamics (SuMD) Insights into the mechanism of action of SARS-CoV-2 main protease inhibitor PF-07321332

The chemical structure of PF-07321332, the first orally available Covid-19 clinical candidate, has recently been revealed by Pfizer. No information has been provided about the interaction pattern between PF-07321332 and its biomolecular counterpart, the SARS-CoV-2 main protease (M^pro). In the present work, we exploited Supervised Molecular Dynamics (SuMD) simulations to elucidate the key features that characterise the interaction between this drug candidate and the protease, emphasising similarities and differences with other structurally related inhibitors such as Boceprevir and PF-07304814. The structural insights provided by SuMD will hopefully be able to inspire the rational discovery of other potent and selective protease inhibitors.

Addressing free fatty acid receptor 1 (FFAR1) activation using supervised molecular dynamics

The free fatty acid receptor 1 (FFAR1, formerly GPR40), is a potential G protein-coupled receptor (GPCR) target for the treatment of type 2 diabetes mellitus (T2DM), as it enhances glucose-dependent insulin secretion upon activation by endogenous long-chain free fatty acids. The presence of two allosterically communicating binding sites and the lack of the conserved GPCR structural motifs challenge the general knowledge of its activation mechanism. To date, four X-ray crystal structures are available for computer-aided drug design. In this study, we employed molecular dynamics (MD) and supervised molecular dynamics (SuMD) to deliver insights into the (un)binding mechanism of the agonist MK-8666, and the allosteric communications between the two experimentally determined FFAR1 binding sites. We found that FFAR1 extracellular loop 2 (ECL2) mediates the binding of the partial agonist MK-8666. Moreover, simulations showed that the agonists MK-8666 and AP8 are reciprocally stabilized and that AP8 influences MK-8666 unbinding from FFAR1.