Activation by substoichiometric inhibition

Molecular Tweezers Block the Functional Pore of a Protein Machine

We present symmetric multivalent tweezers as the first class of supramolecular elements designed to cover and functionally block a protein pore. As a model, we chose the enzyme p97, a hexameric AAA-ATPase that unfolds or segregates substrate proteins by threading them through a pore and channel at the center of the symmetric p97 hexamer fueled by ATP hydrolysis. In a rational design approach, we developed a new class of p97 inhibitors, guided by molecular modeling. These dock onto lysine residues at the entry of the pore via appropriately positioned molecular tweezers. Ligand binding was accompanied by induction of fluorescence of the built-in binding sensitive luminophores which served as a sensor for affinity determination. We further confirmed specific interaction with p97 as well as concomitant inhibition of ATPase activity and protein substrate unfolding using an array of biophysical methods and state-of-the art biochemical assays. Specific binding was also validated by mutagenesis, demonstrating that inhibition of p97 function was mediated by blocking the pore entrance. Especially C3-symmetric multivalent tweezers potently inhibited ATPase activity and protein substrate processing consistent with the symmetry of the docking site. Our data independently confirm substrate threading as a mechanism for protein unfolding by p97 and highlight multivalent tweezers as a supramolecular strategy to target pores in various proteins. Since p97 and related protein machines are vital for protein quality control and cell survival, the new pore binders may open a new approach to combat diseases and be employed in drug discovery.

Mechanism of allosteric activation in human mitochondrial ClpP protease

Human ClpP protease contributes to mitochondrial protein quality control by degrading misfolded proteins. ClpP is overexpressed in cancers such as acute myeloid leukemia (AML), where its inhibition leads to the accumulation of damaged respiratory chain subunits and cell death. Conversely, hyperactivating ClpP with small-molecule activators, such as the recently discovered ONC201, disrupts mitochondrial protein degradation and impairs respiration in cancer cells. Despite its critical role in human health, the mechanism underlying the structural and functional properties of human ClpP remains elusive. Notably, human ClpP is paradoxically activated by active-site inhibitors. All available structures of human ClpP published to date are in the inactive compact or compressed states, surprisingly even when ClpP is bound to an activator molecule such as ONC201. Here, we present structures of human mitochondrial ClpP in the active extended state, including a pair of structures where ClpP is bound to an active-site inhibitor. We demonstrate that amino acid substitutions in the handle region (A192E and E196R) recreate a conserved salt bridge found in bacterial ClpP, stabilizing the extended active state and significantly enhancing ClpP activity. We elucidate the ClpP activation mechanism, highlighting a hormetic effect where substoichiometric inhibitor binding triggers an allosteric transition that drives ClpP into its active extended state. Our findings link the conformational dynamics of ClpP to its catalytic function and provide high-resolution structures for the rational design of potent and specific ClpP inhibitors, with implications for targeting AML and other disorders with ClpP involvement.

Structural Insights into Bortezomib-Induced Activation of the Caseinolytic Chaperone-Protease System in Mycobacterium tuberculosis

The caseinolytic protease (Clp) system has recently emerged as a promising anti-tuberculosis target. The anti-cancer drug bortezomib exhibits potent anti-mycobacterial activity and binds to Mycobacterium tuberculosis (Mtb) Clp protease complexes. We determine cryo-EM structures of Mtb ClpP1P2, ClpC1P1P2 and ClpXP1P2 complexes bound to bortezomib in different conformations. Structural and biochemical data indicate that sub-stoichiometric binding by bortezomib to the protease active sites orthosterically activates the MtbClpP1P2 complex. Bortezomib activation of MtbClpP1P2 induces structural changes promoting the recruitment of the chaperone-unfoldases, MtbClpC1 or MtbClpX, facilitating holoenzyme formation. The structures of the MtbClpC1P1P2 holoenzyme indicate that MtbClpC1 motion, induced by ATP rebinding at the MtbClpC1 spiral seam, translocates the substrate. In the MtbClpXP1P2 holoenzyme structure, we identify a specialized substrate channel gating mechanism involving the MtbClpX pore-2 loop and MtbClpP2 N-terminal domains. Our results provide insights into the intricate regulation of the Mtb Clp system and suggest that bortezomib can disrupt this regulation by sub-stoichiometric binding at the Mtb Clp protease sites.

Importance of integrin transmembrane helical interactions for antagonistic versus agonistic ligand behavior: Consequences for medical applications

Integrins are well-characterized receptors involved in cell adhesion and signaling. With six approved drugs, they are recognized as valuable therapeutic targets. Here, we explore potential activation mechanisms that may clarify the agonist versus antagonist behavior of integrin ligands. The reorganization of the transmembrane domain (TMD) in the integrin receptor, forming homooligomers within focal adhesions, could be key to the understanding of the agonistic properties of integrin ligands at substoichiometric concentrations. This has significant implications for medical applications. While we focus on the RGD peptide-recognizing integrin subfamily, we propose that these mechanistic insights may also apply to other integrin subtypes. For application of integrin ligands in medicine it is essential to consider this mechanism and its consequences for affinity and bioavailability.

Synthesis and biological characterization of a 17β hydroxysteroid dehydrogenase type 10 (17β-HSD10) inhibitor

Alzheimer's disease (AD) is estimated to affect over 55 million people across the world. Small molecule treatment options are limited to symptom management with no impact on disease progression. The need for new protein targets and small molecule hit compounds is unmet and urgent. Hydroxysteroid 17-β dehydrogenase type 10 (17β-HSD10) is a mitochondrial enzyme known to bind amyloid beta, a hallmark of AD, and potentiate its toxicity to neurons. Identification of small molecules capable of interacting with 17β-HSD10 may drive drug discovery efforts for AD. The screening compound BCC0100281 (1), was previously identified as an inhibitor of 17β-HSD10. Herein we report the first synthetic access to the hit compound following a convergent pathway starting from simple heterocyclic building blocks. The compound was found to be toxic to 'neuron-like' cells, specifically those of neuroblastoma origin, providing a potential hit compound for cancer drug discovery, wherein the protein is known to be overexpressed. However, assay of synthetic intermediates identified novel scaffolds with effect to rescue amyloid beta-induced cytotoxicity, showcasing the power of organic synthesis and medicinal chemistry to optimize hit compounds.

Rational correction of pathogenic conformational defects in HTRA1

Loss-of-function mutations in the homotrimeric serine protease HTRA1 cause cerebral vasculopathy. Here, we establish independent approaches to achieve the functional correction of trimer assembly defects. Focusing on the prototypical R274Q mutation, we identify an HTRA1 variant that promotes trimer formation thus restoring enzymatic activity in vitro. Genetic experiments in Htra1R274Q mice further demonstrate that expression of this protein-based corrector in trans is sufficient to stabilize HtrA1-R274Q and restore the proteomic signature of the brain vasculature. An alternative approach employs supramolecular chemical ligands that shift the monomer-trimer equilibrium towards proteolytically active trimers. Moreover, we identify a peptidic ligand that activates HTRA1 monomers. Our findings open perspectives for tailored protein repair strategies.

Heritable epigenetic changes are constrained by the dynamics of regulatory architectures

Interacting molecules create regulatory architectures that can persist despite turnover of molecules. Although epigenetic changes occur within the context of such architectures, there is limited understanding of how they can influence the heritability of changes. Here, I develop criteria for the heritability of regulatory architectures and use quantitative simulations of interacting regulators parsed as entities, their sensors, and the sensed properties to analyze how architectures influence heritable epigenetic changes. Information contained in regulatory architectures grows rapidly with the number of interacting molecules and its transmission requires positive feedback loops. While these architectures can recover after many epigenetic perturbations, some resulting changes can become permanently heritable. Architectures that are otherwise unstable can become heritable through periodic interactions with external regulators, which suggests that mortal somatic lineages with cells that reproducibly interact with the immortal germ lineage could make a wider variety of architectures heritable. Differential inhibition of the positive feedback loops that transmit regulatory architectures across generations can explain the gene-specific differences in heritable RNA silencing observed in the nematode Caenorhabditis elegans. More broadly, these results provide a foundation for analyzing the inheritance of epigenetic changes within the context of the regulatory architectures implemented using diverse molecules in different living systems.

A model for stimulation of enzyme activity by a competitive inhibitor based on the interaction of terazosin and phosphoglycerate kinase 1

The drug terazosin (TZ) binds to and can enhance the activity of the glycolytic enzyme phosphoglycerate kinase 1 (PGK1) and can increase ATP levels. That finding prompted studies of TZ in Parkinson's disease (PD) in which decreased neuronal energy metabolism is a hallmark feature. TZ was neuroprotective in cell-based and animal PD models and in large epidemiological studies of humans. However, how TZ might increase PGK1 activity has remained a perplexing question because structural data revealed that the site of TZ binding to PGK1 overlaps with the site of substrate binding, predicting that TZ would competitively inhibit activity. Functional data also indicate that TZ is a competitive inhibitor. To explore the paradoxical observation of a competitive inhibitor increasing enzyme activity under some conditions, we developed a mass action model of TZ and PGK1 interactions using published data on PGK1 kinetics and the effect of varying TZ concentrations. The model indicated that TZ-binding introduces a bypass pathway that accelerates product release. At low concentrations, TZ binding circumvents slow product release and increases the rate of enzymatic phosphotransfer. However, at high concentrations, TZ inhibits PGK1 activity. The model explains stimulation of enzyme activity by a competitive inhibitor and the biphasic dose-response relationship for TZ and PGK1 activity. By providing a plausible mechanism for interactions between TZ and PGK1, these findings may aid development of TZ or other agents as potential therapeutics for neurodegenerative diseases. The results may also have implications for agents that interact with the active site of other enzymes.

HtrAs are essential for the survival of the haloarchaeon Natrinema gari J7-2 in response to heat, high salinity, and toxic substances

Bacterial and eukaryotic HtrAs can act as an extracytoplasmic protein quality control (PQC) system to help cells survive in stress conditions, but the functions of archaeal HtrAs remain unknown. Particularly, haloarchaea route most secretory proteins to the Tat pathway, enabling them to fold properly in well-controlled cytoplasm with cytosolic PQC systems before secretion. It is unclear whether HtrAs are required for haloarchaeal survival and stress response. The haloarchaeon Natrinema gari J7-2 encodes three Tat signal peptide-bearing HtrAs (NgHtrA, NgHtrB, and NgHtrC), and the signal peptides of NgHtrA and NgHtrC contain a lipobox. Here, the in vitro analysis reveals that the three HtrAs show different profiles of temperature-, salinity-, and metal ion-dependent proteolytic activities and could exhibit chaperone-like activities to prevent the aggregation of reduced lysozyme when their proteolytic activities are inhibited at low temperatures or the active site is disrupted. The gene deletion and complementation assays reveal that NgHtrA and NgHtrC are essential for the survival of strain J7-2 at elevated temperature and/or high salinity and contribute to the resistance of this haloarchaeon to zinc and inhibitory substances generated from tryptone. Mutational analysis shows that the lipobox mediates membrane anchoring of NgHtrA or NgHtrC, and both the membrane-anchored and free extracellular forms of the two enzymes are involved in the stress resistance of strain J7-2, depending on the stress conditions. Deletion of the gene encoding NgHtrB in strain J7-2 causes no obvious growth defect, but NgHtrB can functionally substitute for NgHtrA or NgHtrC under some conditions.IMPORTANCEHtrA-mediated protein quality control plays an important role in the removal of aberrant proteins in the extracytoplasmic space of living cells, and the action mechanisms of HtrAs have been extensively studied in bacteria and eukaryotes; however, information about the function of archaeal HtrAs is scarce. Our results demonstrate that three HtrAs of the haloarchaeon Natrinema gari J7-2 possess both proteolytic and chaperone-like activities, confirming that the bifunctional nature of HtrAs is conserved across all three domains of life. Moreover, we found that NgHtrA and NgHtrC are essential for the survival of strain J7-2 under stress conditions, while NgHtrB can serve as a substitute for the other two HtrAs under certain circumstances. This study provides the first biochemical and genetic evidence of the importance of HtrAs for the survival of haloarchaea in response to stresses.

Small-molecule TIP60 inhibitors enhance regulatory T cell induction through TIP60-P300 acetylation crosstalk

Foxp3 acetylation is essential to regulatory T (Treg) cell stability and function, but pharmacologically increasing it remains an unmet challenge. Here, we report that small-molecule compounds that inhibit TIP60, an acetyltransferase known to acetylate Foxp3, unexpectedly increase Foxp3 acetylation and Treg induction. Utilizing a dual experimental/computational approach combined with a newly developed FRET-based methodology compatible with flow cytometry to measure Foxp3 acetylation, we unraveled the mechanism of action of these small-molecule compounds in murine and human Treg induction cell cultures. We demonstrate that at low-mid concentrations they activate TIP60 to acetylate P300, a different acetyltransferase, which in turn increases Foxp3 acetylation, thereby enhancing Treg cell induction. These results reveal a potential therapeutic target relevant to autoimmunity and transplant.

Heritable epigenetic changes are constrained by the dynamics of regulatory architectures

Interacting molecules create regulatory architectures that can persist despite turnover of molecules. Although epigenetic changes occur within the context of such architectures, there is limited understanding of how they can influence the heritability of changes. Here I develop criteria for the heritability of regulatory architectures and use quantitative simulations of interacting regulators parsed as entities, their sensors and the sensed properties to analyze how architectures influence heritable epigenetic changes. Information contained in regulatory architectures grows rapidly with the number of interacting molecules and its transmission requires positive feedback loops. While these architectures can recover after many epigenetic perturbations, some resulting changes can become permanently heritable. Such stable changes can (1) alter steady-state levels while preserving the architecture, (2) induce different architectures that persist for many generations, or (3) collapse the entire architecture. Architectures that are otherwise unstable can become heritable through periodic interactions with external regulators, which suggests that the evolution of mortal somatic lineages with cells that reproducibly interact with the immortal germ lineage could make a wider variety of regulatory architectures heritable. Differential inhibition of the positive feedback loops that transmit regulatory architectures across generations can explain the gene-specific differences in heritable RNA silencing observed in the nematode C. elegans, which range from permanent silencing to recovery from silencing within a few generations and subsequent resistance to silencing. More broadly, these results provide a foundation for analyzing the inheritance of epigenetic changes within the context of the regulatory architectures implemented using diverse molecules in different living systems.

Activation of the integrated stress response by inhibitors of its kinases

Phosphorylation of the translation initiation factor eIF2α to initiate the integrated stress response (ISR) is a vital signalling event. Protein kinases activating the ISR, including PERK and GCN2, have attracted considerable attention for drug development. Here we find that the widely used ATP-competitive inhibitors of PERK, GSK2656157, GSK2606414 and AMG44, inhibit PERK in the nanomolar range, but surprisingly activate the ISR via GCN2 at micromolar concentrations. Similarly, a PKR inhibitor, C16, also activates GCN2. Conversely, GCN2 inhibitor A92 silences its target but induces the ISR via PERK. These findings are pivotal for understanding ISR biology and its therapeutic manipulations because most preclinical studies used these inhibitors at micromolar concentrations. Reconstitution of ISR activation with recombinant proteins demonstrates that PERK and PKR inhibitors directly activate dimeric GCN2, following a Gaussian activation-inhibition curve, with activation driven by allosterically increasing GCN2 affinity for ATP. The tyrosine kinase inhibitors Neratinib and Dovitinib also activate GCN2 by increasing affinity of GCN2 for ATP. Thus, the mechanism uncovered here might be broadly relevant to ATP-competitive inhibitors and perhaps to other kinases.

10q26 - The enigma in age-related macular degeneration

Despite comprehensive research efforts over the last decades, the pathomechanisms of age-related macular degeneration (AMD) remain far from being understood. Large-scale genome wide association studies (GWAS) were able to provide a defined set of genetic aberrations which contribute to disease risk, with the strongest contributors mapping to distinct regions on chromosome 1 and 10. While the chromosome 1 locus comprises factors of the complement system with well-known functions, the role of the 10q26-locus in AMD-pathophysiology remains enigmatic. 10q26 harbors a cluster of three functional genes, namely PLEKHA1, ARMS2 and HTRA1, with most of the AMD-associated genetic variants mapping to the latter two genes. High linkage disequilibrium between ARMS2 and HTRA1 has kept association studies from reliably defining the risk-causing gene for long and only very recently the genetic risk region has been narrowed to ARMS2, suggesting that this is the true AMD gene at this locus. However, genetic associations alone do not suffice to prove causality and one or more of the 14 SNPs on this haplotype may be involved in long-range control of gene expression, leaving HTRA1 and PLEKHA1 still suspects in the pathogenic pathway. Both, ARMS2 and HTRA1 have been linked to extracellular matrix homeostasis, yet their exact molecular function as well as their role in AMD pathogenesis remains to be uncovered. The transcriptional regulation of the 10q26 locus adds an additional level of complexity, given, that gene-regulatory as well as epigenetic alterations may influence expression levels from 10q26 in diseased individuals. Here, we provide a comprehensive overview on the 10q26 locus and its three gene products on various levels of biological complexity and discuss current and future research strategies to shed light on one of the remaining enigmatic spots in the AMD landscape.

Purification and Characterization of a DegP-Type Protease from the Marine Bacterium Cobetia amphilecti KMM 296

A new member of the DegP-type periplasmic serine endoproteases of the S1C family from the marine bacterium Cobetia amphilecti KMM 296 (CamSP) was expressed in Escherichia coli cells. The calculated molecular weight, number of amino acids, and isoelectric point (pI) of the mature protein CamSP are 69.957 kDa, 666, and 4.84, respectively. The proteolytic activity of the purified recombinant protease CamSP was 2369.4 and 1550.9 U/mg with the use of 1% bovine serum albumin (BSA) and casein as the substrates, respectively. The enzyme CamSP exhibited maximum activity at pH 6.0-6.2, while it was stable over a wide pH range from 5.8 to 8.5. The optimal temperature for the CamSP protease activity was 50 °C. The enzyme required NaCl or KCl at concentrations of 0.3 and 0.5 M, respectively, for its maximum activity. The Michaelis constant (K_m) and V_max for BSA were determined to be 41.7 µg/mL and 0.036 µg/mL min^-1, respectively. The metal ions Zn²⁺, Cu²⁺, Mn²⁺, Li²⁺, Mg²⁺, and Ca²⁺ slightly activated CamSP, while the addition of CoCl₂ to the incubation mixture resulted in a twofold increase in its protease activity. Ethanol, isopropanol, glycerol, and Triton-X-100 increased the activity of CamSP from two- to four-times. The protease CamSP effectively degraded the wheat flour proteins but had no proteolytic activity towards soybean, corn, and the synthetic substrates, α-benzoyl-Arg-p-nitroanilide (BAPNA) and N-Succinyl-L-alanyl-L-alanyl-L-prolyl-L-phenylalanine 4-nitroanilide (SAPNA).

A Single Residue within the MCR-1 Protein Confers Anticipatory Resilience

The envelope stress response (ESR) of Gram-negative enteric bacteria senses fluctuations in nutrient availability and environmental changes to avert damage and promote survival. It has a protective role toward antimicrobials, but direct interactions between ESR components and antibiotic resistance genes have not been demonstrated. Here, we report interactions between a central regulator of ESR viz., the two-component signal transduction system CpxRA (conjugative pilus expression), and the recently described mobile colistin resistance protein (MCR-1). Purified MCR-1 is specifically cleaved within its highly conserved periplasmic bridge element, which links its N-terminal transmembrane domain with the C-terminal active-site periplasmic domain, by the CpxRA-regulated serine endoprotease DegP. Recombinant strains harboring cleavage site mutations in MCR-1 are either protease resistant or degradation susceptible, with widely differing consequences for colistin resistance. Transfer of the gene encoding a degradation-susceptible mutant to strains that lack either DegP or its regulator CpxRA restores expression and colistin resistance. MCR-1 production in Escherichia coli imposes growth restriction in strains lacking either DegP or CpxRA, effects that are reversed by transactive expression of DegP. Excipient allosteric activation of the DegP protease specifically inhibits growth of isolates carrying mcr-1 plasmids. As CpxRA directly senses acidification, growth of strains at moderately low pH dramatically increases both MCR-1-dependent phosphoethanolamine (PEA) modification of lipid A and colistin resistance levels. Strains expressing MCR-1 are also more resistant to antimicrobial peptides and bile acids. Thus, a single residue external to its active site induces ESR activity to confer resilience in MCR-1-expressing strains to commonly encountered environmental stimuli, such as changes in acidity and antimicrobial peptides. Targeted activation of the nonessential protease DegP can lead to the elimination of transferable colistin resistance in Gram-negative bacteria. IMPORTANCE The global presence of transferable mcr genes in a wide range of Gram-negative bacteria from clinical, veterinary, food, and aquaculture environments is disconcerting. Its success as a transmissible resistance factor remains enigmatic, because its expression imposes fitness costs and imparts only moderate levels of colistin resistance. Here, we show that MCR-1 triggers regulatory components of the envelope stress response, a system that senses fluctuations in nutrient availability and environmental changes, to promote bacterial survival in low pH environments. We identify a single residue within a highly conserved structural element of mcr-1 distal to its catalytic site that modulates resistance activity and triggers the ESR. Using mutational analysis, quantitative lipid A profiling and biochemical assays, we determined that growth in low pH environments dramatically increases colistin resistance levels and promotes resistance to bile acids and antimicrobial peptides. We exploited these findings to develop a targeted approach that eliminates mcr-1 and its plasmid carriers.

It's ok to be outnumbered - sub-stoichiometric modulation of homomeric protein complexes

An arsenal of molecular tools with increasingly diversified mechanisms of action is being developed by the scientific community to enable biological interrogation and pharmaceutical modulation of targets and pathways of ever increasing complexity. While most small molecules interact with the target of interest in a 1 : 1 relationship, a noteworthy number of recent examples were reported to bind in a sub-stoichiometric manner to a homomeric protein complex. This approach requires molecular understanding of the physiologically relevant protein assemblies and in-depth characterization of the compound's mechanism of action. The recent literature examples summarized here were selected to illustrate methods used to identify and characterize molecules with such mechanisms. The concept of one small molecule targeting a homomeric protein assembly is not new but the subject deserves renewed inspection in light of emerging technologies and increasingly diverse target biology, to ensure relevant in vitro systems are used and valuable compounds with potentially novel sub-stoichiometric mechanisms of action aren't overlooked.

Allosteric inhibition of HTRA1 activity by a conformational lock mechanism to treat age-related macular degeneration

The trimeric serine protease HTRA1 is a genetic risk factor associated with geographic atrophy (GA), a currently untreatable form of age-related macular degeneration. Here, we describe the allosteric inhibition mechanism of HTRA1 by a clinical Fab fragment, currently being evaluated for GA treatment. Using cryo-EM, X-ray crystallography and biochemical assays we identify the exposed LoopA of HTRA1 as the sole Fab epitope, which is approximately 30 Å away from the active site. The cryo-EM structure of the HTRA1:Fab complex in combination with molecular dynamics simulations revealed that Fab binding to LoopA locks HTRA1 in a non-competent conformational state, incapable of supporting catalysis. Moreover, grafting the HTRA1-LoopA epitope onto HTRA2 and HTRA3 transferred the allosteric inhibition mechanism. This suggests a conserved conformational lock mechanism across the HTRA family and a critical role of LoopA for catalysis, which was supported by the reduced activity of HTRA1-3 upon LoopA deletion or perturbation. This study reveals the long-range inhibition mechanism of the clinical Fab and identifies an essential function of the exposed LoopA for activity of HTRA family proteases.

Substoichiometric action of long noncoding RNAs

Low expression levels and stoichiometric imbalances of long noncoding RNAs (lncRNAs) are often used as evidence for their probable lack of function or for limiting the scope of their potential influence. Recent advances in our understanding of the substoichiometric functions of lncRNAs challenge these notions and suggest routes through which unabundant lncRNAs can affect cellular functions and gene regulatory networks.

An allosteric HTRA1-calpain 2 complex with restricted activation profile

SignificanceClassic serine proteases are synthesized as inactive precursors that are proteolytically processed, resulting in irreversible activation. We report an alternative and reversible mechanism of activation that is executed by an inactive protease. This mechanism involves a protein complex between the serine protease HTRA1 and the cysteine protease calpain 2. Surprisingly, activation is restricted as it improves the proteolysis of soluble tau protein but not the dissociation and degradation of its amyloid fibrils, a task that free HTRA1 is efficiently performing. These data exemplify a challenge for protein quality control proteases in the clearing of pathogenic fibrils and suggest a potential for unexpected side effects of chemical modulators targeting PDZ or other domains located at a distance to the active site.

Computational Design of Novel Allosteric Inhibitors for Plasmodium falciparum DegP

The serine protease, DegP exhibits proteolytic and chaperone activities, essential for cellular protein quality control and normal cell development in eukaryotes. The P. falciparum DegP is essential for the parasite survival and required to combat the oscillating thermal stress conditions during the infection, protein quality checks and protein homeostasis in the extra-cytoplasmic compartments, thereby establishing it as a potential target for drug development against malaria. Previous studies have shown that diisopropyl fluorophosphate (DFP) and the peptide SPMFKGV inhibit E. coli DegP protease activity. To identify novel potential inhibitors specific to PfDegP allosteric and the catalytic binding sites, we performed a high throughput in silico screening using Malaria Box, Pathogen Box, Maybridge library, ChEMBL library and the library of FDA approved compounds. The screening helped identify five best binders that showed high affinity to PfDegP allosteric (T0873, T2823, T2801, RJC02337, CD00811) and the catalytic binding site (T0078L, T1524, T2328, BTB11534 and 552691). Further, molecular dynamics simulation analysis revealed RJC02337, BTB11534 as the best hits forming a stable complex. WaterMap and electrostatic complementarity were used to evaluate the novel bio-isosteric chemotypes of RJC02337, that led to the identification of 231 chemotypes that exhibited better binding affinity. Further analysis of the top 5 chemotypes, based on better binding affinity, revealed that the addition of electron donors like nitrogen and sulphur to the side chains of butanoate group are more favoured than the backbone of butanoate group. In a nutshell, the present study helps identify novel, potent and Plasmodium specific inhibitors, using high throughput in silico screening and bio-isosteric replacement, which may be experimentally validated.