Fragment-based lead discovery to identify novel inhibitors that target the ATP binding site of pyruvate dehydrogenase kinases

Enhanced inhibitory activity of compounds containing purine scaffolds compared to protein kinase CK2α considering crystalline water

We recently reported novel purine-based CK2α inhibitors using the solvent ordering-based method as virtual screening. Among these, the X-ray crystal structure of a complex with CK2α was determined. The results showed that the crystalline water molecules observed in many previously reported complex structures of CK2α and its inhibitors had been eliminated. We then proposed a structure-based drug design. Since the removal of water molecules would be detrimental to inhibitor binding, new groups of compounds were designed by changing the position of the carboxy group located at the point where a water molecule would be present so as not to eliminate it. Compounds with (E)-2-carboxyethenyl and 3-carboxyphenyl substituted at the 2-position on the purine scaffold showed much higher inhibitory potency than 4-carboxyphenyl derivatives. Furthermore, in the presence of a 4-fluorophenyl group at the 9-position on the purine scaffold, the inhibitory activity of the 3-carboxyphenyl derivative against CK2α was 0.18 μM, a 167-fold improvement compared to the 4-carboxyphenyl derivative. The strategy of leaving crystalline water can significantly increase inhibitory activity.

Targeting the PDK/PDH axis to reverse metabolic abnormalities by structure-based virtual screening with in vitro and in vivo experiments

In humans and animals, the pyruvate dehydrogenase kinase (PDK) family proteins (PDKs 1-4) are excessively activated in metabolic disorders such as obesity, diabetes, and cancer, inhibiting the activity of pyruvate dehydrogenase (PDH) which plays a crucial role in energy and fatty acid metabolism and impairing its function. Intervention and regulation of PDH activity have become important research approaches for the treatment of various metabolic disorders. In this study, a small molecule (g25) targeting PDKs and activating PDH, was identified through multi-level computational screening methods. In vivo and in vitro experiments have shown that g25 activated the activity of PDH and reduced plasma lactate and triglyceride level. Besides, g25 significantly decreased hepatic fat deposition in a diet-induced obesity mouse model. Furthermore, g25 enhanced the tumor-inhibiting activity of cisplatin when used in combination. Molecular dynamics simulations and in vitro kinase assay also revealed the specificity of g25 towards PDK2. Overall, these findings emphasize the importance of targeting the PDK/PDH axis to regulate PDH enzyme activity in the treatment of metabolic disorders, providing directions for future related research. This study provides a possible lead compound for the PDK/PDH axis related diseases and offers insights into the regulatory mechanisms of this pathway in diseases.

Thermal Titration Molecular Dynamics (TTMD): Not Your Usual Post-Docking Refinement

Molecular docking is one of the most widely used computational approaches in the field of rational drug design, thanks to its favorable balance between the rapidity of execution and the accuracy of provided results. Although very efficient in exploring the conformational degrees of freedom available to the ligand, docking programs can sometimes suffer from inaccurate scoring and ranking of generated poses. To address this issue, several post-docking filters and refinement protocols have been proposed throughout the years, including pharmacophore models and molecular dynamics simulations. In this work, we present the first application of Thermal Titration Molecular Dynamics (TTMD), a recently developed method for the qualitative estimation of protein-ligand unbinding kinetics, to the refinement of docking results. TTMD evaluates the conservation of the native binding mode throughout a series of molecular dynamics simulations performed at progressively increasing temperatures with a scoring function based on protein-ligand interaction fingerprints. The protocol was successfully applied to retrieve the native-like binding pose among a set of decoy poses of drug-like ligands generated on four different pharmaceutically relevant biological targets, including casein kinase 1δ, casein kinase 2, pyruvate dehydrogenase kinase 2, and SARS-CoV-2 main protease.

Structural studies identify angiotensin II receptor blocker-like compounds as branched-chain ketoacid dehydrogenase kinase inhibitors

The mammalian mitochondrial branched-chain ketoacid dehydrogenase (BCKD) complex is a multienzyme complex involved in the catabolism of branched-chain amino acids. BCKD is regulated by the BCKD kinase, or BCKDK, which binds to the E2 subunit of BCKD, phosphorylates its E1 subunit, and inhibits enzymatic activity. Inhibition of the BCKD complex results in increased levels of branched-chain amino acids and branched-chain ketoacids, and this buildup has been associated with heart failure, type 2 diabetes mellitus, and nonalcoholic fatty liver disease. To find BCKDK inhibitors for potential treatment of these diseases, we performed both NMR and virtual fragment screening and identified tetrazole-bearing fragments that bind BCKDK at multiple sites. Through structure-based virtual screening expanding from these fragments, the angiotensin receptor blocker class antihypertension drugs and angiotensin receptor blocker-like compounds were discovered to be potent BCKDK inhibitors, suggesting potential new avenues for heart failure treatment combining BCKDK inhibition and antihypertension.

Fragment-to-Lead Medicinal Chemistry Publications in 2021

This Perspective is the seventh in an annual series that summarizes successful Fragment-to-Lead (F2L) case studies published in a given year. A tabulated summary of relevant articles published in 2021 is provided, and features such as target class, screening methods, and ligand efficiency are discussed, both for the 2021 examples and for the combined examples over the years 2015-2021. In addition, trends and new developments in the field are summarized. In particular, the use of structural information in fragment-based drug discovery is discussed.

Qualitative Estimation of Protein-Ligand Complex Stability through Thermal Titration Molecular Dynamics Simulations

The prediction of ligand efficacy has long been linked to thermodynamic properties such as the equilibrium dissociation constant, which considers both the association and the dissociation rates of a defined protein-ligand complex. In the last 15 years, there has been a paradigm shift, with an increased interest in the determination of kinetic properties such as the drug-target residence time since they better correlate with ligand efficacy compared to other parameters. In this article, we present thermal titration molecular dynamics (TTMD), an alternative computational method that combines a series of molecular dynamics simulations performed at progressively increasing temperatures with a scoring function based on protein-ligand interaction fingerprints for the qualitative estimation of protein-ligand-binding stability. The protocol has been applied to four different pharmaceutically relevant test cases, including protein kinase CK1δ, protein kinase CK2, pyruvate dehydrogenase kinase 2, and SARS-CoV-2 main protease, on a variety of ligands with different sizes, structures, and experimentally determined affinity values. In all four cases, TTMD was successfully able to distinguish between high-affinity compounds (low nanomolar range) and low-affinity ones (micromolar), proving to be a useful screening tool for the prioritization of compounds in a drug discovery campaign.

Structure-based drug design of novel and highly potent pyruvate dehydrogenase kinase inhibitors

Pyruvate dehydrogenase kinases (PDHKs) are fascinating drug targets for numerous diseases, including diabetes and cancers. In this report, we describe the result of our structure-based drug design from tricyclic lead compounds that led to the discovery of highly potent PDHK2 and PDHK4 dual inhibitors in enzymatic assay. The C3-position of the tricyclic core was explored, and the PDHK2 X-ray structure with a representative compound revealed a novel ATP lid conformation in which the phenyl ring of Phe326 mediated the interaction of the Arg258 sidechain and the compound. Compounds with amide linkers were designed to release the ATP lid by forming an intramolecular pi-pi interaction, and these compounds showed single-digit nM IC₅₀ values in an enzymatic assay. We also explored the C4-position of the tricyclic core to reproduce the interaction observed with the C3-position substitution, and the pyrrolidine compound showed the same level of IC₅₀ values. By optimizing an interaction with the Asn255 sidechain through a docking simulation, compounds with 2-carboxy pyrrole moiety also showed single-digit nM IC₅₀ values without having a cation-pi interaction with the Arg258 sidechain.