From Nonspecific DNA–Protein Encounter Complexes to the Prediction of DNA–Protein Interactions

©2009 Gao, Skolnick. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.doi:10.1371/journal.pcbi.1000341DNA–protein interactions are involved in many essential biological activities. Because there is no simple mapping code between DNA base pairs and protein amino acids, the prediction of DNA–protein interactions is a challenging problem. Here, we present a novel computational approach for predicting DNA-binding protein residues and DNA–protein interaction modes without knowing its specific DNA target sequence. Given the structure of a DNA-binding protein, the method first generates an ensemble of complex structures obtained by rigid-body docking with a nonspecific canonical B-DNA. Representative models are subsequently selected through clustering and ranking by their DNA–protein interfacial energy. Analysis of these encounter complex models suggests that the recognition sites for specific DNA binding are usually favorable interaction sites for the nonspecific DNA probe and that nonspecific DNA–protein interaction modes exhibit some similarity to specific DNA–protein binding modes. Although the method requires as input the knowledge that the protein binds DNA, in benchmark tests, it achieves better performance in identifying DNA-binding sites than three previously established methods, which are based on sophisticated machine-learning techniques. We further apply our method to protein structures predicted through modeling and demonstrate that our method performs satisfactorily on protein models whose root-mean-square Ca deviation from native is up to 5 Å from their native structures. This study provides valuable structural insights into how a specific DNA-binding protein interacts with a nonspecific DNA sequence. The similarity between the specific DNA–protein interaction mode and nonspecific interaction modes may reflect an important sampling step in search of its specific DNA targets by a DNA-binding protein

Investigation of the Interaction between the Large and Small Subunits of Potato ADP-Glucose Pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase), a key allosteric enzyme involved in higher plant starch biosynthesis, is composed of pairs of large (LS) and small subunits (SS). Current evidence indicates that the two subunit types play distinct roles in enzyme function. Recently the heterotetrameric structure of potato AGPase has been modeled. In the current study, we have applied the molecular mechanics generalized born surface area (MM-GBSA) method and identified critical amino acids of the potato AGPase LS and SS subunits that interact with each other during the native heterotetrameric structure formation. We have further shown the role of the LS amino acids in subunit-subunit interaction by yeast two-hybrid, bacterial complementation assay and native gel. Comparison of the computational results with the experiments has indicated that the backbone energy contribution (rather than the side chain energies) of the interface residues is more important in identifying critical residues. We have found that lateral interaction of the LS-SS is much stronger than the longitudinal one, and it is mainly mediated by hydrophobic interactions. This study will not only enhance our understanding of the interaction between the SS and the LS of AGPase, but will also enable us to engineer proteins to obtain better assembled variants of AGPase which can be used for the improvement of plant yield

Advanced Genomic Data Mining

As data banks increase their size, one of the current challenges in bioinformatics is to be able to query them in a sensible way. Information is contained in differen

Detection of a MicroRNA Signal in an In Vivo Expression Set of mRNAs

Background. microRNAs (miRNAs) are approximately 21 nucleotide non-coding transcripts capable of regulating gene expression. The most widely studied mechanism of regulation involves binding of a miRNA to the target mRNA. As a result, translation of the target mRNA is inhibited and the mRNA may be destabilized. The inhibitory effects of miRNAs have been linked to diverse cellular processes including malignant proliferation, apoptosis, development, differentiation, and metabolic processes. We asked whether endogenous fluctuations in a set of mRNA and miRNA profiles contain correlated changes that are statistically distinguishable from the many other fluctuations in the data set. Methodology/Principal Findings. RNA was extracted from 12 human primary brain tumor biopsies. These samples were used to determine genome-wide mRN

Ortho2ExpressMatrix—a web server that interprets cross-species gene expression data by gene family information

<p>Abstract</p> <p>Background</p> <p>The study of gene families is pivotal for the understanding of gene evolution across different organisms and such phylogenetic background is often used to infer biochemical functions of genes. Modern high-throughput experiments offer the possibility to analyze the entire transcriptome of an organism; however, it is often difficult to deduct functional information from that data.</p> <p>Results</p> <p>To improve functional interpretation of gene expression we introduce Ortho2ExpressMatrix, a novel tool that integrates complex gene family information, computed from sequence similarity, with comparative gene expression profiles of two pre-selected biological objects: gene families are displayed with two-dimensional matrices. Parameters of the tool are object type (two organisms, two individuals, two tissues, etc.), type of computational gene family inference, experimental meta-data, microarray platform, gene annotation level and genome build. Family information in Ortho2ExpressMatrix bases on computationally different protein family approaches such as EnsemblCompara, InParanoid, SYSTERS and Ensembl Family. Currently, respective all-against-all associations are available for five species: human, mouse, worm, fruit fly and yeast. Additionally, microRNA expression can be examined with respect to miRBase or TargetScan families. The visualization, which is typical for Ortho2ExpressMatrix, is performed as matrix view that displays functional traits of genes (differential expression) as well as sequence similarity of protein family members (BLAST e-values) in colour codes. Such translations are intended to facilitate the user's perception of the research object.</p> <p>Conclusions</p> <p>Ortho2ExpressMatrix integrates gene family information with genome-wide expression data in order to enhance functional interpretation of high-throughput analyses on diseases, environmental factors, or genetic modification or compound treatment experiments. The tool explores differential gene expression in the light of orthology, paralogy and structure of gene families up to the point of ambiguity analyses. Results can be used for filtering and prioritization in functional genomic, biomedical and systems biology applications. The web server is freely accessible at <url>http://bioinf-data.charite.de/o2em/cgi-bin/o2em.pl</url>.</p

Low potency toxins reveal dense interaction networks in metabolism

Background The chemicals of metabolism are constructed of a small set of atoms and bonds. This may be because chemical structures outside the chemical space in which life operates are incompatible with biochemistry, or because mechanisms to make or utilize such excluded structures has not evolved. In this paper I address the extent to which biochemistry is restricted to a small fraction of the chemical space of possible chemicals, a restricted subset that I call Biochemical Space. I explore evidence that this restriction is at least in part due to selection again specific structures, and suggest a mechanism by which this occurs. Results Chemicals that contain structures that our outside Biochemical Space (UnBiological groups) are more likely to be toxic to a wide range of organisms, even though they have no specifically toxic groups and no obvious mechanism of toxicity. This correlation of UnBiological with toxicity is stronger for low potency (millimolar) toxins. I relate this to the observation that most chemicals interact with many biological structures at low millimolar toxicity. I hypothesise that life has to select its components not only to have a specific set of functions but also to avoid interactions with all the other components of life that might degrade their function. Conclusions The chemistry of life has to form a dense, self-consistent network of chemical structures, and cannot easily be arbitrarily extended. The toxicity of arbitrary chemicals is a reflection of the disruption to that network occasioned by trying to insert a chemical into it without also selecting all the other components to tolerate that chemical. This suggests new ways to test for the toxicity of chemicals, and that engineering organisms to make high concentrations of materials such as chemical precursors or fuels may require more substantial engineering than just of the synthetic pathways involved

Peripheral and central arterial pressure and its relationship to vascular target organ damage in carotid artery, retina and arterial stiffness. Development and validation of a tool. The Vaso risk study

<p>Abstract</p> <p>Background</p> <p>Ambulatory blood pressure monitoring (ABPM) shows a better correlation to target organ damage and cardiovascular morbidity-mortality than office blood pressure. A loss of arterial elasticity and an increase in carotid artery intima-media thickness (IMT) has been associated with increased cardiovascular morbidity-mortality. Tools have been developed that allow estimation of the retinal arteriovenous index but not all studies coincide and there are contradictory results in relation to the evolution of the arteriosclerotic lesions and the caliber of the retinal vessels. The purpose of this study is to analyze the relationship between peripheral and central arterial pressure (clinic and ambulatory) and vascular structure and function as evaluated by the carotid artery intima-media thickness, retina arteriovenous index, pulse wave velocity (PWV) and ankle-brachial index in patients with and without type 2 diabetes. In turn, software is developed and validated for measuring retinal vessel thickness and automatically estimating the arteriovenous index.</p> <p>Methods/Design</p> <p>A cross-sectional study involving a control group will be made, with a posterior 4-year follow-up period in primary care. The study patients will be type 2 diabetics, with a control group of non-diabetic individuals. Consecutive sampling will be used to include 300 patients between 34-75 years of age and no previous cardiovascular disease, one-half being assigned to each group. Main measurements: age, gender, height, weight and abdominal circumference. Lipids, creatinine, microalbuminuria, blood glucose, HbA1c, blood insulin, high sensitivity C-reactive protein and endothelial dysfunction markers. Clinic and ambulatory blood pressure monitoring. Carotid ultrasound to evaluate IMT, and retinography to evaluate the arteriovenous index. ECG to assess left ventricle hypertrophy, ankle-brachial index, and pulse wave analysis (PWA) and pulse wave velocity (PWV) with the Sphigmocor System.</p> <p>Discussion</p> <p>We hope to obtain information on the correlation of different ABPM-derived parameters and PWA to organ target damage - particularly vascular structure and function evaluated from the IMT and PWV - and endothelial dysfunction in patients with and without type 2 diabetes. We also hope to demonstrate the usefulness of the instrument developed for the automated evaluation of retinal vascularization in the early detection of alterations in vascular structure and function and in the prognosis of middle-term cardiovascular morbidity.</p> <p>Trial Registration</p> <p>Clinical Trials.gov Identifier: <a href="http://www.clinicaltrials.gov/ct2/show/NCT01325064">NCT01325064</a></p

Ceacam1 separates graft-versus-host-disease from graft-versus-tumor activity after experimental allogeneic bone marrow transplantation.

BACKGROUND: Allogeneic bone marrow transplantation (allo-BMT) is a potentially curative therapy for a variety of hematologic diseases, but benefits, including graft-versus-tumor (GVT) activity are limited by graft-versus-host-disease (GVHD). Carcinoembryonic antigen related cell adhesion molecule 1 (Ceacam1) is a transmembrane glycoprotein found on epithelium, T cells, and many tumors. It regulates a variety of physiologic and pathological processes such as tumor biology, leukocyte activation, and energy homeostasis. Previous studies suggest that Ceacam1 negatively regulates inflammation in inflammatory bowel disease models. METHODS: We studied Ceacam1 as a regulator of GVHD and GVT after allogeneic bone marrow transplantation (allo-BMT) in mouse models. In vivo, Ceacam1(-/-) T cells caused increased GVHD mortality and GVHD of the colon, and greater numbers of donor T cells were positive for activation markers (CD25(hi), CD62L(lo)). Additionally, Ceacam1(-/-) CD8 T cells had greater expression of the gut-trafficking integrin α(4)β(7), though both CD4 and CD8 T cells were found increased numbers in the gut post-transplant. Ceacam1(-/-) recipients also experienced increased GVHD mortality and GVHD of the colon, and alloreactive T cells displayed increased activation. Additionally, Ceacam1(-/-) mice had increased mortality and decreased numbers of regenerating small intestinal crypts upon radiation exposure. Conversely, Ceacam1-overexpressing T cells caused attenuated target-organ and systemic GVHD, which correlated with decreased donor T cell numbers in target tissues, and mortality. Finally, graft-versus-tumor survival in a Ceacam1(+) lymphoma model was improved in animals receiving Ceacam1(-/-) vs. control T cells. CONCLUSIONS: We conclude that Ceacam1 regulates T cell activation, GVHD target organ damage, and numbers of donor T cells in lymphoid organs and GVHD target tissues. In recipients of allo-BMT, Ceacam1 may also regulate tissue radiosensitivity. Because of its expression on both the donor graft and host tissues, this suggests that targeting Ceacam1 may represent a potent strategy for the regulation of GVHD and GVT after allogeneic transplantation

Comparative Analysis of Serine/Arginine-Rich Proteins across 27 Eukaryotes: Insights into Sub-Family Classification and Extent of Alternative Splicing

Alternative splicing (AS) of pre-mRNA is a fundamental molecular process that generates diversity in the transcriptome and proteome of eukaryotic organisms. SR proteins, a family of splicing regulators with one or two RNA recognition motifs (RRMs) at the N-terminus and an arg/ser-rich domain at the C-terminus, function in both constitutive and alternative splicing. We identified SR proteins in 27 eukaryotic species, which include plants, animals, fungi and “basal” eukaryotes that lie outside of these lineages. Using RNA recognition motifs (RRMs) as a phylogenetic marker, we classified 272 SR genes into robust sub-families. The SR gene family can be split into five major groupings, which can be further separated into 11 distinct sub-families. Most flowering plants have double or nearly double the number of SR genes found in vertebrates. The majority of plant SR genes are under purifying selection. Moreover, in all paralogous SR genes in Arabidopsis, rice, soybean and maize, one of the two paralogs is preferentially expressed throughout plant development. We also assessed the extent of AS in SR genes based on a splice graph approach (http://combi.cs.colostate.edu/as/gmap_SRgenes). AS of SR genes is a widespread phenomenon throughout multiple lineages, with alternative 3′ or 5′ splicing events being the most prominent type of event. However, plant-enriched sub-families have 57%–88% of their SR genes experiencing some type of AS compared to the 40%–54% seen in other sub-families. The SR gene family is pervasive throughout multiple eukaryotic lineages, conserved in sequence and domain organization, but differs in gene number across lineages with an abundance of SR genes in flowering plants. The higher number of alternatively spliced SR genes in plants emphasizes the importance of AS in generating splice variants in these organisms