Molecular Phenotyping Combines Molecular Information, Biological Relevance, and Patient Data to Improve Productivity of Early Drug Discovery.

Today, novel therapeutics are identified in an environment which is intrinsically different from the clinical context in which they are ultimately evaluated. Using molecular phenotyping and an in vitro model of diabetic cardiomyopathy, we show that by quantifying pathway reporter gene expression, molecular phenotyping can cluster compounds based on pathway profiles and dissect associations between pathway activities and disease phenotypes simultaneously. Molecular phenotyping was applicable to compounds with a range of binding specificities and triaged false positives derived from high-content screening assays. The technique identified a class of calcium-signaling modulators that can reverse disease-regulated pathways and phenotypes, which was validated by structurally distinct compounds of relevant classes. Our results advocate for application of molecular phenotyping in early drug discovery, promoting biological relevance as a key selection criterion early in the drug development cascade.

Identification of Drug-Drug Interactions In Vitro: A Case Study Evaluating the Effects of Sofosbuvir and Amiodarone on hiPSC-Derived Cardiomyocytes.

Drug-drug interactions pose a difficult drug safety problem, given the increasing number of individuals taking multiple medications and the relative complexity of assessing the potential for interactions. For example, sofosbuvir-based drug treatments have significantly advanced care for hepatitis C virus-infected patients, yet recent reports suggest interactions with amiodarone may cause severe symptomatic bradycardia and thus limit an otherwise extremely effective treatment. Here, we evaluated the ability of human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) to recapitulate the interaction between sofosbuvir and amiodarone in vitro, and more generally assessed the feasibility of hiPSC-CMs as a model system for drug-drug interactions. Sofosbuvir alone had negligible effects on cardiomyocyte electrophysiology, whereas the sofosbuvir-amiodarone combination produced dose-dependent effects beyond that of amiodarone alone. By comparison, GS-331007, the primary circulating metabolite of sofosbuvir, had no effect alone or in combination with amiodarone. Further mechanistic studies revealed that the sofosbuvir-amiodarone combination disrupted intracellular calcium (Ca2+) handling and cellular electrophysiology at pharmacologically relevant concentrations, and mechanical activity at supra-pharmacological (30x Cmax) concentrations. These effects were independent of the common mechanisms of direct ion channel block and P-glycoprotein activity. These results support hiPSC-CMs as a comprehensive, yet scalable model system for the identification and evaluation of cardioactive pharmacodynamic drug-drug interactions.

A bioluminescent assay for measuring glucose uptake.

Identifying activators and inhibitors of glucose uptake is critical for both diabetes management and anticancer therapy. To facilitate such studies, easy-to-use nonradioactive assays are desired. Here we describe a bioluminescent glucose uptake assay for measuring glucose transport in cells. The assay is based on the uptake of 2-deoxyglucose and the enzymatic detection of the 2-deoxyglucose-6-phosphate that accumulates. Uptake can be measured from a variety of cell types, it can be inhibited by known glucose transporter inhibitors, and the bioluminescent assay yields similar results when compared with the radioactive method. With HCT 116 cells, glucose uptake can be detected in as little as 5000 cells and remains linear up to 50,000 cells with signal-to-background values ranging from 5 to 45. The assay can be used to screen for glucose transporter inhibitors, or by multiplexing with viability readouts, changes in glucose uptake can be differentiated from overall effects on cell health. The assay also can provide a relevant end point for measuring insulin sensitivity. With adipocytes and myotubes, insulin-dependent increases in glucose uptake have been measured with 10- and 2-fold assay windows, respectively. Significant assay signals of 2-fold or more have also been measured with human induced pluripotent stem cell (iPSC)-derived cardiomyocytes and skeletal myoblasts.

Disease modeling and phenotypic drug screening for diabetic cardiomyopathy using human induced pluripotent stem cells.

Diabetic cardiomyopathy is a complication of type 2 diabetes, with known contributions of lifestyle and genetics. We develop environmentally and genetically driven in vitro models of the condition using human-induced-pluripotent-stem-cell-derived cardiomyocytes. First, we mimic diabetic clinical chemistry to induce a phenotypic surrogate of diabetic cardiomyopathy, observing structural and functional disarray. Next, we consider genetic effects by deriving cardiomyocytes from two diabetic patients with variable disease progression. The cardiomyopathic phenotype is recapitulated in the patient-specific cells basally, with a severity dependent on their original clinical status. These models are incorporated into successive levels of a screening platform, identifying drugs that preserve cardiomyocyte phenotype in vitro during diabetic stress. In this work, we present a patient-specific induced pluripotent stem cell (iPSC) model of a complex metabolic condition, showing the power of this technique for discovery and testing of therapeutic strategies for a disease with ever-increasing clinical significance.

Phenotypic screening with human iPS cell-derived cardiomyocytes: HTS-compatible assays for interrogating cardiac hypertrophy.

A major hurdle for cardiovascular disease researchers has been the lack of robust and physiologically relevant cell-based assays for drug discovery. Derivation of cardiomyocytes from human-induced pluripotent stem (iPS) cells at high purity, quality, and quantity enables the development of relevant models of human cardiac disease with source material that meets the demands of high-throughput screening (HTS). Here we demonstrate the utility of iPS cell-derived cardiomyocytes as an in vitro model of cardiac hypertrophy. Exposure of cardiomyocytes to endothelin 1 (ET-1) leads to reactivation of fetal genes, increased cell size, and robust expression of B-type natriuretic peptide (BNP). Using this system, we developed a suite of assays focused on BNP detection, most notably a high-content imaging-based assay designed for phenotypic screening. Miniaturization of this assay to a 384-well format enabled the profiling of a small set of tool compounds known to modulate the hypertrophic response. The assays described here provide consistent and reliable results and have the potential to increase our understanding of the many mechanisms underlying this complex cardiac condition. Moreover, the HTS-compatible workflow allows for the incorporation of human biology into early phases of drug discovery and development.

Systematic comparison of three genomic enrichment methods for massively parallel DNA sequencing.

Massively parallel DNA sequencing technologies have greatly increased our ability to generate large amounts of sequencing data at a rapid pace. Several methods have been developed to enrich for genomic regions of interest for targeted sequencing. We have compared three of these methods: Molecular Inversion Probes (MIP), Solution Hybrid Selection (SHS), and Microarray-based Genomic Selection (MGS). Using HapMap DNA samples, we compared each of these methods with respect to their ability to capture an identical set of exons and evolutionarily conserved regions associated with 528 genes (2.61 Mb). For sequence analysis, we developed and used a novel Bayesian genotype-assigning algorithm, Most Probable Genotype (MPG). All three capture methods were effective, but sensitivities (percentage of targeted bases associated with high-quality genotypes) varied for an equivalent amount of pass-filtered sequence: for example, 70% (MIP), 84% (SHS), and 91% (MGS) for 400 Mb. In contrast, all methods yielded similar accuracies of >99.84% when compared to Infinium 1M SNP BeadChip-derived genotypes and >99.998% when compared to 30-fold coverage whole-genome shotgun sequencing data. We also observed a low false-positive rate with all three methods; of the heterozygous positions identified by each of the capture methods, >99.57% agreed with 1M SNP BeadChip, and >98.840% agreed with the whole-genome shotgun data. In addition, we successfully piloted the genomic enrichment of a set of 12 pooled samples via the MGS method using molecular bar codes. We find that these three genomic enrichment methods are highly accurate and practical, with sensitivities comparable to that of 30-fold coverage whole-genome shotgun data.

Repeat subtraction-mediated sequence capture from a complex genome.

Sequence capture technologies, pioneered in mammalian genomes, enable the resequencing of targeted genomic regions. Most capture protocols require blocking DNA, the production of which in large quantities can prove challenging. A blocker-free, two-stage capture protocol was developed using NimbleGen arrays. The first capture depletes the library of repetitive sequences, while the second enriches for target loci. This strategy was used to resequence non-repetitive portions of an approximately 2.2 Mb chromosomal interval and a set of 43 genes dispersed in the 2.3 Gb maize genome. This approach achieved approximately 1800-3000-fold enrichment and 80-98% coverage of targeted bases. More than 2500 SNPs were identified in target genes. Low rates of false-positive SNP predictions were obtained, even in the presence of captured paralogous sequences. Importantly, it was possible to recover novel sequences from non-reference alleles. The ability to design novel repeat-subtraction and target capture arrays makes this technology accessible in any species.

Quantitative trait loci associated with lettuce seed germination under different temperature and light environments.

Temperature and light are primary environmental cues affecting seed germination. To elucidate the genetic architecture underlying lettuce (Lactuca sativa L.) seed germination under different environmental conditions, an F8 recombinant inbred line population consisting of 131 families was phenotyped for final germination and germination rate. Seeds were imbibed in water at 20 degrees C under continuous red light (20-Rc), 20 degrees C continuous dark (20-Dc), 31.5 degrees C continuous red light (31.5-Rc), 31.5 degrees C continuous dark (31.5-Dc), or 20 degrees C far-red light for 24 h followed by continuous dark (20-FRc-Dc). Thirty-eight quantitative trait loci (QTL) were identified from two seed maturation environments: 10 for final germination and 28 for germination rate. The amount of variation attributed to an individual QTL ranged from 9.3% to 17.2% and from 5.6% to 26.2% for final germination and germination rate, respectively. Path analysis indicated that factors affecting germination under 31.5-Rc or 31.5-Dc are largely the same, and these appear to differ from those employed under 20-FRc-Dc. QTL and path analysis support the notion of common and unique factors for germination under diverse temperature and light regimes. A highly significant effect of the seed maturation environment on subsequent germination capacity under environmental stress was observed.

Effect of destrin mutations on the gene expression profile in vivo.

Remodeling of the actin cytoskeleton through actin dynamics (assembly and disassembly of filamentous actin) is known to be essential for numerous basic biological processes. In addition, recent studies have provided evidence that actin dynamics participate in the control of gene expression. A spontaneous mouse mutant, corneal disease 1 (corn1), is deficient for a regulator of actin dynamics, destrin (DSTN, also known as ADF), which causes epithelial hyperproliferation and neovascularization in the cornea. Dstn(corn1) mice exhibit an actin dynamics defect in the corneal epithelial cells, offering an in vivo model to investigate cellular mechanisms affected by the Dstn mutation and resultant actin dynamics abnormalities. To examine the effect of the Dstn(corn1) mutation on the gene expression profile, we performed a microarray analysis using the cornea from Dstn(corn1) and wild-type mice. A dramatic alteration of the gene expression profile was observed in the Dstn(corn1) cornea, with 1,226 annotated genes differentially expressed. Functional annotation of these genes revealed that the most significantly enriched functional categories are associated with actin and/or cytoskeleton. Among genes that belong to these categories, a considerable number of serum response factor target genes were found, indicating the possible existence of an actin-SRF pathway of transcriptional regulation in vivo. A comparative study using an allelic mutant strain with milder corneal phenotypes suggested that the level of filamentous actin may correlate with the level of gene expression changes. Our study shows that Dstn mutations and resultant actin dynamics abnormalities have a strong impact on the gene expression profile in vivo.

Genetic modification of the schisis phenotype in a mouse model of X-linked retinoschisis.

X-linked retinoschisis (XLRS) is an inherited form of macular degeneration that is caused by mutations in the retinoschisin (RS1) gene. In addition to macular degeneration, other major characteristics of XLRS include splitting of the retina (schisis) and impaired synaptic transmission as indicated by a reduction in the electroretinogram b-wave. It has been known that patients carrying RS1 mutations show a broad range of phenotypic variability. Interestingly, phenotypic variation is observed even among family members with the same RS1 mutation, suggesting the existence of genetic or environmental factors that contribute to the severity of XLRS. However, in the human population, the cause of phenotypic variability and the contribution of genetic modifiers for this relatively rare disease are difficult to study and poorly understood. In this study, using a mouse model for XLRS, we show that genetic factors can contribute to the severity of the retinoschisis phenotype. We report evidence of a major genetic modifier of Rs1, which affects the disease severity in these animals. A quantitative trait locus (QTL), named modifier of Rs1 1 (Mor1), is mapped on chromosome (Chr) 7. When homozygous, the Mor1 allele from the inbred mouse strain AKR/J diminishes the severity of the schisis phenotype in Rs1(tmgc1)/Y male and Rs1(tmgc1)/Rs1(tmgc1) female mice. We also show that the penetrance of the disease phenotype is affected by additional genetic factor(s). Our study suggests that multiple genetic modifiers could potentially be responsible for the phenotypic variation in human XLRS.

The Arabidopsis aleurone layer responds to nitric oxide, gibberellin, and abscisic acid and is sufficient and necessary for seed dormancy.

Seed dormancy is a common phase of the plant life cycle, and several parts of the seed can contribute to dormancy. Whole seeds, seeds lacking the testa, embryos, and isolated aleurone layers of Arabidopsis (Arabidopsis thaliana) were used in experiments designed to identify components of the Arabidopsis seed that contribute to seed dormancy and to learn more about how dormancy and germination are regulated in this species. The aleurone layer was found to be the primary determinant of seed dormancy. Embryos from dormant seeds, however, had a lesser growth potential than those from nondormant seeds. Arabidopsis aleurone cells were examined by light and electron microscopy, and cell ultrastructure was similar to that of cereal aleurone cells. Arabidopsis aleurone cells responded to nitric oxide (NO), gibberellin (GA), and abscisic acid, with NO being upstream of GA in a signaling pathway that leads to vacuolation of protein storage vacuoles and abscisic acid inhibiting vacuolation. Molecular changes that occurred in embryos and aleurone layers prior to germination were measured, and these data show that both the aleurone layer and the embryo expressed the NO-associated gene AtNOS1, but only the embryo expressed genes for the GA biosynthetic enzyme GA3 oxidase.