Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Computational prediction”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48Linked to original sources

Attraction within the membrane. Forces behind transmembrane protein folding and supramolecular complex assembly.

Biological membranes are fascinating two-dimensional microenvironments that exhibit unique solvent behaviours due to their varying lipid composition. Although many important bioenergetic and signalling events involve the transient or permanent assembly of membrane protein complexes, the characterization of the thermodynamic and kinetic properties behind this assembly is just beginning. In particular, the molecular forces that govern protein association within these structures remain poorly understood. An understanding of the docking of transmembrane proteins to supramolecular complexes, which will make possible the development of predictive computational tools, will require detailed knowledge of interaction forces at the atomistic or residue level. Here, I review current data on supramolecular complexes in membrane environments and make a tentative comparison between assembly processes in membranes and those driven by the hydrophobic effect in water. This comparison suggests that, in addition to being controlled by specific characteristics of the lipid molecules themselves, molecular assembly in the membrane milieu also depends more generally on the entropy of the lipid fraction.

Amino Acid Substitution↗

The possible value of synthetic peptides in the diagnosis and therapy of myocarditis and dilated cardiomyopathy.

The adenine nucleotide translocator (ANT) and myosin have been shown to be major autoantigens in myocarditis and dilated cardiomyopathy. We studied the use of synthetic peptides, with sequences derived from ANT and myosin, as antigens in screening tests for autoantibodies in myocarditis (MC) and dilated cardiomyopathy (DCM) and as absorbents for specific elimination of autoantibodies from the sera of patients. Using computer prediction of the secondary structure of the ANT and myosin we identified two sequences of the ANT and three sequences of myosin as possible main antigenic determinants. Using overlapping synthetic peptides and antibodies against them, the antigenicity of the selected determinants was shown. Of 72 sera from patients with MC or DCM 45 (62.5%) bound to the peptides derived from ANT, 32 (44.4%) reacted with the sequences from myosin, in contrast to healthy controls. Using the peptides from the ANT or myosin immobilized on thiopropyl-sepharose, more than 95% of the autoantibodies could be removed specifically from the positive sera. The results demonstrate the usefulness of synthetic peptides as antigens in antibody screening tests in MC and DCM and offers a new approach to the therapy of inflammatory heart disease by specific elimination of autoantibodies.

Autoantibodies↗

Disruption of a novel MFS transporter gene, DIRC2, by a familial renal cell carcinoma-associated t(2;3)(q35;q21).

Previously, we described a family with a significantly increased predisposition for renal cell cancer co-segregating with a t(2;3)(q35;q21) chromosomal translocation. Several primary tumors of the clear cell type from different family members were analyzed at a molecular level. Loss of the derivative chromosome 3 was consistently found. In addition, different somatic Von Hippel Lindau (VHL) gene mutations were observed in most of the tumors analyzed, even within the same patient. Based on these results a multistep tumorigenesis model was proposed in which (non-disjunctional) loss of the derivative chromosome 3 represents an early event and somatic mutation of the VHL gene represents a late event related to tumor progression. More recently, however, we noted that these two anomalies were absent in at least one early-stage tumor sample that we tested. Similar results were obtained in another family with renal cell cancer and t(3;6)(q12;q15), thus suggesting that another genetic event may precede these two oncogenetic steps. We speculate that deregulation of a gene(s) located at or near the translocation breakpoint may act as such. In order to identify such genes, a detailed physical map encompassing the 3q21 breakpoint region was constructed. Through a subsequent positional cloning effort we found that this breakpoint targets a hitherto unidentified gene, designated DIRC2 (disrupted in renal cancer 2). Computer predictions of the putative DIRC2 protein showed significant homology to different members of the major facilitator superfamily (MFS) of transporters. Based on additional DIRC2 expression and mutation analyses, we propose that the observed gene disruption may result in haplo-insufficiency and, through this mechanism, in the onset of tumor growth.

Amino Acid Sequence↗

ALS2, a novel guanine nucleotide exchange factor for the small GTPase Rab5, is implicated in endosomal dynamics.

ALS2 mutations account for a number of recessive motor neuron diseases including forms of amyotrophic lateral sclerosis, primary lateral sclerosis and hereditary spastic paraplegia. Although computational predictions suggest that ALS2 encodes a protein containing multiple guanine nucleotide exchange factor (GEF) domains [RCC1-like domain (RLD), the Dbl homology and pleckstrin homology (DH/PH), and the vacuolar protein sorting 9 (VPS9)], the functions of the ALS2 protein have not been revealed as yet. Here we show that the ALS2 protein specifically binds to small GTPase Rab5 and functions as a GEF for Rab5. Ectopically expressed ALS2 protein localizes with Rab5 and early endosome antigen-1 (EEA1) onto early endosomal compartments and stimulates the enlargement of endosomes in cultured cortical neurons. The carboxy-terminus of ALS2 protein carrying a VPS9 domain mediates not only the activation of Rab5 via a guanine-nucleotide exchanging reaction but also the endosomal localization of the ALS2 protein, while the amino-terminal half containing RLD acts suppressive in its membranous localization. Further, the DH/PH domain in the middle portion of ALS2 protein enhances the VPS9 domain-mediated endosome fusions. Taken together, the ALS2 protein as a novel Rab5-GEF, ALS2rab5GEF seems to be implicated in the endosomal dynamics in vivo. Notably, a feature common to eight reported ALS2 mutations among motor neuron diseases is the loss of VPS9 domain, resulting in the failure of Rab5 activation. Thus, a perturbation of endosomal dynamics caused by loss of ALS2 rab5GEF activity might underlie neuronal dysfunction and degeneration in a number of motor neuron diseases.

Endosomes↗

Classifying the estrogen receptor status of breast cancers by expression profiles reveals a poor prognosis subpopulation exhibiting high expression of the ERBB2 receptor.

Recent work using expression profiling to computationally predict the estrogen receptor (ER) status of breast tumors has revealed that certain tumors are characterized by a high prediction uncertainty ('low-confidence'). We analyzed these 'low-confidence' tumors and determined that their 'uncertain' prediction status arises as a result of widespread perturbations in multiple genes whose expression is important for ER subtype discrimination. Patients with 'low-confidence' ER+ tumors exhibited a significantly worse overall survival (P=0.03) and shorter time to distant metastasis (P=0.004) compared with their 'high-confidence' ER+ counterparts, indicating that the 'high-' and 'low-confidence' binary distinction is clinically meaningful. We then discovered that elevated expression of the ERBB2 receptor is significantly correlated with a breast tumor exhibiting a 'low-confidence' prediction, and this association was subsequently validated across multiple independently derived breast cancer expression datasets employing a variety of different array technologies and patient populations. Although ERBB2 signaling has been proposed to inhibit the transcriptional activity of ER, a large proportion of the perturbed genes in the 'low-confidence'/ERBB2+ samples are not known to be estrogen responsive, and a recently described bioinformatic algorithm (DEREF) was used to demonstrate the absence of potential estrogen-response elements (EREs) in their promoters. We propose that a significant portion of ERBB2's effects on ER+ breast tumors may involve ER-independent mechanisms of gene activation, which may contribute to the clinically aggressive behavior of the 'low-confidence' breast tumor subtype.

Algorithms↗

Response of established human breast tumors to vaccination with mammaglobin-A cDNA.

BACKGROUND: A novel breast cancer-associated antigen, mammaglobin-A, is expressed in 80% of primary breast tumors. The characterization of immune responses against this highly expressed breast cancer-specific antigen would be of value in the development of new therapeutic strategies for breast cancer. METHODS: We developed an in vivo model using human leukocyte antigen-A*0201/human CD8+ (HLA-A2+/hCD8+) double-transgenic mice to define the epitopes and to study the level of protection acquired by mammaglobin-A cDNA vaccination toward mammaglobin-A+/HLA-A2+ breast cancer cell lines. Mammaglobin-A epitopes were identified using an HLA class I peptide binding prediction computer program, and their activity was verified using gamma interferon ELISPOT and cytotoxicity assays. RESULTS: We identified seven mammaglobin-A-derived candidate epitopes that bind the HLA-A*0201 molecule (Mam-A2.1-7). CD8+ cytotoxic T lymphocytes (CTLs) from HLA-A2+/hCD8+ mice reacted to the Mam-A2.1 (amino acids [aa] 83-92, LIYDSSLCDL), Mam-A2.2 (aa 2-10, KLLMVLMLA), Mam-A2.4 (aa 66-74, FLNQTDETL), and Mam-A2.6 (aa 32-40, MQLIYDSSL) epitopes. CD8+ CTLs from breast cancer patients also recognized a similar epitope pattern as did those in the HLA-A2+/hCD8 mice and reacted to the Mam-A2.1, Mam-A2.2, Mam-A2.3, Mam-A2.4, and Mam-A2.7 epitopes. Passive transfer of mammaglobin-A-reactive CTLs into SCID (severe combined immunodeficient) beige mice with actively growing mammaglobin-A+ tumors resulted in statistically significant regression (P<.001) in the growth of the tumors. CONCLUSIONS: The HLA-A2+/hCD8+ mouse represents a valuable animal model to characterize the HLA-A*0201-restricted CD8+ CTL immune response to mammaglobin-A in vivo, and the data reported here demonstrate the immunotherapeutic potential of mammaglobin-A for the treatment and/or prevention of breast cancer.

Animals↗

Disulfide-bond reshuffling in the evolution of an ape placental ribonuclease.

Disulfide bonds play important roles in the folding and stability of proteins and are evolutionarily conserved. A classic example is RNase A (also known as bovine pancreatic ribonuclease), which contains 4 conserved disulfide bonds among 8 cysteines. However, human RNase 8, a paralog of RNase A uniquely expressed in the placenta, has lost one of the conserved cysteines but gained another, when compared with RNase 8 of various monkeys and with RNase A. We here show that both the loss and gain of the cysteines in human RNase 8 occurred in the common ancestor of African great apes (humans, chimps, and gorillas) 7-13 MYA. Computational predictions suggest changes of disulfide bonding by these cysteine substitutions. Site-directed mutagenesis indicates that if the ribonucleolytic activity is essential for RNase 8's function, the gain of the cysteine must have preceded the loss. Human RNase 8 represents one of the first examples in which the presumable evolutionary change of a disulfide bond involves 1 loss and 1 gain of cysteine, instead of 2 losses or 2 gains. Our results provide the foundation for detailed analysis toward understanding the impact of disulfide-bond reshuffling on the structure, function, and evolution of proteins in general and human RNase 8 in particular.

Amino Acid Sequence↗

Combinations of genetic changes in the human cAMP-responsive element modulator gene: a clue towards understanding some forms of male infertility?

The cAMP-responsive element modulator (CREM) gene plays a pivotal role in the mouse spermatogenesis, but its role in the human infertility has not been fully established. We performed a mutation screening in 13 Slovenian men with round spermatid arrest and in six controls. Eleven genetic changes have been identified in the human CREM gene, three novel single-nucleotide polymorphisms [within the promoters P1, P3 and intervening sequence 1 (IVS1)], one insertion (IVS2) and one non-sense mutation (exon gamma). Some infertile patients seem to accumulate potentially harmful genetic changes. We identified a patient with no CREM immunoreactive protein that was homozygous for the nucleotide changes in all promoters, IVS 1, 2, 6, and was heterozygous for the mutation in exon gamma. Interestingly, insertion in IVS2 (IVS2-58_55insT) results in a four-fold decrease in binding of nuclear proteins. Computer predictions suggested the presence of a potential novel CREM promoter, however, random amplification of cDNA ends from the human testis cDNA library was not successful in confirming a novel transcription start site of the CREM gene. Screening of a larger number of patients and controls is required to elucidate whether the observed combinations of genetic changes in the CREM gene can explain some forms of male infertility.

Adult↗

Conformational analysis of the 5' leader and the gag initiation site of Mo-MuLV RNA and allosteric transitions induced by dimerization.

Dimerization of genomic RNA is a key step in the retroviral life cycle and has been postulated to be involved in the regulation of translation, encapsidation and reverse transcription. Here, we have derived a secondary structure model of nucleotides upstream from psi and of the gag initiation region of Mo-MuLV RNA in monomeric and dimeric forms, using chemical probing, sequence comparison and computer prediction. The 5' domain is extensively base-paired and interactions take place between U5 and 5' leader sequences. The U5-PBS subdomain can fold in two mutually exclusive conformations: a very stable and extended helical structure (E form) in which 17 of the 18 nucleotides of the PBS are paired, or an irregular three-branch structure (B form) in which 10 nucleotides of the PBS are paired. The dimeric RNA adopts the B conformation. The monomeric RNA can switch from the E to the B conformation by a thermal treatment. If the E to B transition is associated to dimerization, it may facilitate annealing of the primer tRNAPro to the PBS by lowering the free energy required for melting the PBS. Furthermore, dimerization induces allosteric rearrangements around the SD site and the gag initiation region.

Allosteric Regulation↗

Analysis of the Escherichia coli genome. IV. DNA sequence of the region from 89.2 to 92.8 minutes.

We present the sequence of 176 kilobases of the Escherichia coli K-12 genome, from katG at 89.2 to an open reading frame (ORF) of unknown function at 92.8 minutes on the genetic map. This brings the total of contiguous sequence from the E. coli genome project to 500 kb (81.5 to 92.8 minutes). This segment contains 134 putative coding genes (ORFs) of which 66 genes were previously identified. Eight new genes--acs, pepE, and nrfB-G--were identified as well as the previously mapped gldA and alr genes. Still, 58 ORFs remain unidentified despite literature and similarity searches. The arrangement of proposed genes relative to possible promoters and terminators suggests 55 potential transcription units. Other features include 13 REP elements, one IRU (ERIC) repeat, 59 computer-predicted bends, 42 Chi sites and one new grey hole. Sixteen signal peptides were found, including those of lamB, btuB, and malE. Two ribosomal RNA loci, rrnB and rrnE, are located in this segment, so we have now sequenced four of the seven E. coli rRNA loci. Comparison of the rRNA loci reveals some differences in the ribosomal structural RNAs which are generally compatible with the proposed secondary structures.

Base Sequence↗

GRSF-1: a poly(A)+ mRNA binding protein which interacts with a conserved G-rich element.

Computer predictions identified similarities to a 14-base G-rich element in numerous mRNAs at a variety of locations. A Northwestern screening strategy was used to obtain a cDNA clone from a HeLa cell library using the G-rich RNA element as a probe. A cellular protein (called GRSF-1), which was encoded by this cDNA, binds RNAs containing the G-rich element. GRSF-1 was distinct from DSEF-1, a nuclear protein we have previously identified that interacts with the G-rich element, based on differences in molecular weight and partial peptide maps, as well as the lack of cross-reactivity with GRSF-1 specific monoclonal antibodies. Using indirect immunofluorescence microscopy, we localized GRSF-1 to the cytoplasm. In vivo UV cross-linking further demonstrated that GRSF-1 was bound to poly(A)+ mRNA in living human cells. Western blot analysis revealed four cytoplasmic proteins which expressed GRSF-1 specific epitopes. GRSF-1 contains three potential RNA recognition motifs and two auxiliary domains. Curiously, the domain organization of GRSF-1 is similar to the RNA binding proteins PUB1, ELAV, HuD, Hel-N1, mcs94-1 and RBP9.

Amino Acid Sequence↗

Analysis of the Escherichia coli genome. V. DNA sequence of the region from 76.0 to 81.5 minutes.

The DNA sequence of a 225.4 kilobase segment of the Escherichia coli K-12 genome is described here, from 76.0 to 81.5 minutes on the genetic map. This brings the total of contiguous sequence from the E.coli genome project to 725.1 kb (76.0 to 92.8 minutes). We found 191 putative coding genes (ORFs) of which 72 genes were previously known, and 110 of which remain unidentified despite literature and similarity searches. Seven new genes--arsE, arsF, arsG, treF, xylR, xylG, and xylH--were identified as well as the previously mapped pit and dctA genes. The arrangement of proposed genes relative to possible promoters and terminators suggests 90 potential transcription units. Other features include 19 REP elements, 95 computer-predicted bends, 50 Chi sites, and one grey hole. Thirty-one putative signal peptides were found, including those of thirteen known membrane or periplasmic proteins. One tRNA gene (proK) and two insertion sequences (IS5 and IS150) are located in this segment. The genes in this region are organized with equal numbers oriented with or against replication.

Alcohol Dehydrogenase↗

Protonatable hairpins are conserved in the 5'-untranslated region of tymovirus RNAs.

The secondary structures of the 5'-untranslated region (5'-UTR) of five different tymoviruses have been determined by structure probing, computer prediction and sequence comparison. Despite large sequence differences, there are remarkable similarities in the secondary structure. In all viruses two or four hairpins are found, most of which contain a symmetrical internal loop consisting of adjacent C-C or C-A mismatches. Since it is known that such mismatches can be protonated and protonated cytosines play an important role in RNA-protein interactions in tymoviral virions, the influence of pH on the conformation of the internal loop was studied. UV melting experiments and 1-dimensional proton NMR at varying pH values and salt concentrations confirm that the hairpins can be protonated under relatively mild conditions. The hairpin found in the 5'-UTR of erysimum latent virus, which has an asymmetrical internal loop consisting of cytosines and uridines, shows comparable behaviour. It is concluded that all tymoviral RNAs contain protonatable hairpins in the 5'-UTR. Binding experiments with empty viral capsids, however, do not yet establish a role in capsid protein binding.

Base Sequence↗

A theoretical approach to select effective antisense oligodeoxyribonucleotides at high statistical probability.

Up to now, out of approximately 20 antisense oligodeoxyribonucleotides (as ODN) selected and tested against a given target gene, only one species shows substantial suppression of target gene expression. In part, this seems to be related to the general assumption that the structures of local target sequences or antisense nucleic acids are unfavorable for efficient annealing. Experimental approaches to find effective as ODN are extremely expensive when including a large number of antisense species and when considering their moderate success. Here, we make use of a systematic alignment of computer-predicted secondary structures of local sequence stretches of the target RNA and of semi-empirical rules to identify favorable local target sequences and, hence, to design more effective as ODN. The intercellular adhesion molecule 1 (ICAM-1) gene was chosen as a target because it had been shown earlier to be sensitive to antisense-mediated gene suppression. By applying the protocol described here, 10 ICAM-1-directed as ODN species were found that showed substantially improved inhibition of target gene expression in the endothelial cell line ECV304 when compared with the most effective published as ODN. Further, 17 out of 34 antisense species (50%) selected on the theoretical basis described here showed significant (>50%) inhibition of ICAM-1 expression in mammalian cells.

Cell Line, Transformed↗

A method for detecting the coding DNA for a protein of known sequence in a collection of chimeric plasmids constructed with cDNAs.

A rapid screening method for recognising plasmids containing copies of cDNA corresponding to proteins of known amino acid sequences, is described. The method is based on the computer prediction of the possible restriction sites in the cloned DNAs. It was tested on a series of proteins with known coding DNA sequences and a series of plasmid cloned cDNAs made from pure chicken globin mRNA.

Amino Acid Sequence↗

On the number of protein-protein interactions in the yeast proteome.

Using two different approaches, we estimated that on average there are about five interacting partners per protein in the proteome of the yeast Saccharomyces cerevisiae. In the first approach, we used a novel method to model sampling overlap by a Bernoulli process, compared the results of two independent yeast two-hybrid interaction screens and tested the robustness of the estimate. The most stable estimate of five interactors per protein was obtained when the three most highly connected nodes in the protein interaction network were removed from the analysis (eight interactors per protein if those nodes were kept). In the second approach, we analysed a published high-confidence subset of putative interaction data obtained from multiple sources, including large-scale two-hybrid screens, complex purifications, synthetic lethals, correlated gene expression, computational predictions and previous annotations. Strikingly, the estimate was again five interactors per protein. These estimates suggest a range of approximately 16,000-26,000 different interaction pairs in the yeast, excluding homotypic interactions. We also discuss the approaches to estimating the rate of homotypic interactions.

Protein Binding↗

A systematic analysis of the silencing effects of an active siRNA at all single-nucleotide mismatched target sites.

The specificity of small interfering RNA (siRNA)-mediated gene silencing is a critical consideration for the application of RNA interference (RNAi). While the discovery of potential off-target effects by siRNAs is of concern, no systematic analysis has been conducted to explore the specificity of RNAi. Here, we present a study where a functionally validated siRNA (siCD46) was examined for silencing specificity on all possible 57 permutated target sites, each carrying a single-nucleotide mutation that would generate a mismatch when paired with siRNA antisense strand. We found that it was not only the position of the mismatched base pair, but also the identity of the nucleotides forming the mismatch that influenced silencing. Surprisingly, mismatches formed between adenine (A) and cytosine (C), in addition to the G:U wobble base pair, were well tolerated and target sites containing such mismatches were silenced almost as efficiently as its fully matched counterpart by siCD46. Northern blots showed that the silencing of fusion genes harboring the mutated target sites involved target mRNA degradation. This study provides direct evidence that the target recognition of siRNA is far more degenerative than previously considered. This finding is instrumental in the understanding of RNAi specificity and may aid the computational prediction of RNA secondary structure.

Base Pair Mismatch↗

Mouse microRNA profiles determined with a new and sensitive cloning method.

MicroRNAs (miRNAs) are noncoding RNA molecules of 21 to 24 nt that regulate the expression of target genes in a post-transcriptional manner. Although evidence indicates that miRNAs play essential roles in embryogenesis, cell differentiation and pathogenesis of human diseases, extensive miRNA profiling in cells or tissues has been hampered by the lack of sensitive cloning methods. Here we describe a highly efficient profiling method, termed miRNA amplification profiling (mRAP), as well as its application both to mouse embryos at various developmental stages and to adult mouse organs. A total of 77,436 Small-RNA species was sequenced, with 11,776 of these sequences found to match previously described miRNAs. With the use of a newly developed computational prediction algorithm, we further identified 229 independent candidates for previously unknown miRNAs. The expression of some of these candidate miRNAs was confirmed by northern blot analysis and whole-mount in situ hybridization. Our data thus indicate that the total number of miRNAs in vertebrates is larger than previously appreciated and that the expression of these molecules is tightly controlled in a tissue- and developmental stage-specific manner.

Animals↗