mkgridXf: Consistent Identification of Plausible Binding Sites Despite the Elusive Nature of Cavities and Grooves in Protein Dynamics

From complex data to clear insights: visualizing molecular dynamics trajectories

Advances in simulations, combined with technological developments in high-performance computing, have made it possible to produce a physically accurate dynamic representation of complex biological systems involving millions to billions of atoms over increasingly long simulation times. The analysis of these computed simulations is crucial, involving the interpretation of structural and dynamic data to gain insights into the underlying biological processes. However, this analysis becomes increasingly challenging due to the complexity of the generated systems with a large number of individual runs, ranging from hundreds to thousands of trajectories. This massive increase in raw simulation data creates additional processing and visualization challenges. Effective visualization techniques play a vital role in facilitating the analysis and interpretation of molecular dynamics simulations. In this paper, we focus mainly on the techniques and tools that can be used for visualization of molecular dynamics simulations, among which we highlight the few approaches used specifically for this purpose, discussing their advantages and limitations, and addressing the future challenges of molecular dynamics visualization.

HARIBOSS: a curated database of RNA-small molecules structures to aid rational drug design

MOTIVATION

RNA molecules are implicated in numerous fundamental biological processes and many human pathologies, such as cancer, neurodegenerative disorders, muscular diseases and bacterial infections. Modulating the mode of action of disease-implicated RNA molecules can lead to the discovery of new therapeutical agents and even address pathologies linked to 'undruggable' protein targets. This modulation can be achieved by direct targeting of RNA with small molecules. As of today, only a few RNA-targeting small molecules are used clinically. One of the main obstacles that have hampered the development of a rational drug design protocol to target RNA with small molecules is the lack of a comprehensive understanding of the molecular mechanisms at the basis of RNA-small molecule (RNA-SM) recognition.

RESULTS

Here, we present Harnessing RIBOnucleic acid-Small molecule Structures (HARIBOSS), a curated collection of RNA-SM structures determined by X-ray crystallography, nuclear magnetic resonance spectroscopy and cryo-electron microscopy. HARIBOSS facilitates the exploration of drug-like compounds known to bind RNA, the analysis of ligands and pockets properties and ultimately the development of in silico strategies to identify RNA-targeting small molecules.

AVAILABILITY AND IMPLEMENTATION

HARIBOSS can be explored via a web interface available at http://hariboss.pasteur.cloud.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A highly sensitive cell-based luciferase assay for high-throughput automated screening of SARS-CoV-2 nsp5/3CLpro inhibitors

Effective drugs against SARS-CoV-2 are urgently needed to treat severe cases of infection and for prophylactic use. The main viral protease (nsp5 or 3CLpro) represents an attractive and possibly broad-spectrum target for drug development as it is essential to the virus life cycle and highly conserved among betacoronaviruses. Sensitive and efficient high-throughput screening methods are key for drug discovery. Here we report the development of a gain-of-signal, highly sensitive cell-based luciferase assay to monitor SARS-CoV-2 nsp5 activity and show that it is suitable for the screening of compounds in a 384-well format. A benefit of miniaturisation and automation is that screening can be performed in parallel on a wild-type and a catalytically inactive nsp5, which improves the selectivity of the assay. We performed molecular docking-based screening on a set of 14,468 compounds from an in-house chemical database, selected 359 candidate nsp5 inhibitors and tested them experimentally. We identified two molecules which show anti-nsp5 activity, both in our cell-based assay and in vitro on purified nsp5 protein, and inhibit SARS-CoV-2 replication in A549-ACE2 cells with EC₅₀ values in the 4–8 μM range. The here described high-throughput-compatible assay will allow the screening of large-scale compound libraries for SARS-CoV-2 nsp5 inhibitors. Moreover, we provide evidence that this assay can be adapted to other coronaviruses and viruses which rely on a viral protease.

A highly sensitive cell-based luciferase assay for high-throughput automated screening of SARS-CoV-2 nsp5/3CLpro inhibitors

Effective drugs against SARS-CoV-2 are urgently needed to treat severe cases of infection and for prophylactic use. The main viral protease (nsp5 or 3CLpro) represents an attractive and possibly broad-spectrum target for drug development as it is essential to the virus life cycle and highly conserved among betacoronaviruses. Sensitive and efficient high-throughput screening methods are key for drug discovery. Here we report the development of a gain-of-signal, highly sensitive cell-based luciferase assay to monitor SARS-CoV-2 nsp5 activity and show that it is suitable for high-throughput screening of compounds in a 384-well format. A benefit of miniaturisation and automation is that screening can be performed in parallel on a wild-type and a catalytically inactive nsp5, which improves the selectivity of the assay. We performed molecular docking-based screening on a set of 14,468 compounds from an in-house chemical database, selected 359 candidate nsp5 inhibitors and tested them experimentally. We identified four molecules, including the broad-spectrum antiviral merimepodib/VX-497, which show anti-nsp5 activity and inhibit SARS-CoV-2 replication in A549-ACE2 cells with IC 50 values in the 4-21 µM range. The here described assay will allow the screening of large-scale compound libraries for SARS-CoV-2 nsp5 inhibitors. Moreover, we provide evidence that this assay can be adapted to other coronaviruses and viruses which rely on a viral protease.

DynaBiS: A hierarchical sampling algorithm to identify flexible binding sites for large ligands and peptides

Knowing the ligand or peptide binding site in proteins is highly important to guide drug discovery, but experimental elucidation of the binding site is difficult. Therefore, various computational approaches have been developed to identify potential binding sites in protein structures. However, protein and ligand flexibility are often neglected in these methods due to efficiency considerations despite the recognition that protein-ligand interactions can be strongly affected by mutual structural adaptations. This is particularly true if the binding site is unknown, as the screening will typically be performed based on an unbound protein structure. Herein we present DynaBiS, a hierarchical sampling algorithm to identify flexible binding sites for a target ligand with explicit consideration of protein and ligand flexibility, inspired by our previously presented flexible docking algorithm DynaDock. DynaBiS applies soft-core potentials between the ligand and the protein, thereby allowing a certain protein-ligand overlap resulting in efficient sampling of conformational adaptation effects. We evaluated DynaBiS and other commonly used binding site identification algorithms against a diverse evaluation set consisting of 26 proteins featuring peptide as well as small ligand binding sites. We show that DynaBiS outperforms the other evaluated methods for the identification of protein binding sites for large and highly flexible ligands such as peptides, both with a holo or apo structure used as input.

Analyzing In Silico the Relationship Between the Activation of the Edema Factor and Its Interaction With Calmodulin

Molecular dynamics (MD) simulations have been recorded on the complex between the edema factor (EF) of Bacilllus anthracis and calmodulin (CaM), starting from a structure with the orthosteric inhibitor adefovir bound in the EF catalytic site. The starting structure has been destabilized by alternately suppressing different co-factors, such as adefovir ligand or ions, revealing several long-distance correlations between the conformation of CaM, the geometry of the CaM/EF interface, the enzymatic site and the overall organization of the complex. An allosteric communication between CaM/EF interface and the EF catalytic site, highlighted by these correlations, was confirmed by several bioinformatics approaches from the literature. A network of hydrogen bonds and stacking interactions extending from the helix V of of CaM, and the residues of the switches A, B and C, and connecting to catalytic site residues, is a plausible candidate for the mediation of allosteric communication. The greatest variability in volume between the different MD conditions was also found for cavities present at the EF/CaM interface and in the EF catalytic site. The similarity between the predictions from literature and the volume variability might introduce the volume variability as new descriptor of allostery.

Protein Binding Pocket Optimization for Virtual High-Throughput Screening (vHTS) Drug Discovery

The virtual high-throughput screening (vHTS) approach has been widely used for large database screening to identify potential lead compounds for drug discovery. Due to its high computational demands, docking that allows receptor flexibility has been a challenging problem for virtual screening. Therefore, the selection of protein target conformations is crucial to produce useful vHTS results. Since only a single protein structure is used to screen large databases in most vHTS studies, the main challenge is to reduce false negative rates in selecting compounds for in vitro tests. False negatives are most likely to occur when using apo structures or homology models of protein targets due to the small volume of the binding pocket formed by incorrect side-chain conformations. Even holo protein structures can exhibit high false negative rates due to ligand-induced fit effects, since the shape of the binding pocket highly depends on its bound ligand. To reduce false negative rates and improve success rates for vHTS in drug discovery, we have developed a new Monte Carlo-based approach that optimizes the binding pocket of protein targets. This newly developed Monte Carlo pocket optimization (MCPO) approach was assessed on several datasets showing promising results. The binding pocket optimization approach could be a useful tool for vHTS-based drug discovery, especially in cases when only apo structures or homology models are available.

Bacterial intracellularly active toxins: Membrane localisation of the active domain

Numerous bacterial toxins exert their activity by inactivating or modulating a specific intracellular host target. For this purpose, these toxins have developed efficient strategies to overcome the different host cell defences including specific binding to cell surface, internalisation, passage through the endosome or plasma membrane, exploiting intracellular trafficking and addressing to intracellular targets. Several intracellularly active toxins deliver an active domain into the cytosol that interacts with a target localised to the inner face of the plasma membrane. Thus, the large clostridial glucosylating toxins (LCGTs) target Rho/Ras-GTPases, certain virulence factors of Gram negative bacteria, Rho-GTPases, while Pasteurella multocida toxin (PMT) targets trimeric G-proteins. Others such as botulinum neurotoxins and tetanus neurotoxin have their substrate on synaptic vesicle membrane. LCGTs, PMT, and certain virulence factors from Vibrio sp. show a particular structure constituted of a four-helix bundle membrane (4HBM) protruding from the catalytic site that specifically binds to the membrane phospholipids and then trap the catalytic domain at the proximity of the membrane anchored substrate. Structural and functional analysis indicate that the 4HBM tip of the Clostridium sordellii lethal toxin (TcsL) from the LCGT family contain two loops forming a cavity that mediates the binding to phospholipids and more specifically to phosphatidylserine.