The importance of the quaternary structure to represent conformational ensembles of the major Mycobacterium tuberculosis drug target

Unveiling G-protein coupled receptor kinase-5 inhibitors for chronic degenerative diseases: Multilayered prioritization employing explainable machine learning-driven multi-class QSAR, ligand-based pharmacophore and free energy-inspired molecular simulation

GRK5 holds a pivotal role in cellular signaling pathways, with its overexpression in cardiomyocytes, neuronal cells, and tumor cells strongly associated with various chronic degenerative diseases, which highlights the urgent need for potential inhibitors. In this study, multiclass classification-based QSAR models were developed using diverse machine learning algorithms. These models were built from curated compounds with experimentally derived GRK5 inhibitory activity. Additionally, a pharmacophore model was constructed using active compounds from the dataset. Among the models, the SVM-based approach proved most effective and was initially used to screen DrugBank compounds within the applicability domain. Compounds showing significant GRK5 inhibitory potential underwent evaluation for key pharmacophoric features. Prospective compounds were subjected to molecular docking to assess binding affinity towards GRK5's key active site amino acid residues. Stability at the binding site was analyzed through 200 ns molecular dynamics simulations. MM-GBSA analysis quantified individual free energy components contributing to the total binding energy with respect to binding site residues. Metadynamics analysis, including PCA, FEL, and PDF, provided crucial insights into conformational changes of both apo and holo forms of GRK5 at defined energy states. The study identifies DB02844 (S-Adenosyl-1,8-Diamino-3-Thiooctane) and DB13155 (Esculin) as promising GRK5 inhibitors, warranting further in vitro and in vivo validation studies.

Exploring the dark side of tertiary and quaternary structure dynamics in MtbFBPaseII

Tuberculosis (TB) is a major global cause of mortality, primarily stemming from latent tuberculosis infection (LTBI). Failure to fully treat LTBI can result in drug-resistant forms of TB. Therefore, it is essential to develop novel drugs with unique mechanisms of action to combat TB effectively. One crucial metabolic pathway in Mycobacterium tuberculosis (Mtb), which contributes to TB infection and persistence, is gluconeogenesis. Within this pathway, the enzyme fructose bisphosphatase (FBPase) plays a significant role and is considered a promising target for drug development. By targeting MtbFBPaseII, a specific class of FBPase, researchers have employed molecular dynamics simulations to identify regions capable of binding new drugs, thereby inhibiting the enzyme's activity and potentially paving the way for the development of effective treatments.Communicated by Ramaswamy H. Sarma.

Discovery of new inhibitors of Mycobacterium tuberculosis EPSP synthase - A computational study

Tuberculosis (TB) is a highly infectious disease caused by the pathogen Mycobacterium tuberculosis (Mtb). EPSP Synthase (MtEPSPS), the enzyme responsible for the sixth step of the shikimate pathway, is a potential target for the development of new drugs for the treatment of TB, as it is essential in mycobacteria but absent in humans. In this work, we performed virtual screening using sets of molecules from two databases and three crystallographic structures of MtEPSPS. The initial hits obtained from molecular docking were filtered based on predicted binding affinity and interactions with binding site residues. Subsequently, molecular dynamics simulations were carried out to analyze the stability of protein-ligand complexes. We have found that MtEPSPS forms stable interactions with several candidates, including already approved pharmaceutical drugs such as Conivaptan and Ribavirin monophosphate. In particular, Conivaptan had the highest estimated binding affinity with the open conformation of the enzyme. The complex formed between MtEPSPS and Ribavirin monophosphate was also energetically stable as shown by RMSD, Rg and FEL analyses, and the ligand was stabilized by hydrogen bonds with important residues of the binding site. The findings reported in this work could serve as the basis of promising scaffolds for the discovery, design, and development of new anti-TB drugs.

Rethinking the MtInhA tertiary and quaternary structure flexibility: a molecular dynamics view

Flexibility and function are related properties in the study of protein dynamics. Flexibility reflects in the conformational potential of proteins and thus in their functionalities. The presence of interactions between protein-ligands and protein-protein complexes, substrates, and environmental changes can alter protein plasticity, acting from the rearrangement of the side chains of amino acids to the folding/unfolding of large structural motifs. To evaluate the effects of the flexibility in protein systems, we defined the enzyme 2-trans-enoyl-ACP (CoA) reductase from Mycobacterium tuberculosis, or MtInhA, as our target system. MtInhA is biologically active as a tetramer in solution; however, computational studies commonly use the monomer justifying the independence of its active sites due to their distances. However, differences in flexibility between tertiary and quaternary structures could present impact on the size of the active site, influencing the drug discovery process. In this study, we investigated the influence of flexibility restrictions in A- and B-loops of the MtInhA in order to suggest a monomeric structure that describes the conformational behavior of the tetrameric system. Overall, we observed that simulations where restrictions were applied to the A- and B-loops present a more similar behavior to the native structure when compared to unrestricted simulations. Therefore, our work presents a monomeric model of MtInhA, which has conformational characteristics of the biologically active structure. Thus, the data obtained in this work can be applied to the MtInhA system for the generation of more reliable flexible models for molecular docking experiments, and also for the performance of longer simulations by molecular dynamics and with a lower computational cost.

Molecular recognition and binding of CcrA from Bacteroides fragilis with cefotaxime and ceftazidime by fluorescence spectra and molecular docking

In search of new super-bacterial inhibitor agents, the recognition and binding mechanism of the B1 subclass MβL CcrA from Bacteroides fragilis with cefotaxime (CTX) and ceftazidime (CAZ) were studied using spectroscopy analysis and molecular docking. The results showed that the fluorescence quenching of CcrA induced by CTX and CAZ were all due to the complex formation, which belonged to static quenching and was forced by hydrogen bonds and Van der Waals forces, despite the greater binding ability of CTX with CcrA than CAZ. Upon recognizing CTX or CAZ, the CcrA opened its binding pocket by the microenvironmental and conformational of three loops changing to promote an induced-fit of the freshly introduced antibiotics. In addition, the whole antibiotic molecule ultimately entered the active pocket of CcrA with its original carbonate replaced by the carboxyl oxygen of the hexatomic ring adjacent to the β-lactam ring in CTX or CAZ, forming a new tetrahedral coordination structure at the Zn2 site. Moreover, the difference in steric hindrance and electrostatic effects of the side chain affected the binding ability of the two antibiotics to the CcrA. This work showed the refined procedures of the antibiotics binding to CcrA and might provide useful information hint for the new strategy of developing the novel and innovative super-bacterial antibiotics.