Anticancer drug discovery using artificial intelligence: an application in pharmacological activity prediction

Hematological cancers are a heterogeneous family of diseases that can be divided into leukemias, lymphomas, and myelomas, often called “liquid tumors”. Since they cannot be surgically removable, chemotherapy represents the mainstay of their treatment. However, it still faces several challenges like drug resistance and low response rate, and the need for new anticancer agents is compelling. The drug discovery process is long-term, costly, and prone to high failure rates. With the rapid expansion of biological and chemical "big data", some computational techniques such as machine learning tools have been increasingly employed to speed up and economize the whole process. Machine learning algorithms can create complex models with the aim to determine the biological activity of compounds against several targets, based on their chemical properties. These models are defined as multi-target Quantitative Structure-Activity Relationship (mt-QSAR) and can be used to virtually screen small and large chemical libraries for the identification of new molecules with anticancer activity. The aim of my Ph.D. project was to employ machine learning techniques to build an mt-QSAR classification model for the prediction of cytotoxic drugs simultaneously active against 43 hematological cancer cell lines. For this purpose, first, I constructed a large and diversified dataset of molecules extracted from the ChEMBL database. Then, I compared the performance of different ML classification algorithms, until Random Forest was identified as the one returning the best predictions. Finally, I used different approaches to maximize the performance of the model, which achieved an accuracy of 88% by correctly classifying 93% of inactive molecules and 72% of active molecules in a validation set. This model was further applied to the virtual screening of a small dataset of molecules tested in our laboratory, where it showed 100% accuracy in correctly classifying all molecules. This result is confirmed by our previous in vitro experiments.

Synthesis and characterization of novel selective NKCC1 inhibitors for the treatment of down syndrome and brain disorders characterized by depolarizing GABAergic transmission.

Proper GABAergic transmission through Cl-permeable GABAA receptors is fundamental for physiological brain development and function. Indeed, defective GABAergic signaling – due to a high NKCC1/KCC2 expression ratio – has been implicated in several neurodevelopmental disorders (e.g., Down syndrome, DS, Autism spectrum disorders, ASD). Interestingly, NKCC1 inhibition by the FDA-approved diuretic drug bumetanide reverts cognitive deficits in the TS65Dn mouse models of DS and core symptoms in other models of brain disorders. However, the required chronic treatment with bumetanide is burdened by its diuretic side effects caused by the antagonization of the kidney Cl importer NKCC2. This may lead to hypokalemia, while jeopardizing drug compliance. Crucially, these issues would be solved by selective NKCC1 inhibitors, thus devoid of the diuretic effect of bumetanide. To this aim, starting from bumetanide’s structure, we applied a ligand-based computational approach to design new molecular entities that we tested in vitro for their capacity to selectively block NKCC1. Extensive synthetic efforts and structure-activity relationships analyses allowed us to improve in vitro potency and overall drug-like properties of the initially identified chemical hits. As a result, we identified a new highly potent NKCC1 inhibitor (ARN23746) that displayed excellent solubility, metabolic stability, and no significant effect on NKCC2 in vitro. Moreover, this novel and selective NKCC1 inhibitor was able to rescue cognitive deficits in DS mice and social/repetitive behaviors in ASD mice, with no diuretic effect and no overt toxicity upon chronic treatment in adult animals. Thus, ARN23746 a selective NKCC1 inhibitor devoid of the diuretic effect – represents a suitable and solid therapeutic strategy for the treatment of Down syndrome and all the brain neurological disorders characterized by depolarizing GABAergic transmission.

Cheaper Faster Drug Development Validated By The Repositioning Of Drugs Against Neglected Tropical Diseases.

There is an urgent need to make drug discovery cheaper and faster. This will enable the development of treatments for diseases currently neglected for economic reasons, such as tropical and orphan diseases, and generally increase the supply of new drugs. Here, we report the Robot Scientist 'Eve' designed to make drug discovery more economical. A Robot Scientist is a laboratory automation system that uses artificial intelligence (AI) techniques to discover scientific knowledge through cycles of experimentation. Eve integrates and automates library-screening, hit-confirmation, and lead generation through cycles of quantitative structure activity relationship learning and testing. Using econometric modelling we demonstrate that the use of AI to select compounds economically outperforms standard drug screening. For further efficiency Eve uses a standardized form of assay to compute Boolean functions of compound properties. These assays can be quickly and cheaply engineered using synthetic biology, enabling more targets to be assayed for a given budget. Eve has repositioned several drugs against specific targets in parasites that cause tropical diseases. One validated discovery is that the anti-cancer compound TNP-470 is a potent inhibitor of dihydrofolate reductase from the malaria-causing parasite Plasmodium vivax.

Nectar Robbery By A Hermit Hummingbird: Association To Floral Phenotype And Its Influence On Flowers And Network Structure.

Interactions between flowers and their visitors span the spectrum from mutualism to antagonism. The literature is rich in studies focusing on mutualism, but nectar robbery has mostly been investigated using phytocentric approaches focused on only a few plant species. To fill this gap, we studied the interactions between a nectar-robbing hermit hummingbird, Phaethornis ruber, and the array of flowers it visits. First, based on a literature review of the interactions involving  P. ruber, we characterized the association of floral larceny to floral phenotype. We then experimentally examined the effects of nectar robbing on nectar standing crop and number of visits of the pollinators to the flowers of Canna paniculata. Finally, we asked whether the incorporation of illegitimate interactions into the analysis affects plant-hummingbird network structure. We identified 97 plant species visited by P. ruber and found that P. ruber engaged in floral larceny in almost 30 % of these species. Nectar robbery was especially common in flowers with longer corolla. In terms of the effect on C. paniculata, the depletion of nectar due to robbery by P. ruber was associated with decreased visitation rates of legitimate pollinators. At the community level, the inclusion of the illegitimate visits of P. ruber resulted in modifications of how modules within the network were organized, notably giving rise to a new module consisting of P. ruber and mostly robbed flowers. However, although illegitimate visits constituted approximately 9 % of all interactions in the network, changes in nestedness, modularity, and network-level specialization were minor. Our results indicate that although a flower robber may have a strong effect on the pollination of a particular plant species, the inclusion of its illegitimate interactions has limited capacity to change overall network structure.

Glicoproteína-P, resistência a múltiplas drogas (MDR) e relação estrutura-atividade de moduladores

Multidrug resistance, MDR is a major obstacle for cancer chemotherapy. MDR can be reversed by drugs that vary in their chemical structure and main biological activity. Many efforts have been done to overcome MDR based on studies of structure-activity relationships and in this review we summarize some aspects of MDR mediated by P-glycoprotein (P-gp), as the most experimentally and clinically tested form of drug resistance. The most significant MDR mechanisms revealed until now are shortly discussed. Physicochemical and structural properties of MDR modulators, measures of the MDR reversal, and QSAR studies are included.

Desenhos de estruturas químicas correlacionam-se com propriedades biológicas: MIA-QSAR

Descriptors in multivariate image analysis applied to quantitative structure-activity relationship (MIA-QSAR) are pixels of bidimensional images of chemical structures (drawings), which were used to model the trichomonicidal activities of a series of benzimidazole derivatives. The MIA-QSAR model showed good predictive ability, with r², q² and r val. ext.² of 0.853, 0.519 and 0.778, respectively, which are comparable to the best values obtained by CoMFA e CoMSIA for the same series. A MIA-based analysis was also performed by using images of alphabetic letters with the corresponding numeric ordering as dependent variables, but no correlation was found, supporting that MIA-QSAR is not arbitrary.

Anti-oxidative and anti-ulcerogenic activity of Ipomoea imperati

Ipomoea imperati (Vahl) Griseb., Convolvulaceae, is used in traditional medicine for the treatment of inflammation, swelling and wounds, as well as to treat pains and stomach problems. This work evaluates the anti-oxidative activity by ESR (Electron Spin Resonance spectroscopy) and the preventive and curative actions of I. imperati in gastric ulcer animal model. Ipomoea imperati (200 mg/kg, p.o.) prevented the formation of gastric lesions in 78% (p<0.05) when compared with the negative control tween 80. Lanzoprazole, prevented in 85% the gastric lesions formation induced by ethanol (p<0.05). Therefore, the oral administration of I. imperati one hour before the ulcerogenic agent prevented the ulcer formation, conserving the citoprotection characteristics of the gastric mucosa and assuring the integrity of gastric glands and gastric fossets. The healing activity of I. imperati (200 mg/kg, p.o.) evaluated in chronic ulcer experiments induced by the acetic acid, was 72% (p<0.05). The positive control, ranitidine, healed 78% of the gastric lesions (p<0.05). The histological analysis confirmed the recovery of the mucosal layer and the muscle mucosal layer harmed by the acetic acid. Experiments in vitro with DPPH (2.2-diphenyl-1-picrylhydrazyl) of anti-oxidative activity demonstrated that I. imperati presents an IC50 of 0.73±0.01 mg/mL.

Vocalizations and comments on the relationships of Hypsiboas ericae (Amphibia, Hylidae)

The vocalizations of Hypsiboas ericae (Caramaschi & Cruz, 2000) are described and new information on the external morphology and osteology of the species are presented. H. ericae presents a bony spine in the prepolex and the individuals can present green or brown dorsal color, as other species of the Hypsiboas pulchellus (Duméril & Bibron, 1841) species group. The vocalizations of H. ericae are similar to the vocalizations of Hypsiboas bischoffi (Boulenger, 1887), Hypsiboas guentheri (Boulenger, 1886), and other species in the H. polytaenius (Cope, 1870 "1869") clade of the H. pulchellus species group, but some osteological aspects are different to those found in the majority of the species of this group.

Rates of antimicrobial resistance in latin america (2004-2007) and in vitro activity of the glycylcycline tigecycline and of other antibiotics

As a part of the Tigecycline Evaluation and Surveillance Trial (T.E.S.T.), Gram-positive and Gram-negative bacterial isolates were collected from 33 centers in Latin America (centers in Argentina, Brazil, Chile, Colombia, Guatemala, Honduras, Jamaica, Mexico, Panama, Puerto Rico, and Venezuela) from January 2004 to September 2007. Argentina and Mexico were the greatest contributors of isolates to this study. Susceptibilities were determined according to Clinical Laboratory Standards Institute guidelines. Resistance levels were high for most key organisms across Latin America: 48.3% of Staphylococcus aureus isolates were methicillin-resistant while 21.4% of Acinetobacter spp. isolates were imipenem-resistant. Extended-spectrum β-lactamase were reported in 36.7% of Klebsiella pneumoniae and 20.8% of E. coli isolates. Tigecycline was the most active agent against Gram-positive isolates. Tigecycline was also highly active against all Gram-negative organisms, with the exception of Pseuodomonas aeruginosa, against which piperacillin-tazobactam was the most active agent tested (79.3% of isolates susceptible). The in vitro activity of tigecycline against both Gram-positive and Gram-negative isolates indicates that it may be an useful tool for the treatment of nosocomial infections, even those caused by organisms that are resistant to other antibacterial agents.

Chemical composition and acetylcholinesterase inhibitory activity of essential oils of Myrceugenia myrcioides(Cambess.) O. Berg and Eugenia riedelianaO. Berg, Myrtaceae

The chemical composition of volatile oils from two Myrtaceae species, Myrceugenia myrcioidesand Eugenia riedeliana, both native from the Brazilian Atlantic Rain Forest, was analyzed by GC-MS. Acetylcholinesterase inhibitory activity was colorimetrically evaluated for these oils. For M. myrcioides, monoterpene hydrocarbons represented the major class in the volatile oil, with α-pinene as the most abundant component and a weak inhibitory activity was observed, whilst for E. riedeliana sesquiterpenes were found in higher amounts, being valerianol the major compound, and this oil presented a strong acetylcholinesterase inhibition.

Epidemiology of leisure, transportation, occupational, and household physical activity: prevalence and associated factors

Background: To estimate the prevalence of and identify factors associated with physical activity in leisure, transportation, occupational, and household settings. Methods: This was a cross-sectional study aimed at investigating living and health conditions among the population of São Paulo, Brazil. Data on 1318 adults aged 18 to 65 years were used. To assess physical activity, the long version of the International Physical Activity Questionnaire was applied. Multivariate analysis was conducted using a hierarchical model. Results: The greatest prevalence of insufficient activity related to transportation (91.7%), followed by leisure (77.5%), occupational (68.9%), and household settings (56.7%). The variables associated with insufficient levels of physical activity in leisure were female sex, older age, low education level, nonwhite skin color, smoking, and self-reported poor health; in occupational settings were female sex, white skin color, high education level, self-reported poor health, nonsmoking, and obesity; in transportation settings were female sex; and in household settings, with male sex, separated, or widowed status and high education level. Conclusion: Physical activity in transportation and leisure settings should be encouraged. This study will serve as a reference point in monitoring different types of physical activities and implementing public physical activity policies in developing countries

Epidemiology of Leisure, Transportation, Occupational, and Household Physical Activity: Prevalence and Associated Factors

Background: To estimate the prevalence of and identify factors associated with physical activity in leisure, transportation, occupational, and household settings. Methods: This was a cross-sectional study aimed at investigating living and health conditions among the population of So Paulo, Brazil. Data on 1318 adults aged 18 to 65 years were used. To assess physical activity, the long version of the International Physical Activity Questionnaire was applied. Multivariate analysis was conducted using a hierarchical model. Results: The greatest prevalence of insufficient activity related to transportation (91.7%), followed by leisure (77.5%), occupational (68.9%), and household settings (56.7%). The variables associated with insufficient levels of physical activity in leisure were female sex, older age, low education level, nonwhite skin color, smoking, and self-reported poor health; in occupational settings were female sex, white skin color, high education level, self-reported poor health, nonsmoking, and obesity; in transportation settings were female sex; and in household settings, with male sex, separated, or widowed status and high education level. Conclusion: Physical activity in transportation and leisure settings should be encouraged. This study will serve as a reference point in monitoring different types of physical activities and implementing public physical activity policies in developing countries.

Chemometric Studies on Natural Products as Potential Inhibitors of the NADH Oxidase from Trypanosoma cruzi Using the VolSurf Approach

Natural products have widespread biological activities, including inhibition of mitochondrial enzyme systems. Some of these activities, for example cytotoxicity, may be the result of alteration of cellular bioenergetics. Based on previous computer-aided drug design (CADD) studies and considering reported data on structure-activity relationships (SAR), an assumption regarding the mechanism of action of natural products against parasitic infections involves the NADH-oxidase inhibition. In this study, chemometric tools, such as: Principal Component Analysis (PCA), Consensus PCA (CPCA), and partial least squares regression (PLS), were applied to a set of forty natural compounds, acting as NADH-oxidase inhibitors. The calculations were performed using the VolSurf+ program. The formalisms employed generated good exploratory and predictive results. The independent variables or descriptors having a hydrophobic profile were strongly correlated to the biological data.

LIFE CYCLE AND POPULATION STRUCTURE OF AEGLA PAULENSIS (DECAPODA: ANOMURA: AEGLIDAE)

We describe growth, longevity, sex ratio, reproductive period, and recruitment of Aegla paulensis from Jaragua Stale Park, Sao Paulo, Brazil (23 degrees 27'27.9 '' S; 46 degrees 45'32.3 '' W). The population was sampled monthly (September 2007 through August 2009) with the aid of traps. Over five thousand individuals were captured, sexed, measured (carapace length = CL) and inspected for reproductive traits (females only), and then released back to the sampling site. The pattern of the reproductive cycle was strongly seasonal (austral mid autumn through late winter), with a single recruitment pulse per year. The obtained von Bertalanffy growth equations were CL = 21.25[1-e(-0.041(t + 1.250))] and CL = 16.52[1-e(-0.049(t + 1.823))] for males and females, respectively. Males (mean CL +/- SD = 11.86 +/- 2.79 mm) attain larger sizes than females (mean CL +/- SD = 10.84 +/- 2.36 mm). Aegla paulensis reproduces twice during an estimated life span of 40.2 months for females and 33.9 months for males. Temporal variation of sex ratio showed a distinctive pattern characterized by a sequence of three distinct periods that repeated from one year to another, and which suggested that a behavioral component influence the proportion of sex in adult specimens sampled with traps during reproductive and non-reproductive periods.

Single-nucleotide polymorphism, linkage disequilibrium and geographic structure in the malaria parasite Plasmodium vivax: prospects for genome-wide association studies

Background: The ideal malaria parasite populations for initial mapping of genomic regions contributing to phenotypes such as drug resistance and virulence, through genome-wide association studies, are those with high genetic diversity, allowing for numerous informative markers, and rare meiotic recombination, allowing for strong linkage disequilibrium (LD) between markers and phenotype-determining loci. However, levels of genetic diversity and LD in field populations of the major human malaria parasite P. vivax remain little characterized. Results: We examined single-nucleotide polymorphisms (SNPs) and LD patterns across a 100-kb chromosome segment of P. vivax in 238 field isolates from areas of low to moderate malaria endemicity in South America and Asia, where LD tends to be more extensive than in holoendemic populations, and in two monkey-adapted strains (Salvador-I, from El Salvador, and Belem, from Brazil). We found varying levels of SNP diversity and LD across populations, with the highest diversity and strongest LD in the area of lowest malaria transmission. We found several clusters of contiguous markers with rare meiotic recombination and characterized a relatively conserved haplotype structure among populations, suggesting the existence of recombination hotspots in the genome region analyzed. Both silent and nonsynonymous SNPs revealed substantial between-population differentiation, which accounted for similar to 40% of the overall genetic diversity observed. Although parasites clustered according to their continental origin, we found evidence for substructure within the Brazilian population of P. vivax. We also explored between-population differentiation patterns revealed by loci putatively affected by natural selection and found marked geographic variation in frequencies of nucleotide substitutions at the pvmdr-1 locus, putatively associated with drug resistance. Conclusion: These findings support the feasibility of genome-wide association studies in carefully selected populations of P. vivax, using relatively low densities of markers, but underscore the risk of false positives caused by population structure at both local and regional levels.