PartsList: a web-based system for dynamically ranking protein folds based on disparate attributes, including whole-genome expression and interaction information

As the number of protein folds is quite limited, a mode of analysis that will be increasingly common in the future, especially with the advent of structural genomics, is to survey and re-survey the finite parts list of folds from an expanding number of perspectives. We have developed a new resource, called PartsList, that lets one dynamically perform these comparative fold surveys. It is available on the web at http://bioinfo.mbb.yale.edu/partslist and http://www.partslist.org. The system is based on the existing fold classifications and functions as a form of companion annotation for them, providing ‘global views’ of many already completed fold surveys. The central idea in the system is that of comparison through ranking; PartsList will rank the approximately 420 folds based on more than 180 attributes. These include: (i) occurrence in a number of completely sequenced genomes (e.g. it will show the most common folds in the worm versus yeast); (ii) occurrence in the structure databank (e.g. most common folds in the PDB); (iii) both absolute and relative gene expression information (e.g. most changing folds in expression over the cell cycle); (iv) protein–protein interactions, based on experimental data in yeast and comprehensive PDB surveys (e.g. most interacting fold); (v) sensitivity to inserted transposons; (vi) the number of functions associated with the fold (e.g. most multi-functional folds); (vii) amino acid composition (e.g. most Cys-rich folds); (viii) protein motions (e.g. most mobile folds); and (ix) the level of similarity based on a comprehensive set of structural alignments (e.g. most structurally variable folds). The integration of whole-genome expression and protein–protein interaction data with structural information is a particularly novel feature of our system. We provide three ways of visualizing the rankings: a profiler emphasizing the progression of high and low ranks across many pre-selected attributes, a dynamic comparer for custom comparisons and a numerical rankings correlator. These allow one to directly compare very different attributes of a fold (e.g. expression level, genome occurrence and maximum motion) in the uniform numerical format of ranks. This uniform framework, in turn, highlights the way that the frequency of many of the attributes falls off with approximate power-law behavior (i.e. according to V–b, for attribute value V and constant exponent b), with a few folds having large values and most having small values.

Simultaneous determination of helical unwinding angles and intrinsic association constants in ligand–DNA complexes: The interaction between DNA and calichearubicin B

We present a helical unwinding assay for reversibly binding DNA ligands that uses closed circular DNA, topoisomerase I (Topo I), and two-dimensional agarose gel electrophoresis. Serially diluted Topo I relaxation reactions at constant DNA/ligand ratio are performed, and the resulting apparent unwinding of the closed circular DNA is used to calculate both ligand unwinding angle (φ) and intrinsic association constant (Ka). Mathematical treatment of apparent unwinding is formally analogous to that of apparent extinction coefficient data for optical binding titrations. Extrapolation to infinite DNA concentration yields the true unwinding angle of a given ligand and its association constant under Topo I relaxation conditions. Thus this assay delivers simultaneous structural and thermodynamic information describing the ligand–DNA complex. The utility of this assay has been demonstrated by using calichearubicin B (CRB), a synthetic hybrid molecule containing the anthraquinone chromophore of (DA) and the carbohydrate domain of calicheamicin γ1I. The unwinding angle for CRB calculated by this method is −5.3 ± 0.5°. Its Ka value is 0.20 × 106 M−1. For comparison, the unwinding angles of ethidium bromide and DA have been independently calculated, and the results are in agreement with canonical values for these compounds. Although a stronger binder to selected sites, CRB is a less potent unwinder than its parent compound DA. The assay requires only small amounts of ligand and offers an attractive option for analysis of DNA binding by synthetic and natural compounds.

Neuronal type information encoded in the basic-helix–loop–helix domain of proneural genes

The proneural genes encode basic-helix–loop–helix (bHLH) proteins and promote the formation of distinct types of sensory organs. In Drosophila, two sets of proneural genes, atonal (ato) and members of the achaete–scute complex (ASC), are required for the formation of chordotonal (ch) organs and external sensory (es) organs, respectively. We assayed the production of sensory organs in transgenic flies expressing chimeric genes of ato and scute (sc), a member of ASC, and found that the information that specifies ch organs resides in the bHLH domain of ato; chimeras containing the b domain of ato and the HLH domain of sc also induced ch organ formation, but to a lesser extent than those containing the bHLH domain of ato. The b domains of ato and sc differ in seven residues. Mutations of these seven residues in the b domain of ato suggest that most or perhaps all of these residues are required for induction of ch organs. None of these seven residues is predicted to contact DNA directly by computer simulation using the structure of the myogenic factor MyoD as a model, implying that interaction of ato with other cofactors is likely to be involved in neuronal type specification.

Structural details of an interaction between cardiolipin and an integral membrane protein

Anionic lipids play a variety of key roles in biomembrane function, including providing the immediate environment for the integral membrane proteins that catalyze photosynthetic and respiratory energy transduction. Little is known about the molecular basis of these lipid–protein interactions. In this study, x-ray crystallography has been used to examine the structural details of an interaction between cardiolipin and the photoreaction center, a key light-driven electron transfer protein complex found in the cytoplasmic membrane of photosynthetic bacteria. X-ray diffraction data collected over the resolution range 30.0–2.1 Å show that binding of the lipid to the protein involves a combination of ionic interactions between the protein and the lipid headgroup and van der Waals interactions between the lipid tails and the electroneutral intramembrane surface of the protein. In the headgroup region, ionic interactions involve polar groups of a number of residues, the protein backbone, and bound water molecules. The lipid tails sit along largely hydrophobic grooves in the irregular surface of the protein. In addition to providing new information on the immediate lipid environment of a key integral membrane protein, this study provides the first, to our knowledge, high-resolution x-ray crystal structure for cardiolipin. The possible significance of this interaction between an integral membrane protein and cardiolipin is considered.

Molecular Interaction Map of the Mammalian Cell Cycle Control and DNA Repair Systems

Eventually to understand the integrated function of the cell cycle regulatory network, we must organize the known interactions in the form of a diagram, map, and/or database. A diagram convention was designed capable of unambiguous representation of networks containing multiprotein complexes, protein modifications, and enzymes that are substrates of other enzymes. To facilitate linkage to a database, each molecular species is symbolically represented only once in each diagram. Molecular species can be located on the map by means of indexed grid coordinates. Each interaction is referenced to an annotation list where pertinent information and references can be found. Parts of the network are grouped into functional subsystems. The map shows how multiprotein complexes could assemble and function at gene promoter sites and at sites of DNA damage. It also portrays the richness of connections between the p53-Mdm2 subsystem and other parts of the network.

BIND—The Biomolecular Interaction Network Database

The Biomolecular Interaction Network Database (BIND; http://binddb.org) is a database designed to store full descriptions of interactions, molecular complexes and pathways. Development of the BIND 2.0 data model has led to the incorporation of virtually all components of molecular mechanisms including interactions between any two molecules composed of proteins, nucleic acids and small molecules. Chemical reactions, photochemical activation and conformational changes can also be described. Everything from small molecule biochemistry to signal transduction is abstracted in such a way that graph theory methods may be applied for data mining. The database can be used to study networks of interactions, to map pathways across taxonomic branches and to generate information for kinetic simulations. BIND anticipates the coming large influx of interaction information from high-throughput proteomics efforts including detailed information about post-translational modifications from mass spectrometry. Version 2.0 of the BIND data model is discussed as well as implementation, content and the open nature of the BIND project. The BIND data specification is available as ASN.1 and XML DTD.

A dynamic assembly of diverse transcription factors integrates activation and cell-type information for interleukin 2 gene regulation.

The interleukin 2 (IL-2) gene is subject to two types of regulation: its expression is T-lymphocyte-specific and it is acutely dependent on specific activation signals. The IL-2 transcriptional apparatus integrates multiple types of biochemical information in determining whether or not the gene will be expressed, using multiple diverse transcription factors that are each optimally activated or inhibited by different signaling pathways. When activation of one or two of these factors is blocked IL-2 expression is completely inhibited. The inability of the other, unaffected factors to work is explained by the striking finding that none of the factors interacts stably with its target site in the IL-2 enhancer unless all the factors are present. Coordinate occupancy of all the sites in the minimal enhancer is apparently maintained by continuous assembly and disassembly cycles that respond to the instantaneous levels of each factor in the nuclear compartment. In addition, the minimal enhancer undergoes specific increases in DNase I accessibility, consistent with dramatic changes in chromatin structure upon activation. Still to be resolved is what interaction(s) conveys T-lineage specificity. In the absence of activating signals, the minimal IL-2 enhancer region in mature T cells is apparently unoccupied, exactly as in non-T lineage cells. However, in a conserved but poorly studied upstream region, we have now mapped several novel sites of DNase I hypersensitivity in vivo that constitutively distinguish IL-2 producer type T cells from cell types that cannot express IL-2. Thus a distinct domain of the IL-2 regulatory sequence may contain sites for competence- or lineage-marking protein contacts.

Earthquake prediction: the interaction of public policy and science.

Earthquake prediction research has searched for both informational phenomena, those that provide information about earthquake hazards useful to the public, and causal phenomena, causally related to the physical processes governing failure on a fault, to improve our understanding of those processes. Neither informational nor causal phenomena are a subset of the other. I propose a classification of potential earthquake predictors of informational, causal, and predictive phenomena, where predictors are causal phenomena that provide more accurate assessments of the earthquake hazard than can be gotten from assuming a random distribution. Achieving higher, more accurate probabilities than a random distribution requires much more information about the precursor than just that it is causally related to the earthquake.

The dimer-dimer interaction surface of the replication terminator protein of Bacillus subtilis and termination of DNA replication.

The replication terminator protein (RTP) of Bacillus subtilis causes polar fork arrest at replication termini by sequence-specific interaction of two dimeric proteins with the terminus sequence. The crystal structure of the RTP protein has been solved, and the structure has already provide valuable clues regarding the structural basis of its function. However, it provides little information as to the surface of the protein involved in dimer-dimer interaction. Using site-directed mutagenesis, we have identified three sites on the protein that appear to mediate the dimer-dimer interaction. Crystallographic analysis of one of the mutant proteins (Y88F) showed that its structure is unaltered when compared to the wild-type protein. The locations of the three sites suggested a model for the dimer-dimer interaction that involves an association between two beta-ribbon motifs. This model is supported by a fourth mutation that was predicted to disrupt the interaction and was shown to do so. Biochemical analyses of these mutants provide compelling evidence that cooperative protein-protein interaction between two dimers of RTP is essential to impose polar blocks to the elongation of both DNA and RNA chains.

Protein-protein interaction: a genetic selection for compensating mutations at the barnase-barstar interface.

Barnase and barstar are trivial names of the extracellular RNase and its intracellular inhibitor produced by Bacillus amyloliquefaciens. Inhibition involves the formation of a very tight one-to-one complex of the two proteins. With the crystallographic solution of the structure of the barnase-barstar complex and the development of methods for measuring the free energy of binding, the pair can be used to study protein-protein recognition in detail. In this report, we describe the isolation of suppressor mutations in barstar that compensate for the loss in interaction energy caused by a mutation in barnase. Our suppressor search is based on in vivo selection for barstar variants that are able to protect host cells against the RNAse activity of those barnase mutants not properly inhibited by wild-type barstar. This approach utilizes a plasmid system in which barnase expression is tightly controlled to keep the mutant barnase gene silent. When expression of barnase is turned on, failure to form a complex between the mutant barnase and barstar has a lethal effect on host cells unless overcome by substitution of the wild-type barstar by a functional suppressor derivative. A set of barstar suppressors has been identified for barnase mutants with substitutions in two amino acid positions (residues 102 and 59), which are critically involved in both RNase activity and barstar binding. The mutations selected as suppressors could not have been predicted on the basis of the known protein structures. The single barstar mutation with the highest information content for inhibition of barnase (H102K) has the substitution Y30W. The reduction in binding caused by the R59E mutation in barnase can be partly reversed by changing Glu-76 of barstar, which forms a salt bridge with the Arg-59 in the wild-type complex, to arginine, thus completing an interchange of the two charges.

The crystal structure of an N-terminal two-domain fragment of vascular cell adhesion molecule 1 (VCAM-1): a cyclic peptide based on the domain 1 C-D loop can inhibit VCAM-1-alpha 4 integrin interaction.

Vascular cell adhesion molecule 1 (VCAM-1) represents a structurally and functionally distinct class of immunoglobulin superfamily molecules that bind leukocyte integrins and are involved in inflammatory and immune functions. X-ray crystallography defines the three-dimensional structure of the N-terminal two-domain fragment that participates in ligand binding. Residues in domain 1 important for ligand binding reside in the C-D loop, which projects markedly from one face of the molecule near the contact between domains 1 and 2. A cyclic peptide that mimics this loop inhibits binding of alpha 4 beta 1 integrin-bearing cells to VCAM-1. These data demonstrate how crystallographic structural information can be used to design a small molecule inhibitor of biological function.