How To Design Selective Ligands for Highly Conserved Binding Sites: A Case Study Using N-Myristoyltransferases as a Model System

Exploring Subsite Selectivity within Plasmodium vivax N-Myristoyltransferase Using Pyrazole-Derived Inhibitors

N-myristoyltransferase (NMT) is a promising antimalarial drug target. Despite biochemical similarities between Plasmodium vivax and human NMTs, our recent research demonstrated that high selectivity is achievable. Herein, we report PvNMT-inhibiting compounds aimed at identifying novel mechanisms of selectivity. Various functional groups are appended to a pyrazole moiety in the inhibitor to target a pocket formed beneath the peptide binding cleft. The inhibitor core group polarity, lipophilicity, and size are also varied to probe the water structure near a channel. Selectivity index values range from 0.8 to 125.3. Cocrystal structures of two selective compounds, determined at 1.97 and 2.43 Å, show that extensions bind the targeted pocket but with different stabilities. A bulky naphthalene moiety introduced into the core binds next to instead of displacing protein-bound waters, causing a shift in the inhibitor position and expanding the binding site. Our structure-activity data provide a conceptual foundation for guiding future inhibitor optimizations.

Experimental structure based drug design (SBDD) applications for anti-leishmanial drugs: A paradigm shift?

Leishmaniasis is a group of neglected tropical diseases caused by at least 20 species of Leishmania protozoa, which are spread by the bite of infected sandflies. There are three main forms of the disease: cutaneous leishmaniasis (CL, the most common), visceral leishmaniasis (VL, also known as kala-azar, the most serious), and mucocutaneous leishmaniasis. One billion people live in areas endemic to leishmaniasis, with an annual estimation of 30,000 new cases of VL and more than 1 million of CL. New treatments for leishmaniasis are an urgent need, as the existing ones are inefficient, toxic, and/or expensive. We have revised the experimental structure-based drug design (SBDD) efforts applied to the discovery of new drugs against leishmaniasis. We have grouped the explored targets according to the metabolic pathways they belong to, and the key achieved advances are highlighted and evaluated. In most cases, SBDD studies follow high-throughput screening campaigns and are secondary to pharmacokinetic optimization, due to the majoritarian belief that there are few validated targets for SBDD in leishmaniasis. However, some SBDD strategies have significantly contributed to new drug candidates against leishmaniasis and a bigger number holds promise for future development.

Identification of and Structural Insights into Hit Compounds Targeting N-Myristoyltransferase for Cryptosporidium Drug Development

Each year, approximately 50,000 children under 5 die as a result of diarrhea caused by Cryptosporidium parvum, a protozoan parasite. There are currently no effective drugs or vaccines available to cure or prevent Cryptosporidium infection, and there are limited tools for identifying and validating targets for drug or vaccine development. We previously reported a high throughput screening (HTS) of a large compound library against Plasmodium N-myristoyltransferase (NMT), a validated drug target in multiple protozoan parasite species. To identify molecules that could be effective against Cryptosporidium, we counter-screened hits from the Plasmodium NMT HTS against Cryptosporidium NMT. We identified two potential hit compounds and validated them against CpNMT to determine if NMT might be an attractive drug target also for Cryptosporidium. We tested the compounds against Cryptosporidium using both cell-based and NMT enzymatic assays. We then determined the crystal structure of CpNMT bound to Myristoyl-Coenzyme A (MyrCoA) and structures of ternary complexes with MyrCoA and the hit compounds to identify the ligand binding modes. The binding site architectures display different conformational states in the presence of the two inhibitors and provide a basis for rational design of selective inhibitors.

Identification of potent and selective N-myristoyltransferase inhibitors of Plasmodium vivax liver stage hypnozoites and schizonts

Drugs targeting multiple stages of the Plasmodium vivax life cycle are needed to reduce the health and economic burdens caused by malaria worldwide. N-myristoyltransferase (NMT) is an essential eukaryotic enzyme and a validated drug target for combating malaria. However, previous PvNMT inhibitors have failed due to their low selectivity over human NMTs. Herein, we apply a structure-guided hybridization approach combining chemical moieties of previously reported NMT inhibitors to develop the next generation of PvNMT inhibitors. A high-resolution crystal structure of PvNMT bound to a representative selective hybrid compound reveals a unique binding site architecture that includes a selective conformation of a key tyrosine residue. The hybridized compounds significantly decrease P. falciparum blood-stage parasite load and consistently exhibit dose-dependent inhibition of P. vivax liver stage schizonts and hypnozoites. Our data demonstrate that hybridized NMT inhibitors can be multistage antimalarials, targeting dormant and developing forms of liver and blood stage.

Molecular Dynamics Simulations of the Human Ecto-5'-Nucleotidase (h-ecto-5'-NT, CD73): Insights into Protein Flexibility and Binding Site Dynamics

Human ecto-5'-nucleotidase (h-ecto-5'-NT, CD73) is a homodimeric Zn2+-binding metallophosphoesterase that hydrolyzes adenosine 5'-monophosphate (5'-AMP) to adenosine and phosphate. h-Ecto-5'-NT is a key enzyme in purinergic signaling pathways and has been recognized as a promising biological target for several diseases, including cancer and inflammatory, infectious, and autoimmune diseases. Despite its importance as a biological target, little is known about h-ecto-5'-NT dynamics, which poses a considerable challenge to the design of inhibitors of this target enzyme. Here, to explore h-ecto-5'-NT flexibility, all-atom unbiased molecular dynamics (MD) simulations were performed. Remarkable differences in the dynamics of the open (catalytically inactive) and closed (catalytically active) conformations of the apo-h-ecto-5'-NT were observed during the simulations, and the nucleotide analogue inhibitor AMPCP was shown to stabilize the protein structure in the closed conformation. Our results suggest that the large and complex domain motion that enables the h-ecto-5'-NT open/closed conformational switch is slow, and therefore, it could not be completely captured within the time scale of our simulations. Nonetheless, we were able to explore the faster dynamics of the h-ecto-5'-NT substrate binding site, which is mainly located at the C-terminal domain and well conserved among the protein's open and closed conformations. Using the TRAPP ("Transient Pockets in Proteins") approach, we identified transient subpockets close to the substrate binding site. Finally, conformational states of the substrate binding site with higher druggability scores than the crystal structure were identified. In summary, our study provides valuable insights into h-ecto-5'-NT structural flexibility, which can guide the structure-based design of novel h-ecto-5'-NT inhibitors.

Inhibitor Trapping in N-Myristoyltransferases as a Mechanism for Drug Potency

Predicting inhibitor potency is critical in drug design and development, yet it has remained one of computational biology's biggest unresolved challenges. Here, we show that in the case of the N-myristoyltransferase (NMT), this problem could be traced to the mechanisms by which the NMT enzyme is inhibited. NMT adopts open or closed conformations necessary for orchestrating the different steps of the catalytic process. The results indicate that the potency of the NMT inhibitors is determined by their ability to stabilize the enzyme conformation in the closed state, and that in this state, the small molecules themselves are trapped and locked inside the structure of the enzyme, creating a significant barrier for their dissociation. By using molecular dynamics simulations, we demonstrate that the conformational stabilization of the protein molecule in its closed form is highly correlated with the ligands activity and can be used to predict their potency. Hence, predicting inhibitor potency in silico might depend on modeling the conformational changes of the protein molecule upon binding of the ligand rather than estimating the changes in free binding energy that arise from their interaction.

Exploring N-myristoyltransferase as a promising drug target against parasitic neglected tropical diseases

Neglected tropical diseases (NTDs) constitute a group of approximately 20 infectious diseases that mainly affect the impoverished population without basic sanitation in tropical countries. These diseases are responsible for many deaths worldwide, costing billions of dollars in public health investment to treat and control these infections. Among them are the diseases caused by protozoa of the Trypanosomatid family, which constitute Trypanosoma cruzi (Chagas disease), Trypanosoma brucei (sleeping sickness), and Leishmaniasis. In addition, there is a classification of other diseases, called the big three, AIDS, tuberculosis, and malaria, which are endemic in countries with tropical conditions. Despite the high mortality rates, there is still a gap in the treatment. The drugs have a high incidence of side effects and protozoan resistance, justifying the investment in developing new alternatives. In fact, the Target-Based Drug Design (TBDD) approach is responsible for identifying several promising compounds, and among the targets explored through this approach, N-myristoyltransferase (NMT) stands out. It is an enzyme related to the co-translational myristoylation of N-terminal glycine in various peptides. The myristoylation process is a co-translation that occurs after removing the initiator methionine. This process regulates the assembly of protein complexes and stability, which justifies its potential as a drug target. In order to propose NMT as a potential target for parasitic diseases, this review will address the entire structure and function of this enzyme and the primary studies demonstrating its promising potential against Leishmaniasis, T. cruzi, T. brucei, and malaria. We hope our information can help researchers worldwide search for potential drugs against these diseases that have been threatening the health of the world's population.

Protein-Based Virtual Screening Tools Applied for RNA-Ligand Docking Identify New Binders of the preQ1-Riboswitch

Targeting RNA with small molecules is an emerging field. While several ligands for different RNA targets are reported, structure-based virtual screenings (VSs) against RNAs are still rare. Here, we elucidated the general capabilities of protein-based docking programs to reproduce native binding modes of small-molecule RNA ligands and to discriminate known binders from decoys by the scoring function. The programs were found to perform similar compared to the RNA-based docking tool rDOCK, and the challenges faced during docking, namely, protomer and tautomer selection, target dynamics, and explicit solvent, do not largely differ from challenges in conventional protein-ligand docking. A prospective VS with the Bacillus subtilis preQ₁-riboswitch aptamer domain performed with FRED, HYBRID, and FlexX followed by microscale thermophoresis assays identified six active compounds out of 23 tested VS hits with potencies between 29.5 nM and 11.0 μM. The hits were selected not solely based on their docking score but for resembling key interactions of the native ligand. Therefore, this study demonstrates the general feasibility to perform structure-based VSs against RNA targets, while at the same time it highlights pitfalls and their potential solutions when executing RNA-ligand docking.

QSAR and Pharmacophore Modeling of Nitrogen Heterocycles as Potent Human N-Myristoyltransferase (Hs-NMT) Inhibitors

N-myristoyltransferase (NMT) is an important eukaryotic monomeric enzyme which has emerged as an attractive target for developing a drug for cancer, leishmaniasis, ischemia-reperfusion injury, malaria, inflammation, etc. In the present work, statistically robust machine leaning models (QSAR (Quantitative Structure–Activity Relationship) approach) for Human NMT (Hs-NMT) inhibitory has been performed for a dataset of 309 Nitrogen heterocycles screened for NMT inhibitory activity. Hundreds of QSAR models were derived. Of these, the model 1 and 2 were chosen as they not only fulfil the recommended values for a good number of validation parameters (e.g., R² = 0.77–0.79, Q²_LMO = 0.75–0.76, CCC_ex = 0.86–0.87, Q²-F³ = 0.74–0.76, etc.) but also provide useful insights into the structural features that sway the Hs-NMT inhibitory activity of Nitrogen heterocycles. That is, they have an acceptable equipoise of descriptive and predictive qualities as per Organisation for Economic Co-operation and Development (OECD) guidelines. The developed QSAR models identified a good number of molecular descriptors like solvent accessible surface area of all atoms having specific partial charge, absolute surface area of Carbon atoms, etc. as important features to be considered in future optimizations. In addition, pharmacophore modeling has been performed to get additional insight into the pharmacophoric features, which provided additional results.