Tools for integrated sequence-structure analysis with UCSF Chimera

Background  

Comparing 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 researcher's 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.     

structureViz: linking Cytoscape and UCSF Chimera

structureViz is a Cytoscape plug-in that links the visualization of biological networks provided by Cytoscape with the visualization and analysis of macromolecular structures and sequences provided by UCSF Chimera. When combined with Cytoscape and Chimera, structureViz provides the first tool that links these two critical aspects of computational analysis in a straightforward manner. structureViz includes commands to open structures in Chimera and align them using Chimera's sequence-structure analysis tools. When a structure is opened, structureViz provides an alternative interface to Chimera: the Cytoscape Molecular Structure Navigator. This interface uses a tree-based paradigm to allow users to select and affect the display of models, chains and residues, mostly through the use of context menus.     

UCSF tomography: an integrated software suite for real-time electron microscopic tomographic data collection, alignment, and reconstruction

A real-time alignment and reconstruction scheme for electron microscopic tomography (EMT) has been developed and integrated within our UCSF tomography data collection software. This newly integrated software suite provides full automation from data collection to real-time reconstruction by which the three-dimensional (3D) reconstructed volume is immediately made available at the end of each data collection. Real-time reconstruction is achieved by calculating a weighted back-projection on a small Linux cluster (five dual-processor compute nodes) concurrently with the UCSF tomography data collection running on the microscope's computer, and using the fiducial-marker free alignment data generated during the data collection process. The real-time reconstructed 3D volume provides users with immediate feedback to fully asses all aspects of the experiment ranging from sample choice, ice thickness, experimental parameters to the quality of specimen preparation. This information can be used to guide subsequent data collections. Access to the reconstruction is especially useful in low-dose cryo EMT where such information is very difficult to obtain due to extraordinary low signal to noise ratio in each 2D image. In our environment, we generally collect 2048 x 2048 pixel images which are subsequently computationally binned four-fold for the on-line reconstruction. Based upon experiments performed with thick and cryo specimens at various CCD magnifications (50000x-80000x), alignment accuracy is sufficient to support this reduced resolution but should be refined before calculating a full resolution reconstruction. The reduced resolution has proven to be quite adequate to assess sample quality, or to screen for the best data set for full-resolution reconstruction, significantly improving both productivity and efficiency of system resources. The total time from start of data collection to a final reconstructed volume (512 x 512 x 256 pixels) is about 50 min for a +/-70 degrees 2k x 2k pixel tilt series acquired at every 1 degrees.     

ModEco: an integrated software package for ecological niche modeling

ModEco is a software package for ecological niche modeling. It integrates a range of niche modeling methods within a geographical information system. ModEco provides a user friendly platform that enables users to explore, analyze, and model species distribution data with relative ease. ModEco has several unique features: 1) it deals with different types of ecological observation data, such as presence and absence data, presence‐only data, and abundance data; 2) it provides a range of models when dealing with presence‐only data, such as presence‐only models, pseudo‐absence models, background vs presence data models, and ensemble models; and 3) it includes relatively comprehensive tools for data visualization, feature selection, and accuracy assessment.     

BioSilico: an integrated metabolic database system

BioSilico is a web-based database system that facilitates the search and analysis of metabolic pathways. Heterogeneous metabolic databases including LIGAND, ENZYME, EcoCyc and MetaCyc are integrated in a systematic way, thereby allowing users to efficiently retrieve the relevant information on enzymes, biochemical compounds and reactions. In addition, it provides well-designed view pages for more detailed summary information. BioSilico is developed as an extensible system with a robust systematic architecture.     

IMP dehydrogenase: structural schizophrenia and an unusual base

Textbooks describe enzymes as relatively rigid templates for the transition state of a chemical reaction, and indeed an enzyme such as chymotrypsin, which catalyzes a relatively simple hydrolysis reaction, is reasonably well described by this model. Inosine monophosphate dehydrogenase (IMPDH) undergoes a remarkable array of conformational transitions in the course of a complicated catalytic cycle, offering a dramatic counterexample to this view. IMPDH displays several other unusual mechanistic features, including an Arg residue that may act as a general base catalyst and a dynamic monovalent cation site. Further, IMPDH appears to be involved in 'moon-lighting' functions that may require additional conformational states. How the balance between conformational states is maintained and how the various conformational states interconvert is only beginning to be understood.     

Pathways database system: an integrated system for biological pathways

MOTIVATION: During the next phase of the Human Genome Project, research will focus on functional studies of attributing functions to genes, their regulatory elements, and other DNA sequences. To facilitate the use of genomic information in such studies, a new modeling perspective is needed to examine and study genome sequences in the context of many kinds of biological information. Pathways are the logical format for modeling and presenting such information in a manner that is familiar to biological researchers. RESULTS: In this paper we present an integrated system, called Pathways Database System, with a set of software tools for modeling, storing, analyzing, visualizing, and querying biological pathways data at different levels of genetic, molecular, biochemical and organismal detail. The novel features of the system include: (a) genomic information integrated with other biological data and presented from a pathway, rather than from the DNA sequence, perspective; (b) design for biologists who are possibly unfamiliar with genomics, but whose research is essential for annotating gene and genome sequences with biological functions; (c) database design, implementation and graphical tools which enable users to visualize pathways data in multiple abstraction levels, and to pose predetermined queries; and (d) an implementation that allows for web(XML)-based dissemination of query outputs (i.e. pathways data) to researchers in the community, giving them control on the use of pathways data. AVAILABILITY: Available on request from the authors.     

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 .     

Maintaining and restoring hydrologic habitat connectivity in mediterranean streams: an integrated modeling framework

Hydrologic alterations designed to provide a stable water supply and to prevent flooding are commonly used in mediterranean-climate river (med-rivers) basins, and these alterations have led to habitat loss and significant declines in aquatic biodiversity. Often the health of freshwater ecosystems depends on maintaining and recovering hydrologic habitat connectivity, which includes structural components related to the physical landscape, functionality of flow dynamics, and an understanding of species habitat requirements for movement, reproduction, and survival. To advance our understanding of hydrologic habitat connectivity and benefits of habitat restoration alternatives we provide: (1) a review of recent perspectives on hydrologic connectivity, including quantitative methods; and (2) a modeling framework to quantify the effects of restoration on hydrologic habitat connectivity. We then illustrate this approach through a case study on lateral hydrologic habitat connectivity that results from channel restoration scenarios using scenarios with different historic and climate-change flows to restore fish floodplain habitat in a med-river, the San Joaquin River, California. Case study results show that in addition to the channel alterations, higher flows are required to recover significant flooded habitat area, especially given reductions in flows expected under climate change. These types of studies will help the planning for restoration of hydrologic habitat connectivity in med-rivers, a critical step for mediterranean species recovery.     

High-content proteomics: fluorescence multiplexing using an integrated,high-sensitivity,multiwavelength charge-coupled device imaging system

The detection of proteins in 2-D gels and their subsequent identification by MS is still the "gold standard" in proteomics. Fluorescent detection has increasingly replaced colorimetric and radiometric detection on gels and blots. The reasons for this are multiple and varied and include higher sensitivity, better quantitation, increased dynamic range, speed, safety and ease of use. Unlike other methods, fluorescent protein detection is also typically very consistent in response from protein to protein and in many cases is compatible with MS methods for protein identification. The superior sensitivity and benefits achieved by fluorescent techniques have spurred the development of instrumentation capable of delivering precise, sensitive, high-resolution image acquisition over a wide variety of excitation and emission wavelengths. This report focuses on applications using the highly sensitive, charge-coupled device based ProXPRESS multilabel imager, readily configurable for image acquisition over a wide variety of wavelengths (380-700 nm and ultraviolet (UV)) using xenon lamp or UV excitation. The ability to simultaneously detect enzyme activities or protein modifications with different color fluorescent probes in addition to total protein amounts (multiplexing) allows the further mining of proteomic data content from a single set of protein samples. To this end, the development of instrumentation that enables a multiplexing strategy will become central to in-depth proteomic studies. The ProXPRESS maximizes the efficiency of experimental strategies that require flexibility and multicolor fluorescence detection.     

HOWDY: an integrated database system for human genome research

HOWDY is an integrated database system for accessing and analyzing human genomic information (http://www-alis.tokyo.jst.go.jp/HOWDY/). HOWDY stores information about relationships between genetic objects and the data extracted from a number of databases. HOWDY consists of an Internet accessible user interface that allows thorough searching of the human genomic databases using the gene symbols and their aliases. It also permits flexible editing of the sequence data. The database can be searched using simple words and the search can be restricted to a specific cytogenetic location. Linear maps displaying markers and genes on contig sequences are available, from which an object can be chosen. Any search starting point identifies all the information matching the query. HOWDY provides a convenient search environment of human genomic data for scientists unsure which database is most appropriate for their search.     

Optimal modeling of an arterial system

The aortic blood flow is described by a set of nonlinear hyperbolic partial differential equations that account for mass and momentum conservation, and nonlinear models for the mechanical properties of the artery. Identification is used for determining the wave speed, arterial taper, and cross section: these parameters reflect the elastic characteristics of the aorta wall and control the pulsatile response. The differential equations were numerically integrated by the Lax-Wendroff scheme of Abarbanel and Goldberg [J. Comput. Phys. 10:1–21 (1972)] that avoids nonlinear oscillations. The Gauss-Newton technique was used for the parameter identification. By reference to reported elocity and pressure input-output pairs, a parameter vector is found such that the distance in the L₂ norm between the predicted outputs and the measured functions is minimal. Calculations of the velocity and pressure waves show excellent compatibility of the model with reported experimental data: starting from arbitrary parameter estimates, which yield grossly distorted waveforms, the error is typically reduced to 7–8%. Introduction of viscoelastic behaviour for the arterial wall in the form of a Volterra integral for the cross section does not lead to significant improvement. Numerical examples are presented which prove the convergence, accuracy, and stability of the algorithm. Emphasis is placed on the computational feasibility of the proposed system identification.     

On an interacting particle system modeling an epidemic

We consider an interacting particle system onZ ^dto model an epidemic. Each site ofZ ^dcan be in either one of three states: empty, healthy or infected. Healthy and infected individuals give birth at different rates to healthy individuals on empty sites. Healthy individuals get infected by infected individuals. Infected and healthy individuals die at different rates. We prove that in dimension 1 and with nearest-neighbor interactions the epidemic may persist forever if and only if the rate at which infected individuals give birth to healthy individuals is high enough. This is in sharp contrast with models analysed by Andjel and Schinazi (1994) and Sato et al. (1994) where infected individuals do not give birth. We also show that some results in the latter reference can be obtained easily and rigorously using probabilistic coupling to the contact process.     

Toward an integrated system of clade names

Although the proposition that higher taxa should correspond to clades is widely accepted, current nomenclature does not distinguish clearly between different clades in nested series. In particular, the same name is often applied to a total clade, its crown clade, and clades originating with various nodes, branches, and apomorphies in between. An integrated system of clade names is described based on categories of clades defined with respect to lineages that have survived to the present time. In this system, the most widely known names are applied to crown clades, the names of total clades are formed by adding a standard prefix to the names of the corresponding crowns, and the names of apomorphy clades describe the specific apomorphies with which they originated. Relative to traditional approaches, this integrated approach to naming clades is both more precise concerning the associations of names with particular clades and more efficient with regard to the cognitive effort required to recognize the names of corresponding crown and total clades. It also seems preferable to five alternatives that could be used to make the same distinctions. The integrated system of clade names has several advantages, including the facilitation of communication among biologists who study distantly related clades, promoting a broader conceptualization of the origins of distinctive clades of extant organisms and emphasizing the continuous nature of evolution.     

GCSDB: an integrated database system for the Georgia Centenarian Study

GCSDB is a web-oriented integrated database system for the Georgia Centenarian Study, a phase III, population-based,multidisciplinary study of centenarians. The Study recruited 244 centenarians and near-centenarians (age 98 and older),80 octogenarians and 400 young controls in Northern Georgia. GCSDB incorporates more than 40 relational tablescontaining data about the participants including demographics, family longevity, physical health, cognition,neuropsychology, mental health, neuropathology, functional capacity, and genetics. The GCSDB web site includesdetailed information about these tables and functions for genetic and other kinds of data analysis. More data andfunctions will be added as the study progresses. GCSDB provides a resource that could be used to identify whatbiological, psychological, and social factors as well as their epistatic interactions help these centenarians achieve longlife.

Availability     

Computational modeling of human coreceptor CCR5 antagonist as a HIV-1 entry inhibitor: using an integrated homology modeling,docking, and membrane molecular dynamics simulation analysis approach

Chemokine receptor 5 (CCR5) is an integral membrane protein that is utilized during human immunodeficiency virus type-1 entry into host cells. CCR5 is a G-protein coupled receptor that contains seven transmembrane (TM) helices. However, the crystal structure of CCR5 has not been reported. A homology model of CCR5 was developed based on the recently reported CXCR4 structure as template. Automated docking of the most potent (14), medium potent (37), and least potent (25) CCR5 antagonists was performed using the CCR5 model. To characterize the mechanism responsible for the interactions between ligands (14, 25, and 37) and CCR5, membrane molecular dynamic (MD) simulations were performed. The position and orientation of ligands (14, 25, and 37) were found to be changed after MD simulations, which demonstrated the ability of this technique to identify binding modes. Furthermore, at the end of simulation, it was found that residues identified by docking were changed and some new residues were introduced in the proximity of ligands. Our results are in line with the majority of previous mutational reports. These results show that hydrophobicity is the determining factor of CCR5 antagonism. In addition, salt bridging and hydrogen bond contacts between ligands (14, 25, and 37) and CCR5 are also crucial for inhibitory activity. The residues newly identified by MD simulation are Ser160, Phe166, Ser180, His181, and Trp190, and so far no site-directed mutagenesis studies have been reported. To determine the contributions made by these residues, additional mutational studies are suggested. We propose a general binding mode for these derivatives based on the MD simulation results of higher (14), medium (37), and lower (25) potent inhibitors. Interestingly, we found some trend for these inhibitors such as, salt bridge interaction between basic nitrogen of ligand and acidic Glu283 seemed necessary for inhibitory activity. Also, two aromatic pockets (pocket I – TM1-3 and pocket II – TM3-6) were linked by the central polar region in TM7, and the simulated inhibitors show important interactions with the Trp86, Tyr89, Tyr108, Phe112, Ile198, Tyr251, Leu255, and Gln280 and Glu283 residues. These results shed light on the usage of MD simulation to identify more stable, optimal binding modes of the inhibitors.     

Columba: an integrated database of proteins,structures, and annotations

Background  

Structural and functional research often requires the computation of sets of protein structures based on certain properties of the proteins, such as sequence features, fold classification, or functional annotation. Compiling such sets using current web resources is tedious because the necessary data are spread over many different databases. To facilitate this task, we have created COLUMBA, an integrated database of annotations of protein structures.