A key prerequisite for the successful application of protein crystallography in drug discovery is to establish a robust crystallization system for a new drug-target protein fast enough to deliver crystal structures when the first inhibitors have been identified in the hit-finding campaign or, at the latest, in the subsequent hit-to-lead process. The first crucial step towards generating well folded proteins with a high likelihood of crystallizing is the identification of suitable truncation variants of the target protein. In some cases an optimal length variant alone is not sufficient to support crystallization and additional surface mutations need to be introduced to obtain suitable crystals. In this contribution, four case studies are presented in which rationally designed surface modifications were key to establishing crystallization conditions for the target proteins (the protein kinases Aurora-C, IRAK4 and BUB1, and the KRAS-SOS1 complex). The design process which led to well diffracting crystals is described and the crystal packing is analysed to understand retrospectively how the specific surface mutations promoted successful crystallization. The presented design approaches are routinely used in our team to support the establishment of robust crystallization systems which enable structure-guided inhibitor optimization for hit-to-lead and lead-optimization projects in pharmaceutical research.

In search of new opportunities to develop Trypanosoma brucei phosphodiesterase B1 (TbrPDEB1) inhibitors that have selectivity over the off-target human PDE4 (hPDE4), different stages of a fragment-growing campaign were studied using a variety of biochemical, structural, thermodynamic, and kinetic binding assays. Remarkable differences in binding kinetics were identified and this kinetic selectivity was explored with computational methods, including molecular dynamics and interaction fingerprint analyses. These studies indicate that a key hydrogen bond between GlnQ.50 and the inhibitors is exposed to a water channel in TbrPDEB1, leading to fast unbinding. This water channel is not present in hPDE4, leading to inhibitors with a longer residence time. The computer-aided drug design protocols were applied to a recently disclosed TbrPDEB1 inhibitor with a different scaffold and our results confirm that shielding this key hydrogen bond through disruption of the water channel represents a viable design strategy to develop more selective inhibitors of TbrPDEB1. Our work shows how computational protocols can be used to understand the contribution of solvent dynamics to inhibitor binding, and our results can be applied in the design of selective inhibitors for homologous PDEs found in related parasites.

Plasmodium falciparum aminoacyl tRNA synthetases (PfaaRSs) are potent antimalarial targets essential for proteome fidelity and overall parasite survival in every stage of the parasite\'s life cycle. So far, some of these proteins have been singly targeted yielding inhibitor compounds that have been limited by incidences of resistance which can be overcome via pan-inhibition strategies. Hence, herein, for the first time, we report the identification and in vitro antiplasmodial validation of Mitomycin (MMC) as a probable pan-inhibitor of class 1a (arginyl(A)-, cysteinyl(C), isoleucyl(I)-, leucyl(L), methionyl(M), and valyl(V)-) PfaaRSs which hypothetically may underlie its previously reported activity on the ribosomal RNA to inhibit protein translation and biosynthesis. We combined multiple in silico structure-based discovery strategies that first helped identify functional and druggable sites that were preferentially targeted by the compound in each of the plasmodial proteins: Ins1-Ins2 domain in Pf-ARS; anticodon binding domain in Pf-CRS; CP1-editing domain in Pf-IRS and Pf-MRS; C-terminal domain in Pf-LRS; and CP-core region in Pf-VRS. Molecular dynamics studies further revealed that MMC allosterically induced changes in the global structures of each protein. Likewise, prominent structural perturbations were caused by the compound across the functional domains of the proteins. More so, MMC induced systematic alterations in the binding of the catalytic nucleotide and amino acid substrates which culminated in the loss of key interactions with key active site residues and ultimate reduction in the nucleotide-binding affinities across all proteins, as deduced from the binding energy calculations. These altogether confirmed that MMC uniformly disrupted the structure of the target proteins and essential substrates. Further, MMC demonstrated IC50 < 5 μM against the Dd2 and 3D7 strains of parasite making it a good starting point for malarial drug development. We believe that findings from our study will be important in the current search for highly effective multi-stage antimalarial drugs.

Over the years, human dihydroorotate dehydrogenase (hDHODH), which is a key player in the de novo pyrimidine-biosynthesis pathway, has been targeted in the treatment of several conditions, including autoimmune disorders and acute myelogenous leukaemia, as well as in host-targeted antiviral therapy. A molecular exploration of its inhibitor-binding behaviours yielded promising candidates for innovative drug design. A detailed description of the enzymatic pharmacophore drove the decoration of well-established inhibitory scaffolds, thus gaining further in vitro and in vivo efficacy. In the present work, using X-ray crystallography, an atypical rearrangement was identified in the binding pose of a potent inhibitor characterized by a polar pyridine-based moiety (compound 18). The crystal structure shows that upon binding compound 18 the dynamics of a protein loop involved in a gating mechanism at the cofactor-binding site is modulated by the presence of three water molecules, thus fine-tuning the polarity/hydrophobicity of the binding pocket. These solvent molecules are engaged in the formation of a hydrogen-bond mesh in which one of them establishes a direct contact with the pyridine moiety of compound 18, thus paving the way for a reappraisal of the inhibition of hDHODH. Using an integrated approach, the thermodynamics of such a modulation is described by means of isothermal titration calorimetry coupled with molecular modelling. These structural insights will guide future drug design to obtain a finer Kd/logD7.4 balance and identify membrane-permeable molecules with a drug-like profile in terms of water solubility.

TIGIT is mainly expressed on T cells and is an inhibitory checkpoint receptor that binds to its ligand PVR in the tumor microenvironment. Anti-TIGIT monoclonal antibodies (mAbs) such as Ociperlimab and Tiragolumab block the TIGIT-PVR interaction and are in clinical development. However, the molecular blockade mechanism of these mAbs remains elusive. Here, we report the crystal structures of TIGIT in complex with Ociperlimab_Fab and Tiragolumab_Fab revealing that both mAbs bind TIGIT with a large steric clash with PVR. Furthermore, several critical epitopic residues are identified. Interestingly, the binding affinity of Ociperlimab toward TIGIT increases approximately 17-fold when lowering the pH from 7.4 to 6.0. Our structure shows a strong electrostatic interaction between ASP103HCDR3 and HIS76TIGIT explaining the pH-responsive mechanism of Ociperlimab. In contrast, Tiragolumab does not show an acidic pH-dependent binding enhancement. Our results provide valuable information that could help to improve the efficacy of therapeutic antibodies for cancer treatment.

The launch of AlphaFold series has brought deep-learning techniques into the molecular structural science. As another crucial problem, structure-based prediction of protein-ligand binding affinity urgently calls for advanced computational techniques. Is deep learning ready to decode this problem? Here we review mainstream structure-based, deep-learning approaches for this problem, focusing on molecular representations, learning architectures and model interpretability. A model taxonomy has been generated. To compensate for the lack of valid comparisons among those models, we realized and evaluated representatives from a uniform basis, with the advantages and shortcomings discussed. This review will potentially benefit structure-based drug discovery and related areas.

The ability to quickly discover reliable hits from screening and rapidly convert them into lead compounds, which can be verified in functional assays, is central to drug discovery. The expedited validation of novel targets and the identification of modulators to advance to preclinical studies can significantly increase drug development success. Our SaXPyTM (\"SAR by X-ray Poses Quickly\") platform, which is applicable to any X-ray crystallography-enabled drug target, couples the established methods of protein X-ray crystallography and fragment-based drug discovery (FBDD) with advanced computational and medicinal chemistry to deliver small molecule modulators or targeted protein degradation ligands in a short timeframe. Our approach, especially for elusive or \"undruggable\" targets, allows for (i) hit generation; (ii) the mapping of protein-ligand interactions; (iii) the assessment of target ligandability; (iv) the discovery of novel and potential allosteric binding sites; and (v) hit-to-lead execution. These advances inform chemical tractability and downstream biology and generate novel intellectual property. We describe here the application of SaXPy in the discovery and development of DNA damage response inhibitors against DNA polymerase eta (Pol η or POLH) and apurinic/apyrimidinic endonuclease 1 (APE1 or APEX1). Notably, our SaXPy platform allowed us to solve the first crystal structures of these proteins bound to small molecules and to discover novel binding sites for each target.

Leishmaniasis is caused by ∼20 species of Leishmania that affects millions in endemic areas. Available therapies are not sufficient to effectively control the disease, cause severe side effects and eventually lead to drug resistance, making the discovery of novel therapeutic molecules an immediate need. Molecular target-based drug discovery, where the target is a defined molecular gene, protein or a mechanism, is a rationale driven approach for novel therapeutics. Humans obtain the essential amino acid such as threonine from dietary sources, while Leishmania synthesize it de-novo. Enzymes of the threonine biosynthesis pathway, including the rate limiting Homoserine kinase (HSK) which converts L-homoserine into ortho-phospho homoserine are thus attractive targets for rationale driven therapy. The absence of HSK in humans and its presence in Leishmania donovani enhances the opportunity to exploit HSK as a molecular target for anti-leishmanials therapeutic development. In this study, we utilize structure-based high throughput drug discovery (SBDD), followed by biochemical validation and identified two potential inhibitors (RH00038 and S02587) from Maybridge chemical library that targets L. donovani HSK. These two inhibitors effectively induced the mortality of Leishmania donovani in both amastigote and promastigote stages, with one of them being specific to parasite and twice as effective as the standard therapeutic molecule.

Nowadays, it is widely accepted that water molecules play a key role in binding a ligand to a molecular target. Neglecting water molecules in the process of molecular recognition was the result of several failures of the structure-based drug discovery campaigns. The application of WaterMap, in particular WaterMap-guided molecular docking, enables the reasonably accurate and quick description of the location and energetics of water molecules at the ligand-protein interface.
In this review, the authors shortly discuss the importance of water in drug design and discovery and provide a brief overview of the computational approaches used to predict the solvent-related effects for the purposes of presenting WaterMap in the context of other available techniques and tools. A concise description of WaterMap concept is followed by the presentation of WaterMap-assisted virtual screening literature published between 2013 and 2023.
In recent years, WaterMap software has been extensively used to support structure-based drug design, in particular structure-based virtual screening. Indeed, it is a useful tool to rescore docking results considering water molecules in the binding pocket. Although WaterMap allows for the consideration of the dynamic behavior of water molecules in the binding site, for best accuracy, its application in conjunction with other techniques such as molecular mechanics-generalized Born surface area of FEP (Free Energy Perturbation) is recommended.

The design of accurate virtual screening tools is an open challenge in drug discovery. Several structure-based methods have been developed at different levels of approximation. Among them, molecular docking is an established technique with high efficiency, but typically low accuracy. Moreover, docking performances are known to be target-dependent, which makes the choice of the docking program and corresponding scoring function critical when approaching a new protein target. To compare the performances of different docking protocols, we developed ChemFlow_py, an automated tool to perform docking and rescoring. Using four protein systems extracted from DUD-E with 100 known active compounds and 3000 decoys per target, we compared the performances of several rescoring strategies including consensus scoring. We found that the average docking results can be improved by consensus ranking, which emphasizes the relevance of consensus scoring when little or no chemical information is available for a given target. ChemFlow_py is a free toolkit to optimize the performances of virtual high-throughput screening (vHTS). The software is publicly available at https://github.com/IFMlab/ChemFlow_py .