Estimating the probability for a protein to have a new fold: A statistical computational model

Structural genomics aims to solve a large number of protein structures that represent the protein space. Currently an exhaustive solution for all structures seems prohibitively expensive, so the challenge is to define a relatively small set of proteins with new, currently unknown folds. This paper presents a method that assigns each protein with a probability of having an unsolved fold. The method makes extensive use of protomap, a sequence-based classification, and scop, a structure-based classification. According to protomap, the protein space encodes the relationship among proteins as a graph whose vertices correspond to 13,354 clusters of proteins. A representative fold for a cluster with at least one solved protein is determined after superposition of all scop (release 1.37) folds onto protomap clusters. Distances within the protomap graph are computed from each representative fold to the neighboring folds. The distribution of these distances is used to create a statistical model for distances among those folds that are already known and those that have yet to be discovered. The distribution of distances for solved/unsolved proteins is significantly different. This difference makes it possible to use Bayes' rule to derive a statistical estimate that any protein has a yet undetermined fold. Proteins that score the highest probability to represent a new fold constitute the target list for structural determination. Our predicted probabilities for unsolved proteins correlate very well with the proportion of new folds among recently solved structures (new scop 1.39 records) that are disjoint from our original training set.

In vitro and in vivo differentiation into B cells, T cells, and myeloid cells of primitive yolk sac hematopoietic precursor cells expanded >100-fold by coculture with a clonal yolk sac endothelial cell line

The yolk sac, first site of hematopoiesis during mammalian development, contains not only hematopoietic stem cells but also the earliest precursors of endothelial cells. We have previously shown that a nonadherent yolk sac cell population (WGA+, density <1.077, AA4.1+) can give rise to B cells, T cells, and myeloid cells both in vitro and in vivo. We now report on the ability of a yolk sac-derived cloned endothelial cell line (C166) to provide a suitable microenvironment for expansion of these early precursor cells. Single day 10 embryonic mouse yolk sac hematopoietic stem cells were expanded >100 fold within 8 days by coculture with irradiated C166 cells. Colony-forming ability was retained for at least three passages in vitro, with retention of the ability to differentiate into T-cell, B-cell, and myeloid lineages. Stem cell properties were maintained by a significant fraction of nonadherent cells in the third passage, although these stem cells expressed a somewhat more mature cell surface phenotype than the initial yolk sac stem cells. When reintroduced into adult allogeneic immunocompromised (scid) hosts, they were able to give rise to all of the leukocyte lineages, including T cells, B cells, and myeloid cells. We conclude that yolk sac endothelial cells can support the stable proliferation of multipotential hematopoietic stem cells, thus generating adequate numbers of cells for study of the mechanisms involved in their subsequent development and differentiation, for in vivo hematopoietic restitution, and for potential use as a vehicle for gene transfer.

Small antibody-like proteins with prescribed ligand specificities derived from the lipocalin fold

We demonstrate that the ligand pocket of a lipocalin from Pieris brassicae, the bilin-binding protein (BBP), can be reshaped by combinatorial protein design such that it recognizes fluorescein, an established immunological hapten. For this purpose 16 residues at the center of the binding site, which is formed by four loops on top of an eight-stranded β-barrel, were subjected to random mutagenesis. Fluorescein-binding BBP variants were then selected from the mutant library by bacterial phage display. Three variants were identified that complex fluorescein with high affinity, exhibiting dissociation constants as low as 35.2 nM. Notably, one of these variants effects almost complete quenching of the ligand fluorescence, similarly as an anti-fluorescein antibody. Detailed ligand-binding studies and site-directed mutagenesis experiments indicated (i) that the molecular recognition of fluorescein is specific and (ii) that charged residues at the center of the pocket are responsible for tight complex formation. Sequence comparison of the BBP variants directed against fluorescein with the wild-type protein and with further variants that were selected against several other ligands revealed that all of the randomized amino acid positions are variable. Hence, a lipocalin can be used for generating molecular pockets with a diversity of shapes. We term this class of engineered proteins “anticalins.” Their one-domain scaffold makes them a promising alternative to antibodies to create a stable receptor protein for a ligand of choice.

How does a β-hairpin fold/unfold? Competition between topology and heterogeneity in a solvable model

We study the competition between topological effects and sequence inhomogeneities in determining the thermodynamics and the un/folding kinetics of a β-hairpin. Our work utilizes a new exactly solvable model that allows for arbitrary configurations of native contacts. In general, the competition between heterogeneity and topology results in a crossover of the dominant transition state. Interestingly, near this crossover, the single reaction coordinate picture can be seriously misleading. Our results also suggest that inferring the folding pathway from unfolding simulations is not always justified.

A knob-associated tandem repeat in maize capable of forming fold-back DNA segments: Are chromosome knobs megatransposons?

A class of tandemly repeated DNA sequences (TR-1) of 350-bp unit length was isolated from the knob DNA of chromosome 9 of Zea mays L. Comparative fluorescence in situ hybridization revealed that TR-1 elements are also present in cytologically detectable knobs on other maize chromosomes in different proportions relative to the previously described 180-bp repeats. At least one knob on chromosome 4 is composed predominantly of the TR-1 repeat. In addition, several small clusters of the TR-1 and 180-bp repeats have been found in different chromosomes, some not located in obvious knob heterochromatin. Variation in restriction fragment fingerprints and copy number of the TR-1 elements was found among maize lines and among maize chromosomes. TR-1 tandem arrays up to 70 kilobases in length can be interspersed with stretches of 180-bp tandem repeat arrays. DNA sequence analysis and restriction mapping of one particular stretch of tandemly arranged TR-1 units indicate that these elements may be organized in the form of fold-back DNA segments. The TR-1 repeat shares two short segments of homology with the 180-bp repeat. The longest of these segments (31 bp; 64% identity) corresponds to the conserved region among 180-bp repeats. The polymorphism and complex structure of knob DNA suggest that, similar to the fold-back DNA-containing giant transposons in Drosophila, maize knob DNA may have some properties of transposable elements.

Reverse engineering the (β/α)8 barrel fold

The (β/α)8 barrel is the most commonly occurring fold among protein catalysts. To lay a groundwork for engineering novel barrel proteins, we investigated the amino acid sequence restrictions at 182 structural positions of the prototypical (β/α)8 barrel enzyme triosephosphate isomerase. Using combinatorial mutagenesis and functional selection, we find that turn sequences, α-helix capping and stop motifs, and residues that pack the interface between β-strands and α-helices are highly mutable. Conversely, any mutation of residues in the central core of the β-barrel, β-strand stop motifs, and a single buried salt bridge between amino acids R189 and D227 substantially reduces catalytic activity. Four positions are effectively immutable: conservative single substitutions at these four positions prevent the mutant protein from complementing a triosephosphate isomerase knockout in Escherichia coli. At 142 of the 182 positions, mutation to at least one amino acid of a seven-letter amino acid alphabet produces a triosephosphate isomerase with wild-type activity. Consequently, it seems likely that (β/α)8 barrel structures can be encoded with a subset of the 20 amino acids. Such simplification would greatly decrease the computational burden of (β/α)8 barrel design.

α10: A determinant of nicotinic cholinergic receptor function in mammalian vestibular and cochlear mechanosensory hair cells

We report the cloning and characterization of rat α10, a previously unidentified member of the nicotinic acetylcholine receptor (nAChR) subunit gene family. The protein encoded by the α10 nAChR subunit gene is most similar to the rat α9 nAChR, and both α9 and α10 subunit genes are transcribed in adult rat mechanosensory hair cells. Injection of Xenopus laevis oocytes with α10 cRNA alone or in pairwise combinations with either α2-α6 or β2-β4 subunit cRNAs yielded no detectable ACh-gated currents. However, coinjection of α9 and α10 cRNAs resulted in the appearance of an unusual nAChR subtype. Compared with homomeric α9 channels, the α9α10 nAChR subtype displays faster and more extensive agonist-mediated desensitization, a distinct current–voltage relationship, and a biphasic response to changes in extracellular Ca2+ ions. The pharmacological profiles of homomeric α9 and heteromeric α9α10 nAChRs are essentially indistinguishable and closely resemble those reported for endogenous cholinergic eceptors found in vertebrate hair cells. Our data suggest that efferent modulation of hair cell function occurs, at least in part, through heteromeric nAChRs assembled from both α9 and α10 subunits.

Growth factor treatment enhances vestibular hair cell renewal and results in improved vestibular function

The vestibules of adult guinea pigs were lesioned with gentamicin and then treated with perilymphatic infusion of either of two growth factor mixtures (i.e., GF I or GF II). GF I contained transforming growth factor α (TGFα), insulin-like growth factor type one (IGF-1), and retinoic acid (RA), whereas GF II contained those three factors and brain-derived neurotrophic factor. Treatment with GF I significantly enhanced vestibular hair cell renewal in ototoxin-damaged utricles and the maturation of stereociliary bundle morphology. The addition of brain-derived neurotrophic factor to the GF II infusion mixture resulted in the return of type 1 vestibular hair cells in ototoxin-damaged cristae, and improved vestibular function. These results suggest that growth factor therapy may be an effective treatment for balance disorders that are the result of hair cell dysfunction and/or loss.

Reversibly locking a protein fold in an active conformation with a disulfide bond: Integrin αL I domains with high affinity and antagonist activity in vivo

The integrin αLβ2 has three different domains in its headpiece that have been suggested to either bind ligand or to regulate ligand binding. One of these, the inserted or I domain, has a fold similar to that of small G proteins. The I domain of the αM and α2 subunits has been crystallized in both open and closed conformations; however, the αL I domain has been crystallized in only the closed conformation. We hypothesized that the αL domain also would have an open conformation, and that this would be the ligand binding conformation. Therefore, we introduced pairs of cysteine residues to form disulfides that would lock the αL I domain in either the open or closed conformation. Locking the I domain open resulted in a 9,000-fold increase in affinity to intercellular adhesion molecule-1 (ICAM-1), which was reversed by disulfide reduction. By contrast, the affinity of the locked closed conformer was similar to wild type. Binding completely depended on Mg2+. Orders of affinity were ICAM-1 > ICAM-2 > ICAM-3. The kon, koff, and KD values for the locked open I domain were within 1.5-fold of values previously determined for the αLβ2 complex, showing that the I domain is sufficient for full affinity binding to ICAM-1. The locked open I domain antagonized αLβ2-dependent adhesion in vitro, lymphocyte homing in vivo, and firm adhesion but not rolling on high endothelial venules. The ability to reversibly lock a protein fold in an active conformation with dramatically increased affinity opens vistas in therapeutics and proteomics.

Five-fold symmetry in crystalline quasicrystal lattices

To demonstrate that crystallographic methods can be applied to index and interpret diffraction patterns from well-ordered quasicrystals that display non-crystallographic 5-fold symmetry, we have characterized the properties of a series of periodic two-dimensional lattices built from pentagons, called Fibonacci pentilings, which resemble aperiodic Penrose tilings. The computed diffraction patterns from periodic pentilings with moderate size unit cells show decagonal symmetry and are virtually indistinguishable from that of the infinite aperiodic pentiling. We identify the vertices and centers of the pentagons forming the pentiling with the positions of transition metal atoms projected on the plane perpendicular to the decagonal axis of quasicrystals whose structure is related to crystalline η phase alloys. The characteristic length scale of the pentiling lattices, evident from the Patterson (autocorrelation) function, is ∼τ2 times the pentagon edge length, where τ is the golden ratio. Within this distance there are a finite number of local atomic motifs whose structure can be crystallographically refined against the experimentally measured diffraction data.

Photoactive yellow protein: A structural prototype for the three-dimensional fold of the PAS domain superfamily

PAS domains are found in diverse proteins throughout all three kingdoms of life, where they apparently function in sensing and signal transduction. Although a wealth of useful sequence and functional information has become recently available, these data have not been integrated into a three-dimensional (3D) framework. The very early evolutionary development and diverse functions of PAS domains have made sequence analysis and modeling of this protein superfamily challenging. Limited sequence similarities between the ∼50-residue PAS repeats and one region of the bacterial blue-light photosensor photoactive yellow protein (PYP), for which ground-state and light-activated crystallographic structures have been determined to high resolution, originally were identified in sequence searches using consensus sequence probes from PAS-containing proteins. Here, we found that by changing a few residues particular to PYP function, the modified PYP sequence probe also could select PAS protein sequences. By mapping a typical ∼150-residue PAS domain sequence onto the entire crystallographic structure of PYP, we show that the PAS sequence similarities and differences are consistent with a shared 3D fold (the PAS/PYP module) with obvious potential for a ligand-binding cavity. Thus, PYP appears to prototypically exhibit all the major structural and functional features characteristic of the PAS domain superfamily: the shared PAS/PYP modular domain fold of ∼125–150 residues, a sensor function often linked to ligand or cofactor (chromophore) binding, and signal transduction capability governed by heterodimeric assembly (to the downstream partner of PYP). This 3D PAS/PYP module provides a structural model to guide experimental testing of hypotheses regarding ligand-binding, dimerization, and signal transduction.