• Fragment-based drug design (FBDD) is revolutionizing the identification and optimization of new drug candidates. • Molecular simulations-based approaches are consolidating alongside with experimental and canonical structure-based (SBDD) techniques. • The implementation of molecular simulations-based techniques affects the entire FBDD pipeline.

The use and misuse of antibiotics has resulted in critical conditions for drug-resistant bacteria emergency, accelerating the development of antimicrobial resistance (AMR). In this context, the co-administration of an antibiotic with a compound able to restore sufficient antibacterial activity may be a successful strategy. In particular, the identification of efflux pump inhibitors (EPIs) holds promise for new antibiotic resistance breakers (ARBs). Indeed, bacterial efflux pumps have a key role in AMR development; for instance, NorA efflux pump contributes to Staphylococcus aureus (S. aureus) resistance against fluoroquinolone antibiotics (e.g., ciprofloxacin) by promoting their active extrusion from the cells. Even though NorA efflux pump is known to be a potential target for EPIs development, the absence of structural information about this protein and the little knowledge available on its mechanism of action have strongly hampered rational drug discovery efforts in this area. In the present work, we investigated at the molecular level the substrate recognition pathway of NorA through a Supervised Molecular Dynamics (SuMD) approach, using a NorA homology model. Specific amino acids were identified as playing a key role in the efflux pump-mediated extrusion of its substrate, paving the way for a deeper understanding of both the mechanisms of action and the inhibition of such efflux pumps.

Fragment-Based Drug Discovery (FBDD) has become, in recent years, a consolidated approach in the drug discovery process, leading to several drug candidates under investigation in clinical trials and some approved drugs. Among these successful applications of the FBDD approach, kinases represent a class of targets where this strategy has demonstrated its real potential with the approved kinase inhibitor Vemurafenib. In the Kinase family, protein kinase CK1 isoform δ (CK1δ) has become a promising target in the treatment of different neurodegenerative diseases such as Alzheimer’s disease, Parkinson’s disease, and amyotrophic lateral sclerosis. In the present work, we set up and applied a computational workflow for the identification of putative fragment binders in large virtual databases. To validate the method, the selected compounds were tested in vitro to assess the CK1δ inhibition.

The number of entries in the Protein Data Bank (PDB) has doubled in the last decade, and it has increased tenfold in the last twenty years. The availability of an ever-growing number of structures is having a huge impact on the Structure-Based Drug Discovery (SBDD), allowing investigation of new targets and giving the possibility to have multiple structures of the same macromolecule in a complex with different ligands. Such a large resource often implies the choice of the most suitable complex for molecular docking calculation, and this task is complicated by the plethora of possible posing and scoring function algorithms available, which may influence the quality of the outcomes. Here, we report a large benchmark performed on the PDBbind database containing more than four thousand entries and seventeen popular docking protocols. We found that, even in protein families wherein docking protocols generally showed acceptable results, certain ligand-protein complexes are poorly reproduced in the self-docking procedure. Such a trend in certain protein families is more pronounced, and this underlines the importance in identification of a suitable protein–ligand conformation coupled to a well-performing docking protocol.

One of the most intriguing findings highlighted from G protein-coupled receptor (GPCR) crystallography is the presence, in many members of class A, of a partially hydrated sodium ion in the middle of the seven transmembrane helices (7TM) bundle. In particular, the human adenosine A_2A receptor (A_2A AR) is the first GPCR in which a monovalent sodium ion was crystallized in a distal site from the canonical orthosteric one, corroborating, from a structural point of view, its role as a negative allosteric modulator. However, the molecular mechanism by which the sodium ion influences the recognition of the A_2A AR agonists is not yet fully understood. In this study, the supervised molecular dynamics (SuMD) technique was exploited to analyse the sodium ion recognition mechanism and how its presence influences the binding of the endogenous agonist adenosine. Due to a higher degree of flexibility of the receptor extracellular (EC) vestibule, we propose the sodium-bound A_2A AR as less efficient in stabilizing the adenosine during the different steps of binding.

Fragment-Based Drug Discovery (FBDD) approaches have gained popularity not only in industry but also in academic research institutes. However, the computational prediction of the binding mode adopted by fragment-like molecules within a protein binding site is still a very challenging task. One of the most crucial aspects of fragment binding is related to the large amounts of bound waters in the targeted binding pocket. The binding affinity of fragments may not be sufficient to displace the bound water molecules. In the present work, we confirmed the importance of the bound water molecules in the correct prediction of the fragment binding mode. Moreover, we investigate whether the use of methods based on explicit solvent molecular dynamics simulations can improve the accuracy of fragment posing. The protein chosen for this study is HSP-90.