PubMed Health⌕ Search

Biomedical subjects

I S Mian

Publications and source records attributed to I S Mian.

At least 19 recordsLinked to original sources

Integrated analysis of transcript profiling and protein sequence data.

Transcript profiling can be used to elucidate the molecular and cellular mechanisms involved in ageing and cancer. A recent study of human gastrointestinal stromal tumours (GISTs) with mutations in the KIT gene, Cancer Res. 61 (2001) 8624 exemplifies a common type of investigation. cDNA microarrays were used to generate measurements for 1987 clones in two types of tissues: 13 KIT mutation-positive GISTs and 6 spindle cell tumours from locations outside the gastrointestinal tract. Statistical problems associated with such two-class, high-dimensional profiling data include simultaneous classification and relevant feature identification, probabilistic clustering and protein sequence family modelling. Here, the GIST data were reexamined using specific solutions to these problems, namely sparse hyperplanes, nai;ve Bayes models and profile hidden Markov models respectively. The integrated analysis of molecular profiling and sequence data highlighted 6 clones that may be of clinical and experimental interest. The protein encoded by one of these putative biomarkers defined a novel protein family present in diverse eucarya. The family may be involved in chromosome segregation and/or stability. One family member is a potential biomarker identified recently from a retrospective analysis of transcript profiles for sporadic breast cancer samples from patients with poor and good prognosis, Signal Process. (in press).

Amino Acid Sequence↗

A model for interphase chromosomes and evaluation of radiation-induced aberrations.

We have developed a theoretical model for evaluating radiation-induced chromosomal exchanges by explicitly taking into account interphase (G(0)/G(1)) chromosome structure, nuclear organization of chromosomes, the production of double-strand breaks (DSBs), and the subsequent rejoinings in a faithful or unfaithful manner. Each of the 46 chromosomes for human lymphocytes (40 chromosomes for mouse lymphocytes) is modeled as a random polymer inside a spherical volume. The chromosome spheres are packed randomly inside a spherical nucleus with an allowed overlap controlled by a parameter Omega. The rejoining of DSBs is determined by a Monte Carlo procedure using a Gaussian proximity function with an interaction range parameter sigma. Values of Omega and sigma have been found which yield calculated results of interchromosomal aberration frequencies that agree with a wide range of experimental data. Our preferred solution is one with an interaction range of 0.5 microm coupled with a relatively small overlap parameter of 0.675 microm, which more or less confirms previous estimates. We have used our model with these parameter values and with resolution or detectability limits to calculate yields of translocations and dicentrics for human lymphocytes exposed to low-LET radiation that agree with experiments in the dose range 0.09 to 4 Gy. Five different experimental data sets have been compared with the theoretical results. Essentially all of the experimental data fall between theoretical curves corresponding to resolution limits of 1 Mbp and 20 Mbp, which may reflect the fact that different investigators use different limits for sensitivity or detectability. Translocation yields for mouse lymphocytes have also been calculated and are in good agreement with experimental data from 1 cGy to 10 cGy. There is also good agreement with recent data on complex aberrations. Our model is expected to be applicable to both low- and high-LET radiation, and we include a sample prediction of the yield of interchromosomal rejoining in the dose range 0.22 Gy to 2 Gy of 1000 MeV/nucleon iron particles. This dose range corresponds to average particle traversals per nucleus ranging from 1.0 to 9.12.

Animals↗

C. elegans clk-2, a gene that limits life span, encodes a telomere length regulator similar to yeast telomere binding protein Tel2p.

An important quest in modern biology is to identify genes involved in aging. Model organisms such as the nematode Caenorhabditis elegans are particularly useful in this regard. The C. elegans genome has been sequenced [1], and single gene mutations that extend adult life span have been identified [2]. Among these longevity-controlling loci are four apparently unrelated genes that belong to the clk family. In mammals, telomere length and structure can influence cellular, and possibly organismal, aging. Here, we show that clk-2 encodes a regulator of telomere length in C. elegans.

Aging↗

Conversion of the bifunctional 8-oxoguanine/beta-delta apurinic/apyrimidinic DNA repair activities of Drosophila ribosomal protein S3 into the human S3 monofunctional beta-elimination catalyst through a single amino acid change.

The Drosophila S3 ribosomal protein has important roles in both protein translation and DNA repair. In regards to the latter activity, it has been shown that S3 contains vigorous N-glycosylase activity for the removal of 8-oxoguanine residues in DNA that leaves baseless sites in their places. Drosophila S3 also possesses an apurinic/apyrimidinic (AP) lyase activity in which the enzyme catalyzes a beta-elimination reaction that cleaves phosphodiester bonds 3' and adjacent to an AP lesion in DNA. In certain situations, this is followed by a delta-elimination reaction that ultimately leads to the formation of a single nucleotide gap in DNA bordered by 5'- and 3'-phosphate groups. The human S3 protein, although 80% identical to its Drosophila homolog and shorter by only two amino acids, has only marginal N-glycosylase activity. Its lyase activity only cleaves AP DNA by a beta-elimination reaction, thus further distinguishing itself from the Drosophila S3 protein in lacking a delta-elimination activity. Using a hidden Markov model analysis based on the crystal structures of several DNA repair proteins, the enzymatic differences between Drosophila and human S3 were suggested by the absence of a conserved glutamine residue in human S3 that usually resides at the cleft of the deduced active site pocket of DNA glycosylases. Here we show that the replacement of the Drosophila glutamine by an alanine residue leads to the complete loss of glycosylase activity. Unexpectedly, the delta-elimination reaction at AP sites was also abrogated by a change in the Drosophila glutamine residue. Thus, a single amino acid change converted the Drosophila activity into one that is similar to that possessed by the human S3 protein. In support of this were experiments executed in vivo that showed that human S3 and the Drosophila site-directed glutamine-changed S3 performed poorly when compared with Drosophila wild-type S3 and its ability to protect a bacterial mutant from the harmful effects of DNA-damaging agents.

Amino Acid Sequence↗

Identification of the yeast cytidine deaminase CDD1 as an orphan C-->U RNA editase.

Yeast co-expressing rat APOBEC-1 and a fragment of human apolipoprotein B (apoB) mRNA assembled functional editosomes and deaminated C6666 to U in a mooring sequence-dependent fashion. The occurrence of APOBEC-1-complementing proteins suggested a naturally occurring mRNA editing mechanism in yeast. Previously, a hidden Markov model identified seven yeast genes encoding proteins possessing putative zinc-dependent deaminase motifs. Here, only CDD1, a cytidine deaminase, is shown to have the capacity to carry out C-->U editing on a reporter mRNA. This is only the second report of a cytidine deaminase that can use mRNA as a substrate. CDD1-dependent editing was growth phase regulated and demonstrated mooring sequence-dependent editing activity. Candidate yeast mRNA substrates were identified based on their homology with the mooring sequence-containing tripartite motif at the editing site of apoB mRNA and their ability to be edited by ectopically expressed APOBEC-1. Naturally occurring yeast mRNAs edited to a significant extent by CDD1 were, however, not detected. We propose that CDD1 be designated an orphan C-->U editase until its native RNA substrate, if any, can be identified and that it be added to the CDAR (cytidine deaminase acting on RNA) family of editing enzymes.

APOBEC-1 Deaminase↗

Identifying marker genes in transcription profiling data using a mixture of feature relevance experts.

Transcription profiling experiments permit the expression levels of many genes to be measured simultaneously. Given profiling data from two types of samples, genes that most distinguish the samples (marker genes) are good candidates for subsequent in-depth experimental studies and developing decision support systems for diagnosis, prognosis, and monitoring. This work proposes a mixture of feature relevance experts as a method for identifying marker genes and illustrates the idea using published data from samples labeled as acute lymphoblastic and myeloid leukemia (ALL, AML). A feature relevance expert implements an algorithm that calculates how well a gene distinguishes samples, reorders genes according to this relevance measure, and uses a supervised learning method [here, support vector machines (SVMs)] to determine the generalization performances of different nested gene subsets. The mixture of three feature relevance experts examined implement two existing and one novel feature relevance measures. For each expert, a gene subset consisting of the top 50 genes distinguished ALL from AML samples as completely as all 7,070 genes. The 125 genes at the union of the top 50s are plausible markers for a prototype decision support system. Chromosomal aberration and other data support the prediction that the three genes at the intersection of the top 50s, cystatin C, azurocidin, and adipsin, are good targets for investigating the basic biology of ALL/AML. The same data were employed to identify markers that distinguish samples based on their labels of T cell/B cell, peripheral blood/bone marrow, and male/female. Selenoprotein W may discriminate T cells from B cells. Results from analysis of transcription profiling data from tumor/nontumor colon adenocarcinoma samples support the general utility of the aforementioned approach. Theoretical issues such as choosing SVM kernels and their parameters, training and evaluating feature relevance experts, and the impact of potentially mislabeled samples on marker identification (feature selection) are discussed.

Acute Disease↗

Analysis of telomerase processivity: mechanistic similarity to HIV-1 reverse transcriptase and role in telomere maintenance.

The key protein subunit of the telomerase complex, known as TERT, possesses a reverse transcriptase (RT)-like domain that is conserved in enzymes encoded by retroviruses and retroelements. Structural and functional analysis of HIV-1 RT suggests that RT processivity is governed, in part, by the conserved motif C, motif E, and a C-terminal domain. Mutations in analogous regions of the yeast TERT were found to have anticipated effects on telomerase processivity in vitro, suggesting a great deal of mechanistic and structural similarity between TERT and retroviral RTs, and a similarity that goes beyond the homologous domain. A close correlation was uncovered between telomerase processivity and telomere length in vivo, suggesting that enzyme processivity is a limiting factor for telomere maintenance.

Catalytic Domain↗

SATB1 cleavage by caspase 6 disrupts PDZ domain-mediated dimerization, causing detachment from chromatin early in T-cell apoptosis.

SATB1 is expressed primarily in thymocytes and orchestrates temporal and spatial expression of a large number of genes in the T-cell lineage. SATB1 binds to the bases of chromatin loop domains in vivo, recognizing a special DNA context with strong base-unpairing propensity. The majority of thymocytes are eliminated by apoptosis due to selection processes in the thymus. We investigated the fate of SATB1 during thymocyte and T-cell apoptosis. Here we show that SATB1 is specifically cleaved by a caspase 6-like protease at amino acid position 254 to produce a 65-kDa major fragment containing both a base-unpairing region (BUR)-binding domain and a homeodomain. We found that this cleavage separates the DNA-binding domains from amino acids 90 to 204, a region which we show to be a dimerization domain. The resulting SATB1 monomer loses its BUR-binding activity, despite containing both its DNA-binding domains, and rapidly dissociates from chromatin in vivo. We found this dimerization region to have sequence similarity to PDZ domains, which have been previously shown to be involved in signaling by conferring protein-protein interactions. SATB1 cleavage during Jurkat T-cell apoptosis induced by an anti-Fas antibody occurs concomitantly with the high-molecular-weight fragmentation of chromatin of ~50-kb fragments. Our results suggest that mechanisms of nuclear degradation early in apoptotic T cells involve efficient removal of SATB1 by disrupting its dimerization and cleavage of genomic DNA into loop domains to ensure rapid and efficient disassembly of higher-order chromatin structure.

Amino Acid Sequence↗

Analysis of molecular profile data using generative and discriminative methods.

A modular framework is proposed for modeling and understanding the relationships between molecular profile data and other domain knowledge using a combination of generative (here, graphical models) and discriminative [Support Vector Machines (SVMs)] methods. As illustration, naive Bayes models, simple graphical models, and SVMs were applied to published transcription profile data for 1,988 genes in 62 colon adenocarcinoma tissue specimens labeled as tumor or nontumor. These unsupervised and supervised learning methods identified three classes or subtypes of specimens, assigned tumor or nontumor labels to new specimens and detected six potentially mislabeled specimens. The probability parameters of the three classes were utilized to develop a novel gene relevance, ranking, and selection method. SVMs trained to discriminate nontumor from tumor specimens using only the 50-200 top-ranked genes had the same or better generalization performance than the full repertoire of 1,988 genes. Approximately 90 marker genes were pinpointed for use in understanding the basic biology of colon adenocarcinoma, defining targets for therapeutic intervention and developing diagnostic tools. These potential markers highlight the importance of tissue biology in the etiology of cancer. Comparative analysis of molecular profile data is proposed as a mechanism for predicting the physiological function of genes in instances when comparative sequence analysis proves uninformative, such as with human and yeast translationally controlled tumour protein. Graphical models and SVMs hold promise as the foundations for developing decision support systems for diagnosis, prognosis, and monitoring as well as inferring biological networks.

Bayes Theorem↗

Integrating naive Bayes models and external knowledge to examine copper and iron homeostasis in S. cerevisiae.

A novel suite of analytical techniques and visualization tools are applied to 78 published transcription profiling experiments monitoring 5,687 Saccharomyces cerevisiae genes in studies examining cell cycle, responses to stress, and diauxic shift. A naive Bayes model discovered and characterized 45 classes of gene profile vectors. An enrichment measure quantified the association between these classes and specific external knowledge defined by four sets of categories to which genes can be assigned: 106 protein functions, 5 stages of the cell cycle, 265 transcription factors, and 16 chromosomal locations. Many of the 38 genes in class 42 are known to play roles in copper and iron homeostasis. The 17 uncharacterized open reading frames in this class may be involved in similar homeostatic processes; human homologs of two of them could be associated with as yet undefined disease states arising from aberrant metal ion regulation. The Met4, Met31, and Met32 transcription factors may play a role in coregulating genes involved in copper and iron metabolism. Extensions of the simple graphical model used for clustering to learning more complex models of genetic networks are discussed.

Bayes Theorem↗

Prediction of structural and functional relationships of Repeat 1 of human interphotoreceptor retinoid-binding protein (IRBP) with other proteins.

PURPOSE: We compared the structure and function of interphotoreceptor retinoid-binding protein (IRBP) related proteins and predicted domain and secondary structure within each repeat of IRBP and its relatives. We tested whether tail specific protease (Tsp), which bears sequence similarity to IRBP Domain B, binds fatty acids or retinoids, and whether IRBP possessed protease activity resembling Tsp's catalytic function. These tests helped us to learn whether the primary sequence similarities of family members extended to higher order structural and functional levels. METHODS: Predictions derived from multiple sequence alignments among IRBP and Tsp family members and secondary structure computer programs were carried out. The first repeat of human IRBP (EcR1) and Tsp were expressed, purified, and tested for binding properties. Tsp was examined for fluorescence enhancement of retinol or 16-anthroyloxy-palmitic acid (16-AP) to test for ligand binding. IRBP was tested for protease activity. RESULTS: Tsp did not exhibit fluorescence enhancement with retinol or 16-AP. IRBP did not exhibit protease activity. The positions of critical residues needed for the ligand binding properties of retinol were predicted. Primary sequence and three-dimensional similarity was found between Domain A of IRBP Repeat 3 and eglin c. CONCLUSIONS: The sequence similarity of Tsp and IRBP raised the possibility that each might share the function of the other protein: IRBP might possess protease activity or Tsp might possess retinoid or fatty acid binding activity. Our studies do not support such a shared function hypothesis, and suggest that the sequence similarity is the result of maintenance of structure. The finding of similarity to eglin c in Domain A suggests the possibility of a tight interaction between Domain A and Domain B, possibly implying the need for Domain A in retinoid-binding, and suggesting that both Domains should be present in testing mutations. The positions of predicted critical amino acids suggest models in which a large binding pocket holds the retinoid or fatty acid ligand. These predictions are tested in a companion paper.

Cluster Analysis↗

Effects of dispersed point substitutions in Repeat 1 of human interphotoreceptor retinoid binding protein (IRBP).

PURPOSE: The purpose of this study was to measure the effects of mutations on the retinol binding capability of human Repeat 1 of interphotoreceptor retinoid-binding protein (IRBP). First, we predicted important functional amino acids by several computer programs. We also noted the lack of shared functions between Tail-specific protease (Tsp) and IRBP, which bear sequence similarity, and this aided in predicting functional residues. We analyzed the effects of point substitutions on the retinol and fatty acid binding properties of Repeat 1 of human IRBP at 25 and 50 degrees C. METHODS: To find residues critical to retinol binding that might affect function, a series of thirteen mutations were created by site-specific mutagenesis between positions 140 and 280 in Repeat 1 of human IRBP. These mutants were expressed, purified, and tested for binding properties. The conformations of the proteins were examined by circular dichroism (CD) scans. RESULTS: Seven of the mutations exhibited reduced binding capacity, and five were not expressed at high enough levels to assess binding activity. Four of the mutants were purified, and their CD scans were very similar to those of Repeat 1. Only one of the mutations did not affect binding, folding, or expression when compare to wild type Repeat 1. CONCLUSIONS: Several IRBP mutants containing point mutations retained native structure but lost retinol binding function. The data suggest that retinol binding is affected by many different amino acid substitutions in or near a binding pocket. That even a single point substitution can profoundly affect binding without affecting overall conformation suggests that much of Domain B (from amino acid positions 80 to 300) is involved with ligand binding. This excludes three previously proposed IRBP-retinol binding mechanisms: (1) retinol binds to a small portion of the protein repeat, (2) retinol can bind to any hydrophobic patch in the protein, and (3) native conformation is not required for retinol binding to the repeat.

Amino Acid Substitution↗

Representing and reasoning about protein families using generative and discriminative methods.

This work addresses the issues of data representation and incorporation of domain knowledge into the design of learning systems for reasoning about protein families. Given the limited expressive capacity of a particular method, a mixture of protein annotation and fold recognition experts, each implementing a different underlying representation, should provide a robust method for assigning sequences to families. These ideas are illustrated using two data-driven learning methods that make use of different prior information and employ independent, yet complementary, projections of a family: hidden Markov models (HMMs) based on a multiple sequence alignment and neural networks (NNs) based on global sequence descriptors of proteins. Examination of seven protein families indicates that combining a generative (HMM) and a discriminative (NN) method is better than either method on its own. Biologically, human 4-hydroxyphenylpyruvic acid dioxygenase, involved in tyrosinemia type 3, is predicted to be structurally and functionally related to the glyoxalase I family.

Amino Acid Sequence↗

AZU-1: a candidate breast tumor suppressor and biomarker for tumor progression.

To identify genes misregulated in the final stages of breast carcinogenesis, we performed differential display to compare the gene expression patterns of the human tumorigenic mammary epithelial cells, HMT-3522-T4-2, with those of their immediate premalignant progenitors, HMT-3522-S2. We identified a novel gene, called anti-zuai-1 (AZU-1), that was abundantly expressed in non- and premalignant cells and tissues but was appreciably reduced in breast tumor cell types and in primary tumors. The AZU-1 gene encodes an acidic 571-amino-acid protein containing at least two structurally distinct domains with potential protein-binding functions: an N-terminal serine and proline-rich domain with a predicted immunoglobulin-like fold and a C-terminal coiled-coil domain. In HMT-3522 cells, the bulk of AZU-1 protein resided in a detergent-extractable cytoplasmic pool and was present at much lower levels in tumorigenic T4-2 cells than in their nonmalignant counterparts. Reversion of the tumorigenic phenotype of T4-2 cells, by means described previously, was accompanied by the up-regulation of AZU-1. In addition, reexpression of AZU-1 in T4-2 cells, using viral vectors, was sufficient to reduce their malignant phenotype substantially, both in culture and in vivo. These results indicate that AZU-1 is a candidate breast tumor suppressor that may exert its effects by promoting correct tissue morphogenesis.

Amino Acid Sequence↗

Identification of functionally important domains in the N-terminal region of telomerase reverse transcriptase.

Telomerase is a ribonucleoprotein reverse transcriptase responsible for the maintenance of one strand of telomere terminal repeats. The key protein subunit of the telomerase complex, known as TERT, possesses reverse transcriptase-like motifs that presumably mediate catalysis. These motifs are located in the C-terminal region of the polypeptide. Hidden Markov model-based sequence analysis revealed in the N-terminal region of all TERTs the presence of four conserved motifs, named GQ, CP, QFP, and T. Point mutation analysis of conserved residues confirmed the functional importance of the GQ motif. In addition, the distinct phenotypes of the GQ mutants suggest that this motif may play at least two distinct functions in telomere maintenance. Deletion analysis indicates that even the most N-terminal nonconserved region of yeast TERT (N region) is required for telomerase function. This N region exhibits a nonspecific nucleic acid binding activity that probably reflects an important physiologic function. Expression studies of various portions of the yeast TERT in Escherichia coli suggest that the N region and the GQ motif together may constitute a stable domain. We propose that all TERTs may have a bipartite organization, with an N-GQ domain connected to the other motifs through a flexible linker.

Amino Acid Sequence↗

Chromatin conformation in living cells: support for a zig-zag model of the 30 nm chromatin fiber.

A new method was used to probe the conformation of chromatin in living mammalian cells. The method employs ionizing radiation and is based on the concept that such radiation induces correlated breaks in DNA strands that are in spatial proximity. Human dermal fibroblasts in G0 phase of the cell cycle and Chinese hamster ovary cells in mitosis were irradiated by X-rays or accelerated ions. Following lysis of the cells, DNA fragments induced by correlated breaks were end-labeled and separated according to size on denaturing polyacrylamide gels. A characteristic peak was obtained for a fragment size of 78 bases, which is the size that corresponds to one turn of DNA around the nucleosome. Additional peaks between 175 and 450 bases reflect the relative position of nearest-neighbor nucleosomes. Theoretical calculations that simulate the indirect and direct effect of radiation on DNA demonstrate that the fragment size distributions are closely related to the chromatin structure model used. Comparison of the experimental data with theoretical results support a zig-zag model of the chromatin fiber rather than a simple helical model. Thus, radiation-induced damage analysis can provide information on chromatin structure in the living cell.

Animals↗

Sequence, structural, functional, and phylogenetic analyses of three glycosidase families.

Glycosidases, which cleave the glycosidic bond between a carbohydrate and another moiety, have been classified into over 63 families. Here, a variety of computational techniques have been employed to examine three families important in normal and abnormal pathology with the aim of developing a framework for future homology modeling, experimental and other studies. Family 1 includes bacterial and archaeal enzymes as well as lactase phlorizin-hydrolase and klotho, glycosidases implicated in disaccharide intolerance II and aging respectively. A statistical model, a hidden Markov model (HMM), for the family 1 glycosidase domain was trained and used as the basis for comparative examination of the conserved and variable sequence and structural features as well as the phylogenetic relationships between family members. Although the structures of four family 1 glycosidases have been determined, this is the first comparative examination of all these enzymes. Aspects that are unique to specific members or subfamilies (substrate binding loops) as well those common to all members (a beta/alpha)8 barrel fold) have been defined. Active site residues in some domains in klotho and lactase-phlorizin hydrolases differ from other members and in one instance may bind but not cleave substrate. The four invariant and most highly conserved residues are not residues implicated in catalysis and/or substrate binding. Of these, a histidine may be involved in transition state stabilization. Glucosylceramidase (family 30) and galactosylceramidase (family 59) are mutated in the lysosomal storage disorders Gaucher disease and Krabbe disease, respectively. HMM-based analysis, structure prediction studies and examination of disease mutations reveal a glycosidase domain common to these two families that also occurs in some bacterial glycosidases. Similarities in the reactions catalyzed by families 30 and 59 are reflected in the presence of a structurally and functionally related (beta/alpha)8 barrel fold related to that in family 1.

Amino Acid Sequence↗