RESUMO
The druggable proteome refers to proteins that can bind to small molecules with appropriate chemical affinity, inducing a favorable clinical response. Predicting druggable proteins through screening and in silico modeling is imperative for drug design. To contribute to this field, we developed an accurate predictive classifier for druggable cancer-driving proteins using amino acid composition descriptors of protein sequences and 13 machine learning linear and non-linear classifiers. The optimal classifier was achieved with the support vector machine method, utilizing 200 tri-amino acid composition descriptors. The high performance of the model is evident from an area under the receiver operating characteristics (AUROC) of 0.975 ± 0.003 and an accuracy of 0.929 ± 0.006 (threefold cross-validation). The machine learning prediction model was enhanced with multi-omics approaches, including the target-disease evidence score, the shortest pathways to cancer hallmarks, structure-based ligandability assessment, unfavorable prognostic protein analysis, and the oncogenic variome. Additionally, we performed a drug repurposing analysis to identify drugs with the highest affinity capable of targeting the best predicted proteins. As a result, we identified 79 key druggable cancer-driving proteins with the highest ligandability, and 23 of them demonstrated unfavorable prognostic significance across 16 TCGA PanCancer types: CDKN2A, BCL10, ACVR1, CASP8, JAG1, TSC1, NBN, PREX2, PPP2R1A, DNM2, VAV1, ASXL1, TPR, HRAS, BUB1B, ATG7, MARK3, SETD2, CCNE1, MUTYH, CDKN2C, RB1, and SMARCA4. Moreover, we prioritized 11 clinically relevant drugs targeting these proteins. This strategy effectively predicts and prioritizes biomarkers, therapeutic targets, and drugs for in-depth studies in clinical trials. Scripts are available at https://github.com/muntisa/machine-learning-for-druggable-proteins .
Assuntos
Inteligência Artificial , Neoplasias , Humanos , Neoplasias/tratamento farmacológico , Neoplasias/genética , Neoplasias/metabolismo , Antineoplásicos/farmacologia , Antineoplásicos/uso terapêutico , Antineoplásicos/química , Aprendizado de Máquina , Proteínas de Neoplasias/metabolismo , Proteínas de Neoplasias/genética , Proteínas de Neoplasias/química , Máquina de Vetores de Suporte , Reposicionamento de Medicamentos/métodos , Biologia Computacional/métodos , MultiômicaRESUMO
Introduction: Coronavirus Disease 2019 (COVID-19) is a severe respiratory illness caused by the RNA virus SARS-CoV-2. Globally, there have been over 759.4 million cases and 6.74 million deaths, while Ecuador has reported more than 1.06 million cases and 35.9 thousand deaths. To describe the COVID-19 pandemic impact and the vaccinations effectiveness in a low-income country like Ecuador, we aim to assess the seroprevalence of IgG and IgM antibodies against SARS-CoV-2 in a sample from healthy blood donors at the Cruz Roja Ecuatoriana. Methods: The present seroprevalence study used a lateral flow immunoassay (LFIA) to detect anti-SARS-CoV-2 IgG and IgM antibodies in months with the highest confirmed case rates (May 2020; January, April 2021; January, February, June, July 2022) and months with the highest vaccination rates (May, June, July, August, December 2021) in Quito, Ecuador. The IgG and IgM seroprevalence were also assessed based on sex, age range, blood type and RhD antigen type. The sample size was 8,159, and sampling was performed based on the availability of each blood type. Results: The results showed an overall IgG and IgM seroprevalence of 47.76% and 3.44%, respectively. There were no differences in IgG and IgM seroprevalences between blood groups and sex, whereas statistical differences were found based on months, age range groups, and RhD antigen type. For instance, the highest IgG seroprevalence was observed in February 2022 and within the 17-26 years age range group, while the highest IgM seroprevalence was in April 2021 and within the 47-56 years age range group. Lastly, only IgG seroprevalence was higher in RhD+ individuals while IgM seroprevalence was similar across RhD types. Discussion: This project contributes to limited data on IgG and IgM antibodies against SARS-CoV-2 in Ecuador. It suggests that herd immunity may have been achieved in the last evaluated months, and highlights a potential link between the RhD antigen type and COVID-19 susceptibility. These findings have implications for public health strategies and vaccine distribution not only in Ecuador but also in regions with similar characteristics.
Assuntos
Anticorpos Antivirais , Doadores de Sangue , COVID-19 , Imunoglobulina G , Imunoglobulina M , SARS-CoV-2 , Humanos , COVID-19/epidemiologia , COVID-19/imunologia , Equador/epidemiologia , Imunoglobulina G/sangue , Estudos Soroepidemiológicos , Imunoglobulina M/sangue , Masculino , SARS-CoV-2/imunologia , Adulto , Doadores de Sangue/estatística & dados numéricos , Anticorpos Antivirais/sangue , Feminino , Pessoa de Meia-Idade , Adolescente , Adulto Jovem , Idoso , PandemiasRESUMO
Introduction: Gastric cancer is one of the most prevalent types of cancer worldwide. The World Health Organization (WHO), the International Agency for Research on Cancer (IARC), and the Global Cancer Statistics (GLOBOCAN) reported an age standardized global incidence rate of 9.2 per 100,000 individuals for gastric cancer in 2022, with a mortality rate of 6.1. Despite considerable progress in precision oncology through the efforts of international consortia, understanding the genomic features and their influence on the effectiveness of anti-cancer treatments across diverse ethnic groups remains essential. Methods: Our study aimed to address this need by conducting integrated in silico analyses to identify actionable genomic alterations in gastric cancer driver genes, assess their impact using deleteriousness scores, and determine allele frequencies across nine global populations: European Finnish, European non-Finnish, Latino, East Asian, South Asian, African, Middle Eastern, Ashkenazi Jewish, and Amish. Furthermore, our goal was to prioritize targeted therapeutic strategies based on pharmacogenomics clinical guidelines, in silico drug prescriptions, and clinical trial data. Results: Our comprehensive analysis examined 275,634 variants within 60 gastric cancer driver genes from 730,947 exome sequences and 76,215 whole-genome sequences from unrelated individuals, identifying 13,542 annotated and predicted oncogenic variants. We prioritized the most prevalent and deleterious oncogenic variants for subsequent pharmacogenomics testing. Additionally, we discovered actionable genomic alterations in the ARID1A, ATM, BCOR, ERBB2, ERBB3, CDKN2A, KIT, PIK3CA, PTEN, NTRK3, TP53, and CDKN2A genes that could enhance the efficacy of anti-cancer therapies, as suggested by in silico drug prescription analyses, reviews of current pharmacogenomics clinical guidelines, and evaluations of phase III and IV clinical trials targeting gastric cancer driver proteins. Discussion: These findings underline the urgency of consolidating efforts to devise effective prevention measures, invest in genomic profiling for underrepresented populations, and ensure the inclusion of ethnic minorities in future clinical trials and cancer research in developed countries.
RESUMO
Purpose: Thiopurine S-methyltransferase (TPMT) is an enzyme that metabolizes purine analogs, agents used in the treatment of acute lymphoblastic leukemia. Improper drug metabolism leads to toxicity in chemotherapy patients and reduces treatment effectiveness. TPMT variants associated with reduced enzymatic activity vary across populations. Therefore, studying these variants in heterogeneous populations, such as Ecuadorians, can help identify molecular causes of deficiency for this enzyme. Methods: We sequenced the entire TPMT coding region in 550 Ecuadorian individuals from Afro-Ecuadorian, Indigenous, Mestizo, and Montubio ethnicities. Moreover, we conducted an ancestry analysis using 46 informative ancestry markers. Results: We identified 8 single nucleotide variants in the coding region of TPMT. The most prevalent alleles were TPMT*3A, TPMT*3B, and TPMT*3C, with frequencies of 0.055, 0.012, and 0.015, respectively. Additionally, we found rare alleles TPMT*4 and TPMT*8 with frequencies of 0.005 and 0.003. Correlating the ancestry proportions with TPMT-deficient genotypes, we observed that the Native American ancestry proportion influenced the distribution of the TPMT*1/TPMT*3A genotype (OR = 5.977, p = 0.002), while the contribution of African ancestral populations was associated with the TPMT*1/TPMT*3C genotype (OR = 9.769, p = 0.003). The rates of TPMT-deficient genotypes observed in Mestizo (f = 0.121) and Indigenous (f = 0.273) groups provide evidence for the influence of Native American ancestry and the prevalence of the TPMT*3A allele. In contrast, although Afro-Ecuadorian groups demonstrate similar deficiency rates (f = 0.160), the genetic factors involved are associated with contributions from African ancestral populations, specifically the prevalent TPMT*3C allele. Conclusion: The distribution of TPMT-deficient variants offers valuable insights into the populations under study, underscoring the necessity for genetic screening strategies to prevent thiopurine toxicity events among Latin American minority groups.
RESUMO
BACKGROUND: Glucose-6-phosphate dehydrogenase deficiency (G6PDd) is an X-linked disorder affecting over 400 million people worldwide. Individuals with molecular variants associated with reduced enzymatic activity are susceptible to oxidative stress in red blood cells, thereby increasing the risk of pathophysiological conditions and toxicity to anti-malarial treatments. Globally, the prevalence of G6PDd varies among populations. Accordingly, this study aims to characterize G6PDd distribution within the Ecuadorian population and to describe the spatial distribution of reported malaria cases. METHODS: Molecular variants associated with G6PDd were genotyped in 581 individuals from Afro-Ecuadorian, Indigenous, Mestizo, and Montubio ethnic groups. Additionally, spatial analysis was conducted to identify significant malaria clusters with high incidence rates across Ecuador, using data collected from 2010 to 2021. RESULTS: The A- c.202G > A and A- c.968T > C variants underpin the genetic basis of G6PDd in the studied population. The overall prevalence of G6PDd was 4.6% in the entire population. However, this frequency increased to 19.2% among Afro-Ecuadorian people. Spatial analysis revealed 12 malaria clusters, primarily located in the north of the country and its Amazon region, with relative risks of infection of 2.02 to 87.88. CONCLUSIONS: The findings of this study hold significant implications for public health interventions, treatment strategies, and targeted efforts to mitigate the burden of malaria in Ecuador. The high prevalence of G6PDd among Afro-Ecuadorian groups in the northern endemic areas necessitates the development of comprehensive malaria eradication strategies tailored to this geographical region.
Assuntos
Deficiência de Glucosefosfato Desidrogenase , Malária , Humanos , Equador/epidemiologia , Eritrócitos , Etnicidade , Deficiência de Glucosefosfato Desidrogenase/epidemiologia , Deficiência de Glucosefosfato Desidrogenase/genética , Malária/epidemiologiaRESUMO
Dihydropyrimidine dehydrogenase is one of the main pharmacological metabolizers of fluoropyrimidines, a group of drugs widely used in clinical oncology. Around 20 to 30% of patients treated with fluoropyrimidines experience severe toxicity caused by a partial or total decrease in enzymatic activity. This decrease is due to molecular variants in the DPYD gene. Their prevalence and allelic frequencies vary considerably worldwide, so their description in heterogeneous groups such as the Ecuadorian population will allow for the description of pharmacogenetic variants and proper characterization of this population. Thus, we genotyped all the molecular variants with a predictive value for DPYD in a total of 410 Ecuadorian individuals belonging to Mestizo, Afro-Ecuadorian, and Indigenous ethnic groups. Moreover, we developed a genetic ancestry analysis using 46 autosomal ancestry informative markers. We determined 20 genetic variations in 5 amplified regions, including 3 novel single nucleotide variants. The allele frequencies for DPYD variants c.1627G>A (*5, rs1801159), c.1129-15T>C (rs56293913), c.1218G>A (rs61622928), rs1337752, rs141050810, rs2786783, rs2811178, and g.97450142G>A (chr1, GRCh38.p13) are significantly related to Native American and African ancestry proportions. In addition, the FST calculated from these variants demonstrates the closeness between Indigenous and Mestizo populations, and evidences genetic divergence between Afro-Ecuadorian groups when compared with Mestizo and Indigenous ethnic groups. In conclusion, the genetic variability in the DPYD gene is related to the genetic component of ancestral populations in different Ecuadorian ethnic groups. The absence and low frequency of variants with predictive value for fluoropyrimidine toxicity such as DPYD *2A, HapB3, and c.2846A>T (prevalent in populations with European ancestry) is consistent with the genetic background found.
RESUMO
Background: It is imperative to identify drugs that allow treating symptoms of severe COVID-19. Respiratory failure is the main cause of death in severe COVID-19 patients, and the host inflammatory response at the lungs remains poorly understood. Methods: Therefore, we retrieved data from post-mortem lungs from COVID-19 patients and performed in-depth in silico analyses of single-nucleus RNA sequencing data, inflammatory protein interactome network, and shortest pathways to physiological phenotypes to reveal potential therapeutic targets and drugs in advanced-stage COVID-19 clinical trials. Results: Herein, we analyzed transcriptomics data of 719 inflammatory response genes across 19 cell types (116,313 nuclei) from lung autopsies. The functional enrichment analysis of the 233 significantly expressed genes showed that the most relevant biological annotations were inflammatory response, innate immune response, cytokine production, interferon production, macrophage activation, blood coagulation, NLRP3 inflammasome complex, and the TLR, JAK-STAT, NF-κB, TNF, oncostatin M signaling pathways. Subsequently, we identified 34 essential inflammatory proteins with both high-confidence protein interactions and shortest pathways to inflammation, cell death, glycolysis, and angiogenesis. Conclusion: We propose three small molecules (baricitinib, eritoran, and montelukast) that can be considered for treating severe COVID-19 symptoms after being thoroughly evaluated in COVID-19 clinical trials.
RESUMO
Batrachochytrium dendrobatidis (Bd) infection is one of the principal causes of amphibian declines worldwide. The presence of Bd has been determined in Gastrotheca riobambae tadpoles that inhabit ponds in Quito's Metropolitan Guangüiltagua Park, Ecuador. This study sought to determine whether these tadpoles are infected and to determine the presence of chytridiomycosis in another frog species, Pristimantis unistrigatus, which also inhabits the park and has different reproductive biology and distinct behavioral habits. We used end-point and real-time PCR techniques to detect and quantify Bd infection. At 1 yr, samples were taken from the skin of P. unistrigatus using swabs and were also taken from the mouthparts of G. riobambae tadpoles. It was found that the two species were infected with a Bd prevalence of 39% (53/135) in G. riobambae tadpoles and 15% (57/382) in P. unistrigatus frogs. The two types of samples (tissue and swabs) from mouthparts showed differences in the zoospores per microliter loads (xÌ=1,376.7±3,450.2 vs. xÌ=285.0±652.3). Moreover, a correlation (r2=0.621) was discovered between the monthly mean maximum temperature of the pond with disease prevalence in G. riobambae tadpoles. Infection levels in the P. unistrigatus population varied significantly over time, and distance to the pond was a determinant factor for infection intensity.
Assuntos
Quitridiomicetos , Micoses , Animais , Anuros , Batrachochytrium , Equador/epidemiologia , Micoses/epidemiologia , Micoses/veterináriaRESUMO
OBJECTIVES: According to demographic history, Ecuador has experienced shifts in its Native American populations caused by European colonization and the African slave trade. The continuous admixture events among Europeans, Native Americans, and Africans occurred differently in each region of the country, producing a stratified population. Thus, the aim of this study was to investigate the level of genetic substructure in the Ecuadorian Mestizo population. MATERIALS AND METHODS: A total of 377 male and 209 female samples were genotyped for two sets of X-chromosomal markers (32 X-Indels and 12 X-STRs). Population analyses performed included Hardy-Weinberg equilibrium tests, LD analysis, PCA, pairwise FST s, and AMOVA. RESULTS: Significant levels of LD were observed between markers separated by distances of less than 1 cM, as well as between markers separated by distances varying from 10.891 to 163.53 cM. Among Ecuadorian regions, Amazonia showed the highest average R2 value. DISCUSSION: When X-chromosomal and autosomal differentiation values were compared, a sex-biased admixture between European men and Native American and African women was revealed, as well as between African men and Native American women. Moreover, a distinct Native American ancestry was discernible in the Amazonian population, in addition to sex-biased gene flow between Amazonia and the Andes and Pacific coast regions. Overall, these results underline the importance of integrating X chromosome information to achieve a more comprehensive view of the genetic and demographic histories of South American admixed populations.
Assuntos
Variação Genética/genética , Genética Populacional/métodos , Indígenas Sul-Americanos/genética , Antropologia Física , Cromossomos Humanos X/genética , Equador , Feminino , Humanos , Mutação INDEL/genética , Desequilíbrio de Ligação/genética , Masculino , Repetições de Microssatélites/genéticaRESUMO
BACKGROUND: The 15q11.1-13.1 duplication, also known as Dup15q syndrome, is a rare congenital disease affecting 1 in 30,000 to 1 in 60,000 children worldwide. This condition is characterized by the presence of at least one extra copy of genetical material within the Prader-Willi/Angelman Critical Region (PWACR) of the referred 15q11.2-q13.1 chromosome. Case Report. Our study presents the clinical and genetical features of the first patient with a denovo 15q11.2 interstitial duplication on the maternal allele (inv Dup15q) that mimics a milder Prader-Willi syndrome probably due to an atypical disruption of the SNHG14 gene. Methylation-specific MLPA analysis has confirmed the presence of a very unlikely duplication that lies between breakpoint 1 (BP1) and the middle of BP2 and BP3 (BP3). This atypical alteration might be linked to the milder patient's clinical phenotype. CONCLUSIONS: This is the first Dup15q patient reported in Ecuador and of the very few in South America. This aberration has never been described in a patient with Dup15q, and the unusual clinical presentation is probably due to the atypical distal breakpoint occurring within the gene SNHG14 which lies between BP2 and BP3 and does not therefore contain the whole PWACR. If the duplication disrupted the gene, then it is possible that it is the cause of, or contributing to, the patient's clinical phenotype.
RESUMO
Background: There is pressing urgency to identify therapeutic targets and drugs that allow treating COVID-19 patients effectively. Methods: We performed in silico analyses of immune system protein interactome network, single-cell RNA sequencing of human tissues, and artificial neural networks to reveal potential therapeutic targets for drug repurposing against COVID-19. Results: We screened 1,584 high-confidence immune system proteins in ACE2 and TMPRSS2 co-expressing cells, finding 25 potential therapeutic targets significantly overexpressed in nasal goblet secretory cells, lung type II pneumocytes, and ileal absorptive enterocytes of patients with several immunopathologies. Then, we performed fully connected deep neural networks to find the best multitask classification model to predict the activity of 10,672 drugs, obtaining several approved drugs, compounds under investigation, and experimental compounds with the highest area under the receiver operating characteristics. Conclusion: After being effectively analyzed in clinical trials, these drugs can be considered for treatment of severe COVID-19 patients. Scripts can be downloaded at https://github.com/muntisa/immuno-drug-repurposing-COVID-19.
RESUMO
Preeclampsia is a hypertensive disorder that occurs during pregnancy. It is a complex disease with unknown pathogenesis and the leading cause of fetal and maternal mortality during pregnancy. Using all drugs currently under clinical trial for preeclampsia, we extracted all their possible targets from the DrugBank and ChEMBL databases and labeled them as "targets". The proteins labeled as "off-targets" were extracted in the same way but while taking all antihypertensive drugs which are inhibitors of ACE and/or angiotensin receptor antagonist as query molecules. Classification models were obtained for each of the 55 total proteins (45 targets and 10 off-targets) using the TPOT pipeline optimization tool. The average accuracy of the models in predicting the external dataset for targets and off-targets was 0.830 and 0.850, respectively. The combinations of models maximizing their virtual screening performance were explored by combining the desirability function and genetic algorithms. The virtual screening performance metrics for the best model were: the Boltzmann-Enhanced Discrimination of ROC (BEDROC)α=160.9 = 0.258, the Enrichment Factor (EF)1% = 31.55 and the Area Under the Accumulation Curve (AUAC) = 0.831. The most relevant targets for preeclampsia were: AR, VDR, SLC6A2, NOS3 and CHRM4, while ABCG2, ERBB2, CES1 and REN led to the most relevant off-targets. A virtual screening of the DrugBank database identified estradiol, estriol, vitamins E and D, lynestrenol, mifrepristone, simvastatin, ambroxol, and some antibiotics and antiparasitics as drugs with potential application in the treatment of preeclampsia.
Assuntos
Anti-Hipertensivos/uso terapêutico , Bases de Dados de Produtos Farmacêuticos , Reposicionamento de Medicamentos , Modelos Cardiovasculares , Pré-Eclâmpsia/tratamento farmacológico , Anti-Hipertensivos/química , Feminino , Humanos , GravidezRESUMO
Osteosarcoma is the most common type of primary malignant bone tumor. Although nowadays 5-year survival rates can reach up to 60-70%, acute complications and late effects of osteosarcoma therapy are two of the limiting factors in treatments. We developed a multi-objective algorithm for the repurposing of new anti-osteosarcoma drugs, based on the modeling of molecules with described activity for HOS, MG63, SAOS2, and U2OS cell lines in the ChEMBL database. Several predictive models were obtained for each cell line and those with accuracy greater than 0.8 were integrated into a desirability function for the final multi-objective model. An exhaustive exploration of model combinations was carried out to obtain the best multi-objective model in virtual screening. For the top 1% of the screened list, the final model showed a BEDROC = 0.562, EF = 27.6, and AUC = 0.653. The repositioning was performed on 2218 molecules described in DrugBank. Within the top-ranked drugs, we found: temsirolimus, paclitaxel, sirolimus, everolimus, and cabazitaxel, which are antineoplastic drugs described in clinical trials for cancer in general. Interestingly, we found several broad-spectrum antibiotics and antiretroviral agents. This powerful model predicts several drugs that should be studied in depth to find new chemotherapy regimens and to propose new strategies for osteosarcoma treatment.
RESUMO
Sarcomas are a group of malignant neoplasms of connective tissue with a different etiology than carcinomas. The efforts to discover new drugs with antisarcoma activity have generated large datasets of multiple preclinical assays with different experimental conditions. For instance, the ChEMBL database contains outcomes of 37,919 different antisarcoma assays with 34,955 different chemical compounds. Furthermore, the experimental conditions reported in this dataset include 157 types of biological activity parameters, 36 drug targets, 43 cell lines, and 17 assay organisms. Considering this information, we propose combining perturbation theory (PT) principles with machine learning (ML) to develop a PTML model to predict antisarcoma compounds. PTML models use one function of reference that measures the probability of a drug being active under certain conditions (protein, cell line, organism, etc.). In this paper, we used a linear discriminant analysis and neural network to train and compare PT and non-PT models. All the explored models have an accuracy of 89.19-95.25% for training and 89.22-95.46% in validation sets. PTML-based strategies have similar accuracy but generate simplest models. Therefore, they may become a versatile tool for predicting antisarcoma compounds.
RESUMO
Wuhan, China was the epicenter of the first zoonotic transmission of the severe acute respiratory syndrome coronavirus clade 2 (SARS-CoV-2) in December 2019 and it is the causative agent of the novel human coronavirus disease 2019 (COVID-19). Almost from the beginning of the COVID-19 outbreak several attempts were made to predict possible drugs capable of inhibiting the virus replication. In the present work a drug repurposing study is performed to identify potential SARS-CoV-2 protease inhibitors. We created a Quantitative Structure-Activity Relationship (QSAR) model based on a machine learning strategy using hundreds of inhibitor molecules of the main protease (Mpro) of the SARS-CoV coronavirus. The QSAR model was used for virtual screening of a large list of drugs from the DrugBank database. The best 20 candidates were then evaluated in-silico against the Mpro of SARS-CoV-2 by using docking and molecular dynamics analyses. Docking was done by using the Gold software, and the free energies of binding were predicted with the MM-PBSA method as implemented in AMBER. Our results indicate that levothyroxine, amobarbital and ABP-700 are the best potential inhibitors of the SARS-CoV-2 virus through their binding to the Mpro enzyme. Five other compounds showed also a negative but small free energy of binding: nikethamide, nifurtimox, rebimastat, apomine and rebastinib.
Assuntos
Antivirais/farmacologia , Tratamento Farmacológico da COVID-19 , Proteases 3C de Coronavírus/antagonistas & inibidores , Descoberta de Drogas/métodos , Reposicionamento de Medicamentos/métodos , Inibidores de Proteases/farmacologia , SARS-CoV-2/enzimologia , Amobarbital/farmacologia , Antivirais/química , Sítios de Ligação , Simulação por Computador , Humanos , Aprendizado de Máquina , Simulação de Acoplamento Molecular , Simulação de Dinâmica Molecular , Pandemias , Inibidores de Proteases/química , Ligação Proteica , Relação Quantitativa Estrutura-Atividade , SARS-CoV-2/efeitos dos fármacos , Bibliotecas de Moléculas Pequenas/química , Software , Termodinâmica , Tiroxina/farmacologiaRESUMO
Breast cancer (BC) is a heterogeneous disease where genomic alterations, protein expression deregulation, signaling pathway alterations, hormone disruption, ethnicity and environmental determinants are involved. Due to the complexity of BC, the prediction of proteins involved in this disease is a trending topic in drug design. This work is proposing accurate prediction classifier for BC proteins using six sets of protein sequence descriptors and 13 machine-learning methods. After using a univariate feature selection for the mix of five descriptor families, the best classifier was obtained using multilayer perceptron method (artificial neural network) and 300 features. The performance of the model is demonstrated by the area under the receiver operating characteristics (AUROC) of 0.980 ± 0.0037, and accuracy of 0.936 ± 0.0056 (3-fold cross-validation). Regarding the prediction of 4,504 cancer-associated proteins using this model, the best ranked cancer immunotherapy proteins related to BC were RPS27, SUPT4H1, CLPSL2, POLR2K, RPL38, AKT3, CDK3, RPS20, RASL11A and UBTD1; the best ranked metastasis driver proteins related to BC were S100A9, DDA1, TXN, PRNP, RPS27, S100A14, S100A7, MAPK1, AGR3 and NDUFA13; and the best ranked RNA-binding proteins related to BC were S100A9, TXN, RPS27L, RPS27, RPS27A, RPL38, MRPL54, PPAN, RPS20 and CSRP1. This powerful model predicts several BC-related proteins that should be deeply studied to find new biomarkers and better therapeutic targets. Scripts can be downloaded at https://github.com/muntisa/neural-networks-for-breast-cancer-proteins.
Assuntos
Biomarcadores Tumorais/genética , Neoplasias da Mama/metabolismo , Regulação Neoplásica da Expressão Gênica , Imunoterapia/métodos , Aprendizado de Máquina , Redes Neurais de Computação , RNA/metabolismo , Neoplasias da Mama/secundário , Neoplasias da Mama/terapia , Feminino , Perfilação da Expressão Gênica , Humanos , Metástase NeoplásicaRESUMO
Breast cancer (BC) is the leading cause of cancer-related death among women and the most commonly diagnosed cancer worldwide. Although in recent years large-scale efforts have focused on identifying new therapeutic targets, a better understanding of BC molecular processes is required. Here we focused on elucidating the molecular hallmarks of BC heterogeneity and the oncogenic mutations involved in precision medicine that remains poorly defined. To fill this gap, we established an OncoOmics strategy that consists of analyzing genomic alterations, signaling pathways, protein-protein interactome network, protein expression, dependency maps in cell lines and patient-derived xenografts in 230 previously prioritized genes to reveal essential genes in breast cancer. As results, the OncoOmics BC essential genes were rationally filtered to 140. mRNA up-regulation was the most prevalent genomic alteration. The most altered signaling pathways were associated with basal-like and Her2-enriched molecular subtypes. RAC1, AKT1, CCND1, PIK3CA, ERBB2, CDH1, MAPK14, TP53, MAPK1, SRC, RAC3, BCL2, CTNNB1, EGFR, CDK2, GRB2, MED1 and GATA3 were essential genes in at least three OncoOmics approaches. Drugs with the highest amount of clinical trials in phases 3 and 4 were paclitaxel, docetaxel, trastuzumab, tamoxifen and doxorubicin. Lastly, we collected ~3,500 somatic and germline oncogenic variants associated with 50 essential genes, which in turn had therapeutic connectivity with 73 drugs. In conclusion, the OncoOmics strategy reveals essential genes capable of accelerating the development of targeted therapies for precision oncology.
Assuntos
Biomarcadores Tumorais/genética , Neoplasias da Mama/genética , Neoplasias da Mama/patologia , Regulação Neoplásica da Expressão Gênica , Genes Essenciais , Mutação , Medicina de Precisão , Animais , Biomarcadores Tumorais/metabolismo , Neoplasias da Mama/metabolismo , Feminino , Redes Reguladoras de Genes , Sequenciamento de Nucleotídeos em Larga Escala , Humanos , Camundongos , Prognóstico , Mapas de Interação de Proteínas , Proteoma , Células Tumorais Cultivadas , Ensaios Antitumorais Modelo de XenoenxertoRESUMO
Osteosarcoma is the most common subtype of primary bone cancer, affecting mostly adolescents. In recent years, several studies have focused on elucidating the molecular mechanisms of this sarcoma; however, its molecular etiology has still not been determined with precision. Therefore, we applied a consensus strategy with the use of several bioinformatics tools to prioritize genes involved in its pathogenesis. Subsequently, we assessed the physical interactions of the previously selected genes and applied a communality analysis to this protein-protein interaction network. The consensus strategy prioritized a total list of 553 genes. Our enrichment analysis validates several studies that describe the signaling pathways PI3K/AKT and MAPK/ERK as pathogenic. The gene ontology described TP53 as a principal signal transducer that chiefly mediates processes associated with cell cycle and DNA damage response It is interesting to note that the communality analysis clusters several members involved in metastasis events, such as MMP2 and MMP9, and genes associated with DNA repair complexes, like ATM, ATR, CHEK1, and RAD51. In this study, we have identified well-known pathogenic genes for osteosarcoma and prioritized genes that need to be further explored.
Assuntos
Neoplasias Ósseas/genética , Neoplasias Ósseas/patologia , Osteossarcoma/genética , Osteossarcoma/patologia , Biologia Computacional/métodos , Consenso , Reparo do DNA/genética , Regulação Neoplásica da Expressão Gênica/genética , Ontologia Genética , Redes Reguladoras de Genes/genética , Humanos , Mapas de Interação de Proteínas/genética , Transdução de Sinais/genéticaRESUMO
Head and neck squamous cell carcinoma (HNSCC) is the sixth most common cancer with the highest incidence worldwide. HNSCC is often diagnosed at advanced stages, incurring significant high mortality and morbidity. The use of saliva, as a noninvasive tool for the diagnosis of cancer, has recently increased. Salivary microRNAs (miRNAs) have emerged as a promising molecular tool for early diagnosis of HNSCC. The aim was to identify the differential expression of salivary miRNAs associated with HNSCC in the high altitude mestizo Ecuadorian population. Using PCR Arrays, miR-122-5p, miR-92a-3p, miR-124-3p, miR-205-5p, and miR-146a-5p were found as the most representative ones. Subsequently, miRNAs expression was confirmed in saliva samples from 108 cases and 108 controls. miR-122-5p, miR-92a-3p, miR-124-3p, and miR-146a-5p showed significant statistical difference between cases and controls with areas under the curve (AUC) of 0.73 (p < 0.001), 0.70 (p < 0.001), 0.71 (p = 0.002), and 0.66 (p = 0.008), respectively. miRNAs were also deregulated in between HNSCC localizations. A differentiated expression of miR-122-5p between oral cancer and oropharynx cancer (AUC of 0.96 p = 0.01) was found: miR-124-3p between larynx and pharynx (AUC = 0.97, p < 0.01) and miR-146a-5p between larynx, oropharynx, and oral cavity (AUC = 0.96, p = 0.01). Moreover, miR-122-5p, miR-124-3p, miR-205-5p, and miR-146a-5p could differentiate between HPV+ and HPV- (p=0.004). Finally, the expression profiles of the five miRNAs were evaluated to discriminate HNSCC patient's tumor stages (TNM 2-4). miR-122-5p differentiates TNM 2 and 3 (p = 0.002, AUC = 0.92), miR-124-3p TNM 2, 3, and 4 (p < 0.001, AUC = 98), miR-146a-5p TNM 2 and 3 (p < 0.001, AUC = 0.97), and miR-92a-3p TNM 3 (p < 0.001, AUC = 0.99). Taken together, these findings show that altered expression of miRNAs could be used as biomarkers for HNSCC diagnosis in the high altitude mestizo Ecuadorian population.
Assuntos
Neoplasias de Cabeça e Pescoço/diagnóstico , Neoplasias de Cabeça e Pescoço/metabolismo , MicroRNAs/metabolismo , Saliva/metabolismo , Carcinoma de Células Escamosas de Cabeça e Pescoço/diagnóstico , Carcinoma de Células Escamosas de Cabeça e Pescoço/metabolismo , Adulto , Idoso , Idoso de 80 Anos ou mais , Altitude , Biomarcadores Tumorais/metabolismo , Estudos de Casos e Controles , Detecção Precoce de Câncer/métodos , Equador , Feminino , Perfilação da Expressão Gênica/métodos , Regulação Neoplásica da Expressão Gênica/fisiologia , Neoplasias de Cabeça e Pescoço/patologia , Humanos , Masculino , Pessoa de Meia-Idade , Estadiamento de Neoplasias , Carcinoma de Células Escamosas de Cabeça e Pescoço/patologiaRESUMO
Consensus strategy was proved to be highly efficient in the recognition of gene-disease association. Therefore, the main objective of this study was to apply theoretical approaches to explore genes and communities directly involved in breast cancer (BC) pathogenesis. We evaluated the consensus between 8 prioritization strategies for the early recognition of pathogenic genes. A communality analysis in the protein-protein interaction (PPi) network of previously selected genes was enriched with gene ontology, metabolic pathways, as well as oncogenomics validation with the OncoPPi and DRIVE projects. The consensus genes were rationally filtered to 1842 genes. The communality analysis showed an enrichment of 14 communities specially connected with ERBB, PI3K-AKT, mTOR, FOXO, p53, HIF-1, VEGF, MAPK and prolactin signaling pathways. Genes with highest ranking were TP53, ESR1, BRCA2, BRCA1 and ERBB2. Genes with highest connectivity degree were TP53, AKT1, SRC, CREBBP and EP300. The connectivity degree allowed to establish a significant correlation between the OncoPPi network and our BC integrated network conformed by 51 genes and 62 PPi. In addition, CCND1, RAD51, CDC42, YAP1 and RPA1 were functional genes with significant sensitivity score in BC cell lines. In conclusion, the consensus strategy identifies both well-known pathogenic genes and prioritized genes that need to be further explored.