An Atlas of Peroxiredoxins Created Using an Active Site Profile-Based Approach to Functionally Relevant Clustering of Proteins.

Peroxiredoxins (Prxs or Prdxs) are a large protein superfamily of antioxidant enzymes that rapidly detoxify damaging peroxides and/or affect signal transduction and, thus, have roles in proliferation, differentiation, and apoptosis. Prx superfamily members are widespread across phylogeny and multiple methods have been developed to classify them. Here we present an updated atlas of the Prx superfamily identified using a novel method called MISST (Multi-level Iterative Sequence Searching Technique). MISST is an iterative search process developed to be both agglomerative, to add sequences containing similar functional site features, and divisive, to split groups when functional site features suggest distinct functionally-relevant clusters. Superfamily members need not be identified initially-MISST begins with a minimal representative set of known structures and searches GenBank iteratively. Further, the method's novelty lies in the manner in which isofunctional groups are selected; rather than use a single or shifting threshold to identify clusters, the groups are deemed isofunctional when they pass a self-identification criterion, such that the group identifies itself and nothing else in a search of GenBank. The method was preliminarily validated on the Prxs, as the Prxs presented challenges of both agglomeration and division. For example, previous sequence analysis clustered the Prx functional families Prx1 and Prx6 into one group. Subsequent expert analysis clearly identified Prx6 as a distinct functionally relevant group. The MISST process distinguishes these two closely related, though functionally distinct, families. Through MISST search iterations, over 38,000 Prx sequences were identified, which the method divided into six isofunctional clusters, consistent with previous expert analysis. The results represent the most complete computational functional analysis of proteins comprising the Prx superfamily. The feasibility of this novel method is demonstrated by the Prx superfamily results, laying the foundation for potential functionally relevant clustering of the universe of protein sequences.

DASP3: identification of protein sequences belonging to functionally relevant groups.

BACKGROUND: Development of automatable processes for clustering proteins into functionally relevant groups is a critical hurdle as an increasing number of sequences are deposited into databases. Experimental function determination is exceptionally time-consuming and can't keep pace with the identification of protein sequences. A tool, DASP (Deacon Active Site Profiler), was previously developed to identify protein sequences with active site similarity to a query set. Development of two iterative, automatable methods for clustering proteins into functionally relevant groups exposed algorithmic limitations to DASP. RESULTS: The accuracy and efficiency of DASP was significantly improved through six algorithmic enhancements implemented in two stages: DASP2 and DASP3. Validation demonstrated DASP3 provides greater score separation between true positives and false positives than earlier versions. In addition, DASP3 shows similar performance to previous versions in clustering protein structures into isofunctional groups (validated against manual curation), but DASP3 gathers and clusters protein sequences into isofunctional groups more efficiently than DASP and DASP2. CONCLUSIONS: DASP algorithmic enhancements resulted in improved efficiency and accuracy of identifying proteins that contain active site features similar to those of the query set. These enhancements provide incremental improvement in structure database searches and initial sequence database searches; however, the enhancements show significant improvement in iterative sequence searches, suggesting DASP3 is an appropriate tool for the iterative processes required for clustering proteins into isofunctional groups.

RRDistMaps: a UCSF Chimera tool for viewing and comparing protein distance maps.

MOTIVATION: Contact maps are a convenient method for the structural biologists to identify structural features through two-dimensional simplification. Binary (yes/no) contact maps with a single cutoff distance can be generalized to show continuous distance ranges. We have developed a UCSF Chimera tool, RRDistMaps, to compute such generalized maps in order to analyze pairwise variations in intramolecular contacts. An interactive utility, RRDistMaps, visualizes conformational changes, both local (e.g. binding-site residues) and global (e.g. hinge motion), between unbound and bound proteins through distance patterns. Users can target residue pairs in RRDistMaps for further navigation in Chimera. The interface contains the unique features of identifying long-range residue motion and aligning sequences to simultaneously compare distance maps. AVAILABILITY AND IMPLEMENTATION: RRDistMaps was developed as part of UCSF Chimera release 1.10, which is freely available at http://rbvi.ucsf.edu/chimera/download.html, and operates on Linux, Windows, and Mac OS. CONTACT: conrad@cgl.ucsf.edu.

cddApp: a Cytoscape app for accessing the NCBI conserved domain database.

MOTIVATION: cddApp is a Cytoscape extension that supports the annotation of protein networks with information about domains and specific functional sites from the National Center for Biotechnology Information's conserved domain database (CDD). CDD information is loaded for nodes annotated with NCBI numbers or UniProt identifiers and (optionally) Protein Data Bank structures. cddApp integrates with the Cytoscape apps structureViz2 and enhancedGraphics. Together, these three apps provide powerful tools to annotate nodes with CDD domain and site information and visualize that information in both network and structural contexts. AVAILABILITY AND IMPLEMENTATION: cddApp is written in Java and freely available for download from the Cytoscape app store (http://apps.cytoscape.org). Documentation is provided at http://www.rbvi.ucsf.edu/cytoscape, and the source is publically available from GitHub http://github.com/RBVI/cddApp.

Enhancing UCSF Chimera through web services.

Integrating access to web services with desktop applications allows for an expanded set of application features, including performing computationally intensive tasks and convenient searches of databases. We describe how we have enhanced UCSF Chimera (http://www.rbvi.ucsf.edu/chimera/), a program for the interactive visualization and analysis of molecular structures and related data, through the addition of several web services (http://www.rbvi.ucsf.edu/chimera/docs/webservices.html). By streamlining access to web services, including the entire job submission, monitoring and retrieval process, Chimera makes it simpler for users to focus on their science projects rather than data manipulation. Chimera uses Opal, a toolkit for wrapping scientific applications as web services, to provide scalable and transparent access to several popular software packages. We illustrate Chimera's use of web services with an example workflow that interleaves use of these services with interactive manipulation of molecular sequences and structures, and we provide an example Python program to demonstrate how easily Opal-based web services can be accessed from within an application. Web server availability: http://webservices.rbvi.ucsf.edu/opal2/dashboard?command=serviceList.

The Structure-Function Linkage Database.

The Structure-Function Linkage Database (SFLD, http://sfld.rbvi.ucsf.edu/) is a manually curated classification resource describing structure-function relationships for functionally diverse enzyme superfamilies. Members of such superfamilies are diverse in their overall reactions yet share a common ancestor and some conserved active site features associated with conserved functional attributes such as a partial reaction. Thus, despite their different functions, members of these superfamilies 'look alike', making them easy to misannotate. To address this complexity and enable rational transfer of functional features to unknowns only for those members for which we have sufficient functional information, we subdivide superfamily members into subgroups using sequence information, and lastly into families, sets of enzymes known to catalyze the same reaction using the same mechanistic strategy. Browsing and searching options in the SFLD provide access to all of these levels. The SFLD offers manually curated as well as automatically classified superfamily sets, both accompanied by search and download options for all hierarchical levels. Additional information includes multiple sequence alignments, tab-separated files of functional and other attributes, and sequence similarity networks. The latter provide a new and intuitively powerful way to visualize functional trends mapped to the context of sequence similarity.

Multidomain Assembler (MDA) Generates Models of Large Multidomain Proteins.

Homology modeling predicts protein structures using known structures of related proteins as templates. We developed MULTIDOMAIN ASSEMBLER (MDA) to address the special problems that arise when modeling proteins with large numbers of domains, such as fibronectin with 30 domains, as well as cases with hundreds of templates. These problems include how to spatially arrange nonoverlapping template structures, and how to get the best template coverage when some sequence regions have hundreds of available structures while other regions have a few distant homologs. MDA automates the tasks of template searching, visualization, and selection followed by multidomain model generation, and is part of the widely used molecular graphics package UCSF CHIMERA (University of California, San Francisco). We demonstrate applications and discuss MDA's benefits and limitations.

IHMCIF: An Extension of the PDBx/mmCIF Data Standard for Integrative Structure Determination Methods.

IHMCIF (github.com/ihmwg/IHMCIF) is a data information framework that supports archiving and disseminating macromolecular structures determined by integrative or hybrid modeling (IHM), and making them Findable, Accessible, Interoperable, and Reusable (FAIR). IHMCIF is an extension of the Protein Data Bank Exchange/macromolecular Crystallographic Information Framework (PDBx/mmCIF) that serves as the framework for the Protein Data Bank (PDB) to archive experimentally determined atomic structures of biological macromolecules and their complexes with one another and small molecule ligands (e.g., enzyme cofactors and drugs). IHMCIF serves as the foundational data standard for the PDB-Dev prototype system, developed for archiving and disseminating integrative structures. It utilizes a flexible data representation to describe integrative structures that span multiple spatiotemporal scales and structural states with definitions for restraints from a variety of experimental methods contributing to integrative structural biology. The IHMCIF extension was created with the benefit of considerable community input and recommendations gathered by the Worldwide Protein Data Bank (wwPDB) Task Force for Integrative or Hybrid Methods (wwpdb.org/task/hybrid). Herein, we describe the development of IHMCIF to support evolving methodologies and ongoing advancements in integrative structural biology. Ultimately, IHMCIF will facilitate the unification of PDB-Dev data and tools with the PDB archive so that integrative structures can be archived and disseminated through PDB.

ModBase, a database of annotated comparative protein structure models, and associated resources.

ModBase (http://salilab.org/modbase) is a database of annotated comparative protein structure models. The models are calculated by ModPipe, an automated modeling pipeline that relies primarily on Modeller for fold assignment, sequence-structure alignment, model building and model assessment (http://salilab.org/modeller/). ModBase currently contains 10,355,444 reliable models for domains in 2,421,920 unique protein sequences. ModBase allows users to update comparative models on demand, and request modeling of additional sequences through an interface to the ModWeb modeling server (http://salilab.org/modweb). ModBase models are available through the ModBase interface as well as the Protein Model Portal (http://www.proteinmodelportal.org/). Recently developed associated resources include the SALIGN server for multiple sequence and structure alignment (http://salilab.org/salign), the ModEval server for predicting the accuracy of protein structure models (http://salilab.org/modeval), the PCSS server for predicting which peptides bind to a given protein (http://salilab.org/pcss) and the FoXS server for calculating and fitting Small Angle X-ray Scattering profiles (http://salilab.org/foxs).

UCSF ChimeraX: Tools for structure building and analysis.

Advances in computational tools for atomic model building are leading to accurate models of large molecular assemblies seen in electron microscopy, often at challenging resolutions of 3-4 Å. We describe new methods in the UCSF ChimeraX molecular modeling package that take advantage of machine-learning structure predictions, provide likelihood-based fitting in maps, and compute per-residue scores to identify modeling errors. Additional model-building tools assist analysis of mutations, post-translational modifications, and interactions with ligands. We present the latest ChimeraX model-building capabilities, including several community-developed extensions. ChimeraX is available free of charge for noncommercial use at https://www.rbvi.ucsf.edu/chimerax.

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.

Improving the quality of protein similarity network clustering algorithms using the network edge weight distribution.

MOTIVATION: Clustering protein sequence data into functionally specific families is a difficult but important problem in biological research. One useful approach for tackling this problem involves representing the sequence dataset as a protein similarity network, and afterwards clustering the network using advanced graph analysis techniques. Although a multitude of such network clustering algorithms have been developed over the past few years, comparing algorithms is often difficult because performance is affected by the specifics of network construction. We investigate an important aspect of network construction used in analyzing protein superfamilies and present a heuristic approach for improving the performance of several algorithms. RESULTS: We analyzed how the performance of network clustering algorithms relates to thresholding the network prior to clustering. Our results, over four different datasets, show how for each input dataset there exists an optimal threshold range over which an algorithm generates its most accurate clustering output. Our results further show how the optimal threshold range correlates with the shape of the edge weight distribution for the input similarity network. We used this correlation to develop an automated threshold selection heuristic in order to most optimally filter a similarity network prior to clustering. This heuristic allows researchers to process their protein datasets with runtime efficient network clustering algorithms without sacrificing the clustering accuracy of the final results. AVAILABILITY: Python code for implementing the automated threshold selection heuristic, together with the datasets used in our analysis, are available at http://www.rbvi.ucsf.edu/Research/cytoscape/threshold_scripts.zip.

Designed divergent evolution of enzyme function.

It is generally believed that proteins with promiscuous functions divergently evolved to acquire higher specificity and activity, and that this process was highly dependent on the ability of proteins to alter their functions with a small number of amino acid substitutions (plasticity). The application of this theory of divergent molecular evolution to promiscuous enzymes may allow us to design enzymes with more specificity and higher activity. Many structural and biochemical analyses have identified the active or binding site residues important for functional plasticity (plasticity residues). To understand how these residues contribute to molecular evolution, and thereby formulate a design methodology, plasticity residues were probed in the active site of the promiscuous sesquiterpene synthase gamma-humulene synthase. Identified plasticity residues were systematically recombined based on a mathematical model in order to construct novel terpene synthases, each catalysing the synthesis of one or a few very different sesquiterpenes. Here we present the construction of seven specific and active synthases that use different reaction pathways to produce the specific and very different products. Creation of these enzymes demonstrates the feasibility of exploiting the underlying evolvability of this scaffold, and provides evidence that rational approaches based on these ideas are useful for enzyme design.

Computational tools for the interactive exploration of proteomic and structural data.

Linking proteomics and structural data is critical to our understanding of cellular processes, and interactive exploration of these complementary data sets can be extremely valuable for developing or confirming hypotheses in silico. However, few computational tools facilitate linking these types of data interactively. In addition, the tools that do exist are neither well understood nor widely used by the proteomics or structural biology communities. We briefly describe several relevant tools, and then, using three scenarios, we present in depth two tools for the integrated exploration of proteomics and structural data.

clusterMaker: a multi-algorithm clustering plugin for Cytoscape.

BACKGROUND: In the post-genomic era, the rapid increase in high-throughput data calls for computational tools capable of integrating data of diverse types and facilitating recognition of biologically meaningful patterns within them. For example, protein-protein interaction data sets have been clustered to identify stable complexes, but scientists lack easily accessible tools to facilitate combined analyses of multiple data sets from different types of experiments. Here we present clusterMaker, a Cytoscape plugin that implements several clustering algorithms and provides network, dendrogram, and heat map views of the results. The Cytoscape network is linked to all of the other views, so that a selection in one is immediately reflected in the others. clusterMaker is the first Cytoscape plugin to implement such a wide variety of clustering algorithms and visualizations, including the only implementations of hierarchical clustering, dendrogram plus heat map visualization (tree view), k-means, k-medoid, SCPS, AutoSOME, and native (Java) MCL. RESULTS: Results are presented in the form of three scenarios of use: analysis of protein expression data using a recently published mouse interactome and a mouse microarray data set of nearly one hundred diverse cell/tissue types; the identification of protein complexes in the yeast Saccharomyces cerevisiae; and the cluster analysis of the vicinal oxygen chelate (VOC) enzyme superfamily. For scenario one, we explore functionally enriched mouse interactomes specific to particular cellular phenotypes and apply fuzzy clustering. For scenario two, we explore the prefoldin complex in detail using both physical and genetic interaction clusters. For scenario three, we explore the possible annotation of a protein as a methylmalonyl-CoA epimerase within the VOC superfamily. Cytoscape session files for all three scenarios are provided in the Additional Files section. CONCLUSIONS: The Cytoscape plugin clusterMaker provides a number of clustering algorithms and visualizations that can be used independently or in combination for analysis and visualization of biological data sets, and for confirming or generating hypotheses about biological function. Several of these visualizations and algorithms are only available to Cytoscape users through the clusterMaker plugin. clusterMaker is available via the Cytoscape plugin manager.

UCSF ChimeraX: Structure visualization for researchers, educators, and developers.

UCSF ChimeraX is the next-generation interactive visualization program from the Resource for Biocomputing, Visualization, and Informatics (RBVI), following UCSF Chimera. ChimeraX brings (a) significant performance and graphics enhancements; (b) new implementations of Chimera's most highly used tools, many with further improvements; (c) several entirely new analysis features; (d) support for new areas such as virtual reality, light-sheet microscopy, and medical imaging data; (e) major ease-of-use advances, including toolbars with icons to perform actions with a single click, basic "undo" capabilities, and more logical and consistent commands; and (f) an app store for researchers to contribute new tools. ChimeraX includes full user documentation and is free for noncommercial use, with downloads available for Windows, Linux, and macOS from https://www.rbvi.ucsf.edu/chimerax.

Genetic variants of human organic anion transporter 4 demonstrate altered transport of endogenous substrates.

Apical reabsorption from the urine has been shown to be important for such processes as the maintenance of critical metabolites in the blood and the excretion of nephrotoxic compounds. The solute carrier (SLC) transporter OAT4 (SLC22A11) is expressed on the apical membrane of renal proximal tubule cells and is known to mediate the transport of a variety of xenobiotic and endogenous organic anions. Functional characterization of genetic variants of apical transporters thought to mediate reabsorption, such as OAT4, may provide insight into the genetic factors influencing the complex pathways involved in drug elimination and metabolite reclamation occurring in the kidney. Naturally occurring genetic variants of OAT4 were identified in public databases and by resequencing DNA samples from 272 individuals comprising 4 distinct ethnic groups. Nine total nonsynonymous variants were identified and functionally assessed using uptake of three radiolabeled substrates. A nonsense variant, R48Stop, and three other variants (R121C, V155G, and V155M) were found at frequencies of at least 2% in an ethnic group specific fashion. The L29P, R48Stop, and H469R variants displayed a complete loss of function, and kinetic analysis identified a reduced V(max) in the common nonsynonymous variants. Plasma membrane levels of OAT4 protein were absent or reduced in the nonfunctional variants, providing a mechanistic reason for the observed loss of function. Characterization of the genetic variants of reabsorptive transporters such as OAT4 is an important step in understanding variability in tubular reabsorption with important implications in innate homeostatic processes and drug disposition.

Role of organic cation transporter 3 (SLC22A3) and its missense variants in the pharmacologic action of metformin.

OBJECTIVES: The goals of this study were to determine the role of organic cation transporter 3 (OCT3) in the pharmacological action of metformin and to identify and functionally characterize genetic variants of OCT3 with respect to the uptake of metformin and monoamines. METHODS: For pharmacological studies, we evaluated metformin-induced activation of AMP-activated protein kinase, a molecular target of metformin. We used quantitative PCR and immunostaining to localize the transporter and isotopic uptake studies in cells transfected with OCT3 and its nonsynonymous genetic variants for functional analyses. RESULTS: Quantitative PCR and immunostaining showed that OCT3 was expressed high on the plasma membrane of skeletal muscle and liver, target tissues for metformin action. Both the OCT inhibitor, cimetidine, and OCT3-specific short hairpin RNA significantly reduced the activating effect of metformin on AMP-activated protein kinase. To identify genetic variants in OCT3, we used recent data from the 1000 Genomes and the Pharmacogenomics of Membrane Transporters projects. Six novel missense variants were identified. In functional assays, using various monoamines and metformin, three variants, T44M (c.131C>T), T400I (c.1199C>T) and V423F (c.1267G>T) showed altered substrate specificity. Notably, in cells expressing T400I and V423F, the uptakes of metformin and catecholamines were significantly reduced, but the uptakes of metformin, 1-methyl-4-phenylpyridinium and histamine by T44M were significantly increased more than 50%. Structural modeling suggested that these two variants may be located in the pore lining (T400) or proximal (V423) membrane-spanning helixes. CONCLUSION: Our study suggests that OCT3 plays a role in the therapeutic action of metformin and that genetic variants of OCT3 may modulate metformin and catecholamine action.

Visualization software for molecular assemblies.

Software for viewing three-dimensional models and maps of viruses, ribosomes, filaments, and other molecular assemblies is advancing on many fronts. New developments include molecular representations that offer better control over level of detail, lighting that improves the perception of depth, and two-dimensional projections that simplify data interpretation. Programmable graphics processors offer quality, speed, and visual effects not previously possible, while 3D printers, haptic interaction devices, and auto-stereo displays show promise in more naturally engaging our senses. Visualization methods are developed by diverse groups of researchers with differing goals: experimental biologists, database developers, computer scientists, and package developers. We survey recent developments and problems faced by the developer community in bringing innovative visualization methods into widespread use.

Functional genetic variation in the basal promoter of the organic cation/carnitine transporters OCTN1 (SLC22A4) and OCTN2 (SLC22A5).

The organic cation/ergothioneine transporter OCTN1 (SLC22A4) and the high-affinity carnitine transporter OCTN2 (SLC22A5), play an important role in the disposition of xenobiotics and endogenous compounds. Here, we analyzed the sequence of the proximal promoter regions of OCTN1 and OCTN2 in four ethnic groups and determined the effects of the identified genetic variants on transcriptional activities and mRNA expression. Six variants were found in the proximal promoter of OCTN1, one of which showed high allele frequency ranging from 13 to 34% in samples from individuals with ancestries in Africa, Europe, China, and Mexico. OCTN1 haplotypes had similar activities as the reference in luciferase reporter assays. For OCTN2, three of the seven variants identified in the proximal promoter showed allele frequencies greater than 29.5% in all populations, with the exception of -207C>G (rs2631367) that was monomorphic in Asian Americans. OCTN2 haplotypes containing -207G, present in all populations, were associated with a gain of function in luciferase reporter assays. Consistent with reporter assays, OCTN2 mRNA expression levels in lymphoblastoid cell lines (LCLs) from gene expression analysis were greater in samples carrying a marker for -207G. This SNP seems to contribute to racial differences in OCTN2 mRNA expression levels in LCLs. Our study with healthy subjects (n = 16) homozygous for either -207C or -207G, showed no appreciable effect of this SNP on carnitine disposition. However, there were significant effects of gender on carnitine plasma levels (p < 0.01). Further in vivo studies of OCTN2 promoter variants on carnitine disposition and variation in drug response are warranted.