RESUMEN
BACKGROUND: Lung cancer is the leading cause of cancer-related death worldwide. Among its subtypes, non-small cell lung cancer (NSCLC) is the most common. Recently, the mitochondrial isoform of the enzyme phosphoenolpyruvate carboxykinase (HsPEPCK-M) was identified as responsible for the metabolic adaptation in the NSCLC allowing tumor growth even under conditions of glucose deficiency. This adaptation is possible due to the role of HsPEPCK-M in gluconeogenesis, converting the oxaloacetate to phosphoenolpyruvate in the presence of GTP, which plays an important role in the energetic support of these tumors. In this context, it was shown that the inhibition or knockdown of this enzyme was able to induce apoptosis in NSCLC under low glucose conditions. PURPOSE: In this study, novel putative inhibitors were proposed for the human PEPCK-M (HsPEPCK-M) based on a computer-aided approach. METHODS: Comparative modeling was used to generate 3D models for HsPEPCK-M. Subsequently, the set of natural compounds of the ZINC database was screened against HsPEPCK-M models using structure-based pharmacophore modeling and molecular docking approaches. The selected compounds were evaluated according to its chemical diversity and clustered based on chemical similarity. RESULTS: The pharmacophore hypotheses, generated based on known PEPCK inhibitors, were able to select 7,124 candidate compounds. These compounds were submitted to molecular docking studies using three conformations of HsPEPCK-M generated by comparative modeling. The aim was to select compounds with high predicted binding affinity for at least one of the conformations of HsPEPCK-M. After molecular docking, 612 molecules were selected as potential inhibitors of HsPEPCK-M. These compounds were clustered according to their structural similarity. Chemical profiling and binding mode analyses of these compounds allowed the proposal of four promising compounds: ZINC01656421, ZINC895296, ZINC00895535 and ZINC02571340. CONCLUSION: These compounds may be considered as potential candidates for HsPEPCK-M inhibitors and may also be used as lead compounds for the development of novel HsPEPCK-M inhibitors.
RESUMEN
BACKGROUND: Human TWIST1 is a highly conserved member of the regulatory basic helix-loop-helix (bHLH) transcription factors. TWIST1 forms homo- or heterodimers with E-box proteins, such as E2A (isoforms E12 and E47), MYOD and HAND2. Haploinsufficiency germ-line mutations of the twist1 gene in humans are the main cause of Saethre-Chotzen syndrome (SCS), which is characterized by limb abnormalities and premature fusion of cranial sutures. Because of the importance of TWIST1 in the regulation of embryonic development and its relationship with SCS, along with the lack of an experimentally solved 3D structure, we performed comparative modeling for the TWIST1 bHLH region arranged into wild-type homodimers and heterodimers with E47. In addition, three mutations that promote DNA binding failure (R118C, S144R and K145E) were studied on the TWIST1 monomer. We also explored the behavior of the mutant forms in aqueous solution using molecular dynamics (MD) simulations, focusing on the structural changes of the wild-type versus mutant dimers. RESULTS: The solvent-accessible surface area of the homodimers was smaller on wild-type dimers, which indicates that the cleft between the monomers remained more open on the mutant homodimers. RMSD and RMSF analyses indicated that mutated dimers presented values that were higher than those for the wild-type dimers. For a more careful investigation, the monomer was subdivided into four regions: basic, helix I, loop and helix II. The basic domain presented a higher flexibility in all of the parameters that were analyzed, and the mutant dimer basic domains presented values that were higher than the wild-type dimers. The essential dynamic analysis also indicated a higher collective motion for the basic domain. CONCLUSIONS: Our results suggest the mutations studied turned the dimers into more unstable structures with a wider cleft, which may be a reason for the loss of DNA binding capacity observed for in vitro circumstances.