Chapter 3: Small Molecules and Disease

“Big” molecules such as proteins and genes still continue to capture the imagination of most biologists, biochemists and bioinformaticians. “Small” molecules, on the other hand, are the molecules that most biologists, biochemists and bioinformaticians prefer to ignore. However, it is becoming increasingly apparent that small molecules such as amino acids, lipids and sugars play a far more important role in all aspects of disease etiology and disease treatment than we realized. This particular chapter focuses on an emerging field of bioinformatics called “chemical bioinformatics” – a discipline that has evolved to help address the blended chemical and molecular biological needs of toxicogenomics, pharmacogenomics, metabolomics and systems biology. In the following pages we will cover several topics related to chemical bioinformatics. First, a brief overview of some of the most important or useful chemical bioinformatic resources will be given. Second, a more detailed overview will be given on those particular resources that allow researchers to connect small molecules to diseases. This section will focus on describing a number of recently developed databases or knowledgebases that explicitly relate small molecules – either as the treatment, symptom or cause – to disease. Finally a short discussion will be provided on newly emerging software tools that exploit these databases as a means to discover new biomarkers or even new treatments for disease.

What to Learn in This Chapter

• The meaning of chemical bioinformatics
• Strengths and limitations of existing chemical bioinformatic databases
• Using databases to learn about the cause and treatment of diseases
• The Small Molecule Pathway Database (SMPDB)
• The Human Metabolome Database (HMDB)
• DrugBank
• The Toxin and Toxin-Target Database (T3DB)
• PolySearch and Metabolite Set Enrichment Analysis

This article is part of the “Translational Bioinformatics” collection for PLOS Computational Biology.

     

靶向小分子创新药物

寻找药物新靶点是全球创新药物研究激烈竞争的焦点."组学"、生物信息学、系统生物学、药物筛选现代检测技术等新理论、新技术的发展使新的筛选模型和评价技术不断取得突破.靶向抗肿瘤药物的开发是靶向小分子创新药物的重点任务,多靶点的抗肿瘤药物开发及新靶点的发现是抗肿瘤药物研发的新趋势.     

Block of Nav1.8 by Small Molecules

Sodium channels are key proteins in regulating neuronal excitability and accumulating data suggest that specific subtypes of voltage-dependent sodium channels are important in signaling various types of pain. Consistent with this theme, Jarvis et al (2007) recently reported the identification of a subtype-selective Nav1.8 blocker that was active in several pre-clinical models of pain. During the course of these studies compounds were also identified that showed large differences in potency when tested on Nav1.8 channels from different species. This Addendum illustrates one of these compounds along with the potency correlation between recombinant and native tetrodotoxin-resistant sodium channels for additional examples. These data show that significant differences can be observed for sodium channel blockers across species and highlight the importance of considering this possibility when searching for new compounds and research tools to probe sodium channel function.     

Effects of Small Molecules on Chaperone-Mediated Autophagy

Autophagy, including macroautophagy (MA), chaperone-mediated autophagy (CMA), crinophagy, pexophagy and microautophagy, are processes by which cells select internal components such as proteins, secretory vesicles, organelles, or foreign bodies, and deliver them to lysosomes for degradation. MA and CMA are activated during conditions of serum withdrawal in cell culture and during short-term (MA) and prolonged (CMA) starvation in organisms. Although MA and CMA are activated under similar conditions, they are regulated by different mechanisms. We used pulse/chase analysis under conditions in which most intracellular proteolysis is due to CMA to test a variety of compounds for effects on CMA. We show that inhibitors of MA such as 3-methyladenine, wortmannin, and LY294002 have no effect on CMA. Protein degradation by MA is sensitive to microtubule inhibitors such as colcemide and vinblastine, but protein degradation by CMA is not. Activators of MA such as rapamycin also have no effect on CMA. We demonstrate that CMA, like MA, is inhibited by protein synthesis inhibitors anisomycin and cycloheximide. CMA is also partially inhibited when the P38 mitogen activated protein kinase is blocked. Finally we demonstrate that the glucose-6-phophate dehydrogenase inhibitor, 6-aminonicotinamide, and heat shock protein of 90 kilodaltons inhibitor, geldanamycin, have the ability to activate CMA.     

Molecular Dynamics Simulations of some Small Organic Molecules

Abstract

Free energy differences between different conformers of D-ribofuranose, L-malic acid and meso-tartaric acid in solution were calculated using Molecular Dynamics simulations. In case of ribose the α → β transition was studied. For the acids attention was focussed on the transitions between the three possible staggered conformers with respect to the central C-C bond. In all cases a thermodynamic integration method was employed to evaluate the free energy difference. The use of an alternative technique, umbrella sampling, for ribose did not give promising results.

It was shown that one needs a fairly accurate picture of the accessible conformational space in case of flexible molecules like the ones considered here before one can determine meaningful free energy differences. Large hysteresis effects between forward and reverse simulated transitions were observed, but contrary to the general belief they are no direct measure of the accuracy of the calculated ΔG values. In all cases the ΔG values resulting from the simulations and from NMR experiments agree within the, considerable, error limits and for the different forms of D-ribose, L-malic acid and L-tartaric acid the relative order of their populations is also correctly reproduced.     

Modulating Molecular Functions of p53 with Small Molecules

Association of the p53 with many human cancers makes it a valuable therapeutic target. Stress-induced molecular interactions of p53 with other effector proteins are immensely intertwined with regulation of its functions in orchestrating a wide array of cellular responses, thereby defying analysis of the underlying molecular mechanisms with conventional molecular and cellular biology methods. Recent discoveries of small molecules that can selectively modulate the molecular interactions of p53 offer promising opportunities to address the challenge of dissecting these complex mechanisms and increase the hope for pharmacological control of p53 for clinical benefits of cancer patients.     

MDM2-P53蛋白质相互作用的小分子抑制物

p53作为重要的抑癌基因已经成为一个治疗癌症重点的突破目标之一。直接调节p53基因或调节P53和MDM2蛋白质相互作用是再激活p53基因的两种重要机制。对于表达野生型P53的癌症设计小分子阻断剂阻断MDM2与P53蛋白相互作用是一个很有前景的治疗癌症的方向。文章主要总结了作为治疗癌症的新方法-MDM2-P53蛋白相互作用小分子抑制物的最新研究进展,其中最新的是人工合成化合物Nutlin-3和MI-219。     

Dynamics of Small Molecules in a Dense Polymer Matrix: Molecular Dynamics Studies

Abstract

The anomalous diffusion regime appearing in the self-diffusion of small molecules in bulk amorphous polymers has been extensively studied by molecular dynamics simulations. A rather long simulation of duration ～ 10 ^?8 s is performed on a polyethylene-like simple polymer model containing either oxygen molecules or helium atoms as a diffusant. Dynamic properties evaluated for these diffusants are the mean-square displacement, the van Hove self correlation function, and the self part of the density autocorrelation. It is first confirmed that the anomalous diffusion regime appears in a few hundred picoseconds for oxygen molecule, while the Einstein relation adopted beyond this regime results in the self-diffusion coefficient of the order of ～ 10^?5 cm²/s. This anomaly is still observed for helium that diffuses much faster than oxygen. In the anomalous diffusion regime, it is found that the correlation functions for the two diffusants show characteristic features and become essentially the same as time is scaled appropriately. These features allow the estimation of the two characteristic spatial scales which are probably dominated by the microstructure of the polymer matrix, namely, the cage size and the distance between adjacent cages. The time dependence of the mean-square displacements of the two diffusants can be well interpreted by these characteristic spatial scales as time is scaled with the self-diffusion coefficients. It is shown that the anomalous diffusion regime arises from the inhomogeneous microstructure of the polymer matrix.     

Amyloid-Precursor-Protein-Lowering Small Molecules for Disease Modifying Therapy of Alzheimer's Disease

Alzheimer''s disease (AD) is the most common form of dementia in the elderly with progressive cognitive decline and memory loss. According to the amyloid-hypothesis, AD is caused by generation and subsequent cerebral deposition of β-amyloid (Aβ). Aβ is generated through sequential cleavage of the transmembrane Amyloid-Precursor-Protein (APP) by two endoproteinases termed beta- and gamma-secretase. Increased APP-expression caused by APP gene dosage effects is a risk factor for the development of AD. Here we carried out a large scale screen for novel compounds aimed at decreasing APP-expression. For this we developed a screening system employing a cell culture model of AD. A total of 10,000 substances selected for their ability of drug-likeness and chemical diversity were tested for their potential to decrease APP-expression resulting in reduced Aβ-levels. Positive compounds were further evaluated for their effect at lower concentrations, absence of cytotoxicity and specificity. The six most promising compounds were characterized and structure function relationships were established. The novel compounds presented here provide valuable information for the development of causal therapies for AD.     

Identification of Small Molecules That Suppress Ricin-Induced Stress-Activated Signaling Pathways

Ricin is a member of the ribosome-inactivating protein (RIP) family of plant and bacterial toxins. In this study we used a high-throughput, cell-based assay to screen more than 118,000 compounds from diverse chemical libraries for molecules that reduced ricin-induced cell death. We describe three compounds, PW66, PW69, and PW72 that at micromolar concentrations significantly delayed ricin-induced cell death. None of the compounds had any demonstrable effect on ricin''s ability to arrest protein synthesis in cells or on ricin''s enzymatic activity as assessed in vitro. Instead, all three compounds appear to function by blocking downstream stress-induced signaling pathways associated with the toxin-mediated apoptosis. PW66 virtually eliminated ricin-induced TNF-α secretion by J774A.1 macrophages and concomitantly blocked activation of the p38 MAPK and JNK signaling pathways. PW72 suppressed ricin-induced TNF-α secretion, but not p38 MAPK and JNK signaling. PW69 suppressed activity of the executioner caspases 3/7 in ricin toxin- and Shiga toxin 2-treated cells. While the actual molecular targets of the three compounds have yet to be identified, these data nevertheless underscore the potential of small molecules to down-regulate inflammatory signaling pathways associated with exposure to the RIP family of toxins.     

Analysis of Group Randomized Trials with Multiple Binary Endpoints and Small Number of Groups

The group randomized trial (GRT) is a common study design to assess the effect of an intervention program aimed at health promotion or disease prevention. In GRTs, groups rather than individuals are randomized into intervention or control arms. Then, responses are measured on individuals within those groups. A number of analytical problems beset GRT designs. The major problem emerges from the likely positive intraclass correlation among observations of individuals within a group. This paper provides an overview of the analytical method for GRT data and applies this method to a randomized cancer prevention trial, where multiple binary primary endpoints were obtained. We develop an index of extra variability to investigate group-specific effects on response. The purpose of the index is to understand the influence of individual groups on evaluating the intervention effect, especially, when a GRT study involves a small number of groups. The multiple endpoints from the GRT design are analyzed using a generalized linear mixed model and the stepdown Bonferroni method of Holm.     

A CFF 91-based Continuum Solvation Model: Solvation Free Energies of Small Organic Molecules and Conformations of the Alanine Dipeptide in Solution

Abstract

For molecular mechanics simulations of solvated molecules, it is important to use a consistent approach for calculating both the force field energy and the solvation free energy. A continuum solvation model based upon the atomic charges provided with the CFF91 force field is derived. The electrostatic component of the solvation free energy is described by the Poisson-Bolzmann equation while the nonpolar comonent of the solvation energy is assumed to be proportional to the solvent accessible surface area of the solute. Solute atomic radii used to describe the interface between the solute and solvent are fitted to reproduce the energies of small organic molecules. Data for 140 compounds are presented and compared to experiment and to the results from the well-characterized quantum mechanical solvation model AM1-SM2. In particular, accurate results are obtained for amino acid neutral analogues (mean unsigned error of 0.3 kcal/mol). The conformational energetics of the solvated alanine dipeptide is discussed.     

Antioxidants: Molecules, medicines, and myths

There is an industry-driven public obsession with antioxidants, which are equated to safe, health-giving molecules to be swallowed as mega-dose supplements or in fortified foods. Sometimes they are good for you, but sometimes they may not be, and pro-oxidants can be better for you in some circumstances. This article re-examines and challenges some basic assumptions in the nutritional antioxidant field.