Ligand 1 H NMR Chemical Shifts as Accurate Reporters for Protein-Ligand Binding Interfaces in Solution.

The availability of high-resolution 3D structural information is crucial for investigating guest-host systems across a wide range of fields. In the context of drug discovery, the information is routinely used to establish and validate structure-activity relationships, grow initial hits from screening campaigns, and to guide molecular docking. For the generation of protein-ligand complex structural information, X-ray crystallography is the experimental method of choice, however, with limited information on protein flexibility. An experimentally verified structural model of the binding interface in the native solution-state would support medicinal chemists in their molecular design decisions. Here we demonstrate that protein-bound ligand 1 H NMR chemical shifts are highly sensitive and accurate probes for the immediate chemical environment of protein-ligand interfaces. By comparing the experimental ligand 1 H chemical shift values with those computed from the X-ray structure using quantum mechanics methodology, we identify significant disagreements for parts of the ligand between the two experimental techniques. We show that quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) ensembles can be used to refine initial X-ray co-crystal structures resulting in a better agreement with experimental 1 H ligand chemical shift values. Overall, our findings highlight the usefulness of ligand 1 H NMR chemical shift information in combination with a QM/MM MD workflow for generating protein-ligand ensembles that accurately reproduce solution structural data.

Ligandability assessment of the C-terminal Rel-homology domain of NFAT1.

Transcription factors are generally considered challenging, if not "undruggable", targets but they promise new therapeutic options due to their fundamental involvement in many diseases. In this study, we aim to assess the ligandability of the C-terminal Rel-homology domain of nuclear factor of activated T cells 1 (NFAT1), a TF implicated in T-cell regulation. Using a combination of experimental and computational approaches, we demonstrate that small molecule fragments can indeed bind to this protein domain. The newly identified binder is the first small molecule binder to NFAT1 validated with biophysical methods and an elucidated binding mode by X-ray crystallography. The reported eutomer/distomer pair provides a strong basis for potential exploration of higher potency binders on the path toward degrader or glue modalities.

Fragment-based discovery of a chemical probe for the PWWP1 domain of NSD3.

Here, we report the fragment-based discovery of BI-9321, a potent, selective and cellular active antagonist of the NSD3-PWWP1 domain. The human NSD3 protein is encoded by the WHSC1L1 gene located in the 8p11-p12 amplicon, frequently amplified in breast and squamous lung cancer. Recently, it was demonstrated that the PWWP1 domain of NSD3 is required for the viability of acute myeloid leukemia cells. To further elucidate the relevance of NSD3 in cancer biology, we developed a chemical probe, BI-9321, targeting the methyl-lysine binding site of the PWWP1 domain with sub-micromolar in vitro activity and cellular target engagement at 1 µM. As a single agent, BI-9321 downregulates Myc messenger RNA expression and reduces proliferation in MOLM-13 cells. This first-in-class chemical probe BI-9321, together with the negative control BI-9466, will greatly facilitate the elucidation of the underexplored biological function of PWWP domains.

Discovery of a small-molecule inhibitor of specific serine residue BAD phosphorylation.

Human BCL-2-associated death promoter (hBAD) is an apoptosis-regulatory protein mediating survival signals to carcinoma cells upon phosphorylation of Ser99, among other residues. Herein, we screened multiple small-molecule databases queried in a Laplacian-modified naive Bayesian-based cheminformatics platform and identified a Petasis reaction product as a site-specific inhibitor for hBAD phosphorylation. Based on apoptotic efficacy against mammary carcinoma cells, N-cyclopentyl-3-((4-(2,3-dichlorophenyl) piperazin-1-yl) (2-hydroxyphenyl) methyl) benzamide (NPB) was identified as a potential lead compound. In vitro biochemical analyses demonstrated that NPB inhibited the phosphorylation of hBAD specifically on Ser99. NPB was observed to exert this effect independently of AKT and other kinase activities despite the demonstration of AKT-mediated BAD-Ser99 phosphorylation. Using a structure-based bioinformatics platform, we observed that NPB exhibited predicted interactions with hBAD in silico and verified the same by direct binding kinetics. NPB reduced phosphorylation of BAD-Ser99 and enhanced caspase 3/7 activity with associated loss of cell viability in various human cancer cell lines derived from mammary, endometrial, ovarian, hepatocellular, colon, prostatic, and pancreatic carcinoma. Furthermore, by use of a xenograft model, it was observed that NPB, as a single agent, markedly diminished BAD phosphorylation in tumor tissue and significantly inhibited tumor growth. Similar doses of NPB utilized in acute toxicity studies in mice did not exhibit significant effects. Hence, we report a site-specific inhibitor of BAD phosphorylation with efficacy in tumor models.

PI by NMR: Probing CH-π Interactions in Protein-Ligand Complexes by NMR Spectroscopy.

While CH-π interactions with target proteins are crucial determinants for the affinity of arguably every drug molecule, no method exists to directly measure the strength of individual CH-π interactions in drug-protein complexes. Herein, we present a fast and reliable methodology called PI (π interactions) by NMR, which can differentiate the strength of protein-ligand CH-π interactions in solution. By combining selective amino-acid side-chain labeling with 1 H-13 C NMR, we are able to identify specific protein protons of side-chains engaged in CH-π interactions with aromatic ring systems of a ligand, based solely on 1 H chemical-shift values of the interacting protein aromatic ring protons. The information encoded in the chemical shifts induced by such interactions serves as a proxy for the strength of each individual CH-π interaction. PI by NMR changes the paradigm by which chemists can optimize the potency of drug candidates: direct determination of individual π interactions rather than averaged measures of all interactions.

Enhanced vulnerability of human proteins towards disease-associated inactivation through divergent evolution.

Human proteins are vulnerable towards disease-associated single amino acid replacements affecting protein stability and function. Interestingly, a few studies have shown that consensus amino acids from mammals or vertebrates can enhance protein stability when incorporated into human proteins. Here, we investigate yet unexplored relationships between the high vulnerability of human proteins towards disease-associated inactivation and recent evolutionary site-specific divergence of stabilizing amino acids. Using phylogenetic, structural and experimental analyses, we show that divergence from the consensus amino acids at several sites during mammalian evolution has caused local protein destabilization in two human proteins linked to disease: cancer-associated NQO1 and alanine:glyoxylate aminotransferase, mutated in primary hyperoxaluria type I. We demonstrate that a single consensus mutation (H80R) acts as a disease suppressor on the most common cancer-associated polymorphism in NQO1 (P187S). The H80R mutation reactivates P187S by enhancing FAD binding affinity through local and dynamic stabilization of its binding site. Furthermore, we show how a second suppressor mutation (E247Q) cooperates with H80R in protecting the P187S polymorphism towards inactivation through long-range allosteric communication within the structural ensemble of the protein. Our results support that recent divergence of consensus amino acids may have occurred with neutral effects on many functional and regulatory traits of wild-type human proteins. However, divergence at certain sites may have increased the propensity of some human proteins towards inactivation due to disease-associated mutations and polymorphisms. Consensus mutations also emerge as a potential strategy to identify structural hot-spots in proteins as targets for pharmacological rescue in loss-of-function genetic diseases.

Understanding Structure-Activity Relationships for Trypanosomal Cysteine Protease Inhibitors by Simulations and Free Energy Calculations.

The protozoan cysteine proteases cruzain in Trypanosoma cruzi and rhodesain in Trypanosoma brucei are therapeutic targets for Chagas disease and Human African Trypanosomiasis (HAT), respectively. A benzimidazole series was previously characterized as potent noncovalent competitive cruzain and rhodesain inhibitors with activity against trypanosomes. Common structure-activity relationships (SAR) trends and structural modifications leading to selectivity against each enzyme were described. However, some of these trends could not be understood based on the reported binding mode of lead compound 1. Therefore, we employed microsecond molecular dynamics simulations and free energy calculations to understand qualitative SAR trends and to quantitatively recapitulate them. Simulations revealed the most stable protein-ligand interactions and provided insights concerning enzyme selectivity. Calculated relative binding free energies of compound 1 analogs exhibited deviations of 1.1 and 2.2 kcal/mol from the experimental values for cruzain and rhodesain, respectively. These data encourage prospective thermodynamic integration (TI) studies to optimize this series and facilitate the prioritization of compounds for synthesis.

Asymmetric arginine dimethylation of RelA provides a repressive mark to modulate TNFα/NF-κB response.

Nuclear factor kappa B (NF-κB) is an inducible transcription factor that plays critical roles in immune and stress responses and is often implicated in pathologies, including chronic inflammation and cancer. Although much has been learned about NF-κB-activating pathways, the specific repression of NF-κB is far less well understood. Here we identified the type I protein arginine methyltransferase 1 (PRMT1) as a restrictive factor controlling TNFα-induced activation of NF-κB. PRMT1 forms a cellular complex with NF-κB through direct interaction with the Rel homology domain of RelA. We demonstrate that PRMT1 methylates RelA at evolutionary conserved R30, located in the DNA-binding L1 loop, which is a critical residue required for DNA binding. Asymmetric R30 dimethylation inhibits the binding of RelA to DNA and represses NF-κB target genes in response to TNFα. Molecular dynamics simulations of the DNA-bound RelA:p50 predicted structural changes in RelA caused by R30 methylation or a mutation that interferes with the stability of the DNA-NF-κB complex. Our findings provide evidence for the asymmetric arginine dimethylation of RelA and unveil a unique mechanism controlling TNFα/NF-κB signaling.

Sulfated Ceria Catalyzed Synthesis of Imidazopyridines and Their Implementation as DNA Minor Groove Binders.

The small molecules that bind to DNA minor groove are considered as potential therapeutic agents to fight against many human diseases. They induce cell death by interfering with transcription, replication and progression of cell cycle. Herein, we report the synthesis of imidazopyridine-3-amines using sulfated ceria catalyst by employing Groebkee-Blackburne-Bienayme reaction. We evaluated the possible antiproliferative and antimycobacterial activity against A549 cells and Mycobacterium tuberculosis, respectively. Among the tested compounds, N-tert-butyl-2-(2-butyl-4-chloro-1H-imidazol-5-yl)-5,7-dimethylimidazo[1,2-a]pyridin-3-amine (4g) was identified as cytotoxic heterocycle and antimycobacterial agent. Molecular docking studies of the imidazopyridine derivatives revealed the consistent positioning in the minor groove with a tight shape fit between receptor and ligands. Therefore, we speculate that new imidazopyridines induce their pharmacological effect by targeting the minor groove of DNA.

Electrostatic recognition in substrate binding to serine proteases.

Serine proteases of the Chymotrypsin family are structurally very similar but have very different substrate preferences. This study investigates a set of 9 different proteases of this family comprising proteases that prefer substrates containing positively charged amino acids, negatively charged amino acids, and uncharged amino acids with varying degree of specificity. Here, we show that differences in electrostatic substrate preferences can be predicted reliably by electrostatic molecular interaction fields employing customized GRID probes. Thus, we are able to directly link protease structures to their electrostatic substrate preferences. Additionally, we present a new metric that measures similarities in substrate preferences focusing only on electrostatics. It efficiently compares these electrostatic substrate preferences between different proteases. This new metric can be interpreted as the electrostatic part of our previously developed substrate similarity metric. Consequently, we suggest, that substrate recognition in terms of electrostatics and shape complementarity are rather orthogonal aspects of substrate recognition. This is in line with a 2-step mechanism of protein-protein recognition suggested in the literature.

Peptidic Macrocycles - Conformational Sampling and Thermodynamic Characterization.

Macrocycles are of considerable interest as highly specific drug candidates, yet they challenge standard conformer generators with their large number of rotatable bonds and conformational restrictions. Here, we present a molecular dynamics-based routine that bypasses current limitations in conformational sampling and extensively profiles the free energy landscape of peptidic macrocycles in solution. We perform accelerated molecular dynamics simulations to capture a diverse conformational ensemble. By applying an energetic cutoff, followed by geometric clustering, we demonstrate the striking robustness and efficiency of the approach in identifying highly populated conformational states of cyclic peptides. The resulting structural and thermodynamic information is benchmarked against interproton distances from NMR experiments and conformational states identified by X-ray crystallography. Using three different model systems of varying size and flexibility, we show that the method reliably reproduces experimentally determined structural ensembles and is capable of identifying key conformational states that include the bioactive conformation. Thus, the described approach is a robust method to generate conformations of peptidic macrocycles and holds promise for structure-based drug design.

Mechanisms Responsible for ω-Pore Currents in Cav Calcium Channel Voltage-Sensing Domains.

Mutations of positively charged amino acids in the S4 transmembrane segment of a voltage-gated ion channel form ion-conducting pathways through the voltage-sensing domain, named ω-current. Here, we used structure modeling and MD simulations to predict pathogenic ω-currents in CaV1.1 and CaV1.3 Ca2+ channels bearing several S4 charge mutations. Our modeling predicts that mutations of CaV1.1-R1 (R528H/G, R897S) or CaV1.1-R2 (R900S, R1239H) linked to hypokalemic periodic paralysis type 1 and of CaV1.3-R3 (R990H) identified in aldosterone-producing adenomas conducts ω-currents in resting state, but not during voltage-sensing domain activation. The mechanism responsible for the ω-current and its amplitude depend on the number of charges in S4, the position of the mutated S4 charge and countercharges, and the nature of the replacing amino acid. Functional characterization validates the modeling prediction showing that CaV1.3-R990H channels conduct ω-currents at hyperpolarizing potentials, but not upon membrane depolarization compared with wild-type channels.

DPPC-cholesterol phase diagram using coarse-grained Molecular Dynamics simulations.

Cholesterol-phospholipid bilayers continue to be the current state of the art in membrane models and serve as representative systems for studying the effect of cholesterol on the cell membrane. As the mixing of different lipid species requires long spatio-temporal scales, coarse-grained models have gained increasing popularity in modeling such membrane systems. In this paper, a systematic study of the MARTINI coarse-grained model for the DPPC-cholesterol binary system has been performed. We construct the phase diagram of DPPC lipid bilayers in the presence of different cholesterol concentrations and at different temperatures using coarse-grained Molecular Dynamics (MD) simulations with the MARTINI force field. The phase diagram based on the condensation effect is directly comparable to available experimental data and demonstrates qualitative agreement over all cholesterol concentrations. Self-assembled bilayers quantitatively reproduce experimental observables, such as lateral diffusion of lipids, electron density, area per lipid and lipid order parameters. The phase diagram of the DPPC-cholesterol binary system also reveals the profound effect of cholesterol on the physical properties of phospholipid bilayers such lipid order, diffusion, and fluidity. Cholesterol induces the liquid-ordered phase, which increases the fluidity of the phospholipid hydrocarbon chains above the gel to liquid-crystalline phase transition temperature and decreases it below the phase transition. The present study suggests that the MARTINI force field can be successfully used to obtain molecular level insights into cholesterol-DPPC model membranes.

Identification of Novel Class of Triazolo-Thiadiazoles as Potent Inhibitors of Human Heparanase and their Anticancer Activity.

BACKGROUND: Expression and activity of heparanase, an endoglycosidase that cleaves heparan sulfate (HS) side chains of proteoglycans, is associated with progression and poor prognosis of many cancers which makes it an attractive drug target in cancer therapeutics. METHODS: In the present work, we report the in vitro screening of a library of 150 small molecules with the scaffold bearing quinolones, oxazines, benzoxazines, isoxazoli(di)nes, pyrimidinones, quinolines, benzoxazines, and 4-thiazolidinones, thiadiazolo[3,2-a]pyrimidin-5-one, 1,2,4-triazolo-1,3,4-thiadiazoles, and azaspiranes against the enzymatic activity of human heparanase. The identified lead compounds were evaluated for their heparanase-inhibiting activity using sulfate [35S] labeled extracellular matrix (ECM) deposited by cultured endothelial cells. Further, anti-invasive efficacy of lead compound was evaluated against hepatocellular carcinoma (HepG2) and Lewis lung carcinoma (LLC) cells. RESULTS: Among the 150 compounds screened, we identified 1,2,4-triazolo-1,3,4-thiadiazoles bearing compounds to possess human heparanase inhibitory activity. Further analysis revealed 2,4-Diiodo-6-(3-phenyl-[1, 2, 4]triazolo[3,4-b][1, 3, 4]thiadiazol-6yl)phenol (DTP) as the most potent inhibitor of heparanase enzymatic activity among the tested compounds. The inhibitory efficacy was demonstrated by a colorimetric assay and further validated by measuring the release of radioactive heparan sulfate degradation fragments from [35S] labeled extracellular matrix. Additionally, lead compound significantly suppressed migration and invasion of LLC and HepG2 cells with IC50 value of ~5 µM. Furthermore, molecular docking analysis revealed a favourable interaction of triazolo-thiadiazole backbone with Asn-224 and Asp-62 of the enzyme. CONCLUSIONS: Overall, we identified biologically active heparanase inhibitor which could serve as a lead structure in developing compounds that target heparanase in cancer.

Binding Pose Flip Explained via Enthalpic and Entropic Contributions.

The anomalous binding modes of five highly similar fragments of TIE2 inhibitors, showing three distinct binding poses, are investigated. We report a quantitative rationalization for the changes in binding pose based on molecular dynamics simulations. We investigated five fragments in complex with the transforming growth factor ß receptor type 1 kinase domain. Analyses of these simulations using Grid Inhomogeneous Solvation Theory (GIST), pKA calculations, and a tool to investigate enthalpic differences upon binding unraveled the various thermodynamic contributions to the different binding modes. While one binding mode flip can be rationalized by steric repulsion, the second binding pose flip revealed a different protonation state for one of the ligands, leading to different enthalpic and entropic contributions to the binding free energy. One binding pose is stabilized by the displacement of entropically unfavored water molecules (binding pose determined by solvation entropy), ligands in the other binding pose are stabilized by strong enthalpic interactions, overcompensating the unfavorable water entropy in this pose (binding pose determined by enthalpic interactions). This analysis elucidates unprecedented details determining the flipping of the binding modes, which can elegantly explain the experimental findings for this system.

Fold stability during endolysosomal acidification is a key factor for allergenicity and immunogenicity of the major birch pollen allergen.

BACKGROUND: The search for intrinsic factors, which account for a protein's capability to act as an allergen, is ongoing. Fold stability has been identified as a molecular feature that affects processing and presentation, thereby influencing an antigen's immunologic properties. OBJECTIVE: We assessed how changes in fold stability modulate the immunogenicity and sensitization capacity of the major birch pollen allergen Bet v 1. METHODS: By exploiting an exhaustive virtual mutation screening, we generated mutants of the prototype allergen Bet v 1 with enhanced thermal and chemical stability and rigidity. Structural changes were analyzed by means of x-ray crystallography, nuclear magnetic resonance, and molecular dynamics simulations. Stability was monitored by using differential scanning calorimetry, circular dichroism, and Fourier transform infrared spectroscopy. Endolysosomal degradation was simulated in vitro by using the microsomal fraction of JAWS II cells, followed by liquid chromatography coupled to mass spectrometry. Immunologic properties were characterized in vitro by using a human T-cell line specific for the immunodominant epitope of Bet v 1 and in vivo in an adjuvant-free BALB/c mouse model. RESULTS: Fold stabilization of Bet v 1 was pH dependent and resulted in resistance to endosomal degradation at a pH of 5 or greater, affecting presentation of the immunodominant T-cell epitope in vitro. These properties translated in vivo into a strong allergy-promoting TH2-type immune response. Efficient TH2 cell activation required both an increased stability at the pH of the early endosome and efficient degradation at lower pH in the late endosomal/lysosomal compartment. CONCLUSIONS: Our data indicate that differential pH-dependent fold stability along endosomal maturation is an essential protein-inherent determinant of allergenicity.

Conformational Flexibility Differentiates Naturally Occurring Bet v 1 Isoforms.

The protein Bet v 1 represents the main cause for allergic reactions to birch pollen in Europe and North America. Structurally homologous isoforms of Bet v 1 can have different properties regarding allergic sensitization and Th2 polarization, most likely due to differential susceptibility to proteolytic cleavage. Using NMR relaxation experiments and molecular dynamics simulations, we demonstrate that the initial proteolytic cleavage sites in two naturally occurring Bet v 1 isoforms, Bet v 1.0101 (Bet v 1a) and Bet v 1.0102 (Bet v 1d), are conformationally flexible. Inaccessible cleavage sites in helices and strands are highly flexible on the microsecond-millisecond time scale, whereas those located in loops display faster nanosecond-microsecond flexibility. The data consistently show that Bet v 1.0102 is more flexible and conformationally heterogeneous than Bet v 1.0101. Moreover, NMR hydrogen-deuterium exchange measurements reveal that the backbone amides in Bet v 1.0102 are significantly more solvent exposed, in agreement with this isoform's higher susceptibility to proteolytic cleavage. The differential conformational flexibility of Bet v 1 isoforms, along with the transient exposure of inaccessible sites to the protein surface, may be linked to proteolytic susceptibility, representing a potential structure-based rationale for the observed differences in Th2 polarization and allergic sensitization.

Analysis of Differential Efficacy and Affinity of GABAA (α1/α2) Selective Modulators.

Selective modulators of the γ-amino butyric acid (GABAA) family of receptors have the potential to treat a range of disease states related to cognition, pain, and anxiety. While the development of various α subunit-selective modulators is currently underway for the treatment of anxiety disorders, a mechanistic understanding of the correlation between their bioactivity and efficacy, based on ligand-target interactions, is currently still lacking. In order to alleviate this situation, in the current study we have analyzed, using ligand- and structure-based methods, a data set of 5440 GABAA modulators. The Spearman correlation (ρ) between binding activity and efficacy of compounds was calculated to be 0.008 and 0.31 against the α1 and α2 subunits of GABA receptor, respectively; in other words, the compounds had little diversity in structure and bioactivity, but they differed significantly in efficacy. Two compounds were selected as a case study for detailed interaction analysis due to the small difference in their structures and affinities (ΔpKi(comp1_α1 - comp2_α1) = 0.45 log units, ΔpKi(comp1_α2 - comp2_α2) = 0 log units) as compared to larger relative efficacies (ΔRE(comp1_α1 - comp2_α1) = 1.03, ΔRE(comp1_α2 - comp2_α2) = 0.21). Docking analysis suggested that His-101 is involved in a characteristic interaction of the α1 receptor with both compounds 1 and 2. Residues such as Phe-77, Thr-142, Asn-60, and Arg-144 of the γ chain of the α1γ2 complex also showed interactions with heterocyclic rings of both compounds 1 and 2, but these interactions were disturbed in the case of α2γ2 complex docking results. Binding pocket stability analysis based on molecular dynamics identified three substitutions in the loop C region of the α2 subunit, namely, G200E, I201T, and V202I, causing a reduction in the flexibility of α2 compared to α1. These amino acids in α2, as compared to α1, were also observed to decrease the vibrational and dihedral entropy and to increase the hydrogen bond content in α2 in the apo state. However, freezing of both α1 and α2 was observed in the ligand-bound state, with an increased number of internal hydrogen bonds and increased entropy. Therefore, we hypothesize that the amino acid differences in the loop C region of α2 are responsible for conformational changes in the protein structure compared to α1, as well as for the binding modes of compounds and hence their functional signaling.

Adamantyl-tethered-biphenylic compounds induce apoptosis in cancer cells by targeting Bcl homologs.

Bcl homologs prominently contribute to apoptotic resistance in cancer cells and serve as molecular targets in treatment of various cancers. Herein, we report the synthesis of biphenyl-adamantane derivatives by a ligand free palladium on carbon based Suzuki reaction using diisopropylamine as a base for the coupling of adamantane based aryl chloride with a variety of aryl boronic acids. Among the biphenyl derivatives synthesized, compound 3'-(adamantan-1-yl)-4'-methoxy-[1,1'-biphenyl]-3-ol (AMB) displayed cytotoxic activity against hepatocellular carcinoma cell lines without significantly affecting the normal cell lines. Further, AMB caused increased accumulation of the HCC cells in subG1 phase, decreased the expression of Bcl-2, Bcl-xL, cyclin D1, caspase-3, survivin and increased the cleavage of PARP in a time-dependent manner. In silico molecular interaction studies between Bcl homologs and AMB showed that the biphenyl scaffold is predicted to form π-π interactions with Phe-101 and Tyr-105 and the adamantyl fragment is predicted to occupy another hydrophobic region in the kink region of the binding groove. In summary, we report on the synthesis and biological characterization of adamantyl-tethered biphenylic compounds that induce apoptosis in tumor cells most likely by targeting Bcl homologs.

Synthesis and in vitro evaluation of hydrazinyl phthalazines against malaria parasite, Plasmodium falciparum.

In this report, we describe the synthesis of 1-(Phthalazin-4-yl)-hydrazine using bronsted acidic ionic liquids and demonstrate their ability to inhibit asexual stage development of human malaria parasite, Plasmodium falciparum. Through computational studies, we short-listed chemical scaffolds with potential binding affinity to an essential parasite protein, dihydroorotate dehydrogenase (DHODH). Further, these compounds were synthesized in the lab and tested against P. falciparum. Several compounds from our library showed inhibitory activity at low micro-molar concentrations with minimal cytotoxic effects. These results indicate the potential of hydralazine derivatives as reference scaffolds to develop novel antimalarials.