Faculty Bibliography

The case of a 76-year-old woman with eosinophilic fasciitis is presented. Reported etiologic associations and treatment options are discussed

Genetic screens in zebrafish (Danio rerio) have isolated mutations in hundreds of genes essential for vertebrate development, physiology, and behavior. We have constructed a genetic linkage map that will facilitate the identification of candidate genes for these mutations and allow comparisons among the genomes of zebrafish and other vertebrates. On this map, we have localized 771 zebrafish genes and expressed sequence tags (ESTs) by scoring single-stranded conformational polymorphisms (SSCPs) in a meiotic mapping panel. Of these sequences, 642 represent previously unmapped genes and ESTs. The mapping panel was comprised of 42 homozygous diploid individuals produced by heat shock treatment of haploid embryos at the one-cell stage (HS diploids). This "doubled haploid" strategy combines the advantages of mapping in haploid and standard diploid systems, because heat shock diploid individuals have only one allele at each locus and can survive to adulthood, enabling a relatively large quantity of genomic DNA to be prepared from each individual in the mapping panel. To integrate this map with others, we also scored 593 previously mapped simple-sequence length polymorphisms (SSLPs) in the mapping panel. This map will accelerate the molecular analysis of zebrafish mutations and facilitate comparative analysis of vertebrate genomes.

After the atomic coordinates themselves, the most important data in a homology model are the spatial reliability estimates associated with each of the atoms (atom annotation). Recent blind homology modeling predictions have demonstrated that principally correct sequence-structure alignments are achievable to sequence identities as low as 25% [Martin, A.C., MacArthur, M.W., Thornton, J.M., 1997. Assessment of comparative modeling in CASP2. Proteins Suppl(1), 14-28]. The locations and extent of spatial deviations in the backbone between correctly aligned homologous protein structures remained very poorly estimated however, and these errors were the cause of errant loop predictions [Abagyan, R., Batalov, S., Cardozo, T., Totrov, M., Webber, J., Zhou, Y., 1997. Homology modeling with internal coordinate mechanics: deformation zone mapping and improvements of models via conformational search. Proteins Suppl(1), 29-37]. In order to derive accurate measures for local backbone deviations, we made a systematic study of static local backbone deviations between homologous pairs of protein structures. We found that 'through space' proximity to gaps and chain termini, local three-dimensional 'density', three-dimensional environment conservation, and B-factor of the template contribute to local deviations in the backbone in addition to local sequence identity. Based on these finding, we have identified the meaningful ranges of values within which each of these parameters correlates with static local backbone deviation and produced a combined scoring function to greatly improve the estimation of local backbone deviations. The optimized function has more than twice the accuracy of local sequence identity or B-factor alone and was validated in a recent blind structure prediction experiment. This method may be used to evaluate the utility of a preliminary homology model for a particular biological investigation (e.g. drug design) or to provide an improved starting point for molecular mechanics loop prediction methods

Pathogenicity islands are chromosomal gene clusters, often located adjacent to tRNA genes, that encode virulence factors present in pathogenic organisms but absent or sporadically found in related non-pathogenic species. The selC tRNA locus is the site of integration of different pathogenicity islands in uropathogenic Escherichia coli, enterohaemorrhagic E. coli and Salmonella enterica. We show here that the selC locus of Shigella flexneri, the aetiological agent of bacterial dysentery, also contains a pathogenicity island. This pathogenicity island, designated SHI-2 (Shigella island 2), occupies 23.8 kb downstream of selC and contains genes encoding the aerobactin iron acquisition siderophore system, colicin V immunity and several novel proteins. Remnants of multiple mobile genetic elements are present in SHI-2. SHI-2-hybridizing sequences were detected in all S. flexneri strains tested and parts of the island were also found in other Shigella species. SHI-2 may allow Shigella survival in stressful environments, such as those encountered during infection

Genetic screens in zebrafish (Danio rerio) have isolated mutations in hundreds of genes with essential functions. To facilitate the identification of candidate genes for these mutations, we have genetically mapped 104 genes and expressed sequence tags by scoring single-strand conformational polymorphisms in a panel of haploid siblings. To integrate this map with existing genetic maps, we also scored 275 previously mapped genes, microsatellites, and sequence-tagged sites in the same haploid panel. Systematic phylogenetic analysis defined likely mammalian orthologs of mapped zebrafish genes, and comparison of map positions in zebrafish and mammals identified significant conservation of synteny. This comparative analysis also identified pairs of zebrafish genes that appear to be orthologous to single mammalian genes, suggesting that these genes arose in a genome duplication that occurred in the teleost lineage after the divergence of fish and mammal ancestors. This comparative map analysis will be useful in predicting the locations of zebrafish genes from mammalian gene maps and in understanding the evolution of the vertebrate genome

Phosphatidylinositol 3-kinase (PI3K) mediates a variety of cellular responses by generating PtdIns(3,4)P2 and PtdIns(3,4,5)P3. These 3-phosphoinositides then function directly as second messengers to activate downstream signaling molecules by binding pleckstrin homology (PH) domains in these signaling molecules. We have established a novel assay in the yeast Saccharomyces cerevisiae to identify proteins that bind PtdIns(3,4)P2 and PtdIns(3,4,5)P3 in vivo which we have called TOPIS (Targets of PI3K Identification System). The assay uses a plasma membrane-targeted Ras to complement a temperature-sensitive CDC25 Ras exchange factor in yeast. Coexpression of PI3K and a fusion protein of activated Ras joined to a PH domain known to bind PtdIns(3,4)P2 (AKT) or PtdIns(3,4,5)P3 (BTK) rescues yeast growth at the non-permissive temperature of 37 degreesC. Using this assay, we have identified several amino acids in the beta1-beta2 region of PH domains that are critical for high affinity binding to PtdIns(3,4)P2 and/or PtdIns(3,4,5)P3, and we have proposed a structural model for how these PH domains might bind PI3K products with high affinity. From these data, we derived a consensus sequence which predicts high-affinity binding to PtdIns(3, 4)P2 and/or PtdIns(3,4,5)P3, and we have identified several new PH domain-containing proteins that bind PI3K products, including Gab1, Dos, myosinX, and Sbf1. Use of this assay to screen for novel cDNAs which rescue yeast at the non-permissive temperature should provide a powerful approach for uncovering additional targets of PI3K

Apoptosis (programmed cell death, PCD) is a characteristic type of cell death in which a regulated cellular response pathway mediated by cysteine proteases of the caspase family and Bcl-2 family proteins results in ordered and non-inflammatory involution of the cell. The CED-4 protein and its recently identified mammalian homologue Apaf-1 are critical but functionally uncharacterized components of the cell death machinery. We present here a three-dimensional molecular model for the central domain of CED-4, its alternatively spliced transcript (CED-41) and Apaf-1. A novel protein family is identified and structure prediction for the family identifies a G-protein-like fold with high reliability. The three-dimensional model provides a potential structural explanation for the alternatively spliced variant as well as the known point mutations in CED-4. Regions of the CED-4 and Apaf-1 sequences which may interact with caspases and the Bcl-2 family are proposed. This new information provides a structural molecular framework for the interaction of CED-4-like proteins with the caspases and the Bcl-2 family in the regulation of apoptosis which is analogous to G-protein mediated interactions in well-defined signal transduction pathways

Five models by homology containing insertions and deletions and ranging from 33% to 48% sequence identity to the known homologue, and one high sequence identity (85%) model were built for the CASP2 meeting. For all five low identity targets: (i) our starting models were improved by the Internal Coordinate Mechanics (ICM) energy optimization, (ii) the refined models were consistently better than those built with the automatic SWISS-MODEL program, and (iii) the refined models differed by less than 2% from the best model submitted, as judged by the residue contact area difference (CAD) measure [Abagyan, R.A., Totrov, M.J. Mol. Biol. 268:678-685, 1997]. The CAD measure is proposed for ranking models built by homology instead of global root-mean-square deviation, which is frequently dominated by insignificant yet large contributions from incorrectly predicted fragments or side chains. We demonstrate that the precise identification of regions of local backbone deviation is an independent and crucial step in the homology modeling procedure after alignment, since aligned fragments can strongly deviate from the template at various distances from the alignment gap or even in the ungapped parts of the alignment. We show that a local alignment score can be used as an indicator of such local deviation. While four short loops of the meeting targets were predicted by database search, the best loop 1 target T0028, for which the correct database fragment was not found, was predicted by Internal Coordinate Mechanics global energy optimization at 1.2 A accuracy. A classification scheme for errors in homology modeling is proposed