scholarly journals Combining Different Docking Engines and Consensus Strategies to Design and Validate Optimized Virtual Screening Protocols for the SARS-CoV-2 3CL Protease

Molecules ◽  
2021 ◽  
Vol 26 (4) ◽  
pp. 797
Author(s):  
Candida Manelfi ◽  
Jonas Gossen ◽  
Silvia Gervasoni ◽  
Carmine Talarico ◽  
Simone Albani ◽  
...  

The 3CL-Protease appears to be a very promising medicinal target to develop anti-SARS-CoV-2 agents. The availability of resolved structures allows structure-based computational approaches to be carried out even though the lack of known inhibitors prevents a proper validation of the performed simulations. The innovative idea of the study is to exploit known inhibitors of SARS-CoV 3CL-Pro as a training set to perform and validate multiple virtual screening campaigns. Docking simulations using four different programs (Fred, Glide, LiGen, and PLANTS) were performed investigating the role of both multiple binding modes (by binding space) and multiple isomers/states (by developing the corresponding isomeric space). The computed docking scores were used to develop consensus models, which allow an in-depth comparison of the resulting performances. On average, the reached performances revealed the different sensitivity to isomeric differences and multiple binding modes between the four docking engines. In detail, Glide and LiGen are the tools that best benefit from isomeric and binding space, respectively, while Fred is the most insensitive program. The obtained results emphasize the fruitful role of combining various docking tools to optimize the predictive performances. Taken together, the performed simulations allowed the rational development of highly performing virtual screening workflows, which could be further optimized by considering different 3CL-Pro structures and, more importantly, by including true SARS-CoV-2 3CL-Pro inhibitors (as learning set) when available.

2020 ◽  
Vol 21 (17) ◽  
pp. 5961
Author(s):  
Angelica Mazzolari ◽  
Silvia Gervasoni ◽  
Alessandro Pedretti ◽  
Laura Fumagalli ◽  
Rosanna Matucci ◽  
...  

Structure-based virtual screening is a truly productive repurposing approach provided that reliable target structures are available. Recent progresses in the structural resolution of the G-Protein Coupled Receptors (GPCRs) render these targets amenable for structure-based repurposing studies. Hence, the present study describes structure-based virtual screening campaigns with a view to repurposing known drugs as potential allosteric (and/or orthosteric) ligands for the hM2 muscarinic subtype which was indeed resolved in complex with an allosteric modulator thus allowing a precise identification of this binding cavity. First, a docking protocol was developed and optimized based on binding space concept and enrichment factor optimization algorithm (EFO) consensus approach by using a purposely collected database including known allosteric modulators. The so-developed consensus models were then utilized to virtually screen the DrugBank database. Based on the computational results, six promising molecules were selected and experimentally tested and four of them revealed interesting affinity data; in particular, dequalinium showed a very impressive allosteric modulation for hM2. Based on these results, a second campaign was focused on bis-cationic derivatives and allowed the identification of other two relevant hM2 ligands. Overall, the study enhances the understanding of the factors governing the hM2 allosteric modulation emphasizing the key role of ligand flexibility as well as of arrangement and delocalization of the positively charged moieties.


1994 ◽  
Vol 116 (9) ◽  
pp. 3697-3708 ◽  
Author(s):  
Luigi Gomez Paloma ◽  
Jarrod A. Smith ◽  
Walter J. Chazin ◽  
K. C. Nicolaou

2020 ◽  
Vol 21 (7) ◽  
pp. 2265 ◽  
Author(s):  
Carmine Talarico ◽  
Silvia Gervasoni ◽  
Candida Manelfi ◽  
Alessandro Pedretti ◽  
Giulio Vistoli ◽  
...  

Background: There is an increasing interest in TRPM8 ligands of medicinal interest, the rational design of which can be nowadays supported by structure-based in silico studies based on the recently resolved TRPM8 structures. Methods: The study involves the generation of a reliable hTRPM8 homology model, the reliability of which was assessed by a 1.0 μs MD simulation which was also used to generate multiple receptor conformations for the following structure-based virtual screening (VS) campaigns; docking simulations utilized different programs and involved all monomers of the selected frames; the so computed docking scores were combined by consensus approaches based on the EFO algorithm. Results: The obtained models revealed very satisfactory performances; LiGen™ provided the best results among the tested docking programs; the combination of docking results from the four monomers elicited a markedly beneficial effect on the computed consensus models. Conclusions: The generated hTRPM8 model appears to be amenable for successful structure-based VS studies; cross-talk modulating effects between interacting monomers on the binding sites can be accounted for by combining docking simulations as performed on all the monomers; this strategy can have general applicability for docking simulations involving quaternary protein structures with multiple identical binding pockets.


2020 ◽  
Author(s):  
Samuel C. Gill ◽  
David Mobley

<div>Sampling multiple binding modes of a ligand in a single molecular dynamics simulation is difficult. A given ligand may have many internal degrees of freedom, along with many different ways it might orient itself a binding site or across several binding sites, all of which might be separated by large energy barriers. We have developed a novel Monte Carlo move called Molecular Darting (MolDarting) to reversibly sample between predefined binding modes of a ligand. Here, we couple this with nonequilibrium candidate Monte Carlo (NCMC) to improve acceptance of moves.</div><div>We apply this technique to a simple dipeptide system, a ligand binding to T4 Lysozyme L99A, and ligand binding to HIV integrase in order to test this new method. We observe significant increases in acceptance compared to uniformly sampling the internal, and rotational/translational degrees of freedom in these systems.</div>


2017 ◽  
Author(s):  
Samuel Gill ◽  
Nathan M. Lim ◽  
Patrick Grinaway ◽  
Ariën S. Rustenburg ◽  
Josh Fass ◽  
...  

<div>Accurately predicting protein-ligand binding is a major goal in computational chemistry, but even the prediction of ligand binding modes in proteins poses major challenges. Here, we focus on solving the binding mode prediction problem for rigid fragments. That is, we focus on computing the dominant placement, conformation, and orientations of a relatively rigid, fragment-like ligand in a receptor, and the populations of the multiple binding modes which may be relevant. This problem is important in its own right, but is even more timely given the recent success of alchemical free energy calculations. Alchemical calculations are increasingly used to predict binding free energies of ligands to receptors. However, the accuracy of these calculations is dependent on proper sampling of the relevant ligand binding modes. Unfortunately, ligand binding modes may often be uncertain, hard to predict, and/or slow to interconvert on simulation timescales, so proper sampling with current techniques can require prohibitively long simulations. We need new methods which dramatically improve sampling of ligand binding modes. Here, we develop and apply a nonequilibrium candidate Monte Carlo (NCMC) method to improve sampling of ligand binding modes.</div><div><br></div><div>In this technique the ligand is rotated and subsequently allowed to relax in its new position through alchemical perturbation before accepting or rejecting the rotation and relaxation as a nonequilibrium Monte Carlo move. When applied to a T4 lysozyme model binding system, this NCMC method shows over two orders of magnitude improvement in binding mode sampling efficiency compared to a brute force molecular dynamics simulation. This is a first step towards applying this methodology to pharmaceutically relevant binding of fragments and, eventually, drug-like molecules. We are making this approach available via our new Binding Modes of Ligands using Enhanced Sampling (BLUES) package which is freely available on GitHub.</div>


2018 ◽  
Author(s):  
Samuel Gill ◽  
Nathan M. Lim ◽  
Patrick Grinaway ◽  
Ariën S. Rustenburg ◽  
Josh Fass ◽  
...  

<div>Accurately predicting protein-ligand binding is a major goal in computational chemistry, but even the prediction of ligand binding modes in proteins poses major challenges. Here, we focus on solving the binding mode prediction problem for rigid fragments. That is, we focus on computing the dominant placement, conformation, and orientations of a relatively rigid, fragment-like ligand in a receptor, and the populations of the multiple binding modes which may be relevant. This problem is important in its own right, but is even more timely given the recent success of alchemical free energy calculations. Alchemical calculations are increasingly used to predict binding free energies of ligands to receptors. However, the accuracy of these calculations is dependent on proper sampling of the relevant ligand binding modes. Unfortunately, ligand binding modes may often be uncertain, hard to predict, and/or slow to interconvert on simulation timescales, so proper sampling with current techniques can require prohibitively long simulations. We need new methods which dramatically improve sampling of ligand binding modes. Here, we develop and apply a nonequilibrium candidate Monte Carlo (NCMC) method to improve sampling of ligand binding modes.</div><div><br></div><div>In this technique the ligand is rotated and subsequently allowed to relax in its new position through alchemical perturbation before accepting or rejecting the rotation and relaxation as a nonequilibrium Monte Carlo move. When applied to a T4 lysozyme model binding system, this NCMC method shows over two orders of magnitude improvement in binding mode sampling efficiency compared to a brute force molecular dynamics simulation. This is a first step towards applying this methodology to pharmaceutically relevant binding of fragments and, eventually, drug-like molecules. We are making this approach available via our new Binding Modes of Ligands using Enhanced Sampling (BLUES) package which is freely available on GitHub.</div>


1993 ◽  
Vol 36 (13) ◽  
pp. 1902-1913 ◽  
Author(s):  
Krystina Plucinska ◽  
Takahiro Kataoka ◽  
Mitsuaki Yodo ◽  
Wayne L. Cody ◽  
J. X. He ◽  
...  

2015 ◽  
Vol 12 (110) ◽  
pp. 20150589 ◽  
Author(s):  
Maria C. Z. Meneghetti ◽  
Ashley J. Hughes ◽  
Timothy R. Rudd ◽  
Helena B. Nader ◽  
Andrew K. Powell ◽  
...  

Heparan sulfate (HS) polysaccharides are ubiquitous components of the cell surface and extracellular matrix of all multicellular animals, whereas heparin is present within mast cells and can be viewed as a more sulfated, tissue-specific, HS variant. HS and heparin regulate biological processes through interactions with a large repertoire of proteins. Owing to these interactions and diverse effects observed during in vitro , ex vivo and in vivo experiments, manifold biological/pharmacological activities have been attributed to them. The properties that have been thought to bestow protein binding and biological activity upon HS and heparin vary from high levels of sequence specificity to a dependence on charge. In contrast to these opposing opinions, we will argue that the evidence supports both a level of redundancy and a degree of selectivity in the structure–activity relationship. The relationship between this apparent redundancy, the multi-dentate nature of heparin and HS polysaccharide chains, their involvement in protein networks and the multiple binding sites on proteins, each possessing different properties, will also be considered. Finally, the role of cations in modulating HS/heparin activity will be reviewed and some of the implications for structure–activity relationships and regulation will be discussed.


Sign in / Sign up

Export Citation Format

Share Document