Genomic Biomarkers for Personalized Medicine in Breast Cancer

Breast cancer is the most common malignant disease in Western women. Historically, breast cancer was perceived as a single disease with various clinicalpathological features, and therefore, “one drug fits all” approaches drove the treatment regimens. The advent of genomics studies has led to a new paradigm in which breast cancer is heterogeneous consisting of different diseases from the same organ site. For example, gene expression profiling analysis revealed that estrogen receptor (ER)-positive and ER negative breast cancer are two distinct diseases with different risk factors, clinical presentations, outcomes, and responses to systemic therapies. Consequently, the new paradigm demands a personalized strategy in cancer medicine, in which the selection of treatment regimens for each cancer patient will largely rely on assessment by predictive biomarkers and study of the anatomical and pathological features of the cancer

Development of Estrogen Receptor Beta Binding Prediction Model Using Large Sets of Chemicals

We developed an ERβ binding prediction model to facilitate identification of chemicals specifically bind ERβ or ERα together with our previously developed ERα binding model. Decision Forest was used to train ERβ binding prediction model based on a large set of compounds obtained from EADB. Model performance was estimated through 1000 iterations of 5-fold cross validations. Prediction confidence was analyzed using predictions from the cross validations. Informative chemical features for ERβ binding were identified through analysis of the frequency data of chemical descriptors used in the models in the 5-fold cross validations. 1000 permutations were conducted to assess the chance correlation. The average accuracy of 5-fold cross validations was 93.14% with a standard deviation of 0.64%. Prediction confidence analysis indicated that the higher the prediction confidence the more accurate the predictions. Permutation testing results revealed that the prediction model is unlikely generated by chance. Eighteen informative descriptors were identified to be important to ERβ binding prediction. Application of the prediction model to the data from ToxCast project yielded very high sensitivity of 90-92%. Our results demonstrated ERβ binding of chemicals could be accurately predicted using the developed model. Coupling with our previously developed ERα prediction model, this model could be expected to facilitate drug development through identification of chemicals that specifically bind ERβ or ERα

\u3ci\u3eIn silico\u3c/i\u3e identification of genetic mutations conferring resistance to acetohydroxyacid synthase inhibitors: A case study of \u3ci\u3eKochia scoparia\u3c/i\u3e

Mutations that confer herbicide resistance are a primary concern for herbicide-based chemical control of invasive plants and are often under-characterized structurally and functionally. As the outcome of selection pressure, resistance mutations usually result from repeated long-term applications of herbicides with the same mode of action and are discovered through extensive field trials. Here we used acetohydroxyacid synthase (AHAS) of Kochia scoparia (KsAHAS) as an example to demonstrate that, given the sequence of a target protein, the impact of genetic mutations on ligand binding could be evaluated and resistance mutations could be identified using a biophysics-based computational approach. Briefly, the 3D structures of wild-type (WT) and mutated KsAHAS-herbicide complexes were constructed by homology modeling, docking and molecular dynamics simulation. The resistance profile of two AHAS-inhibiting herbicides, tribenuron methyl and thifensulfuron methyl, was obtained by estimating their binding affinity with 29 KsAHAS (1 WT and 28 mutated) using 6 molecular mechanical (MM) and 18 hybrid quantum mechanical/molecular mechanical (QM/MM) methods in combination with three structure sampling strategies. By comparing predicted resistance with experimentally determined resistance in the 29 biotypes of K. scoparia field populations, we identified the best method (i.e., MM-PBSA with single structure) out of all tested methods for the herbicide-KsAHAS system, which exhibited the highest accuracy (up to 100%) in discerning mutations conferring resistance or susceptibility to the two AHAS inhibitors. Our results suggest that the in silico approach has the potential to be widely adopted for assessing mutation-endowed herbicide resistance on a case-by-case basis

Development of a Nicotinic Acetylcholine Receptor nAChR α7 Binding Activity Prediction Model

Despite the well-known adverse health effects associated with tobacco use, addiction to nicotine found in tobacco products causes difficulty in quitting among users. Nicotinic acetylcholine receptors (nAChRs) are the physiological targets of nicotine and facilitate addiction to tobacco products. The nAChR-α7 subtype plays an important role in addiction; therefore, predicting the binding activity of tobacco constituents to nAChR-α7 is an important component for assessing addictive potential of tobacco constituents. We developed an α7 binding activity prediction model based on a large training data set of 843 chemicals with human α7 binding activity data extracted from PubChem and ChEMBL. The model was tested using 1215 chemicals with rat α7 binding activity data from the same databases. Based on the competitive docking results, the docking scores were partitioned to the key residues that play important roles in the receptor−ligand binding. A decision forest was used to train the human α7 binding activity prediction model based on the partition of docking scores. Five-fold cross validations were conducted to estimate the performance of the decision forest models. The developed model was used to predict the potential human α7 binding activity for 5275 tobacco constituents. The human α7 binding activity data for 84 of the 5275 tobacco constituents were experimentally measured to confirm and empirically validate the prediction results. The prediction accuracy, sensitivity, and specificity were 64.3, 40.0, and 81.6%, respectively. The developed prediction model of human α7 may be a useful tool for high-throughput screening of potential addictive tobacco constituents

Evaluation of gene expression data generated from expired Affymetrix GeneChip® microarrays using MAQC reference RNA samples

BACKGROUND: The Affymetrix GeneChip(®) system is a commonly used platform for microarray analysis but the technology is inherently expensive. Unfortunately, changes in experimental planning and execution, such as the unavailability of previously anticipated samples or a shift in research focus, may render significant numbers of pre-purchased GeneChip(®) microarrays unprocessed before their manufacturer’s expiration dates. Researchers and microarray core facilities wonder whether expired microarrays are still useful for gene expression analysis. In addition, it was not clear whether the two human reference RNA samples established by the MAQC project in 2005 still maintained their transcriptome integrity over a period of four years. Experiments were conducted to answer these questions. RESULTS: Microarray data were generated in 2009 in three replicates for each of the two MAQC samples with either expired Affymetrix U133A or unexpired U133Plus2 microarrays. These results were compared with data obtained in 2005 on the U133Plus2 microarray. The percentage of overlap between the lists of differentially expressed genes (DEGs) from U133Plus2 microarray data generated in 2009 and in 2005 was 97.44%. While there was some degree of fold change compression in the expired U133A microarrays, the percentage of overlap between the lists of DEGs from the expired and unexpired microarrays was as high as 96.99%. Moreover, the microarray data generated using the expired U133A microarrays in 2009 were highly concordant with microarray and TaqMan(®) data generated by the MAQC project in 2005. CONCLUSIONS: Our results demonstrated that microarray data generated using U133A microarrays, which were more than four years past the manufacturer’s expiration date, were highly specific and consistent with those from unexpired microarrays in identifying DEGs despite some appreciable fold change compression and decrease in sensitivity. Our data also suggested that the MAQC reference RNA samples, stored at -80°C, were stable over a time frame of at least four years

Toward predictive models for drug-induced liver injury in humans: are we there yet?

Drug-induced liver injury (DILI) is a frequent cause for the termination of drug development programs and a leading reason of drug withdrawal from the marketplace. Unfortunately, the current preclinical testing strategies, including the regulatory-required animal toxicity studies or simple in vitro tests, are insufficiently powered to predict DILI in patients reliably. Notably, the limited predictive power of such testing strategies is mostly attributed to the complex nature of DILI, a poor understanding of its mechanism, a scarcity of human hepatotoxicity data and inadequate bioinformatics capabilities. With the advent of high-content screening assays, toxicogenomics and bioinformatics, multiple end points can be studied simultaneously to improve prediction of clinically relevant DILIs. This review focuses on the current state of efforts in developing predictive models from diverse data sources for potential use in detecting human hepatotoxicity, and also aims to provide perspectives on how to further improve DILI prediction