Computational Biology and Language

Current scientific research is characterized by increasing specialization, accumulating knowledge at a high speed due to parallel advances in a multitude of sub-disciplines. Recent estimates suggest that human knowledge doubles every two to three years – and with the advances in information and communication technologies, this wide body of scientific knowledge is available to anyone, anywhere, anytime. This may also be referred to as ambient intelligence – an environment characterized by plentiful and available knowledge. The bottleneck in utilizing this knowledge for specific applications is not accessing but assimilating the information and transforming it to suit the needs for a specific application. The increasingly specialized areas of scientific research often have the common goal of converting data into insight allowing the identification of solutions to scientific problems. Due to this common goal, there are strong parallels between different areas of applications that can be exploited and used to cross-fertilize different disciplines. For example, the same fundamental statistical methods are used extensively in speech and language processing, in materials science applications, in visual processing and in biomedicine. Each sub-discipline has found its own specialized methodologies making these statistical methods successful to the given application. The unification of specialized areas is possible because many different problems can share strong analogies, making the theories developed for one problem applicable to other areas of research. It is the goal of this paper to demonstrate the utility of merging two disparate areas of applications to advance scientific research. The merging process requires cross-disciplinary collaboration to allow maximal exploitation of advances in one sub-discipline for that of another. We will demonstrate this general concept with the specific example of merging language technologies and computational biology.

Familial Diarrhea Syndrome Caused by an Activating GUCY2C Mutation

BACKGROUND Familial diarrhea disorders are, in most cases, severe and caused by recessive mutations. We describe the cause of a novel dominant disease in 32 members of a Norwegian family. The affected members have chronic diarrhea that is of early onset, is relatively mild, and is associated with increased susceptibility to inflammatory bowel disease, small-bowel obstruction, and esophagitis. METHODS We used linkage analysis, based on arrays with single-nucleotide polymorphisms, to identify a candidate region on chromosome 12 and then sequenced GUCY2C, encoding guanylate cyclase C (GC-C), an intestinal receptor for bacterial heat-stable enterotoxins. We performed exome sequencing of the entire candidate region from three affected family members, to exclude the possibility that mutations in genes other than GUCY2C could cause or contribute to susceptibility to the disease. We carried out functional studies of mutant GC-C using HEK293T cells. RESULTS We identified a heterozygous missense mutation (c.2519G -> T) in GUCY2C in all affected family members and observed no other rare variants in the exons of genes in the candidate region. Exposure of the mutant receptor to its ligands resulted in markedly increased production of cyclic guanosine monophosphate (cGMP). This may cause hyperactivation of the cystic fibrosis transmembrane regulator (CFTR), leading to increased chloride and water secretion from the enterocytes, and may thus explain the chronic diarrhea in the affected family members. CONCLUSIONS Increased GC-C signaling disturbs normal bowel function and appears to have a proinflammatory effect, either through increased chloride secretion or additional effects of elevated cellular cGMP. Further investigation of the relevance of genetic variants affecting the GC-C-CFTR pathway to conditions such as Crohn's disease is warranted. (Funded by Helse Vest Western Norway Regional Health Authority] and the Department of Science and Technology, Government of India.)

Mechanism of Fragility at BCL2 Gene Minor Breakpoint Cluster Region during t(14;18) Chromosomal Translocation

The t(14;18) translocation in follicular lymphoma is one of the most common chromosomal translocations. Breaks in chromosome 18 are localized at the 3'-UTR of BCL2 gene or downstream and are mainly clustered in either the major breakpoint region or the minor breakpoint cluster region (mcr). The recombination activating gene (RAG) complex induces breaks at IgH locus of chromosome 14, whereas the mechanism of fragility at BCL2 mcr remains unclear. Here, for the first time, we show that RAGs can nick mcr; however, the mechanism is unique. Three independent nicks of equal efficiency are generated, when both Mg2+ and Mn2+ are present, unlike a single nick during V(D)J recombination. Further, we demonstrate that RAG binding and nicking at the mcr are independent of nonamer, whereas a CCACCTCT motif plays a critical role in its fragility, as shown by sequential mutagenesis. More importantly, we recapitulate the BCL2 mcr translocation and find that mcr can undergo synapsis with a standard recombination signal sequence within the cells, in a RAG-dependent manner. Further, mutation to the CCACCTCT motif abolishes recombination within the cells, indicating its vital role. Hence, our data suggest a novel, physiologically relevant, nonamer-independent mechanism of RAG nicking at mcr, which may be important for generation of chromosomal translocations in humans.

Mutations in Drosophila Myosin Rod Cause Defects in Myofibril Assembly

The roles of myosin during muscle contraction are well studied, but how different domains of this protein are involved in myofibril assembly in vivo is far less understood. The indirect flight muscles (IFMs) of Drosophila melanogaster provide a good model for understanding muscle development and function in vivo. We show that two missense mutations in the rod region of the myosin heavy-chain gene, Mhc, give rise to IFM defects and abnormal myofibrils. These defects likely result from thick filament abnormalities that manifest during early sarcomere development or later by hypercontraction. The thick filament defects are accompanied by marked reduction in accumulation of flightin, a myosin binding protein, and its phosphorylated forms, which are required to stabilise thick filaments. We investigated with purified rod fragments whether the mutations affect the coiled-coil structure, rod aggregate size or rod stability. No significant changes in these parameters were detected, except for rod thermodynamic stability in one mutation. Molecular dynamics simulations suggest that these mutations may produce localised rod instabilities. We conclude that the aberrant myofibrils are a result of thick filament defects, but that these in vivo effects cannot be detected in vitro using the biophysical techniques employed. The in vivo investigation of these mutant phenotypes in IFM development and function provides a useful platform for studying myosin rod and thick filament formation generically, with application to the aetiology of human myosin rod myopathies. (C) 2012 Elsevier Ltd. All rights reserved.

Functional characterization of UvrD helicases from Haemophilus influenzae and Helicobacter pylori

Haemophilus influenzae and Helicobacter pylori are major bacterial pathogens that face high levels of genotoxic stress within their host. UvrD, a ubiquitous bacterial helicase that plays important roles in multiple DNA metabolic pathways, is essential for genome stability and might, therefore, be crucial in bacterial physiology and pathogenesis. In this study, the functional characterization of UvrD helicase from Haemophilus influenzae and Helicobacter pylori is reported. UvrD from Haemophilus influenzae (HiUvrD) and Helicobacter pylori (HpUvrD) exhibit strong single-stranded DNA-specific ATPase and 3'5' helicase activities. Mutation of highly conserved arginine (R288) in HiUvrD and glutamate (E206) in HpUvrD abrogated their activities. Both the proteins were able to bind and unwind a variety of DNA structures including duplexes with strand discontinuities and branches, three- and four-way junctions that underpin their role in DNA replication, repair and recombination. HiUvrD required a minimum of 12 nucleotides, whereas HpUvrD preferred 20 or more nucleotides of 3'-single-stranded DNA tail for efficient unwinding of duplex DNA. Interestingly, HpUvrD was able to hydrolyze and utilize GTP for its helicase activity although not as effectively as ATP, which has not been reported to date for UvrD characterized from other organisms. HiUvrD and HpUvrD were found to exist predominantly as monomers in solution together with multimeric forms. Noticeably, deletion of distal C-terminal 48 amino acid residues disrupted the oligomerization of HiUvrD, whereas deletion of 63 amino acids from C-terminus of HpUvrD had no effect on its oligomerization. This study presents the characteristic features and comparative analysis of Haemophilus influenzae and Helicobacter pylori UvrD, and constitutes the basis for understanding the role of UvrD in the biology and virulence of these pathogens.

Recent Developments in the Chemistry and Biology of G-Quadruplexes with Reference to the DNA Groove Binders

DNA is the chemotherapeutic target for treating diseases of genetic origin. Besides well-known double-helical structures (A, B, Z, parallel stranded-DNA etc.), DNA is capable of forming several multi-stranded structures (triplex, tetraplex, i-motif etc.) which have unique biological significance. The G-rich 3'-ends of chromosomes, called telomeres, are synthesized by telomerase, a ribonucleoprotein, and over-expression of telomerase is associated with cancer. The activity of telomerase is suppressed if the G-rich region is folded into the four stranded structures, called G-quadruplexes (G4-DNAs) using small synthetic ligands. Thus design and synthesis of new G4-DNA ligands is an attractive strategy to combat cancer. G4-DNA forming sequences are also prevalent in other genomic regions of biological significance including promoter regions of several oncogenes. Effective gene regulation may be achieved by inducing a G4-DNA structure within the G-rich promoter sequences. To date, several G4-DNA stabilizing ligands are known. DNA groove binders interact with the duplex B-DNA through the grooves (major and minor groove) in a sequence-specific manner. Some of the groove binders are known to stabilize the G4-DNA. However, this is a relatively under explored field of research. In this review, we focus on the recent advances in the understanding of the G4-DNA structures, particularly made from the human telomeric DNA stretches. We summarize the results of various investigations of the interaction of various organic ligands with the G4-DNA while highlighting the importance of groove binder-G4-DNA interactions.

IDH1 Mutations in Diffusely Infiltrating Astrocytomas Grade Specificity, Association With Protein Expression, and Clinical Relevance

IDH1 mutations are frequent genetic alterations in low-grade diffuse gliomas and secondary glioblastoma (GBM). To validate mutation frequency, IDH1 gene at codon 132 was sequenced in 74 diffusely infiltrating astrocytomas: diffuse astrocytoma (DA; World Health Organization WHO] grade II), anaplastic astrocytoma (AA; WHO grade III), and GBM (WHO grade IV). All cases were immunostained with IDH1-R132H monoclonal antibody. Mutational status was correlated with mutant protein expression, patient age, duration of symptoms, and prognosis of patients with GBM. We detected 31 (41.9%) heterozygous IDH1 mutations resulting in arginine-to-histidine substitution (R132H;CGT-CAT). All 12 DAs (100%), 13 of 14 AAs (92.9%), and 6 of 48 GBMs (12.5%) (5/6 83.3%] secondary, and 1/42 2.4%] primary) harbored IDH1 mutations. The correlation between mutational status and protein expression was significant (P < .001). IDH1 mutation status, though not associated with prognosis of patients with GBM, showed significant association with younger age and longer duration of symptoms in the whole cohort (P < .001). Our study validates IDH1 mutant protein expression across various grades of astrocytoma, and demonstrates a high incidence of IDH1 mutations in DA, AA, and secondary GBM.

Stochastic Simulations Suggest that HIV-1 Survives Close to Its Error Threshold

The use of mutagenic drugs to drive HIV-1 past its error threshold presents a novel intervention strategy, as suggested by the quasispecies theory, that may be less susceptible to failure via viral mutation-induced emergence of drug resistance than current strategies. The error threshold of HIV-1, mu(c), however, is not known. Application of the quasispecies theory to determine mu(c) poses significant challenges: Whereas the quasispecies theory considers the asexual reproduction of an infinitely large population of haploid individuals, HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We performed population genetics-based stochastic simulations of the within-host evolution of HIV-1 and estimated the structure of the HIV-1 quasispecies and mu(c). We found that with small mutation rates, the quasispecies was dominated by genomes with few mutations. Upon increasing the mutation rate, a sharp error catastrophe occurred where the quasispecies became delocalized in sequence space. Using parameter values that quantitatively captured data of viral diversification in HIV-1 patients, we estimated mu(c) to be 7 x 10(-5) -1 x 10(-4) substitutions/site/replication, similar to 2-6 fold higher than the natural mutation rate of HIV-1, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs. The latter estimate was weakly dependent on the within-host effective population size of HIV-1. With large population sizes and in the absence of recombination, our simulations converged to the quasispecies theory, bridging the gap between quasispecies theory and population genetics-based approaches to describing HIV-1 evolution. Further, mu(c) increased with the recombination rate, rendering HIV-1 less susceptible to error catastrophe, thus elucidating an added benefit of recombination to HIV-1. Our estimate of mu(c) may serve as a quantitative guideline for the use of mutagenic drugs against HIV-1.

Site-specific N-Linked Glycosylation of Receptor Guanylyl Cyclase C Regulates Ligand Binding, Ligand-mediated Activation and Interaction with Vesicular Integral Membrane Protein 36, VIP36

Guanylyl cyclase C (GC-C) is a multidomain, membrane-associated receptor guanylyl cyclase. GC-C is primarily expressed in the gastrointestinal tract, where it mediates fluid-ion homeostasis, intestinal inflammation, and cell proliferation in a cGMP-dependent manner, following activation by its ligands guanylin, uroguanylin, or the heat-stable enterotoxin peptide (ST). GC-C is also expressed in neurons, where it plays a role in satiation and attention deficiency/hyperactive behavior. GC-C is glycosylated in the extracellular domain, and differentially glycosylated forms that are resident in the endoplasmic reticulum (130 kDa) and the plasma membrane (145 kDa) bind the ST peptide with equal affinity. When glycosylation of human GC-C was prevented, either by pharmacological intervention or by mutation of all of the 10 predicted glycosylation sites, ST binding and surface localization was abolished. Systematic mutagenesis of each of the 10 sites of glycosylation in GC-C, either singly or in combination, identified two sites that were critical for ligand binding and two that regulated ST-mediated activation. We also show that GC-C is the first identified receptor client of the lectin chaperone vesicular integral membrane protein, VIP36. Interaction with VIP36 is dependent on glycosylation at the same sites that allow GC-C to fold and bind ligand. Because glycosylation of proteins is altered in many diseases and in a tissue-dependent manner, the activity and/or glycan-mediated interactions of GC-C may have a crucial role to play in its functions in different cell types.

Biochemical Properties of MutT2 Proteins from Mycobacterium tuberculosis and M. smegmatis and Their Contrasting Antimutator Roles in Escherichia coli

Mycobacterium tuberculosis, the causative agent of tuberculosis, is at increased risk of accumulating damaged guanine nucleotides such as 8-oxo-dGTP and 8-oxo-GTP because of its residency in the oxidative environment of the host macrophages. By hydrolyzing the oxidized guanine nucleotides before their incorporation into nucleic acids, MutT proteins play a critical role in allowing organisms to avoid their deleterious effects. Mycobacteria possess several MutT proteins. Here, we purified recombinant M. tuberculosis MutT2 (MtuMutT2) and M. smegmatis MutT2 (MsmMutT2) proteins from M. tuberculosis (a slow grower) and M. smegmatis (fast growing model mycobacteria), respectively, for their biochemical characterization. Distinct from the Escherichia coli MutT, which hydrolyzes 8-oxo-dGTP and 8-oxo-GTP, the mycobacterial proteins hydrolyze not only 8-oxo-dGTP and 8-oxo-GTP but also dCTP and 5-methyl-dCTP. Determination of kinetic parameters (K-m and V-max) revealed that while MtuMutT2 hydrolyzes dCTP nearly four times better than it does 8-oxo-dGTP, MsmMutT2 hydrolyzes them nearly equally. Also, MsmMutT2 is about 14 times more efficient than MtuMutT2 in its catalytic activity of hydrolyzing 8-oxo-dGTP. Consistent with these observations, MsmMutT2 but not MtuMutT2 rescues E. coli for MutT deficiency by decreasing both the mutation frequency and A-to-C mutations (a hallmark of MutT deficiency). We discuss these findings in the context of the physiological significance of MutT proteins.

Human la protein interaction with GCAC near the initiator AUG enhances hepatitis C virus RNA replication by promoting linkage between 5 ` and 3 ` untranslated regions

Human La protein has been implicated in facilitating the internal initiation of translation as well as replication of hepatitis C virus (HCV) RNA. Previously, we demonstrated that La interacts with the HCV internal ribosome entry site (IRES) around the GCAC motif near the initiator AUG within stem-loop IV by its RNA recognition motif (RRM) (residues 112 to 184) and influences HCV translation. In this study, we have deciphered the role of this interaction in HCV replication in a hepatocellular carcinoma cell culture system. We incorporated mutation of the GCAC motif in an HCV monocistronic subgenomic replicon and a pJFH1 construct which altered the binding of La and checked HCV RNA replication by reverse transcriptase PCR (RT-PCR). The mutation drastically affected HCV replication. Furthermore, to address whether the decrease in replication is a consequence of translation inhibition or not, we incorporated the same mutation into a bicistronic replicon and observed a substantial decrease in HCV RNA levels. Interestingly, La overexpression rescued this inhibition of replication. More importantly, we observed that the mutation reduced the association between La and NS5B. The effect of the GCAC mutation on the translation-to-replication switch, which is regulated by the interplay between NS3 and La, was further investigated. Additionally, our analyses of point mutations in the GCAC motif revealed distinct roles of each nucleotide in HCV replication and translation. Finally, we showed that a specific interaction of the GCAC motif with human La protein is crucial for linking 5' and 3' ends of the HCV genome. Taken together, our results demonstrate the mechanism of regulation of HCV replication by interaction of the cis-acting element GCAC within the HCV IRES with human La protein.

Heat and SDS insensitive NDK dimers are largely stabilised by hydrophobic interaction to form functional hexamer in Mycobacterium smegmatis

The primary structure and function of nucleoside diphosphate kinase (NDK), a substrate non-specific enzyme involved in the maintenance of nucleotide pools is also implicated to play pivotal roles in many other cellular processes. NDK is conserved from bacteria to human and forms a homotetramer or hexamer to exhibit its biological activity. However, the nature of the functional oligomeric form of the enzyme differs among different organisms. The functional form of NDKs from many bacterial systems, including that of the human pathogen, Mycobacterium tuberculosis (MtuNDK), is a hexamer, although some bacterial NDKs are tetrameric in nature. The present study addresses the oligomeric property of MsmNDK and how a dimer, the basic subunit of a functional hexamer, is stabilized by hydrogen bonds and hydrophobic interactions. Homology modeling was generated using the three-dimensional structure of MtuNDK as a template; the residues interacting at the monomer-monomer interface of MsmNDK were mapped. Using recombinant enzymes of wild type, catalytically inactive mutant, and monomer-monomer interactive mutants of MsmNDK, the stability of the dimer was verified under heat, SDS, low pH, and methanol. The predicted residues (Gln17, Ser24 and Glu27) were engaged in dimer formation, however the mutated proteins retained the ATPase and GTPase activity even after introducing single (MsmNDK- Q17A, MsmNDK-E27A, and MsmNDK-E27Q) and double (MsmNDK-E27A/Q17A) mutation. However, the monomer monomer interaction could be abolished using methanol, indicating the stabilization of the monomer-monomer interaction by hydrophobic interaction.

Mutation analysis of the SLC4A11 gene in Indian families with congenital hereditary endothelial dystrophy 2 and a review of the literature

Purpose: Congenital hereditary endothelial dystrophy 2 (CHED2) is an autosomal recessive disorder caused by mutations in the solute carrier family 4, sodium borate transporter, member 11 (SLC4A11) gene. The purpose of this study was to identify the genetic cause of CHED2 in six Indian families and catalog all known mutations in the SLC4A11 gene. Methods: Peripheral blood samples were collected from individuals of the families with CHED2 and used in genomic DNA isolation. PCR primers were used to amplify the entire coding region including intron-exon junctions of SLC4A11. Amplicons were subsequently sequenced to identify the mutations. Results: DNA sequence analysis of the six families identified four novel (viz., p.Thr262Ile, p.Gly417Arg, p.Cys611Arg, and p.His724Asp) mutations and one known p.Arg869His homozygous mutation in the SLC4A11 gene. The mutation p.Gly417Arg was identified in two families. Conclusions: This study increases the mutation spectrum of the SLC4A11 gene. A review of the literature showed that the total number of mutations in the SLC4A11 gene described to date is 78. Most of the mutations are missense, followed by insertions-deletions. The present study will be helpful in genetic diagnosis of the families reported here.

G-Quadruplex Structures Formed at the HOX11 Breakpoint Region Contribute to Its Fragility during t(10;14) Translocation in T-Cell Leukemia

The t(10;14) translocation involving the HOX11 gene is found in several T-cell leukemia patients. Previous efforts to determine the causes of HOX11 fragility were not successful. The role of non-B DNA structures is increasingly becoming an important cause of genomic instability. In the present study, bioinformatics analysis revealed two G-quadruplex-forming motifs at the HOX11 breakpoint cluster. Gel shift assays showed formation of both intra- and intermolecular G-quadruplexes, the latter being more predominant. The structure formation was dependent on four stretches of guanines, as revealed by mutagenesis. Circular dichroism analysis identified parallel conformations for both quadruplexes. The non-B DNA structure could block polymerization during replication on a plasmid, resulting in consistent K K+-dependent pause sites, which were abolished upon mutation of G-motifs, thereby demonstrating the role of the stretches of guanines even on double-stranded DNA. Extrachromosomal assays showed that the G-quadruplex motifs could block transcription, leading to reduced expression of green fluorescent protein (GFP) within cells. More importantly, sodium bisulfite modification assay showed the single-stranded character at regions I and II of HOX11 in the genome. Thus, our findings suggest the occurrence of G-quadruplex structures at the HOX11 breakpoint region, which could explain its fragility during the t(10;14) translocation.