Transcriptomic variation of pharmacogenes in multiple human tissues and lymphoblastoid cell lines.

Variation in the expression level and activity of genes involved in drug disposition and action ('pharmacogenes') can affect drug response and toxicity, especially when in tissues of pharmacological importance. Previous studies have relied primarily on microarrays to understand gene expression differences, or have focused on a single tissue or small number of samples. The goal of this study was to use RNA-sequencing (RNA-seq) to determine the expression levels and alternative splicing of 389 Pharmacogenomics Research Network pharmacogenes across four tissues (liver, kidney, heart and adipose) and lymphoblastoid cell lines, which are used widely in pharmacogenomics studies. Analysis of RNA-seq data from 139 different individuals across the 5 tissues (20-45 individuals per tissue type) revealed substantial variation in both expression levels and splicing across samples and tissue types. Comparison with GTEx data yielded a consistent picture. This in-depth exploration also revealed 183 splicing events in pharmacogenes that were previously not annotated. Overall, this study serves as a rich resource for the research community to inform biomarker and drug discovery and use.

Genetic and epigenetic regulation of the organic cation transporter 3, SLC22A3.

Human organic cation transporter 3 (OCT3 and SLC22A3) mediates the uptake of many important endogenous amines and basic drugs in a variety of tissues. OCT3 is identified as one of the important risk loci for prostate cancer, and is markedly underexpressed in aggressive prostate cancers. The goal of this study was to identify genetic and epigenetic factors in the promoter region that influence the expression level of OCT3. Haplotypes that contained the common variants, g.-81G>delGA (rs60515630) (minor allele frequency 11.5% in African American) and g.-2G>A (rs555754) (minor allele frequency>30% in all ethnic groups) showed significant increases in luciferase reporter activities and exhibited stronger transcription factor-binding affinity than the haplotypes that contained the major alleles. Consistent with the reporter assays, OCT3 messenger RNA expression levels were significantly higher in Asian (P<0.001) and Caucasian (P<0.05) liver samples from individuals who were homozygous for g.-2A/A in comparison with those homozygous for the g.-2G/G allele. Studies revealed that the methylation level in the basal promoter region of OCT3 was associated with OCT3 expression level and tumorigenesis capability in various prostate cancer cell lines. The methylation level of the OCT3 promoter was higher in 62% of prostate tumor samples compared with matched normal samples. Our studies demonstrate that genetic polymorphisms in the proximal promoter region of OCT3 alter the transcription rate of the gene and may be associated with altered expression levels of OCT3 in human liver. Aberrant methylation contributes to the reduced expression of OCT3 in prostate cancer.

A common 5'-UTR variant in MATE2-K is associated with poor response to metformin.

Multidrug and toxin extrusion 2 (MATE2-K (SLC47A2)), a polyspecific organic cation exporter, facilitates the renal elimination of the antidiabetes drug metformin. In this study, we characterized genetic variants of MATE2-K, determined their association with metformin response, and elucidated their impact by means of a comparative protein structure model. Four nonsynonymous variants and four variants in the MATE2-K basal promoter region were identified from ethnically diverse populations. Two nonsynonymous variants-c.485C>T and c.1177G>A-were shown to be associated with significantly lower metformin uptake and reduction in protein expression levels. MATE2-K basal promoter haplotypes containing the most common variant, g.-130G>A (>26% allele frequency), were associated with a significant increase in luciferase activities and reduced binding to the transcriptional repressor myeloid zinc finger 1 (MZF-1). Patients with diabetes who were homozygous for g.-130A had a significantly poorer response to metformin treatment, assessed as relative change in glycated hemoglobin (HbA1c) (-0.027 (-0.076, 0.033)), as compared with carriers of the reference allele, g.-130G (-0.15 (-0.17, -0.13)) (P=0.002). Our study showed that MATE2-K plays a role in the antidiabetes response to metformin.

Functional characterization of liver enhancers that regulate drug-associated transporters.

Little is known about how genetic variations in enhancers influence drug response. In this study, we investigated whether nucleotide variations in enhancers that regulate drug transporters can alter their expression levels. Using comparative genomics and liver-specific transcription factor binding site (TFBS) analyses, we identified evolutionary conserved regions (ECRs) surrounding nine liver membrane transporters that interact with commonly used pharmaceuticals. The top 50 ECRs were screened for enhancer activity in vivo, of which five--located around ABCB11, SLC10A1, SLCO1B1, SLCO1A2, and SLC47A1--exhibited significant enhancer activity. Common variants identified in a large ethnically diverse cohort (n = 272) were assayed for differential enhancer activity, and three variants were found to have significant effects on reporter activity as compared with the reference allele. In addition, one variant was associated with reduced SLCO1A2 mRNA expression levels in human liver tissues, and another was associated with increased methotrexate (MTX) clearance in patients. This work provides a general model for the rapid characterization of liver enhancers and identifies associations between enhancer variants and drug response.

Functional analysis of genetic variants in the human concentrative nucleoside transporter 3 (CNT3; SLC28A3).

The human concentrative nucleoside transporter, CNT3 (SLC28A3), plays an important role in mediating the cellular entry of a broad array of physiological nucleosides and synthetic anticancer nucleoside analog drugs. As a first step toward understanding the genetic basis for interindividual differences in the disposition and response to antileukemic nucleoside analogs, we examined the genetic and functional diversity of CNT3. In all, 56 variable sites in the exons and flanking intronic region of SLC28A3 were identified in a collection of 270 DNA samples from US populations (80 African-Americans, 80 European-Americans, 60 Asian-Americans, and 50 Mexican-Americans). Of the 16 coding region variants, 12 had not been previously reported. Also, 10 resulted in amino-acid changes and three of these had total allele frequencies of >/=1%. Nucleotide diversity (pi) at nonsynonymous and synonymous sites was estimated to be 1.81 x 10(4) and 18.13 x 10(4), respectively, suggesting that SLC28A3 is under negative selection. All nonsynonymous variants, constructed by site-directed mutagenesis and expressed in Xenopus laevis oocytes, transported purine and pyrimidine model substrates, except for c. 1099G>A (p. Gly367Arg). This rare variant alters an evolutionarily conserved site in the putative substrate recognition domain of CNT3. The presence of three additional evolutionarily conserved glycine residues in the vicinity of p. Gly367Arg that are also conserved in human paralogs suggest that these glycine residues are critical in the function of the concentrative nucleoside transporter family. The genetic analysis and functional characterization of CNT3 variants suggest that this transporter does not tolerate nonsynonymous changes and is important for human fitness.

Integrated tools for structural and sequence alignment and analysis.

We have developed new computational methods for displaying and analyzing members of protein superfamilies. These methods (MinRMS, AlignPlot and MSFviewer) integrate sequence and structural information and are implemented as separate but cooperating programs to our Chimera molecular modeling system. Integration of multiple sequence alignment information and three-dimensional structural representations enable researchers to generate hypotheses about the sequence-structure relationship. Structural superpositions can be generated and easily tuned to identify similarities around important characteristics such as active sites or ligand binding sites. Information related to the release of Chimera, MinRMS, AlignPlot and MSFviewer can be obtained at http:¿www.cgl.ucsf.edu/chimera.

The object technology framework: an object-oriented interface to molecular data and its application to collagen.

We describe the Object Technology Framework (OTF) software system developed at the University of California, San Francisco Computer Graphics Laboratory for creating C+2 classes that facilitate rapid biomolecular application development and the application of the OTF to collagen modeling. C+2 class libraries for accessing and manipulating data from standard scientific data sources can be generated using the program genlib and its class library toolkit Molecule, thereby facilitating development of new applications. Use of the OTF for generating ideal collagen structural models (gencollagen) is described. The source code for the OTF is freely available at http:/(/)www.cgl.ucsf.edu/off/ to interested application developers.

Molecules to maps: tools for visualization and interaction in support of computational biology.

The volume of data produced by genome projects, X-ray crystallography, NMR spectroscopy, and electron and confocal microscopy present the bioinformatics community with new challenges for analyzing, understanding, and exchanging this data. At the 1998 Pacific Symposium on Biocomputing, a track entitled 'Molecules to Maps: Tools for Visualization and Interaction in Computational Biology' provided tool developers and users with the opportunity to discuss advances in tools and techniques to assist scientists in evaluating, absorbing, navigating, and correlating this sea of information, through visualization and user interaction. In this paper we present these advances and discuss some of the challenges that remain to be solved.

Annotating PDB files with scene information.

We have implemented extensions to the Brookhaven Protein Data Bank (PDB) file format for incorporating scene information such as viewing parameters, additional molecular information (e.g., van der Waals radii and atom colors), and user-defined graphics. These extensions were made in conformance with the PDB standard and provide sufficient information to render the scene in various styles such as space-filling images and ribbon diagrams. For the past 5 years these extensions have been used in the MidasPlus molecular modeling system and have proved both powerful and sufficient for generating complex molecular images. We propose that the extensions to the PDB presented here be adopted by the molecular modeling community for incorporation into visualization programs.

Automated site-directed drug design using molecular lattices.

Receptor-based drug design is predicated on the knowledge of the structure of a target receptor and the principles of molecule recognition. The objective is to produce a wide diversity of structures that are sterically and electrostatically complementary to a specified receptor site. Many drug-receptor interactions are controlled by a few key receptor groups. This observation leads to a design approach in which one focuses on chemical fragments that putatively interact with the key receptor groups. There then remains the difficult task of joining the fragments into molecular structures that match the spatial patterns of recognition forces in the receptor site. In this paper, we describe a new modeling program, BUILDER, that combines database searching techniques and structure generation algorithms within an interactive graphics modeling environment (MidasPlus). A novel tool for process communication (delegate) is introduced and examples of its use are given. To demonstrate the functionality of the package and its ability to produce novel structures, we examine the active site of HIV-1 protease.

Conic: a fast renderer for space-filling molecules with shadows.

We present an algorithm for generating images of molecules represented as a set of intersecting opaque spheres. Both perspective and shadows are computed to provide realistic visual cues. Compared to existing programs for generating similar images, our algorithm is both more accurate and several times faster. We present in detail the mathematics used in picture generation, along with examples of the computed images.

Molecular structure and dynamics of tricyclic antidepressant drugs.

The molecular structure, electrostatic potentials and dynamics of imipramine, chlorimipramine, amitriptyline and nortriptyline were examined by computer graphics, molecular mechanical energy calculations and molecular dynamics simulations, using the AMBER all atom force field. Starting coordinates for amitriptyline and nortriptyline were generated by model building from the crystal structure of imipramine. The structures were refined by molecular mechanical energy minimization, and used as starting points for molecular dynamics simulations in vacuo and in aqueous solution. The simulations demonstrated considerable flexibility of the molecules, both in the side chain and in the ring system, where the angle between the phenyl rings varied between 90 degrees and 168 degrees. The most frequently observed conformations of imipramine, chlorimipramine and nortriptyline during the simulations had the side chain folded above one of the phenyl rings, while amitriptyline showed both folded and extended side chain conformations during the simulations. The results may provide increased understanding of the molecular recognition and specificity of tricyclic antidepressant drugs in interaction with neurotransmitter receptor molecules.

An affordable approach to interactive desktop molecular modeling.

The amazing revolution in computer hardware performance and cost reduction has yet to be carried over to computer software. In fact, application software today is often more expensive and less reliable than the hardware. New enhancements in software development techniques, such as object oriented programming and interactive graphics based user interface design, finally may be having a significant impact on the time-to-market and reliability of these application programs. We discuss our experiences using one such set of software development tools available on the NeXT workstation and describe the effort required to port our MidasPlus molecular modeling package to the NeXT work-station.

A real-time malleable molecular surface.

We describe a method for generating a molecular surface using a parametric patch representation. Unlike previous methods, this algorithm generates a parametric patch surface which is smooth and G continuous and manipulable in real-time. Crucial to our approach is the creation of a net of approximately equilateral triangles from which we generate the control points used as the basis for describing the surface. We present in detail the method used for generating the triangular net and accompanying control points, along with examples of the resulting surfaces.

Van der waals surfaces in molecular modeling: implementation with real-time computer graphics.

A method is described for generating van der Waals molecular surfaces with a real-time interactive calligraphic color display system. These surfaces maintain their proper representation during bond rotation and global transformations, and an interior atom removal method yields a comprehensible picture of the molecular surface for large molecules. Both algorithms are faster than previous methods. This combination provides a powerful tool for real-time interactive molecular modeling.

Computer graphics in drug design: molecular modeling of thyroid hormone-prealbumin interactions.

Computer graphics modeling of the thyroxine-prealbumin complex provides a detailed picture of the interactions between thyroxine and prealbumin. A wide variety of thyroid hormone analogue-prealbumin complexes were modeled by calculating the molecular surfaces of the analogues and the prealbumin hormone-binding site. Analogues with high binding affinity were observed to fill more of the hormone-binding site than low-affinity analogues. These surface models described many aspects of the hormone-protein interaction which were not obvious using simple wire models and led us to develop a model which accounts for thyroid hormone-prealbumin structure-activity relationships and ultimately to predict and measure the relative binding affinities of four previously untested thyroid hormone analogues to prealbumin.

Real-time color graphics in studies of molecular interactions.

Studies of the structures and interactions of large biological molecules require both coordinate data and three-dimensional visualization. Orthodox molecular models often bear a tenuous relationship to the coordinate data. In contrast, computer graphics requires that the display directly and accurately represent the data, and storage of modified configurations and recovery of original structures are simple. Software has been developed that allows real-time display of color line and surface displays of several interacting molecules, while quantitatively monitoring the stereochemistry.