Conrad C. Huang

UCSF CHIMERA-A VISUALIZATION SYSTEM FOR EXPLORATORY RESEARCH AND ANALYSIS

The design, implementation, and capabilities of an extensible visualization system, UCSF Chimera, are discussed. Chimera is segmented into a core that provides basic services and visualization, and extensions that provide most higher level functionality. This architecture ensures that the extension mechanism satisfies the demands of outside developers who wish to incorporate new features. Two unusual extensions are presented: Multiscale, which adds the ability to visualize large‐scale molecular assemblies such as viral coats, and Collaboratory, which allows researchers to share a Chimera session interactively despite being at separate locales. Other extensions include Multalign Viewer, for showing multiple sequence alignments and associated structures; ViewDock, for screening docked ligand orientations; Movie, for replaying molecular dynamics trajectories; and Volume Viewer, for display and analysis of volumetric data. A discussion of the usage of Chimera in real‐world situations is given, along with anticipated future directions. Chimera includes full user documentation, is free to academic and nonprofit users, and is available for Microsoft Windows, Linux, Apple Mac OS X, SGI IRIX, and HP Tru64 Unix from http://www.cgl.ucsf.edu/chimera/.

The MIDAS display system

Abstract The Molecular Interactive Display and Simulation (MIDAS) system is designed to display and manipulate large macromolecules, such as proteins and nucleic acids. Several ancillary programs allow for such features as computing the surface of a molecule, selecting an active site region within a molecule, and computing electrostatic charge potentials. At the core of MIDAS is a hierarchical database system, designed specifically for macromolecules, that is both compact in its storage requirements and fast in its data access.

Tools for integrated sequence-structure analysis with UCSF Chimera

BackgroundComparing related structures and viewing the structures in the context of sequence alignments are important tasks in protein structure-function research. While many programs exist for individual aspects of such work, there is a need for interactive visualization tools that: (a) provide a deep integration of sequence and structure, far beyond mapping where a sequence region falls in the structure and vice versa; (b) facilitate changing data of one type based on the other (for example, using only sequence-conserved residues to match structures, or adjusting a sequence alignment based on spatial fit); (c) can be used with a researchers own data, including arbitrary sequence alignments and annotations, closely or distantly related sets of proteins, etc.; and (d) interoperate with each other and with a full complement of molecular graphics features. We describe enhancements to UCSF Chimera to achieve these goals.ResultsThe molecular graphics program UCSF Chimera includes a suite of tools for interactive analyses of sequences and structures. Structures automatically associate with sequences in imported alignments, allowing many kinds of crosstalk. A novel method is provided to superimpose structures in the absence of a pre-existing sequence alignment. The method uses both sequence and secondary structure, and can match even structures with very low sequence identity. Another tool constructs structure-based sequence alignments from superpositions of two or more proteins. Chimera is designed to be extensible, and mechanisms for incorporating user-specific data without Chimera code development are also provided.ConclusionThe tools described here apply to many problems involving comparison and analysis of protein structures and their sequences. Chimera includes complete documentation and is intended for use by a wide range of scientists, not just those in the computational disciplines. UCSF Chimera is free for non-commercial use and is available for Microsoft Windows, Apple Mac OS X, Linux, and other platforms from http://www.cgl.ucsf.edu/chimera.

BayGenomics: a resource of insertional mutations in mouse embryonic stem cells

The BayGenomics gene-trap resource (http://baygenomics.ucsf.edu) provides researchers with access to thousands of mouse embryonic stem (ES) cell lines harboring characterized insertional mutations in both known and novel genes. Each cell line contains an insertional mutation in a specific gene. The identity of the gene that has been interrupted can be determined from a DNA sequence tag. Approximately 75% of our cell lines contain insertional mutations in known mouse genes or genes that share strong sequence similarities with genes that have been identified in other organisms. These cell lines readily transmit the mutation to the germline of mice and many mutant lines of mice have already been generated from this resource. BayGenomics provides facile access to our entire database, including sequence tags for each mutant ES cell line, through the World Wide Web. Investigators can browse our resource, search for specific entries, download any portion of our database and BLAST sequences of interest against our entire set of cell line sequence tags. They can then obtain the mutant ES cell line for the purpose of generating knockout mice.

Evolutionary conservation predicts function of variants of the human organic cation transporter, OCT1

The organic cation transporter, OCT1, is a major hepatic transporter that mediates the uptake of many organic cations from the blood into the liver where the compounds may be metabolized or secreted into the bile. Because OCT1 interacts with a variety of structurally diverse organic cations, including clinically used drugs as well as toxic substances (e.g., N-methylpyridinium, MPP+), it is an important determinant of systemic exposure to many xenobiotics. To understand the genetic basis of extensive interindividual differences in xenobiotic disposition, we functionally characterized 15 protein-altering variants of the human liver organic cation transporter, OCT1, in Xenopus oocytes. All variants that reduced or eliminated function (OCT1-R61C, OCT1-P341L, OCT1-G220V, OCT1-G401S, and OCT1-G465R) altered evolutionarily conserved amino acid residues. In general, variants with decreased function had amino acid substitutions that resulted in more radical chemical changes (higher Grantham values) and were less evolutionarily favorable (lower blosum62 values) than variants that maintained function. A variant with increased function (OCT1-S14F) changed an amino acid residue such that the human protein matched the consensus of the OCT1 mammalian orthologs. Our results indicate that changes at evolutionarily conserved positions of OCT1 are strong predictors of decreased function and suggest that a combination of evolutionary conservation and chemical change might be a stronger predictor of function.

Natural variation in human membrane transporter genes reveals evolutionary and functional constraints

Membrane transporters maintain cellular and organismal homeostasis by importing nutrients and exporting toxic compounds. Transporters also play a crucial role in drug response, serving as drug targets and setting drug levels. As part of a pharmacogenetics project, we screened exons and flanking intronic regions for variation in a set of 24 membrane transporter genes (96 kb; 57% coding) in 247 DNA samples from ethnically diverse populations. We identified 680 single nucleotide polymorphisms (SNPs), of which 175 were synonymous and 155 caused amino acid changes, and 29 small insertions and deletions. Amino acid diversity (πNS) in transmembrane domains (TMDs) was significantly lower than in loop domains, suggesting that TMDs have special functional constraints. This difference was especially striking in the ATP-binding cassette superfamily and did not parallel evolutionary conservation: there was little variation in the TMDs, even in evolutionarily unconserved residues. We used allele frequency distribution to evaluate different scoring systems (Grantham, blosum62, SIFT, and evolutionarily conserved/evolutionarily unconserved) for their ability to predict which SNPs affect function. Our underlying assumption was that alleles that are functionally deleterious will be selected against and thus under represented at high frequencies and over represented at low frequencies. We found that evolutionary conservation of orthologous sequences, as assessed by evolutionarily conserved/evolutionarily unconserved and SIFT, was the best predictor of allele frequency distribution and hence of function. European Americans had an excess of high frequency alleles in comparison to African Americans, consistent with a historic bottleneck. In addition, African Americans exhibited a much higher frequency of population specific medium-frequency alleles than did European Americans.

UCSF Chimera, MODELLER, and IMP: an integrated modeling system.

Structural modeling of macromolecular complexes greatly benefits from interactive visualization capabilities. Here we present the integration of several modeling tools into UCSF Chimera. These include comparative modeling by MODELLER, simultaneous fitting of multiple components into electron microscopy density maps by IMP MultiFit, computing of small-angle X-ray scattering profiles and fitting of the corresponding experimental profile by IMP FoXS, and assessment of amino acid sidechain conformations based on rotamer probabilities and local interactions by Chimera.

Polymorphisms in a human kidney xenobiotic transporter, OCT2, exhibit altered function

The completion of the Human Genome Project and the development of high-throughput polymorphism identification methods have allowed researchers to carry out full genetic analyses of many clinically relevant genes. However, few studies have combined genetic analysis with in vitro phenotyping to better understand the relationship between genetic variation and protein function. Many transporters in the kidney are thought to play key roles in defense against a variety of foreign substances. The goal of this study was to understand the relationship between variation in a gene encoding a major renal xenobiotic transporter, OCT2, and transporter function. We report a comprehensive genetic analysis and functional characterization of variants of OCT2. Twenty-eight variable sites in the OCT2 gene were identified in a collection of 247 ethnically diverse DNA samples. Eight caused non-synonymous amino acid changes, of which four were present at >/= 1% in an ethnic population. All four of these altered transporter function assayed in Xenopus laevis oocytes. Analysis of nucleotide diversity (pi) revealed a higher prevalence of synonymous (pi = 22.4 x 10-4) versus non-synonymous (pi = 2.1 x 10-4) changes in OCT2 than in other genes. In addition, the non-synonymous sites had a significant tendency to exhibit more skewed allele frequencies (more negative Tajimas D-values) compared to synonymous sites. The population-genetic analysis, together with the functional characterization, suggests that selection has acted against amino acid changes in OCT2. This selection may be due to a necessary role of OCT2 in the renal elimination of endogenous amines or xenobiotics, including environmental toxins, neurotoxic amines and therapeutic drugs.

Interaction of Methotrexate with Organic-Anion Transporting Polypeptide 1A2 and Its Genetic Variants

Methotrexate (MTX) is used in patients with malignant and autoimmune diseases. This drug is primarily excreted unchanged in the urine, and its net excretion occurs via active secretory and reabsorptive processes. We characterized the interaction of MTX with human organic-anion transporting polypeptide transporter (OATP) 1A2, which is expressed in tissues important for MTX disposition and toxicity, such as the intestine, kidney, liver, and endothelial cells of the blood-brain barrier. In Xenopus laevis oocytes expressing OATP1A2, the uptake of the model substrate, estrone-3-sulfate (ES), was enhanced 30-fold compared with uninjected oocytes. MTX uptake in oocytes expressing OATP1A2 was saturable (Km = 457 ± 118 μM; Vmax = 17.5 ± 4.9 pmol/oocyte/60 min) and sensitive to extracellular pH. That is, acidic pHs stimulated MTX uptake by as much as 7-fold. Seven novel protein-altering variants were identified in 270 ethnically diverse DNA samples. Four protein-altering variants in OATP1A2 exhibited altered transport of ES and/or MTX. The common variant, protein reference sequence (p.) Ile13Thr, was hyperfunctional for ES and MTX and showed a 2-fold increase in the Vmax for ES. The common variant, p. Glu172Asp, exhibited reduced maximal transport capacity for ES and MTX. p. Arg168Cys was hypofunctional, and p. Asn277DEL was nonfunctional. Because of its expression on the apical membrane of the distal tubule and in tissues relevant to MTX disposition and toxicity, these findings suggest that OATP1A2 may play a role in active tubular reabsorption of MTX and in MTX-induced toxicities. Furthermore, genetic variation in OATP1A2 may contribute to variation in MTX disposition and response.

The International Gene Trap Consortium Website: a portal to all publicly available gene trap cell lines in mouse

Gene trapping is a method of generating murine embryonic stem (ES) cell lines containing insertional mutations in known and novel genes. A number of international groups have used this approach to create sizeable public cell line repositories available to the scientific community for the generation of mutant mouse strains. The major gene trapping groups worldwide have recently joined together to centralize access to all publicly available gene trap lines by developing a user-oriented Website for the International Gene Trap Consortium (IGTC). This collaboration provides an impressive public informatics resource comprising ∼45 000 well-characterized ES cell lines which currently represent ∼40% of known mouse genes, all freely available for the creation of knockout mice on a non-collaborative basis. To standardize annotation and provide high confidence data for gene trap lines, a rigorous identification and annotation pipeline has been developed combining genomic localization and transcript alignment of gene trap sequence tags to identify trapped loci. This information is stored in a new bioinformatics database accessible through the IGTC Website interface. The IGTC Website () allows users to browse and search the database for trapped genes, BLAST sequences against gene trap sequence tags, and view trapped genes within biological pathways. In addition, IGTC data have been integrated into major genome browsers and bioinformatics sites to provide users with outside portals for viewing this data. The development of the IGTC Website marks a major advance by providing the research community with the data and tools necessary to effectively use public gene trap resources for the large-scale characterization of mammalian gene function.