Structure-Based Virtual Screening of High-Affinity ATP-Competitive Inhibitors Against Human Lemur Tyrosine Kinase-3 (LMTK3) Domain: A Novel Therapeutic Target for Breast Cancer

Insights to Design New Drugs against Human African Trypanosomiasis Targeting Rhodesain using Covalent Docking, Molecular Dynamics Simulations, and MM-PBSA Calculations

BACKGROUND

Neglected tropical diseases (NTDs) are parasitic and bacterial diseases that affect approximately 149 countries, mainly the poor population without basic sanitation. Among these, Human African Trypanosomiasis (HAT), known as sleeping sickness, shows alarming data, with treatment based on suramin and pentamidine in the initial phase and melarsoprol and eflornithine in the chronic phase. Thus, to discover new drugs, several studies point to rhodesain as a promising drug target due to the function of protein degradation and intracellular transport of proteins between the insect and host cells and is present in all cycle phases of the parasite.

METHODS

Here, based on the previous studies by Nascimento et al. (2021) [5], that show the main rhodesain inhibitors development in the last decade, molecular docking and dynamics were applied in these inhibitors datasets to reveal crucial information that can be into drug design.

RESULTS

Also, our findings using MD simulations and MM-PBSA calculations confirmed Gly19, Gly23, Gly65, Asp161, and Trp184, showing high binding energy (ΔGbind between -72.782 to -124.477 kJ.mol-1). In addition, Van der Waals interactions have a better contribution (-140,930 to -96,988 kJ.mol-1) than electrostatic forces (-43,270 to -6,854 kJ.mol-1), indicating Van der Waals interactions are the leading forces in forming and maintaining ligand-rhodesain complexes. Thus, conventional and covalent docking was employed and highlighted the presence of Michael acceptors in the ligands in a peptidomimetics scaffold, and interaction with Gly19, Gly23, Gly65, Asp161, and Trp184 is essential to the inhibiting activity. Furthermore, the Dynamic Cross-Correlation Maps (DCCM) show more correlated movements for all complexes than the free rhodesain and strong interactions in the regions of the aforementioned residues. Principal Component Analysis (PCA) demonstrates complex stability corroborating with RMSF and RMSD.

CONCLUSION

This study can provide valuable insights that can guide researchers worldwide to discover a new promising drug against HAT.

Data mining and molecular dynamics analysis to detect HIV-1 reverse transcriptase RNase H activity inhibitor

HIV-1 is a deadly virus that affects millions of people worldwide. In this study, we aimed to inhibit viral replication by targeting one of the HIV-1 proteins and identifying a new drug candidate. We used data mining and molecular dynamics methods on HIV-1 genomes. Based on MAUVE analysis, we selected the RNase H activity of the reverse transcriptase (R.T) enzyme as a potential target due to its low mutation rate and high conservation level. We screened about 94,000 small molecule inhibitors by virtual screening. We validated the hit compounds' stability and binding free energy through molecular dynamics simulations and MM/PBSA. Phomoarcherin B, known for its anticancer properties, emerged as the best candidate and showed potential as an HIV-1 reverse transcriptase RNase H activity inhibitor. This study presents a new target and drug candidate for HIV-1 treatment. However, in vitro and in vivo tests are required. Also, the effect of RNase H activity on viral replication and the interaction of Phomoarcherin B with other HIV-1 proteins should be investigated.

The influence of N-alkyl chains in benzoyl-thiourea derivatives on urease inhibition: Soil studies and biophysical and theoretical investigations on the mechanism of interaction

Ureases are enzymes produced by fungi, plants, and bacteria associated with agricultural and clinical problems. The urea hydrolysis in NH₃ and CO₂ leads to the loss of N-urea fertilizers in soils and changes the human stomach microenvironment, favoring the colonization of H. pylori. In this sense, it is necessary to evaluate potential enzyme inhibitors to mitigate the effects of their activities and respond to scientific and market demands to produce fertilizers with enhanced efficiency. Thus, biophysical and theoretical studies were carried out to evaluate the influence of the N-alkyl chain in benzoyl-thiourea derivatives on urease enzyme inhibition. A screening based on IC₅₀, binding constants, and theoretical studies demonstrated that BTU1 without the N-alkyl chain (R = H) was more active than other compounds, so the magnitude of the interaction was determined as BTU1 > BTU2 > BTU3 > BTU4 > BTU5, corresponding to progressively increased chain length. Thus, BTU1 was selected for interaction and soil application essays. The binding constants (K_b) for the supramolecular urease-BTU1 complex ranged from 7.95 to 5.71 × 10³ M^-1 at different temperatures (22, 30, and 38 °C), indicating that the preferential forces responsible for the stabilization of the complex are hydrogen bonds and van der Waals forces (ΔH = -15.84 kJ mol^-1 and ΔS = -36.61 J mol^-1 K^-1). Theoretical and experimental results (thermodynamics, synchronous fluorescence, and competition assay) agree and indicate that BTU1 is a mixed inhibitor. Finally, urease inhibition was evaluated in the four soil samples, where BTU1 was as efficient as NBPT (based on ANOVA two-way and Tukey test with 95% confidence), with an average inhibition of 20% of urease activity. Thus, the biophysics and theoretical studies are strategies for evaluating potential inhibitors and showed that increasing the N-alkyl chain in benzoyl-thiourea derivatives did not favor urease inhibition.

ACW-02 an Acridine Triazolidine Derivative Presents Antileishmanial Activity Mediated by DNA Interaction and Immunomodulation

The present study proposed the synthesis of a novel acridine derivative not yet described in the literature, chemical characterization by NMR, MS, and IR, followed by investigations of its antileishmanial potential. In vitro assays were performed to assess its antileishmanial activity against L. amazonensis strains and cytotoxicity against macrophages through MTT assay and annexin V-FITC/PI, and the ability to perform an immunomodulatory action using CBA. To investigate possible molecular targets, its interaction with DNA in vitro and in silico targets were evaluated. As results, the compound showed good antileishmanial activity, with IC50 of 6.57 (amastigotes) and 94.97 (promastigotes) µg mL-1, associated with non-cytotoxicity to macrophages (CC50 > 256.00 µg mL-1). When assessed by flow cytometry, 99.8% of macrophages remained viable. The compound induced an antileishmanial effect in infected macrophages and altered TNF-α, IL-10 and IL-6 expression, suggesting a slight immunomodulatory activity. DNA assay showed an interaction with the minor grooves due to the hyperchromic effect of 47.53% and Kb 1.17 × 106 M-1, and was sustained by docking studies. Molecular dynamics simulations and MM-PBSA calculations propose cysteine protease B as a possible target. Therefore, this study demonstrates that the new compound is a promising molecule and contributes as a model for future works.

Repurposing FDA-approved Drugs Targeting SARS-CoV2 3CLpro: a study by applying Virtual Screening, Molecular Dynamics, MM-PBSA Calculations and Covalent Docking

Background:

Since the end of 2019, the etiologic agent SAR-CoV-2 responsible for one of the most significant epidemics in history has caused severe global economic, social, and health damages. The drug repurposing approach and application of Structure-based Drug Discovery (SBDD) using in silico techniques are increasingly frequent, leading to the identification of several molecules that may represent promising potential.

Method:

In this context, here we use in silico methods of virtual screening (VS), pharmacophore modeling (PM), and fragment-based drug design (FBDD), in addition to molecular dynamics (MD), molecular mechanics/Poisson-Boltzmann surface area (MM -PBSA) calculations, and covalent docking (CD) for the identification of potential treatments against SARS-CoV-2. We initially validated the docking protocol followed by VS in 1,613 FDA-approved drugs obtained from the ZINC database. Thus, we identified 15 top hits, of which three of them were selected for further simulations. In parallel, for the compounds with a fit score value ≤ of 30, we performed the FBDD protocol, where we designed 12 compounds

Result:

By applying a PM protocol in the ZINC database, we identified three promising drug candidates. Then, the 9 top hits were evaluated in simulations of MD, MM-PBSA, and CD. Subsequently, MD showed that all identified hits showed stability at the active site without significant changes in the protein's structural integrity, as evidenced by the RMSD, RMSF, Rg, SASA graphics. They also showed interactions with the catalytic dyad (His41 and Cys145) and other essential residues for activity (Glu166 and Gln189) and high affinity for MM-PBSA, with possible covalent inhibition mechanism.

Conclution:

Finally, our protocol helped identify potential compounds wherein ZINC896717 (Zafirlukast), ZINC1546066 (Erlotinib), and ZINC1554274 (Rilpivirine) were more promising and could be explored in vitro, in vivo, and clinical trials to prove their potential as antiviral agents.

The multifaceted role of lemur tyrosine kinase 3 in health and disease

In the last decade, LMTK3 (lemur tyrosine kinase 3) has emerged as an important player in breast cancer, contributing to the advancement of disease and the acquisition of resistance to therapy through a strikingly complex set of mechanisms. Although the knowledge of its physiological function is largely limited to receptor trafficking in neurons, there is mounting evidence that LMTK3 promotes oncogenesis in a wide variety of cancers. Recent studies have broadened our understanding of LMTK3 and demonstrated its importance in numerous signalling pathways, culminating in the identification of a potent and selective LMTK3 inhibitor. Here, we review the roles of LMTK3 in health and disease and discuss how this research may be used to develop novel therapeutics to advance cancer treatment.

Computer-aided design of 1,4-naphthoquinone-based inhibitors targeting cruzain and rhodesain cysteine proteases

Chagas disease and Human African Trypanosomiasis (HAT) are caused by Trypanosoma cruzi and T. brucei parasites, respectively. Cruzain (CRZ) and Rhodesain (RhD) are cysteine proteases that share 70% of identity and play vital functions in these parasites. These macromolecules represent promising targets for designing new inhibitors. In this context, 26 CRZ and 5 RhD 3D-structures were evaluated by molecular redocking to identify the most accurate one to be utilized as a target. Posteriorly, a virtual screening of a library containing 120 small natural and nature-based compounds was performed on both of them. In total, 14 naphthoquinone-based analogs were identified, synthesized, and biologically evaluated. In total, five compounds were active against RhD, being three of them also active on CRZ. A derivative of 1,4-naphthoquinonepyridin-2-ylsulfonamide was found to be the most active molecule, exhibiting IC₅₀ values of 6.3 and 1.8 µM for CRZ and RhD, respectively. Dynamic simulations at 100 ns demonstrated good stability and do not alter the targets' structures. MM-PBSA calculations revealed that it presents a higher affinity for RhD (-25.3 Kcal mol^-1) than CRZ, in which van der Waals interactions were more relevant. A mechanistic hypothesis (via C3-Michael-addition reaction) involving a covalent mode of inhibition for this compound towards RhD was investigated by covalent molecular docking and DFT B3LYP/6-31 + G* calculations, exhibiting a low activation energy (ΔG^‡) and providing a stable product (ΔG), with values of 7.78 and - 39.72 Kcal mol^-1, respectively; similar to data found in the literature. Nevertheless, a reversibility assay by dilution revealed that JN-11 is a time-dependent and reversible inhibitor. Finally, this study applies modern computer-aided techniques to identify promising inhibitors from a well-known chemical class of natural products. Then, this work could inspire other future studies in the field, being useful for designing potent naphthoquinones as RhD inhibitors.

Comparative study on the interaction of oxyresveratrol and piceatannol with trypsin and lysozyme: binding ability, activity and stability

Natural polyphenols showing a variety of beneficial effects will interact with multiple proteases after administration. The interactions of oxyresveratrol and piceatannol with trypsin and lysozyme were investigated using fluorescence spectroscopy, UV-vis absorption spectroscopy, circular dichroism spectroscopy, differential scanning calorimetry and molecular docking. Fluorescence quenching results and UV-vis absorption difference spectra revealed that the quenching process was a static mode initiated by ground-state complex formation. The different binding ability of oxyresveratrol and piceatannol with trypsin and lysozyme was discussed based on their different molecular structures. Moreover, the major driving force for the binding process was elucidated as hydrogen bonding and van der Waals forces by the negative enthalpy and entropy changes. Synchronous fluorescence, three-dimensional fluorescence and circular dichroism spectral analysis suggested that the binding of oxyresveratrol and piceatannol to trypsin and lysozyme induced some microenvironmental and conformational changes of the two enzymes. The thermal stability of the enzymes in the presence of polyphenols was studied based on the change in melting temperature by differential scanning calorimetry. The above experimental results were validated by the protein-ligand docking studies which showed the location of the two ligands in the enzymes and the surrounding amino acid residues. Furthermore, enzyme activity assays indicated that the enzymatic activity of trypsin and lysozyme was inhibited by oxyresveratrol and piceatannol. The effect of trypsin and lysozyme on the antioxidant activity and stability of oxyresveratrol and piceatannol was also investigated. In conclusion, the comparative study on the interaction of oxyresveratrol and piceatannol with trypsin and lysozyme showed that the positions of hydroxyl groups of the polyphenols had an important influence on their interaction with enzymes and their antioxidant activity and stability as well as the enzyme activities. The obtained results are expected to provide a theoretical basis for the application of polyphenols in functional foods and pharmaceuticals.

The lncRNA MIR2052HG regulates ERα levels and aromatase inhibitor resistance through LMTK3 by recruiting EGR1

Background

Our previous genome-wide association study using the MA.27 aromatase inhibitors adjuvant trial identified SNPs in the long noncoding RNA MIR2052HG associated with breast cancer-free interval. MIR2052HG maintained ERα both by promoting AKT/FOXO3-mediated ESR1 transcription and by limiting ubiquitin-mediated ERα degradation. Our goal was to further elucidate MIR2052HG’s mechanism of action.

Methods

RNA-binding protein immunoprecipitation assays were performed to demonstrate that the transcription factor, early growth response protein 1 (EGR1), worked together with MIR2052HG to regulate that lemur tyrosine kinase-3 (LMTK3) transcription in MCF7/AC1 and CAMA-1 cells. The location of EGR1 on the LMTK3 gene locus was mapped using chromatin immunoprecipitation assays. The co-localization of MIR2052HG RNA and the LMTK3 gene locus was determined using RNA-DNA dual fluorescent in situ hybridization. Single-nucleotide polymorphisms (SNP) effects were evaluated using a panel of human lymphoblastoid cell lines.

Results

MIR2052HG depletion in breast cancer cells results in a decrease in LMTK3 expression and cell growth. Mechanistically, MIR2052HG interacts with EGR1 and facilitates its recruitment to the LMTK3 promoter. LMTK3 sustains ERα levels by reducing protein kinase C (PKC) activity, resulting in increased ESR1 transcription mediated through AKT/FOXO3 and reduced ERα degradation mediated by the PKC/MEK/ERK/RSK1 pathway. MIR2052HG regulated LMTK3 in a SNP- and aromatase inhibitor-dependent fashion: the variant SNP increased EGR1 binding to LMTK3 promoter in response to androstenedione, relative to wild-type genotype, a pattern that can be reversed by aromatase inhibitor treatment. Finally, LMTK3 overexpression abolished the effect of MIR2052HG on PKC activity and ERα levels.

Conclusions

Our findings support a model in which the MIR2052HG regulates LMTK3 via EGR1, and LMTK3 regulates ERα stability via the PKC/MEK/ERK/RSK1 axis. These results reveal a direct role of MIR2052HG in LMTK3 regulation and raise the possibilities of targeting MIR2052HG or LMTK3 in ERα-positive breast cancer.

Electronic supplementary material

The online version of this article (10.1186/s13058-019-1130-3) contains supplementary material, which is available to authorized users.