Discovery of novel allosteric effectors based on the predicted allosteric sites for Escherichia coli D-3-phosphoglycerate dehydrogenase

A Comprehensive Survey of Prospective Structure-Based Virtual Screening for Early Drug Discovery in the Past Fifteen Years

Structure-based virtual screening (SBVS), also known as molecular docking, has been increasingly applied to discover small-molecule ligands based on the protein structures in the early stage of drug discovery. In this review, we comprehensively surveyed the prospective applications of molecular docking judged by solid experimental validations in the literature over the past fifteen years. Herein, we systematically analyzed the novelty of the targets and the docking hits, practical protocols of docking screening, and the following experimental validations. Among the 419 case studies we reviewed, most virtual screenings were carried out on widely studied targets, and only 22% were on less-explored new targets. Regarding docking software, GLIDE is the most popular one used in molecular docking, while the DOCK 3 series showed a strong capacity for large-scale virtual screening. Besides, the majority of identified hits are promising in structural novelty and one-quarter of the hits showed better potency than 1 μM, indicating that the primary advantage of SBVS is to discover new chemotypes rather than highly potent compounds. Furthermore, in most studies, only in vitro bioassays were carried out to validate the docking hits, which might limit the further characterization and development of the identified active compounds. Finally, several successful stories of SBVS with extensive experimental validations have been highlighted, which provide unique insights into future SBVS drug discovery campaigns.

Uncovering the Dominant Motion Modes of Allosteric Regulation Improves Allosteric Site Prediction

Allostery is an important mechanism that biological systems use to regulate function at a distant site. Allosteric drugs have attracted much attention in recent years due to their high specificity and the possibility of overcoming existing drug-resistant mutations. However, the discovery of allosteric drugs remains challenging as allosteric regulation mechanisms are not clearly understood and allosteric sites cannot be accurately predicted. In this study, we analyzed the dominant modes that determine motion correlations between allosteric and orthosteric sites using the Gaussian network model and found that motion correlations between allosteric and orthosteric sites are dominated by either fast or slow vibrational modes. This dependence of modes results from the relative locations of the two sites and local secondary structures. Based on these analyses, we developed CorrSite2.0 to predict allosteric sites by taking the maximum of the Z-scores calculated from using either slow or fast modes. The prediction accuracy of CorrSite2.0 outperformed other commonly used allosteric site prediction methods with prediction accuracy over 90.0%. Our study uncovers the relationship of protein structure, dynamics, and allosteric regulation and demonstrates that using the dominant motion modes greatly improves allosteric site prediction accuracy. CorrSite2.0 has been integrated into the CavityPlus web server, which can be accessed at http://www.pkumdl.cn/cavityplus. CorrSite2.0 provides a powerful and user-friendly tool for allosteric drug and protein design.

Diversified amino acid-mediated allosteric regulation of phosphoglycerate dehydrogenase for serine biosynthesis in land plants

The phosphorylated pathway of serine biosynthesis is initiated with 3-phosphoglycerate dehydrogenase (PGDH). The liverwort Marchantia polymorpha possesses an amino acid-sensitive MpPGDH which is inhibited by l-serine and activated by five proteinogenic amino acids, while the eudicot Arabidopsis thaliana has amino acid-sensitive AtPGDH1 and AtPGDH3 as well as amino acid-insensitive AtPGDH2. In this study, we analyzed PGDH isozymes of the representative land plants: the monocot Oryza sativa (OsPGDH1–3), basal angiosperm Amborella trichopoda (AmtriPGDH1–2), and moss Physcomitrium (Physcomitrella) patens (PpPGDH1–4). We demonstrated that OsPGDH1, AmtriPGDH1, PpPGDH1, and PpPGDH3 were amino acid-sensitive, whereas OsPGDH2, OsPGDH3, AmtriPGDH2, PpPGDH2, and PpPGDH4 were either sensitive to only some of the six effector amino acids or insensitive to all effectors. This indicates that PGDH sensitivity to effectors has been diversified among isozymes and that the land plant species examined, except for M. polymorpha, possess different isozyme types in terms of regulation. Phylogenetic analysis suggested that the different sensitivities convergently evolved in the bryophyte and angiosperm lineages. Site-directed mutagenesis of AtPGDH1 revealed that Asp⁵³⁸ and Asn⁵⁵⁶ residues in the ACT domain are involved in allosteric regulation by the effectors. These findings provide insight into the evolution of PGDH isozymes, highlighting the functional diversification of allosteric regulation in land plants.

Validating the inhibitory effects of d- and l-serine on the enzyme activity of d-3-phosphoglycerate dehydrogenases that are purified from Pseudomonas aeruginosa, Escherichia coli and human colon

Background

We previously demonstrated that the serA gene is associated with bacterial pathogenicity, including bacterial penetration through the Caco-2 cell monolayers, bacterial motility, bacterial adherence, and fly mortality. l-Serine is known to inhibit the d-3-phosphoglycerate dehydrogenase (PGDH) activity of the SerA protein, and it significantly reduced the bacterial pathogenicity as described above. We also demonstrated that in a PGDH assay using crude extracts isolated from overnight cultures of E. coli overexpressing the P. aeruginosa serA gene, l-serine inhibited the PGDH activity of the SerA protein. The basal PGDH activity of the negative control strain was high, presumably due to contamination of unknown proteins in the crude extracts. Therefore, to further confirm the direct inhibition of PGDH activity of P. aeruginosa SerA by l-serine, we purified and characterized the PGDH from P. aeruginosa and compared it with the previously characterized PGDHs from E. coli, and the human colon as controls.

Results

Optimum pH and ionic strength of the purified PGDHs were different depending on the three species; optimal activity of P. aeruginosa PGDH was at pH 7.5 with 50-100 mM Tris-HCl, E. coli PGDH was at pH 8.5 with 100-200 mM Tris-HCl, and human PGDH was at pH 9.0 with 100-200 mM Tris-HCl. The addition of l-serine reduced the activity of PGDH from P. aeruginosa and E. coli, but not the PGDH from human colon. The median inhibitory concentration (IC₅₀) of l-serine was 630 μM for P. aeruginosa and 250 μM for E. coli, while IC₅₀ of d-serine was much higher than that of l-serine; 76 mM in P. aeruginosa PGDH and 45 mM in E. coli PGDH.

Conclusions

These results suggest that l-serine significantly repressed P. aeruginosa pathogenicity through direct inhibition of the PGDH activity, but was not able to inhibit the human PGDH activity. Oral administration of l-serine to compromised hosts might interfere with bacterial translocation and prevent gut-derived sepsis caused by P. aeruginosa through inhibition of the function of the serA gene product.

Correlation Between Allosteric and Orthosteric Sites

Correlation between an allosteric site and its orthosteric site refers to the phenomenon that perturbations like ligand binding, mutation, or posttranslational modifications at the allosteric site leverage variation in the orthosteric site. Understanding this kind of correlation not only helps to disclose how information is transmitted in allosteric regulation but also provides clues for allosteric drug discovery. This chapter starts with an overview of correlation studies on allosteric and orthosteric sites and then introduces recent progress in evolutionary and simulation-based dynamic studies. Discussions and perspectives on future directions are also given.

D-3-Phosphoglycerate Dehydrogenase

l-Serine is the immediate precursor of d-serine, a major agonist of the N-methyl-d-aspartate (NMDA) receptor. l-Serine is a pivotal amino acid since it serves as a precursor to a large number of essential metabolites besides d-serine. In all non-photosynthetic organisms, including mammals, a major source of l-serine is the phosphorylated pathway of l-serine biosynthesis. The pathway consists of three enzymes, d-3-phosphoglycerate dehydrogenase (PGDH), phosphoserine amino transferase (PSAT), and l-phosphoserine phosphatase (PSP). PGDH catalyzes the first step in the pathway by converting d-3-phosphoglycerate (PGA), an intermediate in glycolysis, to phosphohydroxypyruvate (PHP) concomitant with the reduction of NAD⁺. In some, but not all organisms, the catalytic activity of PGDH can be regulated by feedback inhibition by l-serine. Three types of PGDH can be distinguished based on their domain structure. Type III PGDHs contain only a nucleotide binding and substrate binding domain. Type II PGDHs contain an additional regulatory domain (ACT domain), and Type I PGDHs contain a fourth domain, termed the ASB domain. There is no consistent pattern of domain content that correlates with organism type, and even when additional domains are present, they are not always functional. PGDH deficiency results in metabolic defects of the nervous system whose systems range from microcephaly at birth, seizures, and psychomotor retardation. Although deficiency of any of the pathway enzymes have similar outcomes, PGDH deficiency is predominant. Dietary or intravenous supplementation with l-serine is effective in controlling seizures but has little effect on psychomotor development. An increase in PGDH levels, due to overexpression, is also associated with a wide array of cancers. In culture, PGDH is required for tumor cell proliferation, but extracellular l-serine is not able to support cell proliferation. This has led to the hypothesis that the pathway is performing some function related to tumor growth other than supplying l-serine. The most well-studied PGDHs are bacterial, primarily from Escherichia coli and Mycobacterium tuberculosis, perhaps because they have been of most interest mechanistically. However, the relatively recent association of PGDH with neuronal defects and human cancers has provoked renewed interest in human PGDH.

Novel Allosteric Activators for Ferroptosis Regulator Glutathione Peroxidase 4

Glutathione peroxidase 4 (GPX4) is essential for cell membrane repair, inflammation suppression, and ferroptosis inhibition. GPX4 upregulation provides unique drug discovery opportunities for inflammation and ferroptosis-related diseases. However, rational design of protein activators is challenging. Until now, no compound has been reported to activate the enzyme activity of GPX4. Here, we identified a potential allosteric site in GPX4 and successfully found eight GPX4 activators using a novel computational strategy and experimental studies. Compound 1 from the virtual screen increased GPX4 activity, suppressed ferroptosis, reduced pro-inflammatory lipid mediator production, and inhibited NF-κB pathway activation. Further chemical synthesis and structure-activity relationship studies revealed seven more activators. The strongest compound, 1d4, increased GPX4 activity to 150% at 20 μM in the cell-free assay and 61 μM in cell extracts. Therefore, we demonstrated that GPX4 can be directly activated using chemical compounds to suppress ferroptosis and inflammation. Meanwhile, the discovery of GPX4 activators verified the possibility of rational design of allosteric activators.

Engineered control of enzyme structural dynamics and function

Enzymes undergo a range of internal motions from local, active site fluctuations to large-scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post-translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus-responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.

A novel allosteric inhibitor that prevents IKKβ activation

I kappa B kinase β (IKKβ) is one of the primary targets to regulate canonical NF-κB activity. The misregulation of NF-κB is associated with various diseases, including chronic inflammation and cancers. Most of the known IKKβ inhibitors target its active form and suffer from poor selectivity. In the present study, we aim to design inhibitors that can bind to the IKKβ inactive form and block its activation. We identified a potential allosteric site between the kinase domain (KD) and ubiquitin-like domain (ULD) of human IKKβ and used it to virtually screen a chemical library for allosteric inhibitors. Among the 133 compounds tested, 16 inhibited NF-κB activity by over 50% at 50 μM in a reporter gene assay. Further quantitative measurements and cytotoxicity study gave one compound 124 (3,4-dichloro-2-ethoxy-N-(2,2,6,6-tetramethylpiperidin-4-yl)benzenesulfonamide) which specifically targets the IKKβ inactive form. In cells, 124 inhibited IκBα phosphorylation and NF-κB transcriptional activity for the reporter gene with an IC₅₀ of 35 μM by decreasing the phosphorylation level of Ser177/181 on IKKβ and blocking its activation upon TNFα stimulation. Molecular dynamics simulations demonstrated that 124 binds to the pocket between KD and ULD in the inactive conformation of IKKβ rather than the active conformation. As the first allosteric inhibitor that prevents IKKβ activation, 124 provides a good starting point for further inhibitor discovery and a probe for IKKβ enzyme cycle and regulatory mechanism study.

Identification of a New Potent Inhibitor Targeting KRAS in Non-small Cell Lung Cancer Cells

KRAS (v-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog) is an oncogenic driver with mutations in 30% of non-small cell lung cancer (NSCLC). However, there is no effective clinical drug even though it has been identified as an oncogene for 30 years. In this study, we identified a small molecule inhibitor compound 0375-0604 targeting KRAS by using molecular docking based virtual screening approach. Compound 0375-0604 had a good binding affinity to KRAS in vitro and exhibited cytotoxicity in oncogenic KRAS expressing NSCLC cell lines. Further mechanism study showed that compound 0375-0604 can block the formation of the complex of guanosine triphosphate (GTP) and KRAS in vitro. In addition, compound 0375-0604 inhibited KRAS downstream signaling pathway RAF/MEK/ERK and RAF/PI3K/AKT. Finally, we also found that this compound can inhibit the cell growth through G2/M cell cycle arrest and induce apoptosis on the NSCLC cell lines harboring KRAS mutation. Therefore, compound 0375-0604 may be considered as a potential KRAS inhibitor for treatment of NSCLC carrying KRAS oncogene.

Novel regulatory mechanism of serine biosynthesis associated with 3-phosphoglycerate dehydrogenase in Arabidopsis thaliana

The proteinogenic amino acid l-serine is a precursor for various essential biomolecules in all organisms. 3-Phosphoglycerate dehydrogenase (PGDH) is the first committed enzyme of the phosphorylated pathway of l-serine biosynthesis, and is regulated by negative feedback from l-serine in bacteria and plants. In the present study, two Arabidopsis PGDH isoforms were inhibited by l-serine but were activated by l-amino acids such as l-homocysteine in vitro. Activation and inhibition by these amino acids was cooperative, suggesting an allosteric mechanism. Moreover, the half maximal effective concentration of l-homocysteine was 2 orders of magnitude lower than that of l-serine, suggesting greater regulatory potency. These are the first data to show that PGDH is activated by various biomolecules and indicate that serine biosynthesis is regulated by multiple pathways.

Rational Design of Selective Allosteric Inhibitors of PHGDH and Serine Synthesis with Anti-tumor Activity

Metabolic reprogramming in cancer cells facilitates growth and proliferation. Increased activity of the serine biosynthetic pathway through the enzyme phosphoglycerate dehydrogenase (PHGDH) contributes to tumorigenesis. With a small substrate and a weak binding cofactor, (NAD⁺), inhibitor development for PHGDH remains challenging. Instead of targeting the PHGDH active site, we computationally identified two potential allosteric sites and virtually screened compounds that can bind to these sites. With subsequent characterization, we successfully identified PHGDH non-NAD⁺-competing allosteric inhibitors that attenuate its enzyme activity, selectively inhibit de novo serine synthesis in cancer cells, and reduce tumor growth in vivo. Our study not only identifies novel allosteric inhibitors for PHGDH to probe its function and potential as a therapeutic target, but also provides a general strategy for the rational design of small-molecule modulators of metabolic enzyme function.

Allosteric sites can be identified based on the residue-residue interaction energy difference

Allosteric drugs act at a distance to regulate protein functions. They have several advantages over conventional orthosteric drugs, including diverse regulation types and fewer side effects. However, the rational design of allosteric ligands remains a challenge, especially when it comes to the identification allosteric binding sites. As the binding of allosteric ligands may induce changes in the pattern of residue-residue interactions, we calculated the residue-residue interaction energies within the allosteric site based on the molecular mechanics generalized Born surface area energy decomposition scheme. Using a dataset of 17 allosteric proteins with structural data for both the apo and the ligand-bound state available, we used conformational ensembles generated by molecular dynamics simulations to compute the differences in the residue-residue interaction energies in known allosteric sites from both states. For all the known sites, distinct interaction energy differences (>25%) were observed. We then used CAVITY, a binding site detection program to identify novel putative allosteric sites in the same proteins. This yielded a total of 31 "druggable binding sites," of which 21 exhibited >25% difference in residue interaction energies, and were hence predicted as novel allosteric sites. Three of the predicted allosteric sites were supported by recent experimental studies. All the predicted sites may serve as novel allosteric sites for allosteric ligand design. Our study provides a computational method for identifying novel allosteric sites for allosteric drug design.

Discovery of Non-ATP-Competitive Inhibitors of Polo-like Kinase 1

Polo-like kinase 1 (Plk1) is an evolutionarily conserved serine/threonine kinase, and its N-terminal kinase domain (KD) controls cell signaling through phosphorylation. Inhibitors of Plk1 are potential anticancer drugs. Most known Plk1 KD inhibitors are ATP-competitive compounds, which may suffer from low selectivity. In this study we discovered novel non-ATP-competitive Plk1 KD inhibitors by virtual screening and experimental studies. Potential binding sites in Plk1 KD were identified by using the protein binding site detection program Cavity. The identified site was subjected to molecular-docking-based virtual screening. The activities of top-ranking compounds were evaluated by in vitro enzyme assay with full-length Plk1 and direct binding assay with Plk1 KD. Several compounds showed inhibitory activity, and the most potent was found to be 3-((2-oxo-2-(thiophen-2-yl)ethyl)thio)-6-(pyridin-3-ylmethyl)-1,2,4-triazin-5(4H)-one (compound 4) with an IC50 value of 13.1 ± 1.7 μm. Our work provides new insight into the design of kinase inhibitors that target non-ATP binding sites.

Binding cavities and druggability of intrinsically disordered proteins

To assess the potential of intrinsically disordered proteins (IDPs) as drug design targets, we have analyzed the ligand-binding cavities of two datasets of IDPs (containing 37 and 16 entries, respectively) and compared their properties with those of conventional ordered (folded) proteins. IDPs were predicted to possess more binding cavity than ordered proteins at similar length, supporting the proposed advantage of IDPs economizing genome and protein resources. The cavity number has a wide distribution within each conformation ensemble for IDPs. The geometries of the cavities of IDPs differ from the cavities of ordered proteins, for example, the cavities of IDPs have larger surface areas and volumes, and are more likely to be composed of a single segment. The druggability of the cavities was examined, and the average druggable probability is estimated to be 9% for IDPs, which is almost twice that for ordered proteins (5%). Some IDPs with druggable cavities that are associated with diseases are listed. The optimism versus obstacles for drug design for IDPs is also briefly discussed.

Comparison of Type 1 D-3-phosphoglycerate dehydrogenases reveals unique regulation in pathogenic Mycobacteria

D-3-phosphoglycerate dehydrogenases (PGDH) from all organisms catalyze the conversion of D-3-phosphoglycerate to phosphohydroxypyruvate as the first step in the biosynthesis of l-serine. This investigation compares the properties of Type 1 PGDHs from seven different species and demonstrates that conserved residues in the ACT and ASB domains of some allow l-serine to act as a feedback inhibitor at low micromolar concentrations. In addition, the serine sensitivity is dependent on the presence of phosphate ions. These residues are most highly conserved among PGDHs from the actinomycetales family, but only certain pathogenic mycobacteria appear to have the full complement of residues required for high sensitivity to serine. These basic residues are also responsible for the presence of dual pH optima in the acidic region that is also phosphate dependent. Analytical ultracentrifugation analysis demonstrates that the dual pH optima do not require changes in oligomeric state. This study also demonstrates that substrate inhibition is a common feature of Type 1 PGDHs and that it is suppressed by phosphate, indicating that phosphate likely interacts at both the catalytic and regulatory sites. The unique features resulting from the complement of basic residues conserved in pathogenic mycobacteria may impart important metabolic advantages to these organisms.

Succinate Overproduction: A Case Study of Computational Strain Design Using a Comprehensive Escherichia coli Kinetic Model

Computational strain-design prediction accuracy has been the focus for many recent efforts through the selective integration of kinetic information into metabolic models. In general, kinetic model prediction quality is determined by the range and scope of genetic and/or environmental perturbations used during parameterization. In this effort, we apply the k-OptForce procedure on a kinetic model of E. coli core metabolism constructed using the Ensemble Modeling (EM) method and parameterized using multiple mutant strains data under aerobic respiration with glucose as the carbon source. Minimal interventions are identified that improve succinate yield under both aerobic and anaerobic conditions to test the fidelity of model predictions under both genetic and environmental perturbations. Under aerobic condition, k-OptForce identifies interventions that match existing experimental strategies while pointing at a number of unexplored flux re-directions such as routing glyoxylate flux through the glycerate metabolism to improve succinate yield. Many of the identified interventions rely on the kinetic descriptions that would not be discoverable by a purely stoichiometric description. In contrast, under fermentative (anaerobic) condition, k-OptForce fails to identify key interventions including up-regulation of anaplerotic reactions and elimination of competitive fermentative products. This is due to the fact that the pathways activated under anaerobic condition were not properly parameterized as only aerobic flux data were used in the model construction. This study shed light on the importance of condition-specific model parameterization and provides insight on how to augment kinetic models so as to correctly respond to multiple environmental perturbations.

Recent computational advances in the identification of allosteric sites in proteins

Allosteric modulators have the potential to fine-tune protein functional activity. Therefore, the targeting of allosteric sites, as a strategy in drug design, is gaining increasing attention. Currently, it is not trivial to find and characterize new allosteric sites by experimental approaches. Alternatively, computational approaches are useful in helping researchers analyze and select potential allosteric sites for drug discovery. Here, we review state-of-the-art computational approaches directed at predicting putative allosteric sites in proteins, along with examples of successes in identifying allosteric sites utilizing these methods. We also discuss the challenges in developing reliable methods for predicting allosteric sites and tactics to resolve demanding tasks.