A unified statistical framework for sequence comparison and structure comparison

We present an approach for assessing the significance of sequence and structure comparisons by using nearly identical statistical formalisms for both sequence and structure. Doing so involves an all-vs.-all comparison of protein domains [taken here from the Structural Classification of Proteins (scop) database] and then fitting a simple distribution function to the observed scores. By using this distribution, we can attach a statistical significance to each comparison score in the form of a P value, the probability that a better score would occur by chance. As expected, we find that the scores for sequence matching follow an extreme-value distribution. The agreement, moreover, between the P values that we derive from this distribution and those reported by standard programs (e.g., blast and fasta validates our approach. Structure comparison scores also follow an extreme-value distribution when the statistics are expressed in terms of a structural alignment score (essentially the sum of reciprocated distances between aligned atoms minus gap penalties). We find that the traditional metric of structural similarity, the rms deviation in atom positions after fitting aligned atoms, follows a different distribution of scores and does not perform as well as the structural alignment score. Comparison of the sequence and structure statistics for pairs of proteins known to be related distantly shows that structural comparison is able to detect approximately twice as many distant relationships as sequence comparison at the same error rate. The comparison also indicates that there are very few pairs with significant similarity in terms of sequence but not structure whereas many pairs have significant similarity in terms of structure but not sequence.

The structure of the two amino-terminal domains of human ICAM-1 suggests how it functions as a rhinovirus receptor and as an LFA-1 integrin ligand

The normal function of human intercellular adhesion molecule-1 (ICAM-1) is to provide adhesion between endothelial cells and leukocytes after injury or stress. ICAM-1 binds to leukocyte function-associated antigen (LFA-1) or macrophage-1 antigen (Mac-1). However, ICAM-1 is also used as a receptor by the major group of human rhinoviruses and is a catalyst for the subsequent viral uncoating during cell entry. The three-dimensional atomic structure of the two amino-terminal domains (D1 and D2) of ICAM-1 has been determined to 2.2-Å resolution and fitted into a cryoelectron microscopy reconstruction of a rhinovirus–ICAM-1 complex. Rhinovirus attachment is confined to the BC, CD, DE, and FG loops of the amino-terminal Ig-like domain (D1) at the end distal to the cellular membrane. The loops are considerably different in structure to those of human ICAM-2 or murine ICAM-1, which do not bind rhinoviruses. There are extensive charge interactions between ICAM-1 and human rhinoviruses, which are mostly conserved in both major and minor receptor groups of rhinoviruses. The interaction of ICAMs with LFA-1 is known to be mediated by a divalent cation bound to the insertion (I)-domain on the α chain of LFA-1 and the carboxyl group of a conserved glutamic acid residue on ICAMs. Domain D1 has been docked with the known structure of the I-domain. The resultant model is consistent with mutational data and provides a structural framework for the adhesion between these molecules.

The underlying pathway structure of biochemical reaction networks

Bioinformatics is yielding extensive, and in some cases complete, genetic and biochemical information about individual cell types and cellular processes, providing the composition of living cells and the molecular structure of its components. These components together perform integrated cellular functions that now need to be analyzed. In particular, the functional definition of biochemical pathways and their role in the context of the whole cell is lacking. In this study, we show how the mass balance constraints that govern the function of biochemical reaction networks lead to the translation of this problem into the realm of linear algebra. The functional capabilities of biochemical reaction networks, and thus the choices that cells can make, are reflected in the null space of their stoichiometric matrix. The null space is spanned by a finite number of basis vectors. We present an algorithm for the synthesis of a set of basis vectors for spanning the null space of the stoichiometric matrix, in which these basis vectors represent the underlying biochemical pathways that are fundamental to the corresponding biochemical reaction network. In other words, all possible flux distributions achievable by a defined set of biochemical reactions are represented by a linear combination of these basis pathways. These basis pathways thus represent the underlying pathway structure of the defined biochemical reaction network. This development is significant from a fundamental and conceptual standpoint because it yields a holistic definition of biochemical pathways in contrast to definitions that have arisen from the historical development of our knowledge about biochemical processes. Additionally, this new conceptual framework will be important in defining, characterizing, and studying biochemical pathways from the rapidly growing information on cellular function.

Crystal structure of a Staphylococcus aureus protein A domain complexed with the Fab fragment of a human IgM antibody: Structural basis for recognition of B-cell receptors and superantigen activity

Staphylococcus aureus produces a virulence factor, protein A (SpA), that contains five homologous Ig-binding domains. The interactions of SpA with the Fab region of membrane-anchored Igs can stimulate a large fraction of B cells, contributing to lymphocyte clonal selection. To understand the molecular basis for this activity, we have solved the crystal structure of the complex between domain D of SpA and the Fab fragment of a human IgM antibody to 2.7-Å resolution. In the complex, helices II and III of domain D interact with the variable region of the Fab heavy chain (VH) through framework residues, without the involvement of the hypervariable regions implicated in antigen recognition. The contact residues are highly conserved in human VH3 antibodies but not in other families. The contact residues from domain D also are conserved among all SpA Ig-binding domains, suggesting that each could bind in a similar manner. Features of this interaction parallel those reported for staphylococcal enterotoxins that are superantigens for many T cells. The structural homology between Ig VH regions and the T-cell receptor Vβ regions facilitates their comparison, and both types of interactions involve lymphocyte receptor surface remote from the antigen binding site. However, T-cell superantigens reportedly interact through hydrogen bonds with T-cell receptor Vβ backbone atoms in a primary sequence-independent manner, whereas SpA relies on a sequence-restricted conformational binding with residue side chains, suggesting that this common bacterial pathogen has adopted distinct molecular recognition strategies for affecting large sets of B and T lymphocytes.

Crystal structure of the breakpoint cluster region-homology domain from phosphoinositide 3-kinase p85α subunit

Proteins such as the product of the breakpoint cluster region, chimaerin, and the Src homology 3-binding protein 3BP1, are GTPase activating proteins (GAPs) for members of the Rho subfamily of small GTP-binding proteins (G proteins or GTPases). A 200-residue region, named the breakpoint cluster region-homology (BH) domain, is responsible for the GAP activity. We describe here the crystal structure of the BH domain from the p85 subunit of phosphatidylinositol 3-kinase at 2.0 Å resolution. The domain is composed of seven helices, having a previously unobserved arrangement. A core of four helices contains most residues that are conserved in the BH family. Their packing suggests the location of a G-protein binding site. This structure of a GAP-like domain for small GTP-binding proteins provides a framework for analyzing the function of this class of molecules.

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

Large vortex-like structure of dipole field in computer models of liquid water and dipole-bridge between biomolecules

We propose a framework to describe the cooperative orientational motions of water molecules in liquid water and around solute molecules in water solutions. From molecular dynamics (MD) simulation a new quantity “site-dipole field” is defined as the averaged orientation of water molecules that pass through each spatial position. In the site-dipole field of bulk water we found large vortex-like structures of more than 10 Å in size. Such coherent patterns persist more than 300 ps although the orientational memory of individual molecules is quickly lost. A 1-ns MD simulation of systems consisting of two amino acids shows that the fluctuations of site-dipole field of solvent are pinned around the amino acids, resulting in a stable dipole-bridge between side-chains of amino acids. The dipole-bridge is significantly formed even for the side-chain separation of 14 Å, which corresponds to five layers of water. The way that dipole-bridge forms sensitively depends on the side-chain orientations and thereby explains the specificity in the solvent-mediated interactions between biomolecules.

A framework for interpreting the leucine-rich repeats of the Listeria internalins

The surface protein InlB of the bacterial pathogen Listeria monocytogenes is required for inducing phagocytosis in various nonphagocytic mammalian cell types in vitro. InlB causes tyrosine phosphorylation of host cell adaptor proteins, activation of phosphoinositide 3-kinase, and rearrangements of the actin cytoskeleton. These events lead to phagocytic uptake of the bacterium by the host cell. InlB belongs to the internalin family of Listeria proteins, which also includes InlA, another surface protein involved in host cell invasion. The internalins are the largest class of bacterial proteins containing leucine-rich repeats (LRR), a motif associated with protein–protein interactions. The LRR motif is found in a functionally diverse array of proteins, including those involved in the plant immune system and in the mammalian innate immune response. Structural and functional interpretations of the sequences of internalin family members are presented in light of the recently determined x-ray crystal structure of the InlB LRR domain.

The solution structure of the Raf-1 cysteine-rich domain: a novel ras and phospholipid binding site.

The Raf-1 protein kinase is the best-characterized downstream effector of activated Ras. Interaction with Ras leads to Raf-1 activation and results in transduction of cell growth and differentiation signals. The details of Raf-1 activation are unclear, but our characterization of a second Ras-binding site in the cysteine-rich domain (CRD) and the involvement of both Ras-binding sites in effective Raf-1-mediated transformation provides insight into the molecular aspects and consequences of Ras-Raf interactions. The Raf-1 CRD is a member of an emerging family of domains, many of which are found within signal transducing proteins. Several contain binding sites for diacylglycerol (or phorbol esters) and phosphatidylserine and are believed to play a role in membrane translocation and enzyme activation. The CRD from Raf-1 does not bind diacylglycerol but interacts with Ras and phosphatidylserine. To investigate the ligand-binding specificities associated with CRDs, we have determined the solution structure of the Raf-1 CRD using heteronuclear multidimensional NMR. We show that there are differences between this structure and the structures of two related domains from protein kinase C (PKC). The differences are confined to regions of the CRDs involved in binding phorbol ester in the PKC domains. Since phosphatidylserine is a common ligand, we expect its binding site to be located in regions where the structures of the Raf-1 and PKC domains are similar. The structure of the Raf-1 CRD represents an example of this family of domains that does not bind diacylglycerol and provides a framework for investigating its interactions with other molecules.

The crystal structure of the immunity protein of colicin E7 suggests a possible colicin-interacting surface.

The immunity protein of colicin E7 (ImmE7) can bind specifically to the DNase-type colicin E7 and inhibit its bactericidal activity. Here we report the 1.8-angstrom crystal structure of the ImmE7 protein. This is the first x-ray structure determined in the superfamily of colicin immunity proteins. The ImmE7 protein consists of four antiparallel alpha-helices, folded in a topology similar to the architecture of a four-helix bundle structure. A region rich in acidic residues is identified. This negatively charged area has the greatest variability within the family of DNase-type immunity proteins; thus, it seems likely that this area is involved in specific binding to colicin. Based on structural, genetic, and kinetic data, we suggest that all the DNase-type immunity proteins, as well as colicins, share a "homologous-structural framework" and that specific interaction between a colicin and its cognate immunity protein relies upon how well these two proteins' charged residues match on the interaction surface, thus leading to specific immunity of the colicin.

Tertiary structure of an amyloid immunoglobulin light chain protein: a proposed model for amyloid fibril formation.

An immunoglobulin light chain protein was isolated from the urine of an individual (BRE) with systemic amyloidosis. Complete amino acid sequence of the variable region of the light chain (VL) protein established it as a kappa I, which when compared with other kappa I amyloid associated proteins had unique residues, including Ile-34, Leu-40, and Tyr-71. To study the tertiary structure, BRE VL was expressed in Escherichia coli by using a PCR product amplified from the patient BRE's bone marrow DNA. The PCR product was ligated into pCZ11, a thermal-inducible replication vector. Recombinant BRE VL was isolated, purified to homogeneity, and crystallized by using ammonium sulfate as the precipitant. Two crystal forms were obtained. In crystal form I the BRE VL kappa domain crystallizes as a dimer with unit cell constants isomorphous to previously published kappa protein structures. Comparison with a nonamyloid VL kappa domain from patient REI, identified significant differences in position of residues in the hypervariable segments plus variations in framework region (FR) segments 40-46 (FR2) and 66-67 (FR3). In addition, positional differences can be seen along the two types of local diads, corresponding to the monomer-monomer and dimer-dimer interfaces. From the packing diagram, a model for the amyloid light chain (AL) fibril is proposed based on a pseudohexagonal spiral structure with a rise of approximately the width of two dimers per 360 degree turn. This spiral structure could be consistent with the dimensions of amyloid fibrils as determined by electron microscopy.

Structure of human T-cell receptors specific for an immunodominant myelin basic protein peptide: positioning of T-cell receptors on HLA-DR2/peptide complexes.

T-cell receptors (TCRs) recognize peptide bound within the relatively conserved structural framework of major histocompatibility complex (MHC) class I or class II molecules but can discriminate between closely related MHC molecules. The structural basis for the specificity of ternary complex formation by the TCR and MHC/peptide complexes was examined for myelin basic protein (MBP)-specific T-cell clones restricted by different DR2 subtypes. Conserved features of this system allowed a model for positioning of the TCR on DR2/peptide complexes to be developed: (i) The DR2 subtypes that presented the immunodominant MBP peptide differed only at a few polymorphic positions of the DR beta chain. (ii) TCR recognition of a polymorphic residue on the helical portion of the DR beta chain (position DR beta 67) was important in determining the MHC restriction. (iii) The TCR variable region (V) alpha 3.1 gene segment was used by all of the T-cell clones. TCR V beta usage was more diverse but correlated with the MHC restriction--i.e., with the polymorphic DR beta chains. (iv) Two clones with conserved TCR alpha chains but different TCR beta chains had a different MHC restriction but a similar peptide specificity. The difference in MHC restriction between these T-cell clones appeared due to recognition of a cluster of polymorphic DR beta-chain residues (DR beta 67-71). MBP-(85-99)-specific TCRs therefore appeared to be positioned on the DR2/peptide complex such that the TCR beta chain contacted the polymorphic DR beta-chain helix while the conserved TCR alpha chain contacted the nonpolymorphic DR alpha chain.